The contour method is a new approach to calculating the non-adiabatic pulsation frequencies of stars. These frequencies can be found by solving for the complex roots of a characteristic equation constructed from the linear non-adiabatic stellar pulsation equations. A complex-root solver requires an initial trial frequency for each non-adiabatic root. A standard method for obtaining initial trial frequencies is to use a star's adiabatic pulsation frequencies, but this method can fail to converge to non-adiabatic roots, especially as the growth and/or damping rate of the pulsations becomes large. The contour method provides an alternative way for obtaining initial trial frequencies that robustly converges to non-adiabatic roots, even for stellar models with extremely non-adiabatic pulsations and thus large growth/damping rates. We describe the contour method implemented in the GYRE stellar pulsation code and use it to calculate the non-adiabatic pulsation frequencies of $10 \, \rm M_\odot$ and $20 \,\rm M_\odot$ $\beta$ Cephei models, and of a $0.9 \, \rm M_\odot$ extreme helium model.

Laplace's tidal equations govern the angular dependence of oscillations in stars when uniform rotation is treated within the so-called traditional approximation. Using a perturbation expansion approach, I derive improved expressions for the eigenvalue associated with these equations, valid in the asymptotic limit of large spin parameter $q$. These expressions have a relative accuracy of order $q^{-3}$ for gravito-inertial modes, and $q^{-1}$ for Rossby and Kelvin modes; the corresponding absolute accuracy is of order $q^{-1}$ for all three mode types. I validate the expressions against numerical calculations, and demonstrate how they can be applied to derive a formula for the periods of Rossby modes.

Context. The large quantity of high-quality asteroseismic data that obtained from space-based photometric missions and the accuracy of the resulting frequencies motivate a careful consideration of the accuracy of computed oscillation frequencies of stellar models, when applied as diagnostics of the model properties. Aims. Based on models of red-giant stars that have been independently calculated using different stellar evolution codes, we investigate the extent to which the differences in the model calculation affect the model oscillation frequencies. Methods. For each of the models, which cover four different masses and different evolution stages on the red-giant branch, we computed full sets of low-degree oscillation frequencies using a single pulsation code and, from these frequencies, typical asteroseismic diagnostics. In addition, we carried out preliminary analyses to relate differences in the oscillation properties to the corresponding model differences. Results. In general, the differences in asteroseismic properties between the different models greatly exceed the observational precision of these properties, in particular for the nonradial modes whose mixed acoustic and gravity-wave character makes them sensitive to the structure of the deep stellar interior. In some cases, identifying these differences led to improvements in the final models presented here and in Paper I; here we illustrate particular examples of this. Conclusions. Further improvements in stellar modelling are required in order fully to utilise the observational accuracy to probe intrinsic limitations in the modelling. However, our analysis of the frequency differences and their relation to stellar internal properties provides a striking illustration of the potential of the mixed modes of red-giant stars for the diagnostics of stellar interiors.

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.

Astronomy, as a field, has long encouraged the development of free, open digital infrastructure (e.g., National Research Council 2010, 2011). Examples range from simple scripts that enable individual scientific research, through software instruments for entire communities, to data reduction pipelines for telescope operations at national facilities. As with the digital infrastructure of our larger society today (e.g., Eghbal 2016), nearly all astronomical research relies on free, open source software (FOSS) written and maintained by a small number of developers. And like the physical infrastructure of roads or bridges, digital infrastructure needs regular upkeep and maintenance (e.g., Eghbal 2016). In astronomy, financial support for maintaining existing digital infrastructure is generally much harder to secure than funding for developing new digital infrastructures that promise new science. Sustaining astronomy's digital infrastructure is a new topic for many, the sustainability challenges are not always widely known...

(Abridged). We introduce the Aarhus Red Giants Challenge, a series of detailed comparisons between widely used stellar evolution and oscillation codes aiming at establishing the minimum level of uncertainties in properties of red giants arising solely from numerical implementations. Using 9 state-of-the-art stellar evolution codes, we defined a set of input physics and physical constants for our calculations and calibrated the convective efficiency to a specific point on the main sequence. We produced evolutionary tracks and stellar structure models at fixed radius along the red-giant branch for masses of 1.0 M$_\odot$, 1.5 M$_\odot$, 2.0 M$_\odot$, and 2.5 M$_\odot$, and compared the predicted stellar properties. Once models have been calibrated on the main sequence we find a residual spread in the predicted effective temperatures across all codes of ~20 K at solar radius and ~30-40 K in the RGB regardless of the considered stellar mass. The predicted ages show variations of 2-5% (increasing with stellar mass) which we track down to differences in the numerical implementation of energy generation. The luminosity of the RGB-bump shows a spread of about 10% for the considered codes, which translates into magnitude differences of ~0.1 mag in the optical V-band. We also compare the predicted [C/N] abundance ratio and found a spread of 0.1 dex or more for all considered masses. Our comparisons show that differences at the level of a few percent still remain in evolutionary calculations of red giants branch stars despite the use of the same input physics. These are mostly due to differences in the energy generation routines and interpolation across opacities, and call for further investigations on these matters in the context of using properties of red giants as benchmarks for astrophysical studies.

The Tayler instability (TI) is a non-axisymmetric linear instability of an axisymmetric toroidal magnetic field in magneto-hydrostatic equilibrium (MHSE). Spruit (1999, 2002) has proposed that in a differentially rotating radiative region of a star, the TI drives a dynamo which generates magnetic fields that can efficiently transport angular momentum; a parameterized version of this dynamo has been implemented in stellar structure and evolution codes and shown to be important for determining interior spin. Numerical simulations, however, have yet to definitively demonstrate the operation of the dynamo. A criterion for the MHSE to develop the TI was derived using fully-compressible magneto-hydrodynamics, while numerical simulations of dynamical processes in stars frequently use an anelastic approximation. This motivates us to derive a new anelastic Tayler instability (anTI) criterion. We find that some MHSE configurations unstable in the fully-compressible case, become stable in the anelastic case. We find and characterize the unstable modes of a simple family of cylindrical MHSE configurations using numerical calculations, and discuss the implications for fully non-linear anelastic simulations.

We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). RSP is a new functionality in MESAstar that models the non-linear radial stellar pulsations that characterize RR Lyrae, Cepheids, and other classes of variable stars. We significantly enhance numerical energy conservation capabilities, including during mass changes. For example, this enables calculations through the He flash that conserve energy to better than 0.001 %. To improve the modeling of rotating stars in MESA, we introduce a new approach to modifying the pressure and temperature equations of stellar structure, and a formulation of the projection effects of gravity darkening. A new scheme for tracking convective boundaries yields reliable values of the convective-core mass, and allows the natural emergence of adiabatic semiconvection regions during both core hydrogen- and helium-burning phases. We quantify the parallel performance of MESA on current generation multicore architectures and demonstrate improvements in the computational efficiency of radiative levitation. We report updates to the equation of state and nuclear reaction physics modules. We briefly discuss the current treatment of fallback in core-collapse supernova models and the thermodynamic evolution of supernova explosions. We close by discussing the new MESA Testhub software infrastructure to enhance source-code development.

KIC 08626021 is a pulsating DB white dwarf of considerable recent interest, and first of its class to be extensively monitored by Kepler for its pulsation properties. Fitting the observed oscillation frequencies of KIC 08626021 to a model can yield insights into its otherwise-hidden internal structure. Template-based white dwarf models choose a luminosity profile where the luminosity is proportional to the enclosed mass, $L_r \propto M_r$, independent of the effective temperature $T_{\rm eff}$. Evolutionary models of young white dwarfs with $T_{\rm eff} \gtrsim$ 25,000 K suggest neutrino emission gives rise to luminosity profiles with $L_r$ $\not\propto$ $M_r$. We explore this contrast by comparing the oscillation frequencies between two nearly identical white dwarf models: one with an enforced $L_r \propto M_r$ luminosity profile and the other with a luminosity profile determined by the star's previous evolution history. We find the low order g-mode frequencies differ by up to $\simeq$ 70 $\mu$Hz over the range of Kepler observations for KIC 08626021. This suggests that by neglecting the proper thermal structure of the star (e.g., accounting for the effect of plasmon neutrino losses), the model frequencies calculated by using an $L_r \propto M_r$ profile may have uncorrected, effectively-random errors at the level of tens of $\mu$Hz. A mean frequency difference of 30 $\mu$Hz, based on linearly extrapolating published results, suggests a template model uncertainty in the fit precision of $\simeq$ 12% in white dwarf mass, $\simeq$ 9% in the radius, and $\simeq$ 3% in the central oxygen mass fraction.

Context. While rotation has a major impact on stellar structure and evolution, its effects are not well understood. Thanks to high- quality and long timebase photometric observations obtained with recent space missions, we are now able to study stellar rotation more precisely. Aims. We aim to constrain radial differential rotation profiles in gamma Doradus (gamma Dor) stars, and to develop new theoretical seismic diagnosis for such stars with rapid and potentially non-uniform rotation. Methods. We derive a new asymptotic description which accounts for the impact of weak differential near-core rotation on gravity- mode period spacings. The theoretical predictions are illustrated from pulsation computations with the code GYRE and compared with observations of gamma Dor stars. When possible, we also derive the surface rotation rates in these stars by detecting and analysing signatures of rotational modulation, and compute the core-to-surface rotation ratios. Results. Stellar rotation has to be strongly differential before its effects on period spacing patterns can be detected, unless multiple period spacing patterns can be compared. Six stars in our sample exhibit a single unexplained period spacing pattern of retrograde modes. We hypothesise that these are Yanai modes. Finally, we find signatures of rotational spot modulation in the photometric data of eight targets. Conclusions. If only one period spacing pattern is detected and analysed for a star, it is difficult to detect differential rotation. A rigidly rotating model will often provide the best solution. Differential rotation can only be detected when multiple period spacing patterns have been found for a single star or its surface rotation rate is known as well. This is the case for eight stars in our sample, revealing surface-to-core rotation ratios between 0.95 and 1.05.

After an optical peak, a classical or recurrent nova settles into a brief (days to years) period of quasi-stable thermonuclear burning in a compact configuration nearly at the white dwarf (WD) radius. During this time, the underlying WD becomes visible as a strong emitter of supersoft X-rays. Observations during this phase have revealed oscillations in the X-ray emission with periods on the order of tens of seconds. A proposed explanation for the source of these oscillations are internal gravity waves excited by nuclear reactions at the base of the hydrogen-burning layer. In this work, we present the first models exhibiting unstable surface $g$-modes with periods similar to oscillation periods found in galactic novae. However, when comparing mode periods of our models to the observed oscillations of several novae, we find that the modes which are excited have periods shorter than that observed.

We update the capabilities of the software instrument Modules for Experiments in Stellar Astrophysics (MESA) and enhance its ease of use and availability. Our new approach to locating convective boundaries is consistent with the physics of convection, and yields reliable values of the convective core mass during both hydrogen and helium burning phases. Stars with $M<8\,{\rm M_\odot}$ become white dwarfs and cool to the point where the electrons are degenerate and the ions are strongly coupled, a realm now available to study with MESA due to improved treatments of element diffusion, latent heat release, and blending of equations of state. Studies of the final fates of massive stars are extended in MESA by our addition of an approximate Riemann solver that captures shocks and conserves energy to high accuracy during dynamic epochs. We also introduce a 1D capability for modeling the effects of Rayleigh-Taylor instabilities that, in combination with the coupling to a public version of the STELLA radiation transfer instrument, creates new avenues for exploring Type II supernovae properties. These capabilities are exhibited with exploratory models of pair-instability supernova, pulsational pair-instability supernova, and the formation of stellar mass black holes. The applicability of MESA is now widened by the capability of importing multi-dimensional hydrodynamic models into MESA. We close by introducing software modules for handling floating point exceptions and stellar model optimization, and four new software tools -- MESAWeb, MESA-Docker, pyMESA, and mesastar.org -- to enhance MESA's education and research impact.

Motivated by recent interest in the phenomenon of waves transport in massive stars, we examine whether the heat-driven gravity (g) modes excited in slowly-pulsating B (SPB) stars can significantly modify the stars' internal rotation. We develop a formalism for the differential torque exerted by g modes, and implement this formalism using the GYRE oscillation code and the MESASTAR stellar evolution code. Focusing first on a $4.21$ $M_\odot$ model, we simulate 1,000 years of stellar evolution under the combined effects of the torque due to a single unstable prograde g mode (with an amplitude chosen on the basis of observational constraints), and diffusive angular momentum transport due to convection, overshooting, and rotational instabilities. We find that the g mode rapidly extracts angular momentum from the surface layers, depositing it deeper in the stellar interior. The angular momentum transport is so efficient that by the end of the simulation the initially non-rotating surface layers are spun in the retrograde direction to $\approx30\%$ of the critical rate. However, the additional inclusion of magnetic stresses in our simulations, almost completely inhibits this spin-up. Expanding our simulations to cover the whole instability strip, we show that the same general behavior is seen in all SPB stars. After providing some caveats to contextualize our results, we hypothesize that the observed slower surface rotation of SPB stars (as compared to other B-type stars) may be the direct consequence of the angular momentum transport that our simulations demonstrate.

Aims. The X-shooter archive of several thousand telluric star spectra was skimmed for Be and Be-shell stars to derive the stellar fundamental parameters and statistical properties, in particular for the less investigated late type Be stars, and the extension of the Be phenomenon into early A stars. Methods. An adapted version of the BCD method is used, utilizing the Balmer discontinuity parameters to determine effective temperature and surface gravity. This method is optimally suited for late B stars. The projected rotational velocity was obtained by profile fitting to the Mg ii lines of the targets, and the spectra were inspected visually for the presence of peculiar features such as the infrared Ca ii triplet or the presence of a double Balmer discontinuity. The Balmer line equivalent widths were measured, but due to uncertainties in determining the photospheric contribution are useful only in a subsample of Be stars for determining the pure emission contribution. Results. A total of 78 Be stars, mostly late type ones, were identified in the X-shooter telluric standard star archive, out of which 48 had not been reported before. The general trend of late type Be stars having more tenuous disks and being less variable than early type ones is confirmed. The relatively large number (48) of relatively bright (V > 8.5) additional Be stars casts some doubt on the statistics of late type Be stars; they are more common than currently thought: The Be/B star fraction may not strongly depend on spectral subtype.

We derive semi-analytic solutions for optically thick, super-Eddington stellar winds, induced by an assumed steady energy addition $\Delta {\dot E}$ concentrated around a near-surface heating radius $R$ in a massive star of central luminosity $L_\ast$. We show that obtaining steady wind solutions requires both that the resulting total luminosity $L_o = L_\ast + \Delta {\dot E}$ exceed the Eddington luminosity, $\Gamma_o \equiv L_o/L_{Edd} > 1$, and that the induced mass loss rate be such that the "photon-tiring" parameter $m \equiv {\dot M} GM/R L_o \le 1-1/\Gamma_o$, ensuring the luminosity is sufficient to overcome the gravitational potential $GM/R$. Our analysis unifies previous super-Eddington wind models that either: (1) assumed a direct radiative flux-driving without accounting for the advection of radiative enthalpy that can become important in such an optically thick flow; or (2) assumed that such super-Eddington outflows are adiabatic, neglecting the effects of the diffusive radiative flux. We show that these distinct models become applicable in the asymptotic limits of small vs. large values of $m \Gamma_o $, respectively. By solving the coupled differential equations for radiative diffusion and wind momentum, we obtain general solutions that effectively bridge the behaviours of these limiting models. Two key scaling results are for the terminal wind speed to escape speed, which is found to vary as $v_\infty^2/v_{esc}^2 = \Gamma_o/(1+m \Gamma_o) -1$, and for the final observed luminosity $L_{ obs}$, which for all allowed steady-solutions with $m < 1 - 1/\Gamma_o$, exceeds the Eddington luminosity, $L_{obs} > L_{Edd}$. Our super-Eddington wind solutions have potential applicability for modeling phases of eruptive mass loss from massive stars, classical novae, and the remnants of stellar mergers.

Recent spectropolarimetric surveys of bright, hot stars have found that ~10% of OB-type stars contain strong (mostly dipolar) surface magnetic fields (~kG). The prominent paradigm describing the interaction between the stellar winds and the surface magnetic field is the magnetically confined wind shock (MCWS) model. In this model, the stellar wind plasma is forced to move along the closed field loops of the magnetic field, colliding at the magnetic equator, and creating a shock. As the shocked material cools radiatively it will emit X-rays. Therefore, X-ray spectroscopy is a key tool in detecting and characterizing the hot wind material confined by the magnetic fields of these stars. Some B-type stars are found to have very short rotational periods. The effects of the rapid rotation on the X-ray production within the magnetosphere have yet to be explored in detail. The added centrifugal force due to rapid rotation is predicted to cause faster wind outflows along the field lines, leading to higher shock temperatures and harder X-rays. However, this is not observed in all rapidly rotating magnetic B-type stars. In order to address this from a theoretical point of view, we use the X-ray Analytical Dynamical Magnetosphere (XADM) model, originally developed for slow rotators, with an implementation of new rapid rotational physics. Using X-ray spectroscopy from ESA's XMM-Newton space telescope, we observed 5 rapidly rotating B-type stars to add to the previous list of observations. Comparing the observed X-ray luminosity and hardness ratio to that predicted by the XADM allows us to determine the role the added centrifugal force plays in the magnetospheric X-ray emission of these stars.

We substantially update the capabilities of the open-source software instrument Modules for Experiments in Stellar Astrophysics (MESA). MESA can now simultaneously evolve an interacting pair of differentially rotating stars undergoing transfer and loss of mass and angular momentum, greatly enhancing the prior ability to model binary evolution. New MESA capabilities in fully coupled calculation of nuclear networks with hundreds of isotopes now allow MESA to accurately simulate advanced burning stages needed to construct supernova progenitor models. Implicit hydrodynamics with shocks can now be treated with MESA, enabling modeling of the entire massive star lifecycle, from pre-main sequence evolution to the onset of core collapse and nucleosynthesis from the resulting explosion. Coupling of the GYRE non-adiabatic pulsation instrument with MESA allows for new explorations of the instability strips for massive stars while also accelerating the astrophysical use of asteroseismology data. We improve treatment of mass accretion, giving more accurate and robust near-surface profiles. A new MESA capability to calculate weak reaction rates "on-the-fly" from input nuclear data allows better simulation of accretion induced collapse of massive white dwarfs and the fate of some massive stars. We discuss the ongoing challenge of chemical diffusion in the strongly coupled plasma regime, and exhibit improvements in MESA that now allow for the simulation of radiative levitation of heavy elements in hot stars. We close by noting that the MESA software infrastructure provides bit-for-bit consistency for all results across all the supported platforms, a profound enabling capability for accelerating MESA's development.

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.

Magnetic B-type stars exhibit photometric variability due to diverse causes, and consequently on a variety of timescales. In this paper we describe interpretation of BRITE photometry and related ground-based observations of 4 magnetic B-type systems: $\epsilon$ Lupi, $\tau$ Sco, a Cen and $\epsilon$ CMa.

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.

We generalize the Rigid-Field Hydrodynamic equations to accommodate arbitrary magnetic field topologies, resulting in a new Arbitrary Rigid-Field hydrodynamic (ARFHD) formalism. We undertake a critical point calculation of the steady-state ARFHD equations with a CAK-type radiative acceleration and determine the effects of a dipole magnetic field on the usual CAK mass-loss rate and velocity structure. Enforcing the proper optically-thin limit for the radiative line-acceleration is found to decrease both the mass-loss and wind acceleration, while rotation boosts both properties. We define optically-thin-correction and rotation parameters to quantify these effects on the global mass-loss rate and develop scaling laws for the surface mass-flux as a function of surface colatitude. These scaling laws are found to agree with previous laws derived from magnetohydrodynamic simulations of magnetospheres. The dipole magnetosphere velocity structure is found to differ from a global beta-velocity law, which contradicts a central assumption of the previously-developed XADM model of X-ray emission from magnetospheres.

Slowly rotating magnetic massive stars develop "dynamical magnetospheres" (DM's), characterized by trapping of stellar wind outflow in closed magnetic loops, shock heating from collision of the upflow from opposite loop footpoints, and subsequent gravitational infall of radiatively cooled material. In 2D and 3D magnetohydrodynamic (MHD) simulations the interplay among these three components is spatially complex and temporally variable, making it difficult to derive observational signatures and discern their overall scaling trends.Within a simplified, steady-state analysis based on overall conservation principles, we present here an "analytic dynamical magnetosphere" (ADM) model that provides explicit formulae for density, temperature and flow speed in each of these three components -- wind outflow, hot post-shock gas, and cooled inflow -- as a function of colatitude and radius within the closed (presumed dipole) field lines of the magnetosphere. We compare these scalings with time-averaged results from MHD simulations, and provide initial examples of application of this ADM model for deriving two key observational diagnostics, namely hydrogen H-alpha emission line profiles from the cooled infall, and X-ray emission from the hot post-shock gas. We conclude with a discussion of key issues and advantages in applying this ADM formalism toward derivation of a broader set of observational diagnostics and scaling trends for massive stars with such dynamical magnetospheres.

KIC 7760680 is so far the richest slowly pulsating B star, by exhibiting 36 consecutive dipole ($\ell=1$) gravity (g-) modes. The monotonically decreasing period spacing of the series, in addition to the local dips in the pattern confirm that KIC 7760680 is a moderate rotator, with clear mode trapping in chemically inhomogeneous layers. We employ the traditional approximation of rotation to incorporate rotational effects on g-mode frequencies. Our detailed forward asteroseismic modelling of this g-mode series reveals that KIC 7760680 is a moderately rotating B star with mass $\sim3.25$ M$_\odot$. By simultaneously matching the slope of the period spacing, and the number of modes in the observed frequency range, we deduce that the equatorial rotation frequency of KIC 7760680 is 0.4805 day$^{-1}$, which is 26\% of its Roche break up frequency. The relative deviation of the model frequencies and those observed is less than one percent. We succeed to tightly constrain the exponentially-decaying convective core overshooting parameter to $f_{\rm ov}\approx0.024\pm0.001$. This means that convective core overshooting can coexist with moderate rotation. Moreover, models with exponentially-decaying overshoot from the core outperform those with the classical step-function overshoot. The best value for extra diffusive mixing in the radiatively stable envelope is confined to $\log D_{\rm ext}\approx0.75\pm0.25$ (with $D_{\rm ext}$ in cm$^2$ sec$^{-1}$), which is notably smaller than theoretical predictions.

Spectroscopic observations of Achernar over the past decades, have shown the photospheric line width, as measured by the rotational parameter $v \sin i$, to vary in correlation with the emission activity. Here we present new observations, covering the most recent activity phase, and further archival data collected from the archives. The $v \sin i$ variation is confirmed. On the basis of the available data it cannot be decided with certainty whether the increased line width precedes the emission activity, i.e. is a signature of the ejection mechanism, or postdates is, which would make it a signature of re-accretion of some of the disk-material. However, the observed evidence leans towards the re-accretion hypothesis. Two further stars showing the effect of variable line width in correlation with emission activity, namely 66 Oph and $\pi$ Aqr, are presented as well.

The MiMeS project is a large-scale, high resolution, sensitive spectropolarimetric investigation of the magnetic properties of O and early B type stars. Initiated in 2008 and completed in 2013, the project was supported by 3 Large Program allocations, as well as various programs initiated by independent PIs and archival resources. Ultimately, over 4800 circularly polarized spectra of 560 O and B stars were collected with the instruments ESPaDOnS at the Canada-France-Hawaii Telescope, Narval at the T\'elescope Bernard Lyot, and HARPSpol at the European Southern Observatory La Silla 3.6m telescope, making MiMeS by far the largest systematic investigation of massive star magnetism ever undertaken. In this paper, the first in a series reporting the general results of the survey, we introduce the scientific motivation and goals, describe the sample of targets, review the instrumentation and observational techniques used, explain the exposure time calculation designed to provide sensitivity to surface dipole fields above approximately 100 G, discuss the polarimetric performance, stability and uncertainty of the instrumentation, and summarize the previous and forthcoming publications.

This work explores a Pulsationally Driven Orbital Mass Ejection (PDOME) model for the launching of Classical Be star disks. Under this model, a combination of rapid rotation and non-radial pulsation modes contribute to placing material into the circumstellar environment. Several varieties of non-radial pulsation modes, characterized by their propagation direction and the relative phase of their velocity and density perturbations, are considered. As well, the orbital stability of material launched by such a mechanism is investigated.

The initial success of the Rigidly Rotating Magnetosphere (RRM) model application to the B2Vp star sigma OriE by Townsend, Owocki & Groote (2005) triggered a renewed era of observational monitoring of this archetypal object. We utilize high-resolution spectropolarimetry and the magnetic Doppler imaging (MDI) technique to simultaneously determine the magnetic configuration, which is predominately dipolar, with a polar strength Bd = 7.3-7.8 kG and a smaller non-axisymmetric quadrupolar contribution, as well as the surface distribution of abundance of He, Fe, C, and Si. We describe a revised RRM model that now accepts an arbitrary surface magnetic field configuration, with the field topology from the MDI models used as input. The resulting synthetic Ha emission and broadband photometric observations generally agree with observations, however, several features are poorly fit. To explore the possibility of a photospheric contribution to the observed photometric variability, the MDI abundance maps were used to compute a synthetic photospheric light curve to determine the effect of the surface inhomogeneities. Including the computed photospheric brightness modulation fails to improve the agreement between the observed and computed photometry. We conclude that the discrepancies cannot be explained as an effect of inhomogeneous surface abundance. Analysis of the UV light variability shows good agreement between observed variability and computed light curves, supporting the accuracy of the photospheric light variation calculation. We thus conclude that significant additional physics is necessary for the RRM model to acceptably reproduce observations of not only sigma Ori E, but also other similar stars with significant stellar wind-magnetic field interactions.

In this paper we report 23 magnetic field measurements of the B3IV star HD 23478: 12 obtained from high resolution Stokes $V$ spectra using the ESPaDOnS (CFHT) and Narval (TBL) spectropolarimeters, and 11 from medium resolution Stokes $V$ spectra obtained with the DimaPol spectropolarimeter (DAO). HD 23478 was one of two rapidly rotating stars identified as potential "centrifugal magnetosphere" hosts based on IR observations from the Apache Point Observatory Galactic Evolution Experiment survey. We derive basic physical properties of this star including its mass ($M=6.1^{+0.8}_{-0.7}\,M_\odot$), effective temperature ($T_{\rm eff}=20\pm2\,$kK), radius ($R=2.7^{+1.6}_{-0.9}\,R_\odot$), and age ($\tau_{\rm age}=3^{+37}_{-1}\,$Myr). We repeatedly detect weakly-variable Zeeman signatures in metal, He and H lines in all our observations corresponding to a longitudinal magnetic field of $\langle B_z\rangle\approx-2.0\,$kG. The rotational period is inferred from Hipparcos photometry ($P_{\rm rot}=1.0498(4)\,$d). Under the assumption of the Oblique Rotator Model, our obsevations yield a surface dipole magnetic field of strength $B_d\geq9.5\,$kG that is approximately aligned with the stellar rotation axis. We confirm the presence of strong and broad H$\alpha$ emission and gauge the volume of this star's centrifugal magnetosphere to be consistent with those of other H$\alpha$ emitting centrifugal magnetosphere stars based on the large inferred Alfv\'en to Kepler radius ratio.

We analyze a high resolution spectropolarimetric dataset collected for the He-weak B3p IV star HR 2949. The Zeeman effect is visible in the circularly polarized component of numerous spectral lines. The longitudinal magnetic field varies between approximately $-650$ and $+150$ G. The polar strength of the surface magnetic dipole is calculated to be 2.4$^{+0.3}_{-0.2}$ kG. The star has strong overabundances of Fe-peak elements, along with extremely strong overabundances of rare-earth elements; however, He, Al, and S are underabundant. This implies that HR 2949 is a chemically peculiar star. Variability is seen in all photospheric lines, likely due to abundance patches as seen in many Ap/Bp stars. Longitudinal magnetic field variations measured from different spectral lines yield different results, likely a consequence of uneven sampling of the photospheric magnetic field by the abundance patches. Analysis of photometric and spectroscopic data for both HR 2949 and its companion star, HR 2948, suggests a revision of HR 2949's fundamental parameters: in particular, it is somewhat larger, hotter, and more luminous than previously believed. There is no evidence of optical or ultraviolet emission originating in HR 2949's magnetosphere, despite its moderately strong magnetic field and relatively rapid rotation; however, when calculated using theoretical and empirical boundaries on the initial rotational velocity, the spindown age is compatible with the stellar age. With the extensive phase coverage presented here, HR 2949 will make an excellent subject for Zeeman Doppler Imaging.

We report magnetic and spectroscopic observations and modeling of the Of?p star CPD -28 2561. Using more than 75 new spectra, we have measured the equivalent width variations and examined the dynamic spectra of photospheric and wind-sensitive spectral lines. A period search results in an unambiguous 73.41 d variability period. High resolution spectropolarimetric data analyzed using Least-Squares Deconvolution yield a Zeeman signature detected in the mean Stokes V profile corresponding to phase 0.5 of the spectral ephemeris. Interpreting the 73.41 d period as the stellar rotational period, we have phased the equivalent widths and inferred longitudinal field measurements. The phased magnetic data exhibit a weak sinusoidal variation, with maximum of about 565 G at phase 0.5, and a minimum of about -335 G at phase 0.0, with extrema approximately in phase with the (double-wave) Halpha equivalent width variation. Modeling of the Halpha equivalent width variation assuming a quasi-3D magnetospheric model produces a unique solution for the ambiguous couplet of inclination and magnetic obliquity angles: (i, beta) or (beta, i)=(35 deg,90 deg). Adopting either geometry, the longitudinal field variation yields a dipole polar intensity Bd=2.6\pm 0.9~kG, consistent with that obtained from direct modelling of the Stokes V profiles. We derive a wind magnetic confinement parameter eta*\simeq 100, leading to an Alfv\'en radius RA\simeq 3-5~R*, and a Kepler radius RK\simeq 20~R*. This supports a physical scenario in which the Halpha emission and other line variability have their origin in an oblique, co-rotating 'dynamical magnetosphere' structure resulting from a magnetically channeled wind. Nevertheless, the details of the formation of spectral lines and their variability within this framework remain generally poorly understood.

With recent advances in asteroseismology it is now possible to peer into the cores of red giants, potentially providing a way to study processes such as nuclear burning and mixing through their imprint as sharp structural variations -- glitches -- in the stellar cores. Here we show how such core glitches can affect the oscillations we observe in red giants. We derive an analytical expression describing the expected frequency pattern in the presence of a glitch. This formulation also accounts for the coupling between acoustic and gravity waves. From an extensive set of canonical stellar models we find glitch-induced variation in the period spacing and inertia of non-radial modes during several phases of red-giant evolution. Significant changes are seen in the appearance of mode amplitude and frequency patterns in asteroseismic diagrams such as the power spectrum and the \'echelle diagram. Interestingly, along the red-giant branch glitch-induced variation occurs only at the luminosity bump, potentially providing a direct seismic indicator of stars in that particular evolution stage. Similarly, we find the variation at only certain post-helium-ignition evolution stages, namely, in the early phases of helium-core burning and at the beginning of helium-shell burning signifying the asymptotic-giant-branch bump. Based on our results, we note that assuming stars to be glitch-free, while they are not, can result in an incorrect estimate of the period spacing. We further note that including diffusion and mixing beyond classical Schwarzschild, could affect the characteristics of the glitches, potentially providing a way to study these physical processes.

Binary systems with an accreting compact object are a unique chance to investigate the strong, clumpy, line-driven winds of early type supergiants by using the compact object's X-rays to probe the wind structure. We analyze the two-component wind of HDE 226868, the O9.7Iab giant companion of the black hole Cyg X-1 using 4.77 Ms of RXTE observations of the system taken over the course of 16 years. Absorption changes strongly over the 5.6 d binary orbit, but also shows a large scatter at a given orbital phase, especially at superior conjunction. The orbital variability is most prominent when the black hole is in the hard X-ray state. Our data are poorer for the intermediate and soft state, but show signs for orbital variability of the absorption column in the intermediate state. We quantitatively compare the data in the hard state to a toy model of a focussed Castor-Abbott-Klein-wind: as it does not incorporate clumping, the model does not describe the observations well. A qualitative comparison to a simplified simulation of clumpy winds with spherical clumps shows good agreement in the distribution of the equivalent hydrogen column density for models with a porosity length on the order of the stellar radius at inferior conjunction; we conjecture that the deviations between data and model at superior conjunction could be either due to lack of a focussed wind component in the model or a more complicated clump structure.

Approximately 10% of B-type stars possess strong magnetic fields, and of these, 25% host centrifugal magnetospheres (CMs) in which the radiative wind, magnetic field, and rotational support interact to form a dense circumstellar plasma visible in a variety of diagnostic lines. In this article we review the basic theory behind CMs, outline current theoretical and observational problems, compare the observational properties of CM host stars to those of classical Be stars, and finally present preliminary results of a population study aimed at clarifying the characteristics of this growing sub-class.

We use 2D MHD simulations to examine the effects of radiative cooling and inverse Compton (IC) cooling on X-ray emission from magnetically confined wind shocks (MCWS) in magnetic massive stars with radiatively driven stellar winds. For the standard dependence of mass loss rate on luminosity $\Mdot \sim L^{1.7} $, the scaling of IC cooling with $L$ and radiative cooling with $\Mdot$ means that IC cooling become formally more important for lower luminosity stars. However, because the sense of the trends is similar, we find the overall effect of including IC cooling is quite modest. More significantly, for stars with high enough mass loss to keep the shocks radiative, the MHD simulations indicate a linear scaling of X-ray luminosity with mass loss rate; but for lower luminosity stars with weak winds, X-ray emission is reduced and softened by a {\em shock retreat} resulting from the larger post-shock cooling length, which within the fixed length of a closed magnetic loop forces the shock back to lower pre-shock wind speeds. A semi-analytic scaling analysis that accounts both for the wind magnetic confinement and this shock retreat yields X-ray luminosities that have a similar scaling trend, but a factor few higher values, compared to time-averages computed from the MHD simulations. The simulation and scaling results here thus provide a good basis for interpreting available X-ray observations from the growing list of massive stars with confirmed large-scale magnetic fields.

The main theoretical problem for the formation of a Keplerian disk around Be stars is how to supply angular momentum from the star to the disk, even more so since Be stars probably rotate somewhat sub-critically. For instance, nonradial pulsation may transport angular momentum to the stellar surface until (part of) this excess supports the disk formation/replenishment. The nearby Be star Achernar is presently building a new disk and offers an excellent opportunity to observe this process from relatively close-up. Spectra from various sources and epochs are scrutinized to identify the salient stellar parameters characterizing the disk life cycle as defined by H\alpha emission. Variable strength of the non-radial pulsation is confirmed, but does not affect the further results. For the first time it is demonstrated that the photospheric line width does vary in a Be star, by as much as \Delta v sin i \lesssim 35kms^{-1}. However, contrary to assumptions in which a photospheric spin-up accumulates during the diskless phase and then is released into the disk as it is fed, the apparent photospheric spin-up is positively correlated with the appearance of H\alpha line emission: The photospheric line widths and circumstellar emission increase together, and the apparent stellar rotation declines to the value at quiescence after the H\alpha line emission becomes undetectable.

It has long been known that rotation can have an appreciable impact on stellar pulsation --- by modifying the usual p and g modes found in the non-rotating case, and by introducing new classes of mode. However, it's only relatively recently that advances in numerical simulations and complementary theoretical treatments have enabled us to model these phenomena in any great detail. In this talk I'll review highlights in this area (the 'Greatest Hits'), before considering the flip side (or the 'B-side', for those of us old enough to remember vinyl records) of the pulsation-rotation interaction: how pulsation can itself influence internal rotation profiles.

We introduce GYRE, a new open-source stellar oscillation code which solves the adiabatic/non-adiabatic pulsation equations using a novel Magnus Multiple Shooting (MMS) numerical scheme. The code has a global error scaling of up to 6th order in the grid spacing, and can therefore achieve high accuracy with few grid points. It is moreover robust and efficiently makes use of multiple processor cores and/or nodes. We present an example calculation using GYRE, and discuss recent work to integrate GYRE into the asteroseismic optimization module of the MESA stellar evolution code.

We present a new oscillation code, GYRE, which solves the stellar pulsation equations (both adiabatic and non-adiabatic) using a novel Magnus Multiple Shooting numerical scheme devised to overcome certain weaknesses of the usual relaxation and shooting schemes appearing in the literature. The code is accurate (up to 6th order in the number of grid points), robust, efficiently makes use of multiple processor cores and/or nodes, and is freely available in source form for use and distribution. We verify the code against analytic solutions and results from other oscillation codes, in all cases finding good agreement. Then, we use the code to explore how the asteroseismic observables of a 1.5 M_sun star change as it evolves through the red-giant bump.

We substantially update the capabilities of the open source software package Modules for Experiments in Stellar Astrophysics (MESA), and its one-dimensional stellar evolution module, MESA Star. Improvements in MESA Star's ability to model the evolution of giant planets now extends its applicability down to masses as low as one-tenth that of Jupiter. The dramatic improvement in asteroseismology enabled by the space-based Kepler and CoRoT missions motivates our full coupling of the ADIPLS adiabatic pulsation code with MESA Star. This also motivates a numerical recasting of the Ledoux criterion that is more easily implemented when many nuclei are present at non-negligible abundances. This impacts the way in which MESA Star calculates semi-convective and thermohaline mixing. We exhibit the evolution of 3-8 Msun stars through the end of core He burning, the onset of He thermal pulses, and arrival on the white dwarf cooling sequence. We implement diffusion of angular momentum and chemical abundances that enable calculations of rotating-star models, which we compare thoroughly with earlier work. We introduce a new treatment of radiation-dominated envelopes that allows the uninterrupted evolution of massive stars to core collapse. This enables the generation of new sets of supernovae, long gamma-ray burst, and pair-instability progenitor models. We substantially modify the way in which MESA Star solves the fully coupled stellar structure and composition equations, and we show how this has improved MESA's performance scaling on multi-core processors. Updates to the modules for equation of state, opacity, nuclear reaction rates, and atmospheric boundary conditions are also provided. We describe the MESA Software Development Kit (SDK) that packages all the required components needed to form a unified and maintained build environment for MESA. [Abridged]

We present results from three weeks' photometric monitoring of the magnetic helium-strong star sigma Ori E using the MOST microsatellite. The star's light curve is dominated by twice-per-rotation eclipse-like dimmings arising when magnetospheric clouds transit across and occult the stellar disk. However, no evidence is found for any abrupt centrifugal breakout of plasma from the magnetosphere, either in the residual flux or in the depths of the light minima. Motivated by this finding we compare the observationally inferred magnetospheric mass against that predicted by a breakout analysis. The large discrepancy between the values leads us to argue that centrifugal breakout does not play a significant role in establishing the magnetospheric mass budget of sigma Ori E.

Some massive stars possess strong magnetic fields that confine plasma in the circumstellar environment. These \textit{magnetospheres} have been studied spectroscopically, photometrically and, more recently, interferometrically. Here we report on the first firm detection of a magnetosphere in continuum linear polarization, as a result of monitoring of $\sigma$\,Ori\,E at the Pico dos Dias Observatory. {The non-zero intrinsic polarization indicates an asymmetric structure, whose minor elongation axis is oriented $150\fdg0$ east of the celestial north.} A modulation of the polarization was observed, with a period of half of the rotation period, which supports the theoretical prediction of the presence of two diametrally opposed, co-rotating blobs of gas. A phase lag of -0.085 was detected between the polarization minimum and the primary minimum of the light curve, suggestive of a complex shape of the plasma clouds. We present a preliminary analysis of the data with the Rigidly Rotating Magnetosphere model, which could not reproduce simultaneously the photometric and polarimetric data. A toy model comprising two spherical co-rotating blobs {joined by a thin disk} proved more successful in reproducing the polarization modulation. {With this model we were able to determine that the total scattering mass of the thin disk is similar to the mass of the blobs ($2M_{\rm b}/M_{\rm d}=1.2$) and that the blobs are rotating counterclockwise on the plane of the sky.} This result shows that polarimetry can provide a diagnostic of the geometry of clouds, which will serve as an important constraint for {improving} the Rigidly Rotating Magnetosphere model.

Building on results from the Magnetism in Massive Stars (MiMeS) project, this paper shows how a two-parameter classification of massive-star magnetospheres in terms of the magnetic wind confinement (which sets the Alfv\'en radius RA) and stellar rotation (which sets the Kepler co-rotation radius RK) provides a useful organisation of both observational signatures and theoretical predictions. We compile the first comprehensive study of inferred and observed values for relevant stellar and magnetic parameters of 64 confirmed magnetic OB stars with Teff > 16 kK. Using these parameters, we locate the stars in the magnetic confinement-rotation diagram, a log-log plot of RK vs. RA. This diagram can be subdivided into regimes of centrifugal magnetospheres (CM), with RA > RK, vs. dynamical magnetospheres (DM), with RK > RA. We show how key observational diagnostics, like the presence and characteristics of Halpha emission, depend on a star's position within the diagram, as well as other parameters, especially the expected wind mass-loss rates. In particular, we identify two distinct populations of magnetic stars with Halpha emission: namely, slowly rotating O-type stars with narrow emission consistent with a DM, and more rapidly rotating B-type stars with broader emission associated with a CM. For O-type stars, the high mass-loss rates are sufficient to accumulate enough material for line emission even within the relatively short free-fall timescale associated with a DM: this high mass-loss rate also leads to a rapid magnetic spindown of the stellar rotation. For the B-type stars, the longer confinement of a CM is required to accumulate sufficient emitting material from their relatively weak winds, which also lead to much longer spindown timescales. [abbreviated]

The spectral and magnetic properties and variability of the B2Vnp emission-line magnetosphere star HR7355 were analyzed. The object rotates at almost 90% of the critical value, meaning it is a magnetic star for which oblateness and gravity darkening effects cannot be ignored any longer. A detailed modeling of the photospheric parameters indicate that the star is significantly cooler than suggested by the B2 spectral type, with T_eff=17500K atypically cool for a star with a helium enriched surface. The spectroscopic variability of helium and metal lines due to the photospheric abundance pattern is far more complex than a largely dipolar, oblique magnetic field of about 11 to 12kG may suggest. Doppler imaging shows that globally the most He enriched areas coincide with the magnetic poles and metal enriched areas with the magnetic equator. While most of the stellar surface is helium enriched with respect to the solar value, some isolated patches are depleted. The stellar wind in the circumstellar environment is governed by the magnetic field, i.e. the stellar magnetosphere is rigidly corotating with the star. The magnetosphere of HR7355 is similar to the well known \sigma Ori E: the gas trapped in the magnetospheric clouds is fairly dense, and at the limit to being optically thick in the hydrogen emission. Apart from a different magnetic obliquity, HR7355 and the more recently identified HR5907 have virtually identical stellar and magnetic parameters.

We present the first fully 3D MHD simulation for magnetic channeling and confinement of a radiatively driven, massive-star wind. The specific parameters are chosen to represent the prototypical slowly rotating magnetic O star \theta^1 Ori C, for which centrifugal and other dynamical effects of rotation are negligible. The computed global structure in latitude and radius resembles that found in previous 2D simulations, with unimpeded outflow along open field lines near the magnetic poles, and a complex equatorial belt of inner wind trapping by closed loops near the stellar surface, giving way to outflow above the Alfv\'{e}n radius. In contrast to this previous 2D work, the 3D simulation described here now also shows how this complex structure fragments in azimuth, forming distinct clumps of closed loop infall within the Alfv\'{e}n radius, transitioning in the outer wind to radial spokes of enhanced density with characteristic azimuthal separation of $15-20 \degr$. Applying these results in a 3D code for line radiative transfer, we show that emission from the associated 3D `dynamical magnetosphere' matches well the observed H\alpha emission seen from \theta^1 Ori C, fitting both its dynamic spectrum over rotational phase, as well as the observed level of cycle to cycle stochastic variation. Comparison with previously developed 2D models for Balmer emission from a dynamical magnetosphere generally confirms that time-averaging over 2D snapshots can be a good proxy for the spatial averaging over 3D azimuthal wind structure. Nevertheless, fully 3D simulations will still be needed to model the emission from magnetospheres with non-dipole field components, such as suggested by asymmetric features seen in the H\alpha equivalent-width curve of \theta^1 Ori C.

The young O9.5 V spectroscopic binary Theta 2 Ori A shows moderately hard X-ray emission and relatively narrow X-ray lines, suggesting that it may be a Magnetically Confined Wind Shock (MCWS) source, similar to its more massive analogue Theta 1 Ori C. X-ray flares occurring near periastron led to the suggestion that the flares are produced via magnetic reconnection as magnetospheres on both components of the Theta 2 Ori A binary interact at closest approach. We use a series of high-resolution spectropolarimetric observations of Theta 2 Ori A to place an upper limit on the magnetic field strength of 135 G (95% credible region). Such a weak dipole field would not produce magnetic confinement, or a large magnetosphere. A sub-pixel analysis of the Chandra ACIS images of Theta 2 Ori A obtained during quiescence and flaring show that the hard, flaring X-rays are offset from the soft, quiescent emission by 0.4 arcsec. If the soft emission is associated with the A1/A2 spectroscopic binary, the offset and position angle of the hard, flaring source place it at the location of the intermediate-mass A3 companion, discovered via speckle interferometry. The spectropolarimetric and X-ray results taken together point to the A3 companion, not the massive A1/A2 binary, as the source of hard, flaring X-ray emission. We also discuss a similar analysis performed for the magnetic Bp star Sigma Ori E. We find a similar origin for its X-ray flaring.

The magnetic O-star HD191612 exhibits strongly variable, cyclic Balmer line emission on a 538-day period. We show here that its variable Halpha emission can be well reproduced by the rotational phase variation of synthetic spectra computed directly from full radiation magneto-hydrodynamical simulations of a magnetically confined wind. In slow rotators such as HD191612, wind material on closed magnetic field loops falls back to the star, but the transient suspension of material within the loops leads to a statistically overdense, low velocity region around the magnetic equator, causing the spectral variations. We contrast such "dynamical magnetospheres" (DMs) with the more steady-state "centrifugal magnetospheres" of stars with rapid rotation, and discuss the prospects of using this DM paradigm to explain periodic line emission from also other non-rapidly rotating magnetic massive stars.

In the framework of the Magnetism in Massive Stars (MiMeS) project, a HARPSpol Large Program at the 3.6m-ESO telescope has recently started to collect high-resolution spectropolarimetric data of a large number of Southern massive OB stars in the field of the Galaxy and in many young clusters and associations. In this Letter, we report on the first discoveries of magnetic fields in two massive stars with HARPSpol - HD 130807 and HD 122451, and confirm the presence of a magnetic field at the surface of HD 105382 that was previously observed with a low spectral resolution device. The longitudinal magnetic field measurements are strongly varying for HD 130807 from $\sim$-100 G to $\sim$700 G. Those of HD 122451 and HD 105382 are less variable with values ranging from $\sim$-40 to -80 G, and from $\sim$-300 to -600 G, respectively. The discovery and confirmation of three new magnetic massive stars, including at least two He-weak stars, is an important contribution to one of the MiMeS objectives: the understanding of origin of magnetic fields in massive stars and their impacts on stellar structure and evolution.

We present a generalised formalism for treating the porosity-associated reduction in continuum opacity that occurs when individual clumps in a stochastic medium become optically thick. We consider geometries resulting in either isotropic or anisotropic effective opacity, and, in addition to an idealised model in which all clumps have the same local overdensity and scale, we also treat an ensemble of clumps with optical depths set by Markovian statistics. This formalism is then applied to the specific case of bound-free absorption of X-rays in hot star winds, a process not directly affected by clumping in the optically thin limit. We find that the Markov model gives surprisingly similar results to those found previously for the single clump model, suggesting that porous opacity is not very sensitive to details of the assumed clump distribution function. Further, an anisotropic effective opacity favours escape of X-rays emitted in the tangential direction (the `venetian blind' effect), resulting in a 'bump' of higher flux close to line centre as compared to profiles computed from isotropic porosity models. We demonstrate how this characteristic line shape may be used to diagnose the clump geometry, and we confirm previous results that for optically thick clumping to significantly influence X-ray line profiles, very large porosity lengths, defined as the mean free path between clumps, are required. Moreover, we present the first X-ray line profiles computed directly from line-driven instability simulations using a 3-D patch method, and find that porosity effects from such models also are very small. This further supports the view that porosity has, at most, a marginal effect on X-ray line diagnostics in O stars, and therefore that these diagnostics do indeed provide a good `clumping insensitive' method for deriving O star mass-loss rates.

We report magnetic and spectroscopic observations and modeling of the Of?p star HD 148937 within the context of the MiMeS LP at the CFHT. Thirty-two high signal-to-noise ratio circularly polarised (Stokes V) spectra and 13 unpolarised (Stokes I) spectra of HD 148937 were acquired in 2009 and 2010. A definite detection of a Stokes V Zeeman signature is obtained in the grand mean of all observations (in both LSD mean profiles and individual spectral lines). The longitudinal magnetic field inferred from the Stokes V LSD profiles is consistently negative, in contrast to the essentially zero field strength measured from the diagnostic null profiles. A period search of equivalent width measurements confirms the previously-reported 7.03 d variability period. The variation of equivalent widths is not strictly periodic: we present evidence for evolution of the amount or distribution of circumstellar plasma. Interpreting the 7.03 d period as the stellar rotational period within the context of the ORM, we have phased the equivalent widths and longitudinal field measurements. The longitudinal field measurements show a weak sinusoidal variation of constant sign, with extrema out of phase with the H{\alpha} variation by about 0.25 cycles. The inferred magnetic configuration confirms the suggestion of Naz\'e et al (2010), who proposed that the weaker variability of HD 148937 as compared to other members of this class is a consequence of the stellar geometry. Based on the derived magnetic properties and published wind characteristics, we find a wind magnetic confinement parameter \eta\ast \simeq 20 and rotation parameter W = 0.12, supporting a picture in which the Halpha emission and other line variability have their origin in an oblique, rigidly rotating magnetospheric structure resulting from a magnetically channeled wind. (Abridged.)

We report the discovery and analysis of a very strong magnetic field in the rapidly rotating early B-type star HR 5907, based on observations obtained as part of the Magnetism in Massive Stars (MiMeS) project. We infer a rotation period of 0.508276 +0.000015/-0.000012 d from photometric and H{\alpha} EW measurements, making this the shortest period, non-degenerate, magnetic massive star known to date. From the comparison of IUE UV and optical spectroscopy with LTE BRUCE/KYLIE models we find a solid-angle integrated, uniform black-body temperature of 17 000 \pm 1000 K, a projected rotational velocity of 290 \pm 10 km/s, an equatorial radius of 3.1 \pm 0.2 R_sun, a stellar mass of 5.5 \pm 0.5 M_sun, and an inclination angle of the rotation axis to our line-of-sight of 70 \pm 10\circ. Our measurements of the longitudinal magnetic field, which vary between -500 and -2000 G, phase coherently with the rotation period and imply a surface dipole field strength of \sim15.7 kG. On the other hand, from fits to mean Least-Squares Deconvolved Stokes V line profiles we infer a dipole field strength of \sim10.4 kG. This disagreement may result from a magnetic configuration more complex than our model, and/or from the non-uniform helium surface abundance distribution. In either case we obtain a magnetic obliquity nearly aligned with the rotation axis ({\beta} = 7+2/-1\circ). Our optical spectroscopy also shows weak variability in carbon, silicon and nitrogen lines. The emission variability in hydrogen Balmer and Paschen lines indicates the presence of a dense, highly structured magnetosphere, interpreted as a centrifugally supported, magnetically confined circumstellar disk.

We have obtained 18 new high-resolution spectropolarimetric observations of the B2Vp star sigma Ori E with both the Narval and ESPaDOnS spectropolarimeters. The aim of these observations is to test, with modern data, the assumptions of the Rigidly Rotating Magnetosphere (RRM) model of Townsend & Owocki (2005), applied to the specific case of sigma Ori E by Townsend et al. (2005). This model includes a substantially offset dipole magnetic field configuration, and approximately reproduces previous observational variations in longitudinal field strength, photometric brightness, and Halpha emission. We analyze new spectroscopy, including H I, He I, C II, Si III and Fe III lines, confirming the diversity of variability in photospheric lines, as well as the double S-wave variation of circumstellar hydrogen. Using the multiline analysis method of Least-Squares Deconvolution (LSD), new, more precise longitudinal magnetic field measurements reveal a substantial variance between the shapes of the observed and RRM model time-varying field. The phase resolved Stokes V profiles of He I 5876 A and 6678 A lines are fit poorly by synthetic profiles computed from the magnetic topology assumed by Townsend et al. (2005). These results challenge the offset dipole field configuration assumed in the application of the RRM model to sigma Ori E, and indicate that future models of its magnetic field should also include complex, higher-order components.

This paper reports high-precision Stokes V spectra of HD 191612 acquired using the ESPaDOnS spectropolarimeter at the Canada-France-Hawaii Telescope, in the context of the Magnetism in Massive stars (MiMeS) Project. Using measurements of the equivalent width of the Halpha line and radial velocities of various metallic lines, we have updated both the spectroscopic and orbital ephemerides of this star. We confirm the presence of a strong magnetic field in the photosphere of HD 191612, and detect its variability. We establish that the longitudinal field varies in a manner consistent with the spectroscopic period of 537.6 d, in an approximately sinusoidal fashion. This demonstrates a firm connection between the magnetic field and the processes responsible for the line and continuum variability. Interpreting the variation of the longitudinal magnetic field within the context of the dipole oblique rotator model we obtain a best-fit surface magnetic field model with obliquity beta=67\pm 5 deg and polar strength Bd=2450\pm 400 G . The inferred magnetic field strength implies an equatorial wind magnetic confinement parameter eta*~50, supporting a picture in which the Halpha emission and photometric variability have their origin in an oblique, rigidly rotating magnetospheric structure resulting from a magnetically channeled wind. This interpretation is supported by our successful Monte Carlo radiative transfer modeling of the photometric variation, which assumes the enhanced plasma densities in the magnetic equatorial plane above the star implied by such a picture. Predictions of the continuum linear polarisation resulting from Thompson scattering from the magnetospheric material indicate that the Stokes Q and U variations are highly sensitive to the magnetospheric geometry, and that expected amplitudes are in the range of current instrumentation. (abridged)

We report the results of a new survey of massive, OB stars throughout the Carina Nebula using the X-ray point source catalog provided by the Chandra Carina Complex Project (CCCP) in conjunction with infrared (IR) photometry from the Two Micron All-Sky Survey and the Spitzer Space Telescope Vela--Carina survey. Mid-IR photometry is relatively unaffected by extinction, hence it provides strong constraints on the luminosities of OB stars, assuming that their association with the Carina Nebula, and hence their distance, is confirmed. We fit model stellar atmospheres to the optical (UBV) and IR spectral energy distributions (SEDs) of 182 OB stars with known spectral types and measure the bolometric luminosity and extinction for each star. We find that the extinction law measured toward the OB stars has two components: Av=1--1.5 mag produced by foreground dust with a ratio of total-to-selective absorption Rv=3.1 plus a contribution from local dust with Rv>4.0 in the Carina molecular clouds that increases as Av increases. Using X-ray emission as a strong indicator of association with Carina, we identify 94 candidate OB stars with Lbol\geq10^4 Lsun by fitting their IR SEDs. If the candidate OB stars are eventually confirmed by follow-up spectroscopic observations, the number of cataloged OB stars in the Carina Nebula will increase by ~50%. Correcting for incompleteness due to OB stars falling below the Lbol cutoff or the CCCP detection limit, these results potentially double the size of the young massive stellar population.

We discuss the recent detection of a strong, organized magnetic field in the bright, broad-line B2V star, HR 5907, using the ESPaDOnS spectropolarimeter on the CFHT as part of the Magnetism in Massive Stars (MiMeS) survey. We find a rotational period of 0.50833 days, making it the fastest-rotating, non-degenerate magnetic star ever detected. Like the previous rapid-rotation record holder HR 7355 (also discovered by MiMeS: Oksala et al. 2010, Rivinius et al. 2010), this star shows emission line variability that is diagnostic of a structured magnetosphere.

Despite a century of remarkable progress in understanding stellar interiors, we know surprisingly little about the inner workings of stars spinning near their critical limit. New interferometric imaging of these so-called ``rapid rotators'' combined with breakthroughs in asteroseismology promise to lift this veil and probe the strongly latitude-dependent photospheric characteristics and even reveal the internal angular momentum distribution of these luminous objects. Here, we report the first high precision photometry on the low-amplitude delta cuti variable star Rasalhague (alpha Oph, A5IV, 2.18 Msun, omega/omega_c~0.88) based on 30 continuous days of monitoring using the MOST satellite. We have identified 57+/-1 distinct pulsation modes above a stochastic granulation spectrum with a cutoff of ~26 cycles per day. Remarkably, we have also discovered that the fast rotation period of 14.5~hours modulates low-frequency modes (1-10 day periods) that we identify as a rich family of g-modes (|m| up to 7). The spacing of the g-modes is surprisingly linear considering Coriolis forces are expected to strongly distort the mode spectrum, suggesting we are seeing prograde ``equatorial Kelvin'' waves (modes l=m). We emphasize the unique aspects of Rasalhague motivating future detailed asteroseismic modeling -- a source with a precisely measured parallax distance, photospheric oblateness, latitude temperature structure, and whose low-mass companion provides an astrometric orbit for precise mass determinations.

I introduce a new code for fast calculation of the Lomb-Scargle periodogram, that leverages the computing power of graphics processing units (GPUs). After establishing a background to the newly emergent field of GPU computing, I discuss the code design and narrate key parts of its source. Benchmarking calculations indicate no significant differences in accuracy compared to an equivalent CPU-based code. However, the differences in performance are pronounced; running on a low-end GPU, the code can match 8 CPU cores, and on a high-end GPU it is faster by a factor approaching thirty. Applications of the code include analysis of long photometric time series obtained by ongoing satellite missions and upcoming ground-based monitoring facilities; and Monte-Carlo simulation of periodogram statistical properties.

This paper presents results obtained from Stokes I and V spectra of the B2Vp star sigma Ori E, observed by both the Narval and ESPaDOnS spectropolarimeters. Using Least- Squares Deconvolution, we investigate the longitudinal magnetic field at the current epoch, including period analysis exploiting current and historical data. sigma Ori E is the prototypical helium-strong star that has been shown to harbor a strong magnetic field, as well as a magnetosphere, consisting of two clouds of plasma forced by magnetic and centrifugal forces to co-rotate with the star on its 1.19 day period. The Rigidly Rotating Magnetosphere (RRM) model of Townsend & Owocki (2005) approximately reproduces the observed variations in longitudinal field strength, photometric brightness, Halpha emission, and various other observables. There are, however, small discrepancies between the observations and model in the photometric light curve, which we propose arise from inhomogeneous chemical abundances on the star's surface. Using Magnetic Doppler Imaging (MDI), future work will attempt to identify the contributions to the photometric variation due to abundance spots and due to circumstellar material.

We report on the detection of a strong, organized magnetic field in the helium-variable early B-type star HR 7355 using spectropolarimetric data obtained with ESPaDOnS on CFHT by the MiMeS large program. We also present results from new V-band differential photometry obtained with the CTIO 0.9m telescope. We investigate the longitudinal field, using a technique called Least-Squares Deconvolution (LSD), and the rotational period of HR 7355. These new observations strongly support the proposal that HR 7355 harbors a structured magnetosphere similar to that in the prototypical helium-strong star, sigma Ori E.

The Magnetism in Massive Stars (MiMeS) Project is a consensus collaboration among many of the foremost international researchers of the physics of hot, massive stars, with the basic aim of understanding the origin, evolution and impact of magnetic fields in these objects. At the time of writing, MiMeS Large Programs have acquired over 950 high-resolution polarised spectra of about 150 individual stars with spectral types from B5-O4, discovering new magnetic fields in a dozen hot, massive stars. The quality of this spectral and magnetic mat\'eriel is very high, and the Collaboration is keen to connect with colleagues capable of exploiting the data in new or unforeseen ways. In this paper we review the structure of the MiMeS observing programs and report the status of observations, data modeling and development of related theory.

For early type magnetic stars slow, at most moderate rotational velocities have been considered an observational fact. The detection of a multi-kilogauss magnetic field in the B2Vpn star with P approx. 0.52d and v sin approx. 300km/s has brought down this narrative. We have obtained more than 100 high-resolution, high-S/N echelle spectra in 2009. These spectra provide the most detailed description of the variability of any He-strong star to date. The circumstellar environment is dominated by a rotationally locked magnetosphere out to several stellar radii, causing hydrogen emission. The photosphere is characterized by surface chemical abundance inhomogeneities, with much stronger amplitudes, at least for helium, than slower rotating stars like \sigma Ori E. The highly complex rotational line profile modulations of metal lines are probably a consequence the equatorial gravity darkening of HR 7355, and thus may offer an independent measurement of the von Zeipel parameter \beta .

A magnetic field and rotational line profile variability (lpv) is found in the He-weak star HR2949. The field measured from metallic lines varies in a clearly non-sinusoidal way, and shows a phase lag relative to the morphologically similar HeI equivalent width variations. The surface abundance patterns are strong and complex, and visible even in the hydrogen lines.

FN CMa is visually double with a separation of about 0.6arcsec. Sixty high-cadence VLT/UVES spectra permit the A and B components to be disentangled, as the relative contribution of each star to the total light entering the spectrograph fluctuates between exposures due to changes in seeing. Component A exhibits rapid line-profile variations, leading us to attribute the photometric variability seen by HIPPARCOS (with a derived P=0.08866d) to this component. From a total of 122 archival and new echelle spectra it is shown that component A is an SB1 binary with an orbital period of 117.55 days. The eccentricity of 0.6 may result in tidal modulation of the pulsation(s) of component Aa.

We discuss the recent detection of a strong, organized magnetic field in the bright, broad-line B2V star, HD 142184, using the ESPaDOnS spectropolarimeter on the CFHT as part of the Magnetism in Massive Stars (MiMeS) survey. We find a rotational period of 0.50833 days, making it the fastest-rotating, non-degenerate magnetic star ever detected. Like the previous rapid-rotation record holder HR 7355 (also discovered by MiMeS: Oksala et al. 2010, Rivinius et al. 2010), this star shows emission line variability that is diagnostic of a structured magnetosphere.

We obtained relatively high signal-to-noise X-ray spectral data of the early massive star tau Sco (B0.2V) with the Suzaku XIS instrument. This source displays several unusual features that motivated our study: (a) redshifted absorption in UV P Cygni lines to approximately +250 km/s suggestive of infalling gas, (b) unusually hard X-ray emission requiring hot plasma at temperatures in excess of 10 MK whereas most massive stars show relatively soft X-rays at a few MK, and (c) a complex photospheric magnetic field of open and closed field lines. In an attempt to understand the hard component better, X-ray data were obtained at six roughly equally spaced phases within the same epoch of tau Sco's 41 day rotation period. The XIS instrument has three operable detectors: XIS1 is back illuminated with sensitivity down to 0.2 keV; XIS0 and XIS2 are front illuminated with sensivitity only down to 0.4 keV and have overall less effective area than XIS1. The XIS0 and XIS3 detectors show relatively little variability. In contrast, there is a 4sigma detection of a 4% drop in the count rate of the XIS1 detector at one rotational phase. In addition, all three detectors show a 3% increase in count rate at the same phase. The most optimistic prediction of X-ray variability allows for a 40% change in the count rate, particularly near phases where we have pointings. Observed modulations in the X-ray light curve on the rotation cycle is an order of magnitude smaller than this, which places new stringent constraints on future modeling of this interesting magnetic massive star.

We present new U-band photometry of the magnetic Helium-strong star Sigma Ori E, obtained over 2004-2009 using the SMARTS 0.9-m telescope at Cerro Tololo Inter-American Observatory. When combined with historical measurements, these data constrain the evolution of the star's 1.19 d rotation period over the past three decades. We are able to rule out a constant period at the p_null = 0.05% level, and instead find that the data are well described (p_null = 99.3%) by a period increasing linearly at a rate of 77 ms per year. This corresponds to a characteristic spin-down time of 1.34 Myr, in good agreement with theoretical predictions based on magnetohydrodynamical simulations of angular momentum loss from magnetic massive stars. We therefore conclude that the observations are consistent with Sigma Ori E undergoing rotational braking due to its magnetized line-driven wind.

We report the detection of a strong, organized magnetic field in the helium-variable early B-type star HR 7355 using spectropolarimetric data obtained with ESPaDOnS on the 3.6-m Canada-France-Hawaii Telescope within the context of the Magnetism in Massive Stars (MiMeS) Large Program. HR 7355 is both the most rapidly rotating known main-sequence magnetic star and the most rapidly rotating helium-strong star, with $v \sin i$ = 300 $\pm$ 15 km s$^{-1}$ and a rotational period of 0.5214404 $\pm$ 0.0000006 days. We have modeled our eight longitudinal magnetic field measurements assuming an oblique dipole magnetic field. Constraining the inclination of the rotation axis to be between $38^{\circ}$ and $86^{\circ}$, we find the magnetic obliquity angle to be between $30^{\circ}$ and $85^{\circ}$, and the polar strength of the magnetic field at the stellar surface to be between 13-17 kG. The photometric light curve constructed from HIPPARCOS archival data and new CTIO measurements shows two minima separated by 0.5 in rotational phase and occurring 0.25 cycles before/after the magnetic extrema. This photometric behavior coupled with previously-reported variable emission of the H$\alpha$ line (which we confirm) strongly supports the proposal that HR 7355 harbors a structured magnetosphere similar to that in the prototypical helium-strong star, $\sigma$ Ori E.

The B2Vn star HR 7355 is found to be a He-rich magnetic star. Spectropolarimetric data were obtained with FORS1 at UT2 on Paranal observatory to measure the disk-averaged longitudinal magnetic field at various phases of the presumed 0.52 d cycle. A variable magnetic field with strengths between B_z = -2200 and +3200G was found, with confidence limits of 100 to 130G. The field topology is that of an oblique dipole, while the star itself is seen about equator-on. In the intensity spectra the HeI-lines show the typical equivalent width variability of He-strong stars, usually attributed to surface abundance spots. The amplitudes of the equivalent width variability of the HeI lines are extraordinarily strong compared to other cases. These results not only put HR 7355 unambiguously among the early-type magnetic stars, but confirm its outstanding nature: With v sin i = 320 km/s the parameter space in which He-strong stars are known to exist has doubled in terms of rotational velocity.

We fit every emission line in the high-resolution Chandra grating spectrum of zeta Pup with an empirical line profile model that accounts for the effects of Doppler broadening and attenuation by the bulk wind. For each of sixteen lines or line complexes that can be reliably measured, we determine a best-fitting fiducial optical depth, tau_* = kappa*Mdot/4{pi}R_{\ast}v_{\infty}, and place confidence limits on this parameter. These sixteen lines include seven that have not previously been reported on in the literature. The extended wavelength range of these lines allows us to infer, for the first time, a clear increase in tau_* with line wavelength, as expected from the wavelength increase of bound-free absorption opacity. The small overall values of tau_*, reflected in the rather modest asymmetry in the line profiles, can moreover all be fit simultaneously by simply assuming a moderate mass-loss rate of 3.5 \pm 0.3 \times 10^{-6} Msun/yr, without any need to invoke porosity effects in the wind. The quoted uncertainty is statistical, but the largest source of uncertainty in the derived mass-loss rate is due to the uncertainty in the elemental abundances of zeta Pup, which affects the continuum opacity of the wind, and which we estimate to be a factor of two. Even so, the mass-loss rate we find is significantly below the most recent smooth-wind H-alpha mass-loss rate determinations for zeta Pup, but is in line with newer determinations that account for small-scale wind clumping. If zeta Pup is representative of other massive stars, these results will have important implications for stellar and galactic evolution.

Shear flow instabilities can profoundly affect the diffusion of momentum in jets, stars, and disks. The Richardson criterion gives a sufficient condition for instability of a shear flow in a stratified medium. The velocity gradient $V'$ can only destabilize a stably stratified medium with squared Brunt-Vaisala frequency $N^2$ if $V'^2/4>N^2$. We find this is no longer true when the medium is a magnetized plasma. We investigate the effect of stable stratification on magnetic field and velocity profiles unstable to magneto-shear instabilities, i.e., instabilities which require the presence of both magnetic field and shear flow. We show that a family of profiles originally studied by Tatsuno & Dorland (2006) remain unstable even when $V'^2/4<N^2$, violating the Richardson criterion. However, not all magnetic fields can result in a violation of the Richardson criterion. We consider a class of flows originally considered by Kent (1968), which are destabilized by a constant magnetic field, and show that they become stable when $V'^2/4<N^2$, as predicted by the Richardson criterion. This suggests that magnetic free energy is required to violate the Richardson criterion. This work implies that the Richardson criterion cannot be used when evaluating the ideal stability of a sheared, stably stratified, and magnetized plasma. We briefly discuss the implications for astrophysical systems.

We present spectropolarimetric radio images of the supernova remnant (SNR) G296.5+10.0 at frequencies near 1.4 GHz, observed with the Australia Telescope Compact Array. By applying rotation measure (RM) synthesis to the data, a pixel-by-pixel map of Faraday rotation has been produced for the entire remnant. We find G296.5+10.0 to have a highly ordered RM structure, with mainly positive RMs (mean RM of +28 rad/m**2) on the eastern side and negative RMs (mean RM of -14 rad/m**2) on the western side, indicating a magnetic field which is directed away from us on one side and toward us on the other. We consider several possible mechanisms for creating the observed RM pattern. Neither Faraday rotation in foreground interstellar gas nor in a homogeneous ambient medium swept up by the SNR shell can easily explain the magnitude and sign of the observed RM pattern. Instead, we propose that the observed RMs are the imprint of an azimuthal magnetic field in the stellar wind of the progenitor star. Specifically, we calculate that a swept-up magnetized wind from a red supergiant can produce RMs of the observed magnitude, while the azimuthal pattern of the magnetic field at large distances from the star naturally produces the anti-symmetric RM pattern observed. Expansion into such a wind can possibly also account for the striking bilateral symmetry of the SNR's radio and X-ray morphologies.

We report the detection of a strong, organised magnetic field in the O9IV star HD 57682, using spectropolarimetric observations obtained with ESPaDOnS at the 3.6-m Canada-France-Hawaii Telescope within the context of the Magnetism in Massive Stars (MiMeS) Large Program. From the fitting of our spectra using NLTE model atmospheres we determined that HD 57682 is a $17^{+19}_{-9}$ M$_{\odot}$ star with a radius of $7.0^{+2.4}_{-1.8}$ R$_\odot$, and a relatively low mass-loss rate of $1.4^{+3.1}_{-0.95}\times10^{-9}$ M$_{\odot}$ yr$^{-1}$. The photospheric absorption lines are narrow, and we use the Fourier transform technique to infer $v\sin i=15\pm3$ km s$^{-1}$. This $v\sin i$ implies a maximum rotational period of 31.5 d, a value qualitatively consistent with the observed variability of the optical absorption and emission lines, as well as the Stokes $V$ profiles and longitudinal field. Using a Bayesian analysis of the velocity-resolved Stokes $V$ profiles to infer the magnetic field characteristics, we tentatively derive a dipole field strength of $1680^{+134}_{-356}$ G. The derived field strength and wind characteristics imply a wind that is strongly confined by the magnetic field.

This is a white paper submitted to the Stars and Stellar Evolution (SSE) Science Frontier Panel (SFP) of the NRC's 2010 Astronomy and Astrophysics Decadal Survey. The white paper is endorsed by the NSF Physics Frontier Center for Magnetic Self-Organization in Laboratory and Astrophysical Plasmas (CMSO).

A new scheme for incorporating radiative cooling in hydrodynamical codes is presented, centered around exact integration of the governing semi-discrete cooling equation. Using benchmark calculations based on the cooling downstream of a radiative shock, I demonstrate that the new scheme outperforms traditional explicit and implicit approaches in terms of accuracy, while remaining competitive in terms of execution speed.

In the subclass of high-mass X-ray binaries known as "microquasars", relativistic hadrons in the jets launched by the compact object can interact with cold protons from the star's radiatively driven wind, producing pions that then quickly decay into gamma rays. Since the resulting gamma-ray emissivity depends on the target density, the detection of rapid variability in microquasars with GLAST and the new generation of Cherenkov imaging arrays could be used to probe the clumped structure of the stellar wind. We show here that the fluctuation in gamma rays can be modeled using a "porosity length" formalism, usually applied to characterize clumping effects. In particular, for a porosity length defined by h=l/f, i.e. as the ratio of the characteristic size l of clumps to their volume filling factor f, we find that the relative fluctuation in gamma-ray emission in a binary with orbital separation a scales as sqrt(h/pi a) in the "thin-jet" limit, and is reduced by a factor 1/sqrt(1 + phi a/(2 l)) for a jet with a finite opening angle phi. For a thin jet and quite moderate porosity length h ~ 0.03 a, this implies a ca. 10 % variation in the gamma-ray emission. Moreover, the illumination of individual large clumps might result in isolated flares, as has been recently observed in some massive gamma-ray binaries.

We examine the angular momentum loss and associated rotational spin-down for magnetic hot stars with a line-driven stellar wind and a rotation-aligned dipole magnetic field. Our analysis here is based on our previous 2-D numerical MHD simulation study that examines the interplay among wind, field, and rotation as a function of two dimensionless parameters, W(=Vrot/Vorb) and 'wind magnetic confinement', $\eta_\ast$ defined below. We compare and contrast the 2-D, time variable angular momentum loss of this dipole model of a hot-star wind with the classical 1-D steady-state analysis by Weber and Davis (WD), who used an idealized monopole field to model the angular momentum loss in the solar wind. Despite the differences, we find that the total angular momentum loss averaged over both solid angle and time follows closely the general WD scaling $\dot {J} \sim \dot {M} \Omega R_A^2$. The key distinction is that for a dipole field Alfv\`en radius $R_A$ is significantly smaller than for the monopole field WD used in their analyses. This leads to a slower stellar spin-down for the dipole field with typical spin-down times of order 1 Myr for several known magnetic massive stars.

We examine the angular momentum loss and associated rotational spindown for magnetic hot stars with a line-driven stellar wind and a rotation-aligned dipole magnetic field. Our analysis here is based on our previous 2-D numerical MHD simulation study that examines the interplay among wind, field, and rotation as a function of two dimensionless parameters, one characterizing the wind magnetic confinement ($\eta_{\ast} \equiv B_{eq}^{2} R_{\ast}^{2}/{\dot M} v_{\infty}$), and the other the ratio ($W \equiv V_{rot}/V_{orb}$) of stellar rotation to critical (orbital) speed. We compare and contrast the 2-D, time variable angular momentum loss of this dipole model of a hot-star wind with the classical 1-D steady-state analysis by Weber and Davis (WD), who used an idealized monopole field to model the angular momentum loss in the solar wind. Despite the differences, we find that the total angular momentum loss ${\dot J}$ averaged over both solid angle and time follows closely the general WD scaling ${\dot J} = (2/3) {\dot M} \Omega R_{A}^{2}$, where ${\dot M}$ is the mass loss rate, $\Omega$ is the stellar angular velocity, and $R_{A}$ is a characteristic Alfv\'{e}n radius. However, a key distinction here is that for a dipole field, this Alfv\'{e}n radius has a strong-field scaling $R_{A}/R_{\ast} \approx \eta_{\ast}^{1/4}$, instead of the scaling $R_{A}/R_{\ast} \sim \sqrt{\eta_{\ast}}$ for a monopole field. This leads to a slower stellar spindown time that in the dipole case scales as $\tau_{spin} = \tau_{mass} 1.5k/\sqrt{\eta_{\ast}}$, where $\tau_{mass} \equiv M/{\dot M}$ is the characteristic mass loss time, and $k$ is the dimensionless factor for stellar moment of inertia. The full numerical scaling relation we cite gives typical spindown times of order 1 Myr for several known magnetic massive stars.

The photometric variations due to magnetically confined material around He-strong stars are investigated within the framework of the Rigidly Rotating Magnetosphere (RRM) model. For dipole field topologies, the model is used to explore how the morphology of light curves evolves in response to changes to the observer inclination, magnetic obliquity, rotation rate and optical depth. The general result is that double-minimum light curves arise when the obliquity and/or inclination are close to 90 degrees; no light variations are seen in the opposite limit; and for intermediate cases, single-minimum light curves occur. These findings are interpreted with the aid of a simple, analytical torus model, paving the way for the development of new photometric-based constraints on the fundamental parameters of He-strong stars. Illustrative applications to five stars in the class are presented.

We present results of a 90 ksec Chandra X-ray observation of the young sigma Orionis cluster (age ~3 Myr) obtained with the High Energy Transmission Grating Spectrometer. We use the high resolution grating spectrum and moderate resolution CCD spectrum of the massive central star sigma Ori AB (O9.5V + B0.5V) to test wind shock theories of X-ray emission and also analyze the high spatial resolution zero-order ACIS-S image of the central cluster region. Chandra detected 42 X-ray sources on the primary CCD (ACIS-S3). All but five have near-IR or optical counterparts and about one-fourth are variable. Notable high-mass stellar detections are sigma Ori AB, the magnetic B star sigma Ori E, and the B5V binary HD 37525. Most of the other detections have properties consistent with lower mass K or M-type stars. We present the first X-ray spectrum of the unusual infrared source IRS1 located 3.3 arc-sec north of sigma Ori AB, which is likely an embedded T Tauri star whose disk/envelope is being photoevaporated by sigma Ori AB. We focus on the radiative wind shock interpretation of the soft luminous X-ray emission from sigma Ori AB, but also consider possible alternatives including magnetically-confined wind shocks and colliding wind shocks. Its emission lines show no significant asymmetries or centroid shifts and are moderately broadened to HWHM ~ 264 km/s, or one-fourth the terminal wind speed. Forbidden lines in He-like ions are formally undetected, implying strong UV suppression. The Mg XI triplet forms in the wind acceleration zone within one stellar radius above the surface. These X-ray properties are consistent in several respects with the predictions of radiative wind shock theory for an optically thin wind, but explaining the narrow line widths presents a challenge to the theory.

We investigate whether angular momentum transport due to unstable pulsation modes can play a significant role in the rotational evolution of massive stars. We find that these modes can redistribute appreciable angular momentum, and moreover trigger shear-instability mixing in the molecular weight gradient zone adjacent to stellar cores, with significant evolutionary impact.

Using archival spectroscopic and photometric data, we searched for massive stars with Balmer-emission consistent with magnetically confined circumstellar material. HR 7355 is a formerly unknown He-strong star showing Balmer emission. At V=6.02 mag, it is one of the brightest objects simultaneously showing anomalous helium absorption and hydrogen emission. Among similar objects, only sigma Ori E has so far been subjected to any systematic analysis of the circumstellar material responsible for the emission. We argue that the double-wave photometric period of 0.52d corresponds to the rotation period. In tandem with the high projected equatorial velocity, v sin i=320 km/s, this short period suggests that HR 7355 is the most rapidly rotating He-strong star known to date; a class that was hitherto expected to host stars with slow to moderate rotation only.

Building upon our previous MHD simulation study of magnetic channeling in radiatively driven stellar winds, we examine here the additional dynamical effects of stellar {\em rotation} in the (still) 2-D axisymmetric case of an aligned dipole surface field. In addition to the magnetic confinement parameter $\eta_{\ast}$ introduced in Paper I, we characterize the stellar rotation in terms of a parameter $W \equiv V_{\rm{rot}}/V_{\rm{orb}}$ (the ratio of the equatorial surface rotation speed to orbital speed), examining specifically models with moderately strong rotation $W =$ 0.25 and 0.5, and comparing these to analogous non-rotating cases. Defining the associated Alfv\'{e}n radius $R_{\rm{A}} \approx \eta_{\ast}^{1/4} \Rstar$ and Kepler corotation radius $R_{\rm{K}} \approx W^{-2/3} \Rstar$, we find rotation effects are weak for models with $R_{\rm{A}} < R_{\rm{K}}$, but can be substantial and even dominant for models with $R_{\rm{A}} \gtwig R_{\rm{K}}$. In particular, by extending our simulations to magnetic confinement parameters (up to $\eta_{\ast} = 1000$) that are well above those ($\eta_{\ast} = 10$) considered in Paper I, we are able to study cases with $R_{\rm{A}} \gg R_{\rm{K}}$; we find that these do indeed show clear formation of the {\em rigid-body} disk predicted in previous analytic models, with however a rather complex, dynamic behavior characterized by both episodes of downward infall and outward breakout that limit the buildup of disk mass. Overall, the results provide an intriguing glimpse into the complex interplay between rotation and magnetic confinement, and form the basis for a full MHD description of the rigid-body disks expected in strongly magnetic Bp stars like $\sigma$ Ori E.

By quantitatively fitting simple emission line profile models that include both atomic opacity and porosity to the Chandra X-ray spectrum of $\zeta$ Pup, we are able to explore the trade-offs between reduced mass-loss rates and wind porosity. We find that reducing the mass-loss rate of $\zeta$ Pup by roughly a factor of four, to 1.5 \times 10^{-6} M_sun/yr, enables simple non-porous wind models to provide good fits to the data. If, on the other hand, we take the literature mass-loss rate of 6 \times 10^{-6} M_sun/yr, then to produce X-ray line profiles that fit the data, extreme porosity lengths -- of $h_{\infty} \approx 3$ Rstar -- are required. Moreover, these porous models do not provide better fits to the data than the non-porous, low optical depth models. Additionally, such huge porosity lengths do not seem realistic in light of 2-D numerical simulations of the wind instability.

We present multiple epochs of H-alpha spectroscopy for 47 members of the open cluster NGC 3766 to investigate the long term variability of its Be stars. Sixteen of the stars in this sample are Be stars, including one new discovery. Of these, we observe an unprecedented 11 Be stars that undergo disk appearances and/or near disappearances in our H-alpha spectra, making this the most variable population of Be stars known to date. NGC 3766 is therefore an excellent location to study the formation mechanism of Be star disks. From blue optical spectra of 38 cluster members and existing Stromgren photometry of the cluster, we also measure rotational velocities, effective temperatures, and polar surface gravities to investigate the physical and evolutionary factors that may contribute to the Be phenomenon. Our analysis also provides improvements to the reddening and distance of NGC 3766, and we find E(B-V) = 0.22 +/- 0.03 and (V-M_V)_0 = 11.6 +/- 0.2, respectively. The Be stars are not associated with a particular stage of main-sequence evolution, but they are a population of rapidly rotating stars with a velocity distribution generally consistent with rotation at 70-80% of the critical velocity, although systematic effects probably underestimate the true rotational velocities so that the rotation is much closer to critical. Our measurements of the changing disk sizes are consistent with the idea that transitory, nonradial pulsations contribute to the formation of these highly variable disks.

In this contribution, I discuss some basic techniques that can be used to simulate radiative transfer through porous media. As specific examples, I consider scattering transfer through a clumped slab, and X-ray emission line formation in a clumped wind.

To what extent can pulsational instabilities resolve the mass-loss problem of massive stars? How important is pulsation in structuring and modulating the winds of these stars? What role does pulsation play in redistributing angular momentum in massive stars? Although I cannot offer answers to these questions, I hope at the very least to explain how they come to be asked.

We introduce a new Rigid-Field Hydrodynamics approach to modeling the magnetospheres of massive stars in the limit of very-strong magnetic fields. Treating the field lines as effectively rigid, we develop hydrodynamical equations describing the 1-dimensional flow along each, subject to pressure, radiative, gravitational, and centrifugal forces. We solve these equations numerically for a large ensemble of field lines, to build up a 3-dimensional time-dependent simulation of a model star with parameters similar to the archetypal Bp star sigma Ori E. Since the flow along each field line can be solved for independently of other field lines, the computational cost of this approach is a fraction of an equivalent magnetohydrodynamical treatment. The simulations confirm many of the predictions of previous analytical and numerical studies. Collisions between wind streams from opposing magnetic hemispheres lead to strong shock heating. The post-shock plasma cools initially via X-ray emission, and eventually accumulates into a warped, rigidly rotating disk defined by the locus of minima of the effective (gravitational plus centrifugal) potential. But a number of novel results also emerge. For field lines extending far from the star, the rapid area divergence enhances the radiative acceleration of the wind, resulting in high shock velocities (up to ~3,000 km/s) and hard X-rays. Moreover, the release of centrifugal potential energy continues to heat the wind plasma after the shocks, up to temperatures around twice those achieved at the shocks themselves. Finally, in some circumstances the cool plasma in the accumulating disk can oscillate about its equilibrium position, possibly due to radiative cooling instabilities in the adjacent post-shock regions.

Gamma-rays can be produced by the interaction of a relativistic jet and the matter of the stellar wind in the subclass of massive X-ray binaries known as "microquasars". The relativistic jet is ejected from the surroundings of the compact object and interacts with cold protons from the stellar wind, producing pions that then quickly decay into gamma-rays. Since the resulting gamma-ray emissivity depends on the target density, the detection of rapid variability in microquasars with GLAST and the new generation of Cherenkov imaging arrays could be used to probe the clumped structure of the stellar wind. In particular, we show here that the relative fluctuation in gamma rays may scale with the square root of the ratio of porosity length to binary separation, \sqrt{h/a}, implying for example a ca. 10 % variation in gamma ray emission for a quite moderate porosity, h/a ~ 0.01.

We present a semi-analytic model for shaping the nebula around eta Carinae that accounts for the simultaneous production of bipolar lobes and an equatorial disk through a rotating surface explosion. Material is launched normal to the surface of an oblate rotating star with an initial kick velocity that scales approximately with the local escape speed. Thereafter, ejecta follow ballistic orbital trajectories, feeling only a central force corresponding to a radiatively reduced gravity. Our model is conceptually similar to the wind-compressed disk model of Bjorkman & Cassinelli, but we modify it to an explosion instead of a steady line-driven wind, we include a rotationally-distorted star, and we treat the dynamics somewhat differently. Continuum-driving avoids the disk inhibition that normally operates in line-driven winds. Our model provides a simple method by which rotating hot stars can simultaneously produce intrinsically bipolar and equatorial mass ejections, without an aspherical environment or magnetic fields. Although motivated by eta Carinae, the model may have generic application to other LBVs, B[e] stars, or SN1987A's nebula. When near-Eddington radiative driving is less influential, our model generalizes to produce bipolar morphologies without disks, as seen in many PNe.

We present results of numerical simulations of wind variability in the magnetic B1 IVe star $\beta$ Cephei. 2D-MHD simulations are used to determine the structure of the wind. From these wind models we calculate line profiles for different aspect angles to simulate rotation. The results are compared with the observed UV wind line profiles.

We examine the stability of l=1 and l=2 g modes in a pair of nitrogen-rich Wolf-Rayet stellar models characterized by differing hydrogen abundances. We find that modes with intermediate radial orders are destabilized by a kappa mechanism operating on an opacity bump at an envelope temperature log T ~ 6.25. This `deep opacity bump' is due primarily to L-shell bound-free transitions of iron. Periods of the unstable modes span ~ 11-21 hr in the model containing some hydrogen, and ~ 3-12 hr in the hydrogen-depleted model. Based on the latter finding, we suggest that self-excited g modes may be the source of the 9.8 hr-periodic variation of WR 123 recently reported by Lefevre et al. (2005).

We present 2D MHD simulations of the radiatively driven outflow from a rotating hot star with a dipole magnetic field aligned with the star's rotation axis. We focus primarily on a model with moderately rapid rotation (half the critical value), and also a large magnetic confinement parameter, $\eta_{\ast} \equiv B_{\ast}^2 R_{\ast}^{2} / \dot{M} V_{\infty} = 600$. The magnetic field channels and torques the wind outflow into an equatorial, rigidly rotating disk extending from near the Kepler corotation radius outwards. Even with fine-tuning at lower magnetic confinement, none of the MHD models produce a stable Keplerian disk. Instead, material below the Kepler radius falls back on to the stellar surface, while the strong centrifugal force on material beyond the corotation escape radius stretches the magnetic loops outwards, leading to episodic breakout of mass when the field reconnects. The associated dissipation of magnetic energy heats material to temperatures of nearly $10^{8}$K, high enough to emit hard (several keV) X-rays. Such \emph{centrifugal mass ejection} represents a novel mechanism for driving magnetic reconnection, and seems a very promising basis for modeling X-ray flares recently observed in rotating magnetic Bp stars like $\sigma$ Ori E.

I review the basic concepts of the Rigidly Rotating Magnetosphere model for the circumstellar plasma distribution around the helium-strong star sigma Ori E. I demonstrate that the model can furnish a good fit to the photometric, spectroscopic and magneticv ariability exhibited by this star, and argue that the variability of other helium-strong stars may be amenable to a similar interpretation.

The stability of retrograde mixed modes in rotating B-type stars is investigated. It is found that these modes are susceptible to kappa-mechanism excitation, due to the iron opacity bump at log T ~ 5.3. The findings are discussed in the context of the pulsation of SPB and Be stars.

We present new UBVRI observations of the magnetic Bp star sigma Ori E. The basic features of the star's lightcurve have not changed since the previous monitoring by Hesser et al. (1977), indicating that the star's magnetosphere has remained stable over the past three decades. Interestingly, we find a rotation period that is slightly longer than in the Hesser et al. (1977) analysis, suggesting possible spindown of the star.

I examine the stability of retrograde mixed modes in rotating B-type stars. These modes can be regarded as a hybrid between the Rossby modes that arise from conservation of vorticity, and the Poincare modes that are gravity waves modified by the Coriolis force. Using a non-adiabatic pulsation code based around the traditional approximation, I find that the modes are unstable in mid- to late-B type stars, due to the same iron-bump opacity mechanism usually associated with SPB and $\beta$ Cep stars. At one half of the critical rotation rate, the instability for $m=1...4$ modes spans the spectral types B4 to A0. Inertial-frame periods of the unstable modes range from 100 days down to a fraction of a day, while normalized growth rates can reach in excess of $10^{-5}$. I discuss the relevance of these findings to SPB and pulsating Be stars, and to the putative Maia class of variable star. I also outline some of the questions raised by this discovery of a wholly-new class of pulsational instability in early-type stars.

We attempt to characterize the observed variability of the magnetic helium-strong star sigma Ori E in terms of a recently developed rigidly rotating magnetosphere model. This model predicts the accumulation of circumstellar plasma in two co-rotating clouds, situated in magnetohydrostatic equilibrium at the intersection between magnetic and rotational equators. We find that the model can reproduce well the periodic modulations observed in the star's light curve, H alpha emission-line profile, and longitudinal field strength, confirming that it furnishes an essentially correct, quantitative description of the star's magnetically controlled circumstellar environment.

We report on four Chandra grating observations of the oblique magnetic rotator theta^1 Ori C (O5.5 V) covering a wide range of viewing angles with respect to the star's 1060 G dipole magnetic field. We employ line-width and centroid analyses to study the dynamics of the X-ray emitting plasma in the circumstellar environment, as well as line-ratio diagnostics to constrain the spatial location, and global spectral modeling to constrain the temperature distribution and abundances of the very hot plasma. We investigate these diagnostics as a function of viewing angle and analyze them in conjunction with new MHD simulations of the magnetically channeled wind shock mechanism on theta^1 Ori C. This model fits all the data surprisingly well, predicting the temperature, luminosity, and occultation of the X-ray emitting plasma with rotation phase.

Recent observations of blue stragglers by De Marco et al. (2004) have revealed continuum deficits on the blue side of the Balmer discontinuity, leading these authors to infer the presence of discs around the stars. This intriguing possibility may throw light on aspects of the mechanisms responsible for at least some of these objects; current theories of blue straggler formation invoke stellar collisions or interacting binaries, both of which appear capable of forming a circumstellar disc. However, by synthesizing photospheric spectra for models of rotating blue stragglers, we demonstrate that the Balmer jump enhancements can be wholly attributed to the influence of oblateness and gravity darkening on the formation of the continuum. Therefore, we are led to conclude that the observations of De Marco et al. can be ascribed a more prosaic explanation, that of rapid stellar rotation arising from the merger/interaction formation process.

This paper explores the effect of rotation on the kappa-mechanism instability of slowly pulsating B stars. A new nonadiabatic code, that adopts the so-called `traditional approximation' to treat the Coriolis force, is used to investigate the influence exerted by rotation over the stability of stellar models covering the mass range 2.5 M_sun <= M_* <= 13.0 M_sun. The principal finding is that, for all modes considered apart from the prograde sectoral class, rotation shifts the kappa-mechanism instability toward higher luminosities and effective temperatures; these shifts are accompanied by broadenings in the extent of instability strips. Such behaviour is traced to the shortening of mode periods under the action of the Coriolis force. Instability strips associated with prograde sectoral modes behave rather differently, being shifted to marginally lower luminosities and effective temperatures under the influence of rotation. The implications of these results are discussed in the context of the observational scarcity of pulsation in B-type stars having significant rotation; various scenarios are explored to explain the apparent dichotomy between theory and observations. Furthermore, the possible significance of the findings to Be stars is briefly examined.

We discuss three different observational diagnostics related to disks around hot stars: absorption line determinations of rotational velocities of Be stars; polarization diagnostics of circumstellar disks; and X-ray line diagnostics of one specific magnetized hot star, $\theta^1$ Ori C. Some common themes that emerge from these studies include (a) the benefits of having a specific physical model as a framework for interpreting diagnostic data; (b) the importance of combining several different types of observational diagnostics of the same objects; and (c) that while there is often the need to reinterpret traditional diagnostics in light of new theoretical advances, there are many new and powerful diagnostics that are, or will soon be, available for the study of disks around hot stars.

We present a semi-analytical approach for modeling circumstellar emission from rotating hot stars with a strong dipole magnetic field tilted at an arbitrary angle to the rotation axis. By assuming the rigid-field limit in which material driven (e.g., in a wind outflow) from the star is forced to remain in strict rigid-body co-rotation, we are able to solve for the effective centrifugal-plus-gravitational potential along each field line, and thereby identify the location of potential minima where material is prone to accumulate. Applying basic scalings for the surface mass flux of a radiatively driven stellar wind, we calculate the circumstellar density distribution that obtains once ejected plasma settles into hydrostatic stratification along field lines. The resulting accumulation surface resembles a rigidly rotating, warped disk, tilted such that its average surface normal lies between the rotation and magnetic axes. Using a simple model of the plasma emissivity, we calculate time-resolved synthetic line spectra for the disk. Initial comparisons show an encouraging level of correspondence with the observed rotational phase variations of Balmer-line emission profiles from magnetic Bp stars like sigma Ori E.

We argue that, in general, observational studies of Be-star rotation have paid insufficient attention to the effects of equatorial gravity darkening. We present new line-profile calculations that emphasize the insensitivity of line width to rotation for fast rotators. Coupled with a critical review of observational procedures, these calculations suggest that the observational parameter v sin(i) may systematically underestimate the true projected equatorial rotation velocity, ve sin(i), by some tens of per cent for rapid rotators. The crucial implication of this work is that Be stars may be rotating much closer to their critical velocities than is generally supposed, bringing a range of new processes into contention for the elusive physical mechanism responsible for the circumstellar disk thought to be central to the Be phenomenon.

We discuss the possibility of using OH and H2O gigamasers to trace the redshift distribution of luminous, dust-obscured, star-forming galaxies. It has long been thought that ultraluminous, interacting galaxies should host gigamasers due to their vast pumping IR luminosity, the large column density of molecules available to populate the maser states and the turbulent motion of the gas in these dynamically complex systems which allows unsaturated maser emission. OH masers may thus be well-suited to the redshift-blind detection of ultraluminous and hyperluminous infrared galaxies (L_FIR >= 10^12 L_sun) such as those uncovered by the SCUBA submillimetre camera. The bandwidth requirement is low, <1 GHz for z=1-10 (lower still if additional redshift constraints are available) and the dual-line 1665-/1667-MHz OH spectral signature can act as a check on the reality of detections.

A qualitative study of the surface trapping of low-frequency non-radial g-modes in rotating early-type stars is undertaken within the Cowling, adiabatic and traditional approximations. It is demonstrated that, at frequencies below a cut-off, waves cannot be fully trapped within the star, and will leak through the surface. Expressions for the cut-off frequency are derived in both the non-rotating and rotating cases; it is found from these expressions that the cut-off frequency increases with the rotation rate for all but prograde sectoral modes. The results are of possible relevance to the 53 Per and SPB classes of variable star, which exhibit pulsation frequencies of the same order of magnitude as the cut-off frequencies found for the stellar model. It is suggested that observations either of an upper limit on variability periods (corresponding to the cut-off), or of line-profile variations due to leaking modes, may permit asteroseismological studies of the outer layers of these systems.