-
On the spatial distribution of electron energy loss due to gyro-cooling in hot-star magnetospheres
Authors:
B. Das,
S. P. Owocki
Abstract:
Hot magnetic stars often exhibit incoherent circularly polarized radio emission thought to arise from gyro-synchrotron emission by energetic electrons trapped in the circumstellar magnetosphere. Theoretical scalings for electron acceleration by magnetic reconnection driven by centrifugal breakout match well the empirical scalings for observed radio luminosity with both the magnetic field strength…
▽ More
Hot magnetic stars often exhibit incoherent circularly polarized radio emission thought to arise from gyro-synchrotron emission by energetic electrons trapped in the circumstellar magnetosphere. Theoretical scalings for electron acceleration by magnetic reconnection driven by centrifugal breakout match well the empirical scalings for observed radio luminosity with both the magnetic field strength and the stellar rotation rate. This paper now examines how energetic electrons introduced near the top of closed magnetic loops are subsequently cooled by the energy loss associated with their gyro-synchrotron radio emission. For sample assumed distributions for energetic electron deposition about the loop apex, we derive the spatial distribution of the radiated energy from such "gyro-cooling". For sub-relativistic electrons, we show explicitly that this is independent of the input energy, but also find that even extensions to the relativistic regime still yield a quite similar spatial distribution. However, cooling by coulomb collisions with even a modest ambient density of thermal electrons can effectively quench the emission from sub-relativistic electrons, indicating that the observed radio emission likely stems from relativistic electrons that are less affected by such collisional cooling. The overall results form an initial basis for computing radio emission spectra in future models that account for such cooling and multimode excitation about the fundamental gyro-frequency. Though motivated in the context of hot-stars, the basic results here could also be applied to gyro-emission in any dipole magnetospheres, including those of ultra-cool dwarfs and even (exo)-planets.
△ Less
Submitted 2 August, 2023;
originally announced August 2023.
-
MOBSTER -- VII. Using light curves to infer magnetic and rotational properties of stars with centrifugal magnetospheres
Authors:
I. D. Berry,
M. E. Shultz,
S. P. Owocki,
A. ud-Doula
Abstract:
Early-type B stars with strong magnetic fields and rapid rotation form centrifugal magnetospheres (CMs), as the relatively weak stellar wind becomes magnetically confined and centrifugally supported above the Kepler co-rotation radius. CM plasma is concentrated at and above the Kepler co-rotation radius at the intersection between the rotation and magnetic field axis. Stellar rotation can cause th…
▽ More
Early-type B stars with strong magnetic fields and rapid rotation form centrifugal magnetospheres (CMs), as the relatively weak stellar wind becomes magnetically confined and centrifugally supported above the Kepler co-rotation radius. CM plasma is concentrated at and above the Kepler co-rotation radius at the intersection between the rotation and magnetic field axis. Stellar rotation can cause these clouds of material to intersect the viewer's line-of-sight, leading to photometric eclipses. However, for stars with strong ($\sim 10\,{\rm kG}$) magnetic fields and rapid rotation, CMs can become optically thick enough for emission to occur via electron scattering. Using high-precision space photometry from a sample of stars with strong H$α$ emission, we apply simulated light curves from the Rigidly Rotating Magnetosphere model to directly infer magnetic and rotational properties of these stars. By comparing the values inferred from photometric modelling to those independently determined by spectropolarimetry, we find that magnetic obliquity angle $β$, viewer inclination $i$ and critical rotation fraction $W$ can be approximately recovered for 3 of the 4 stars studied here. However, there are large discrepancies between the optical depth at the Kepler radius $τ_{\rm K}$ expected from magnetometry, and the values required to match the observations. We show that $τ_{\rm K}$ of order unity is needed to reasonably match the light curve morphology of our sample stars.
△ Less
Submitted 6 June, 2023;
originally announced June 2023.
-
3D MHD models of the centrifugal magnetosphere from a massive star with an oblique dipole field
Authors:
Asif ud-Doula,
Stanley P. Owocki,
Christopher Russell,
Marc Gagne,
Simon Daley-Yates
Abstract:
We present results from new self-consistent 3D MHD simulations of the magnetospheres from massive stars with a dipole magnetic axis that has a non-zero obliquity angle ($β$) to the star's rotation axis. As an initial direct application, we compare the global structure of co-rotating disks for nearly aligned ($β=5^o$) versus half-oblique ($β=45^o$) models, both with moderately rapid rotation (…
▽ More
We present results from new self-consistent 3D MHD simulations of the magnetospheres from massive stars with a dipole magnetic axis that has a non-zero obliquity angle ($β$) to the star's rotation axis. As an initial direct application, we compare the global structure of co-rotating disks for nearly aligned ($β=5^o$) versus half-oblique ($β=45^o$) models, both with moderately rapid rotation ($\sim$ 0.5 critical). We find that accumulation surfaces broadly resemble the forms predicted by the analytic Rigidly Rotating Magnetosphere (RRM) model, but the mass buildup to near the critical level for centrifugal breakout against magnetic confinement distorts the field from the imposed initial dipole. This leads to an associated warping of the accumulation surface toward the rotational equator, with the highest density concentrated in {\em wings} centered on the intersection between the magnetic and rotational equators. These MHD models can be used to synthesize rotational modulation of photometric absorption and H$α$ emission for a direct comparison with observations.
△ Less
Submitted 27 January, 2023;
originally announced January 2023.
-
Discovery and origin of the radio emission from the multiple stellar system KQVel
Authors:
P. Leto,
L. M. Oskinova,
C. S. Buemi,
M. E. Shultz,
F. Cavallaro,
C. Trigilio,
G. Umana,
L. Fossati,
I. Pillitteri,
J. Krticka,
R. Ignace,
C. Bordiu,
F. Bufano,
G. Catanzaro,
L. Cerrigone,
M. Giarrusso,
A. Ingallinera,
S. Loru,
S. P. Owocki,
K. A. Postnov,
S. Riggi,
J. Robrade,
F. Leone
Abstract:
KQVel is a binary system composed of a slowly rotating magnetic Ap star with a companion of unknown nature. In this paper, we report the detection of its radio emission. We conducted a multi-frequency radio campaign using the ATCA interferometer (band-names: 16cm, 4cm, and 15mm). The target was detected in all bands. The most obvious explanation for the radio emission is that it originates in the…
▽ More
KQVel is a binary system composed of a slowly rotating magnetic Ap star with a companion of unknown nature. In this paper, we report the detection of its radio emission. We conducted a multi-frequency radio campaign using the ATCA interferometer (band-names: 16cm, 4cm, and 15mm). The target was detected in all bands. The most obvious explanation for the radio emission is that it originates in the magnetosphere of the Ap star, but this is shown unfeasible. The known stellar parameters of the Ap star enable us to exploit the scaling relationship for non-thermal gyro-synchrotron emission from early-type magnetic stars. This is a general relation demonstrating how radio emission from stars with centrifugal magnetospheres is supported by rotation. Using KQVel's parameters the predicted radio luminosity is more than five orders of magnitudes lower than the measured one. The extremely long rotation period rules out the Ap star as the source of the observed radio emission. Other possible explanations for the radio emission from KQVel, involving its unknown companion, have been explored. A scenario that matches the observed features (i.e. radio luminosity and spectrum, correlation to X-rays) is a hierarchical stellar system, where the possible companion of the magnetic star is a close binary (possibly of RSCVn type) with at least one magnetically active late-type star. To be compatible with the total mass of the system, the last scenario places strong constraints on the orbital inclination of the KQVel stellar system.
△ Less
Submitted 28 July, 2022;
originally announced July 2022.
-
Method and new tabulations for flux-weighted line-opacity and radiation line-force in supersonic media
Authors:
L. G. Poniatowski,
N. D. Kee,
J. O. Sundqvist,
F. A. Driessen,
N. Moens,
S. P. Owocki,
K. G. Gayley,
L. Decin,
A. de Koter,
H. Sana
Abstract:
In accelerating and supersonic media, the interaction of photons with spectral lines can be of ultimate importance. However, fully accounting for such line forces currently can only be done by specialised codes in 1-D steady-state flows. More general cases and higher dimensions require alternative approaches. We presented a comprehensive and fast method for computing the radiation line-force using…
▽ More
In accelerating and supersonic media, the interaction of photons with spectral lines can be of ultimate importance. However, fully accounting for such line forces currently can only be done by specialised codes in 1-D steady-state flows. More general cases and higher dimensions require alternative approaches. We presented a comprehensive and fast method for computing the radiation line-force using tables of spectral line-strength distribution parameters, which can be applied in arbitrary (multi-D, time-dependent) simulations, including those accounting for the line-deshadowing instability, to compute the appropriate opacities. We assumed local thermodynamic equilibrium to compute a flux-weighted line opacity from $>4$ million spectral lines. We derived the spectral line strength and tabulated the corresponding line-distribution parameters for a range of input densities $ρ\in[10^{-20},10^{-10}]gcm^{-3}$ and temperatures $T\in[10^4,10^{4.7}]K$. We found that the variation of the line distribution parameters plays an essential role in setting the wind dynamics in our models. In our benchmark study, we also found a good overall agreement between the O-star mass-loss rates of our models and those derived from steady-state studies using more detailed radiative transfer. Our models reinforce that self-consistent variation of the line-distribution parameters is important for the dynamics of line-driven flows. Within a well-calibrated O-star regime, our results support the proposed methodology. In practice, utilising the provided tables, yielded a factor $>100$ speed-up in computational time compared to specialised 1-D model-atmosphere codes of line-driven winds, which constitutes an important step towards efficient multi-D simulations. We conclude that our method and tables are ready to be exploited in various radiation-hydrodynamic simulations where the line force is important.
△ Less
Submitted 21 April, 2022;
originally announced April 2022.
-
Centrifugal breakout reconnection as the electron acceleration mechanism powering the radio magnetospheres of early-type stars
Authors:
Stanley P Owocki,
Matt E. Shultz,
Asif ud-Doula,
Poonam Chandra,
Barnali Das,
Paulo Leto
Abstract:
Magnetic B-stars often exhibit circularly polarized radio emission thought to arise from gyrosynchrotron emission by energetic electrons trapped in the circumstellar magnetosphere. Recent empirical analyses show that the onset and strength of the observed radio emission scale with both the magnetic field strength and the stellar rotation rate. This challenges the existing paradigm that the energet…
▽ More
Magnetic B-stars often exhibit circularly polarized radio emission thought to arise from gyrosynchrotron emission by energetic electrons trapped in the circumstellar magnetosphere. Recent empirical analyses show that the onset and strength of the observed radio emission scale with both the magnetic field strength and the stellar rotation rate. This challenges the existing paradigm that the energetic electrons are accelerated in the current sheet between opposite-polarity field lines in the outer regions of magnetised stellar winds, which includes no role for stellar rotation. Building on recent success in explaining a similar rotation-field dependence of H$α$ line emission in terms of a model in which magnetospheric density is regulated by centrifugal breakout (CBO), we examine here the potential role of the associated CBO-driven magnetic reconnection in accelerating the electrons that emit the observed gyrosynchrotron radio. We show in particular that the theoretical scalings for energy production by CBO reconnection match well the empirical trends for observed radio luminosity, with a suitably small, nearly constant conversion efficiency $ε\approx 10^{-8}$. We summarize the distinct advantages of our CBO scalings over previous associations with an electromotive force, and discuss the potential implications of CBO processes for X-rays and other observed characteristics of rotating magnetic B-stars with centrifugal magnetospheres.
△ Less
Submitted 10 February, 2022;
originally announced February 2022.
-
Electron Scattering Emission in the Light Curves of Stars with Centrifugal Magnetospheres
Authors:
I. D. Berry,
S. P. Owocki,
M. E. Shultz,
A. ud-Doula
Abstract:
Strongly magnetic, rapidly rotating B-type stars with relatively weak winds form centrifugal magnetospheres (CMs), as the stellar wind becomes magnetically confined above the Kepler co-rotation radius. Approximating the magnetic field as a dipole tilted by an angle $β$ with respect to the rotation axis, the CM plasma is concentrated in clouds at and above the Kepler radius along the intersection o…
▽ More
Strongly magnetic, rapidly rotating B-type stars with relatively weak winds form centrifugal magnetospheres (CMs), as the stellar wind becomes magnetically confined above the Kepler co-rotation radius. Approximating the magnetic field as a dipole tilted by an angle $β$ with respect to the rotation axis, the CM plasma is concentrated in clouds at and above the Kepler radius along the intersection of the rotational and magnetic equatorial planes. Stellar rotation can bring such clouds in front of the stellar disk, leading to absorption of order 0.1 magnitude ($\sim 10 \%$ of continuum flux). However some stars with prominent CMs, such as $σ$ Ori E, show an emission bump in addition to absorption dips, which has been so far unexplained. We show that emission can occur from electron scattering toward the observer when CM clouds are projected off the stellar limb. Using the Rigidly Rotating Magnetosphere model, modified with a centrifugal breakout density scaling, we present a model grid of photometric light curves spanning parameter space in observer inclination angle $i$, magnetic obliquity angle $β$, critical rotation fraction $W$, and optical depth at the Kepler radius $τ_{\text{K}}$. We show that $τ_{\text{K}}$ of order unity can produce emission bumps of the magnitude $\sim 0.05$ seen in $σ$ Ori E. We discuss the implications for modeling the light curves of CM stars, as well as future work for applying the radiative transfer model developed here to 3D MHD simulations of CMs.
△ Less
Submitted 1 February, 2022;
originally announced February 2022.
-
MOBSTER -- VI. The crucial influence of rotation on the radio magnetospheres of hot stars
Authors:
M. E. Shultz,
S. P. Owocki,
A. ud-Doula,
A. Biswas,
D. Bohlender,
P. Chandra,
B. Das,
A. David-Uraz,
V. Khalack,
O. Kochukhov,
J. D. Landstreet,
P. Leto,
D. Monin,
C. Neiner,
Th. Rivinius,
G. A. Wade
Abstract:
Numerous magnetic hot stars exhibit gyrosynchrotron radio emission. The source electrons were previously thought to be accelerated to relativistic velocities in the current sheet formed in the middle magnetosphere by the wind opening magnetic field lines. However, a lack of dependence of radio luminosity on the wind power, and a strong dependence on rotation, has recently challenged this paradigm.…
▽ More
Numerous magnetic hot stars exhibit gyrosynchrotron radio emission. The source electrons were previously thought to be accelerated to relativistic velocities in the current sheet formed in the middle magnetosphere by the wind opening magnetic field lines. However, a lack of dependence of radio luminosity on the wind power, and a strong dependence on rotation, has recently challenged this paradigm. We have collected all radio measurements of magnetic early-type stars available in the literature. When constraints on the magnetic field and/or the rotational period are not available, we have determined these using previously unpublished spectropolarimetric and photometric data. The result is the largest sample of magnetic stars with radio observations that has yet been analyzed: 131 stars with rotational and magnetic constraints, of which 50 are radio-bright. We confirm an obvious dependence of gyrosynchrotron radiation on rotation, and furthermore find that accounting for rotation neatly separates stars with and without detected radio emission. There is a close correlation between H$α$ emission strength and radio luminosity. These factors suggest that radio emission may be explained by the same mechanism responsible for H$α$ emission from centrifugal magnetospheres, i.e. centrifugal breakout (CBO), however, whereas the H$α$-emitting magnetosphere probes the cool plasma before breakout, radio emission is a consequence of electrons accelerated in centrifugally-driven magnetic reconnection.
△ Less
Submitted 14 January, 2022;
originally announced January 2022.
-
Ultraviolet Spectropolarimetry with Polstar: Clumping and Mass-loss Rate Corrections
Authors:
Ken Gayley,
Jorick S. Vink,
Asif ud-Doula,
Alexandre David-Uraz,
Richard Ignace,
Raman Prinja,
Nicole St-Louis,
Sylvia Ekström,
Yaël Nazé,
Tomer Shenar,
Paul A. Scowen,
Natallia Sudnik,
Stan P. Owocki,
Jon O. Sundqvist,
Florian A. Driessen,
Levin Hennicker
Abstract:
The most massive stars are thought to lose a significant fraction of their mass in a steady wind during the main-sequence and blue supergiant phases. This in turn sets the stage for their further evolution and eventual supernova, with consequences for ISM energization and chemical enrichment. Understanding these processes requires accurate observational constraints on the mass-loss rates of the mo…
▽ More
The most massive stars are thought to lose a significant fraction of their mass in a steady wind during the main-sequence and blue supergiant phases. This in turn sets the stage for their further evolution and eventual supernova, with consequences for ISM energization and chemical enrichment. Understanding these processes requires accurate observational constraints on the mass-loss rates of the most luminous stars, which can also be used to test theories of stellar wind generation. In the past, mass-loss rates have been characterized via collisional emission processes such as H$α$ and free-free radio emission, but these so-called "density squared" diagnostics require correction in the presence of widespread clumping. Recent observational and theoretical evidence points to the likelihood of a ubiquitously high level of such clumping in hot-star winds, but quantifying its effects requires a deeper understanding of the complex dynamics of radiatively driven winds. Furthermore, large-scale structures arising from surface anisotropies and propagating throughout the wind can further complicate the picture by introducing further density enhancements, affecting mass-loss diagnostics. Time series spectroscopy of UV resonance lines with high resolution and high signal-to-noise are required to better understand this complex dynamics, and help correct "density squared" diagnostics of mass-loss rates. The proposed Polstar mission easily provides the necessary resolution at the sound-speed scale of 20 km s$^{-1}$, on three dozen bright targets with signal-to noise an order of magnitude above that of the celebrated IUE MEGA campaign, via continuous observations that track structures advecting through the wind in real time.
△ Less
Submitted 12 December, 2021; v1 submitted 22 November, 2021;
originally announced November 2021.
-
Detection of an extremely strong magnetic field in the double-degenerate binary merger product HD 144941
Authors:
M. E. Shultz,
O. Kochukhov,
J. Labadie-Bartz,
A. David-Uraz,
S. P. Owocki
Abstract:
HD 144941 is an extreme He (EHe) star, a rare class of subdwarf OB star formed from the merger of two white dwarf (WD) stars. Uniquely amongst EHe stars, its light curve has been reported to be modulated entirely by rotation, suggesting the presence of a magnetic field. Here we report the first high-resolution spectropolarimetric observations of HD 144941, in which we detect an extremely strong ma…
▽ More
HD 144941 is an extreme He (EHe) star, a rare class of subdwarf OB star formed from the merger of two white dwarf (WD) stars. Uniquely amongst EHe stars, its light curve has been reported to be modulated entirely by rotation, suggesting the presence of a magnetic field. Here we report the first high-resolution spectropolarimetric observations of HD 144941, in which we detect an extremely strong magnetic field both in circular polarization (with a line-of-sight magnetic field averaged over the stellar disk $\langle B_z \rangle \sim -8$ kG) and in Zeeman splitting of spectral lines (yielding a magnetic modulus of $\langle B \rangle \sim 17$ kG). We also report for the first time weak H$α$ emission consistent with an origin an a Centrifugal Magnetosphere (CM). HD 144941's atmospheric parameters could be consistent with either a subdwarf or a main sequence (MS) star, and its surface abundances are neither similar to other EHe stars nor to He-strong magnetic stars. However, its H$α$ emission properties can only be reproduced if its mass is around 1 M$_\odot$, indicating that it must be a post-MS object. Since there is no indication of binarity, it is unlikely to be a stripped star, and was therefore most likely produced in a WD merger. HD 144941 is therefore further evidence that mergers are a viable pathway for the generation of fossil magnetic fields.
△ Less
Submitted 23 July, 2021;
originally announced July 2021.
-
Dynamically inflated wind models of classical Wolf-Rayet stars
Authors:
L. G. Poniatowski,
J. O. Sundqvist,
N. D. Kee,
S. P. Owocki,
P. Marchant,
L. Decin,
A. de Koter,
L. Mahy,
H. Sana
Abstract:
Vigorous mass loss in the classical Wolf-Rayet (WR) phase is important for the late evolution and final fate of massive stars. We develop spherically symmetric time-dependent and steady-state hydrodynamical models of the radiation-driven wind outflows and associated mass loss from classical WR stars. The simulations are based on combining the opacities typically used in static stellar structure an…
▽ More
Vigorous mass loss in the classical Wolf-Rayet (WR) phase is important for the late evolution and final fate of massive stars. We develop spherically symmetric time-dependent and steady-state hydrodynamical models of the radiation-driven wind outflows and associated mass loss from classical WR stars. The simulations are based on combining the opacities typically used in static stellar structure and evolution models with a simple parametrised form for the enhanced line-opacity expected within a supersonic outflow. Our simulations reveal high mass-loss rates initiated in deep and hot optically thick layers around T\approx 200kK. The resulting velocity structure is non-monotonic and can be separated into three phases: i) an initial acceleration to supersonic speeds ii) stagnation and even deceleration, and iii) an outer region of rapid re-acceleration. The characteristic structures seen in converged steady-state simulations agree well with the outflow properties of our time-dependent models. By directly comparing our dynamic simulations to corresponding hydrostatic models, we demonstrate explicitly that the need to invoke extra energy transport in convectively inefficient regions of stellar structure and evolution models is merely an artefact of enforcing a hydrostatic outer boundary. Moreover, the "dynamically inflated" inner regions of our simulations provide a natural explanation for the often-found mismatch between predicted hydrostatic WR radii and those inferred from spectroscopy. Finally, we contrast our simulations with alternative recent WR wind models based on co-moving frame radiative transfer for computing the radiation force. Since CMF transfer currently cannot handle non-monotonic velocity fields, the characteristic deceleration regions found here are avoided in such simulations by invoking an ad-hoc very high degree of clumping.
△ Less
Submitted 10 December, 2020;
originally announced December 2020.
-
Simulating the formation of $η$ Carinae's surrounding nebula through unstable triple evolution and stellar merger-induced eruption
Authors:
Ryosuke Hirai,
Philipp Podsiadlowski,
Stanley P. Owocki,
Fabian R. N. Schneider,
Nathan Smith
Abstract:
$η$ Carinae is an extraordinary massive star famous for its 19th century Great Eruption and the surrounding Homunculus nebula ejected in that event. The cause of this eruption has been the centre of a long-standing mystery. Recent observations, including light-echo spectra of the eruption, suggest that it most likely resulted from a stellar merger in an unstable triple system. Here we present a de…
▽ More
$η$ Carinae is an extraordinary massive star famous for its 19th century Great Eruption and the surrounding Homunculus nebula ejected in that event. The cause of this eruption has been the centre of a long-standing mystery. Recent observations, including light-echo spectra of the eruption, suggest that it most likely resulted from a stellar merger in an unstable triple system. Here we present a detailed set of theoretical calculations for this scenario; from the dynamics of unstable triple systems and the mass ejection from close binary encounters, to the mass outflow from the eruption caused by the stellar merger and the post-merger wind phase. In our model the bipolar post-merger wind is the primary agent for creating the Homunculus, as it sweeps up external eruption ejecta into a thin shell. Our simulations reproduce many of the key aspects of the shape and kinematics of both the Homunculus nebula and its complex surrounding structure, providing strong support for the merger-in-a-triple scenario.
△ Less
Submitted 24 February, 2021; v1 submitted 24 November, 2020;
originally announced November 2020.
-
Line-drag damping of Alfvén waves in radiatively driven winds of magnetic massive stars
Authors:
F. A. Driessen,
N. D. Kee,
J. O. Sundqvist,
S. P. Owocki
Abstract:
Line-driven stellar winds from massive (OB) stars are subject to a strong line-deshadowing instability. Recently, spectropolarimetric surveys have collected ample evidence that a subset of Galactic massive stars hosts strong surface magnetic fields. We investigate here the propagation and stability of magneto-radiative waves in such a magnetised, line-driven wind. Our analytic, linear stability an…
▽ More
Line-driven stellar winds from massive (OB) stars are subject to a strong line-deshadowing instability. Recently, spectropolarimetric surveys have collected ample evidence that a subset of Galactic massive stars hosts strong surface magnetic fields. We investigate here the propagation and stability of magneto-radiative waves in such a magnetised, line-driven wind. Our analytic, linear stability analysis includes line-scattering from the stellar radiation, and accounts for both radial and non-radial perturbations. We establish a bridging law for arbitrary perturbation wavelength after which we analyse separately the long- and short-wavelength limits. While long-wavelength radiative and magnetic waves are found to be completely decoupled, a key result is that short-wavelength, radially propagating Alfvén waves couple to the scattered radiation field and are strongly damped due to the line-drag effect. This damping of magnetic waves in a scattering-line-driven flow could have important effects on regulating the non-linear wind dynamics, and so might also have strong influence on observational diagnostics of the wind structure and clumping of magnetic line-driven winds.
△ Less
Submitted 12 October, 2020;
originally announced October 2020.
-
How the breakout-limited mass in B-star centrifugal magnetospheres controls their circumstellar H-alpha emission
Authors:
Stanley P. Owocki,
Matt E. Shultz,
Asif ud-Doula,
Jon O. Sundqvist,
Richard H. D. Townsend,
Steven R. Cranmer
Abstract:
Strongly magnetic B-type stars with moderately rapid rotation form `centrifugal magnetospheres' (CMs), from the magnetic trapping of stellar wind material in a region above the Kepler co-rotation radius. A longstanding question is whether the eventual loss of such trapped material occurs from gradual drift and/or diffusive leakage, or through sporadic `{\em centrifugal break out}' (CBO) events, wh…
▽ More
Strongly magnetic B-type stars with moderately rapid rotation form `centrifugal magnetospheres' (CMs), from the magnetic trapping of stellar wind material in a region above the Kepler co-rotation radius. A longstanding question is whether the eventual loss of such trapped material occurs from gradual drift and/or diffusive leakage, or through sporadic `{\em centrifugal break out}' (CBO) events, wherein magnetic tension can no longer contain the built-up mass. We argue here that recent empirical results for Balmer-$α$ emission from such B-star CMs strongly favor the CBO mechanism. Most notably, the fact that the onset of such emission depends mainly on the field strength at the Kepler radius, and is largely {\em independent} of the stellar luminosity, strongly disfavors any drift/diffusion process, for which the net mass balance would depend on the luminosity-dependent wind feeding rate. In contrast, we show that in a CBO model the {\em maximum confined mass} in the magnetosphere is independent of this wind feeding rate, and has a dependence on field strength and Kepler radius that naturally explains the empirical scalings for the onset of H$α$ emission, its associated equivalent width, and even its line profile shapes. However, the general lack of observed Balmer emission in late-B and A-type stars could still be attributed to a residual level of diffusive or drift leakage that does not allow their much weaker winds to fill their CMs to the breakout level needed for such emission; alternatively this might result from a transition to a metal-ion wind that lacks the requisite Hydrogen.
△ Less
Submitted 25 September, 2020;
originally announced September 2020.
-
The effects of surface fossil magnetic fields on massive star evolution: II. Implementation of magnetic braking in MESA and implications for the evolution of surface rotation in OB stars
Authors:
Z. Keszthelyi,
G. Meynet,
M. E. Shultz,
A. David-Uraz,
A. ud-Doula,
R. H. D. Townsend,
G. A. Wade,
C. Georgy,
V. Petit,
S. P. Owocki
Abstract:
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 magneti…
▽ More
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$, $Ω/Ω_{\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$, $Ω/Ω_{\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.
△ Less
Submitted 21 January, 2020; v1 submitted 17 January, 2020;
originally announced January 2020.
-
Magnetic OB[A] Stars with TESS: probing their Evolutionary and Rotational properties (MOBSTER) - I. First-light observations of known magnetic B and A stars
Authors:
A. David-Uraz,
C. Neiner,
J. Sikora,
D. M. Bowman,
V. Petit,
S. Chowdhury,
G. Handler,
M. Pergeorelis,
M. Cantiello,
D. H. Cohen,
C. Erba,
Z. Keszthelyi,
V. Khalack,
O. Kobzar,
O. Kochukhov,
J. Labadie-Bartz,
C. C. Lovekin,
R. MacInnis,
S. P. Owocki,
H. Pablo,
M. E. Shultz,
A. ud-Doula,
G. A. Wade,
the MOBSTER Collaboration
Abstract:
In this paper we introduce the MOBSTER collaboration and lay out its scientific goals. We present first results based on the analysis of nineteen previously known magnetic O, B and A stars observed in 2-minute cadence in sectors 1 and 2 of the Transiting Exoplanet Survey Satellite (TESS) mission. We derive precise rotational periods from the newly obtained light curves and compare them to previous…
▽ More
In this paper we introduce the MOBSTER collaboration and lay out its scientific goals. We present first results based on the analysis of nineteen previously known magnetic O, B and A stars observed in 2-minute cadence in sectors 1 and 2 of the Transiting Exoplanet Survey Satellite (TESS) mission. We derive precise rotational periods from the newly obtained light curves and compare them to previously published values. We also discuss the overall photometric phenomenology of the known magnetic massive and intermediate-mass stars and propose an observational strategy to augment this population by taking advantage of the high-quality observations produced by TESS.
△ Less
Submitted 25 April, 2019;
originally announced April 2019.
-
Hydrodynamical simulations and similarity relations for eruptive mass loss from massive stars
Authors:
Stanley P. Owocki,
Ryo Hirai,
Philipp Podsiadlowski,
Fabian Schneider
Abstract:
Motivated by the eruptive mass loss inferred from Luminous Blue Variable (LBV) stars, we present 1D hydrodynamical simulations of the response from sudden energy injection into the interior of a very massive ($100 \, M_\odot$) star. For a fiducial case with total energy addition set to a factor $f=0.5$ of the net stellar binding energy, and applied within the stellar envelope, we detail the dynami…
▽ More
Motivated by the eruptive mass loss inferred from Luminous Blue Variable (LBV) stars, we present 1D hydrodynamical simulations of the response from sudden energy injection into the interior of a very massive ($100 \, M_\odot$) star. For a fiducial case with total energy addition set to a factor $f=0.5$ of the net stellar binding energy, and applied within the stellar envelope, we detail the dynamical response that leads to ejection of the outermost $7.2 \, M_\odot$. We find that the ejecta's variations in time $t$ and radius $r$ for the velocity $v$, density $ρ$, and temperature $T$ are quite well fit by similarity forms in the variable $r/t \approx v$. Specifically the scaled density follows a simple exponential decline $ρt^{3} \sim \exp (-r/v_{\rm o} t)$. This `exponential similarity' leads to analytic scaling relations for total ejecta mass $ΔM$ and kinetic energy $ΔK$ that agree well with the hydrodynamical simulations, with the specific-energy-averaged speed related to the exponential scale speed $v_{\rm o}$ through ${\bar v} \equiv \sqrt{2 ΔK/ΔM} = \sqrt{12} \, v_{\rm o}$, and a value comparable to the star's surface escape speed, $v_{\rm esc}$. Models with energy added in the core develop a surface shock breakout that propels an initial, higher-speed ejecta ($>$5000km s$^{-1}$), but the bulk of the ejected material still follows the same exponential similarity scalings with ${\bar v} \approx v_{\rm esc}$. A broader parameter study examines how the ejected mass and energy depends on the energy-addition factor $f$, for three distinct model series that locate the added energy in either the core, envelope, or near-surface. We conclude by discussing the relevance of these results for understanding LBV outbursts and other eruptive phenomena, such as failed supernovae and pulsational pair instability events.
△ Less
Submitted 17 February, 2019;
originally announced February 2019.
-
Extreme resonance line profile variations in the ultraviolet spectra of NGC 1624-2: probing the giant magnetosphere of the most strongly magnetized known O-type star
Authors:
A. David-Uraz,
C. Erba,
V. Petit,
A. W. Fullerton,
F. Martins,
N. R. Walborn,
R. MacInnis,
R. H. Barbá,
D. H. Cohen,
J. Maíz Apellániz,
Y. Nazé,
S. P. Owocki,
J. O. Sundqvist,
A. ud-Doula,
G. A. Wade
Abstract:
In this paper, we present high-resolution HST/COS observations of the extreme magnetic O star NGC 1624-2. These represent the first ultraviolet spectra of this archetypal object. We examine the variability of its wind-sensitive resonance lines, comparing it to that of other known magnetic O stars. In particular, the observed variations in the profiles of the CIV and SiIV doublets between low state…
▽ More
In this paper, we present high-resolution HST/COS observations of the extreme magnetic O star NGC 1624-2. These represent the first ultraviolet spectra of this archetypal object. We examine the variability of its wind-sensitive resonance lines, comparing it to that of other known magnetic O stars. In particular, the observed variations in the profiles of the CIV and SiIV doublets between low state and high state are the largest observed in any magnetic O-type star, consistent with the expected properties of NGC 1624-2's magnetosphere. We also observe a redshifted absorption component in the low state, a feature not seen in most stars. We present preliminary modelling efforts based on the Analytic Dynamical Magnetosphere (ADM) formalism, demonstrating the necessity of using non-spherically symmetric models to determine wind/magnetospheric properties of magnetic O stars.
△ Less
Submitted 25 November, 2018;
originally announced November 2018.
-
Radiation Transport Through Super-Eddington Stellar Winds
Authors:
Joyce A. Guzik,
Chris Fryer,
Todd J. Urbatsch,
Stanley P. Owocki
Abstract:
We present results of simulations to assess the feasibility of modeling outflows from massive stars using the Los Alamos 3-D radiation hydrodynamics code Cassio developed for inertial confinement fusion (ICF) applications. We find that a 1-D stellar envelope simulation relaxes into hydrostatic equilibrium using computing resources that would make the simulation tractable in 2-D. We summarize next…
▽ More
We present results of simulations to assess the feasibility of modeling outflows from massive stars using the Los Alamos 3-D radiation hydrodynamics code Cassio developed for inertial confinement fusion (ICF) applications. We find that a 1-D stellar envelope simulation relaxes into hydrostatic equilibrium using computing resources that would make the simulation tractable in 2-D. We summarize next steps to include more physics fidelity and model the response to a large and abrupt energy deposition at the base of the envelope.
△ Less
Submitted 2 June, 2018;
originally announced June 2018.
-
Characterizing the turbulent porosity of stellar-wind structure generated by the line-deshadowing instability
Authors:
Stanley P. Owocki,
Jon O. Sundqvist
Abstract:
We analyze recent 2D simulations of the nonlinear evolution of the line-deshadowing instability (LDI) in hot-star winds, to quantify how the associated highly clumped density structure can lead to a "turbulent porosity" reduction in continuum absorption and/or scattering. The basic method is to examine the statistical variations of mass column as a function of path length, and fit these to analyti…
▽ More
We analyze recent 2D simulations of the nonlinear evolution of the line-deshadowing instability (LDI) in hot-star winds, to quantify how the associated highly clumped density structure can lead to a "turbulent porosity" reduction in continuum absorption and/or scattering. The basic method is to examine the statistical variations of mass column as a function of path length, and fit these to analytic forms that lead to simple statistical scalings for the associated mean extinction. A key result is that one can characterize porosity effects on continuum transport in terms of a single "turbulent porosity length", found here to scale as $H \approx (f_{\rm cl} - 1) a$, where $f_{\rm cl} \equiv \left < ρ^2 \right >/\left < ρ\right >^2$ is the clumping factor in density $ρ$, and $a$ is the density autocorrelation length. For continuum absorption or scattering in an optically thick layer, we find the associated effective reduction in opacity scales as $\sim 1/\sqrt{1+τ_{\rm H}}$, where $τ_{\rm H} \equiv κρH$ is the local optical thickness of this porosity length. For these LDI simulations, the inferred porosity lengths are small, only about a couple percent of the stellar radius, $H \approx 0.02 R_\ast$. For continuum processes like bound-free absorption of X-rays that are only marginally optically thick throughout the full stellar wind, this implies $τ_{\rm H} \ll 1$, and thus that LDI-generated porosity should have little effect on X-ray transport in such winds. The formalism developed here could however be important for understanding the porous regulation of continuum-driven, super-Eddington outflows from luminous blue variables.
△ Less
Submitted 9 December, 2017;
originally announced December 2017.
-
Diffusion-plus-drift models for the mass leakage from centrifugal magnetospheres of magnetic hot-stars
Authors:
Stanley P. Owocki,
Steven R. Cranmer
Abstract:
In the subset of luminous, early-type stars with strong, large-scale magnetic fields and moderate to rapid rotation, material from the star's radiatively driven stellar wind outflow becomes trapped by closed magnetic loops, forming a centrifugally supported, co-rotating magnetosphere. We present here a semi-analytic analysis of how this quasi-steady accumulation of wind mass can be balanced by los…
▽ More
In the subset of luminous, early-type stars with strong, large-scale magnetic fields and moderate to rapid rotation, material from the star's radiatively driven stellar wind outflow becomes trapped by closed magnetic loops, forming a centrifugally supported, co-rotating magnetosphere. We present here a semi-analytic analysis of how this quasi-steady accumulation of wind mass can be balanced by losses associated with a combination of an outward, centrifugally driven drift in the region beyond the Kepler co-rotation radius, and an inward/outward diffusion near this radius. We thereby derive scaling relations for the equilibrium spatial distribution of mass, and the associated emission measure for observational diagnostics like Balmer line emission. We discuss the potential application of these relations for interpreting surveys of the emission line diagnostics for OB stars with centrifugally supported magnetospheres. For a specific model of turbulent field-line-wandering rooted in surface motions associated with the iron opacity bump, we estimate values for the associated diffusion and drift coefficients.
△ Less
Submitted 15 November, 2017;
originally announced November 2017.
-
2D wind clumping in hot, massive stars from hydrodynamical line-driven instability simulations using a pseudo-planar approach
Authors:
J. O. Sundqvist,
S. P. Owocki,
J. Puls
Abstract:
Context: Clumping in the radiation-driven winds of hot, massive stars arises naturally due to the strong, intrinsic instability of line-driving (the `LDI'). But LDI wind models have so far mostly been limited to 1D, mainly because of severe computational challenges regarding calculation of the multi-dimensional radiation force. Aims: To simulate and examine the dynamics and multi-dimensional natur…
▽ More
Context: Clumping in the radiation-driven winds of hot, massive stars arises naturally due to the strong, intrinsic instability of line-driving (the `LDI'). But LDI wind models have so far mostly been limited to 1D, mainly because of severe computational challenges regarding calculation of the multi-dimensional radiation force. Aims: To simulate and examine the dynamics and multi-dimensional nature of wind structure resulting from the LDI. Methods: We introduce a `pseudo-planar', `box-in-a-wind' method that allows us to efficiently compute the line-force in the radial and lateral directions, and then use this approach to carry out 2D radiation-hydrodynamical simulations of the time-dependent wind. Results: Our 2D simulations show that the LDI first manifests itself by mimicking the typical shell-structure seen in 1D models, but how these shells then quickly break up into complex 2D density and velocity structures, characterized by small-scale density `clumps' embedded in larger regions of fast and rarefied gas. Key results of the simulations are that density-variations in the well-developed wind statistically are quite isotropic and that characteristic length-scales are small; a typical clump size is ~0.01R at 2R, thus resulting also in rather low typical clump-masses ~10^17 g. Overall, our results agree well with the theoretical expectation that the characteristic scale for LDI-generated wind-structure is of order the Sobolev length. We further confirm some earlier results that lateral `filling-in' of radially compressed gas leads to somewhat lower clumping factors in 2D simulations than in comparable 1D models. We conclude by discussing an extension of our method toward rotating LDI wind models that exhibit an intriguing combination of large- and small-scale structure extending down to the wind base.
△ Less
Submitted 21 October, 2017;
originally announced October 2017.
-
On the optically-thick winds of Wolf-Rayet stars
Authors:
G. Gräfener,
S. P. Owocki,
L. Grassitelli,
N. Langer
Abstract:
(abridged) The strong winds of Wolf-Rayet (WR) stars are important for the mechanical and chemical feedback of the most massive stars and determine whether they end their lives as neutron stars or black holes. In this work we investigate theoretically the mass-loss properties of H-free WR stars of the nitrogen sequence (WN stars). We connect stellar structure models for He stars with wind models f…
▽ More
(abridged) The strong winds of Wolf-Rayet (WR) stars are important for the mechanical and chemical feedback of the most massive stars and determine whether they end their lives as neutron stars or black holes. In this work we investigate theoretically the mass-loss properties of H-free WR stars of the nitrogen sequence (WN stars). We connect stellar structure models for He stars with wind models for optically-thick winds and assess how both types of models can simultaneously fulfill their respective sonic-point conditions. Fixing the outer wind law and terminal wind velocity, we obtain unique solutions for the mass-loss rates of optically-thick, radiatively-driven winds of WR stars in the phase of core He-burning. The resulting mass-loss relations as a function of stellar parameters, agree well with previous empirical relations. Furthermore, we encounter stellar mass limits below which no continuous solutions exist. While these mass limits agree with observations of WR stars in the Galaxy, they are in conflict with observations in the LMC. While our results confirm in particular the slope of oft-used empirical mass-loss relations, they imply that only part of the observed WN population can be understood in the framework of the standard assumptions of a smooth transonic flow and compact stellar core. This means that alternative approaches, such as a clumped and inflated wind structure, or deviations from the diffusion limit at the sonic point may have to be invoked. Qualitatively, the existence of mass limits for the formation of WR-type winds may be relevant for the non-detection of low-mass WR stars in binary systems, which are believed to be progenitors of Type Ib/c supernovae. The sonic-point conditions derived in this work may provide a possibility to include optically-thick winds in stellar evolution models in a more physically motivated form than in current models.
△ Less
Submitted 12 October, 2017;
originally announced October 2017.
-
Super-Eddington stellar winds: unifying radiative-enthalpy vs. flux-driven models
Authors:
Stanley P. Owocki,
Richard H. D. Townsend,
Eliot Quataert
Abstract:
We derive semi-analytic solutions for optically thick, super-Eddington stellar winds, induced by an assumed steady energy addition $Δ{\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 + Δ{\dot E}$ exceed the Eddington luminosity…
▽ More
We derive semi-analytic solutions for optically thick, super-Eddington stellar winds, induced by an assumed steady energy addition $Δ{\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 + Δ{\dot E}$ exceed the Eddington luminosity, $Γ_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/Γ_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 Γ_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 = Γ_o/(1+m Γ_o) -1$, and for the final observed luminosity $L_{ obs}$, which for all allowed steady-solutions with $m < 1 - 1/Γ_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.
△ Less
Submitted 25 August, 2017;
originally announced August 2017.
-
Investigating the origin of cyclical wind variability in hot, massive stars - II. Hydrodynamical simulations of co-rotating interaction regions using realistic spot parameters for the O giant $ξ$ Persei
Authors:
A. David-Uraz,
S. P. Owocki,
G. A. Wade,
J. O. Sundqvist,
N. D. Kee
Abstract:
OB stars exhibit various types of spectral variability historically associated with wind structures, including the apparently ubiquitous discrete absorption components (DACs). These features have been proposed to be caused either by magnetic fields or non-radial pulsations. In this second paper of this series, we revisit the canonical phenomenological hydrodynamical modelling used to explain the f…
▽ More
OB stars exhibit various types of spectral variability historically associated with wind structures, including the apparently ubiquitous discrete absorption components (DACs). These features have been proposed to be caused either by magnetic fields or non-radial pulsations. In this second paper of this series, we revisit the canonical phenomenological hydrodynamical modelling used to explain the formation of DACs by taking into account modern observations and more realistic theoretical predictions. Using constraints on putative bright spots located on the surface of the O giant $ξ$ Persei derived from high precision space-based broadband optical photometry obtained with the Microvariability and Oscillations of STars (MOST) space telescope, we generate two-dimensional hydrodynamical simulations of co-rotating interaction regions in its wind. We then compute synthetic ultraviolet (UV) resonance line profiles using Sobolev Exact Integration and compare them with historical timeseries obtained by the International Ultraviolet Explorer (IUE) to evaluate if the observed behaviour of $ξ$ Persei's DACs is reproduced. Testing three different models of spot size and strength, we find that the classical pattern of variability can be successfully reproduced for two of them: the model with the smallest spots yields absorption features that are incompatible with observations. Furthermore, we test the effect of the radial dependence of ionization levels on line driving, but cannot conclusively assess the importance of this factor. In conclusion, this study self-consistently links optical photometry and UV spectroscopy, paving the way to a better understanding of cyclical wind variability in massive stars in the context of the bright spot paradigm.
△ Less
Submitted 12 June, 2017;
originally announced June 2017.
-
New Insights into the Puzzling P-Cygni Profiles of Magnetic Massive Stars
Authors:
Christiana Erba,
Alexandre David-Uraz,
Veronique Petit,
Stanley P. Owocki
Abstract:
Magnetic massive stars comprise approximately 10% of the total OB star population. Modern spectropolarimetry shows these stars host strong, stable, large-scale, often nearly dipolar surface magnetic fields of 1 kG or more. These global magnetic fields trap and deflect outflowing stellar wind material, forming an anisotropic magnetosphere that can be probed with wind-sensitive UV resonance lines. R…
▽ More
Magnetic massive stars comprise approximately 10% of the total OB star population. Modern spectropolarimetry shows these stars host strong, stable, large-scale, often nearly dipolar surface magnetic fields of 1 kG or more. These global magnetic fields trap and deflect outflowing stellar wind material, forming an anisotropic magnetosphere that can be probed with wind-sensitive UV resonance lines. Recent HST UV spectra of NGC 1624-2, the most magnetic O star observed to date, show atypically unsaturated P-Cygni profiles in the Civ resonant doublet, as well as a distinct variation with rotational phase. We examine the effect of non-radial, magnetically-channeled wind outflow on P-Cygni line formation, using a Sobolev Exact Integration (SEI) approach for direct comparison with HST UV spectra of NGC 1624-2. We demonstrate that the addition of a magnetic field desaturates the absorption trough of the P-Cygni profiles, but further efforts are needed to fully account for the observed line profile variation. Our study thus provides a first step toward a broader understanding of how strong magnetic fields affect mass loss diagnostics from UV lines.
△ Less
Submitted 27 February, 2017;
originally announced February 2017.
-
Investigating the Magnetospheres of Rapidly Rotating B-type Stars
Authors:
C. L. Fletcher,
V. Petit,
Y. Naze,
G. A. Wade,
R. H. Townsend,
S. P. Owocki,
D. H. Cohen,
A. David-Uraz,
M. Shultz
Abstract:
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 lo…
▽ More
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.
△ Less
Submitted 21 February, 2017;
originally announced February 2017.
-
Magnetic massive stars as progenitors of "heavy" stellar-mass black holes
Authors:
V. Petit,
Z. Keszthelyi,
R. MacInnis,
D. H. Cohen,
R. H. D. Townsend,
G. A. Wade,
S. L. Thomas,
S. P. Owocki,
J. Puls,
J. A. ud-Doula
Abstract:
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 w…
▽ More
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.
△ Less
Submitted 27 November, 2016;
originally announced November 2016.
-
An `Analytic Dynamical Magnetosphere' formalism for X-ray and optical emission from slowly rotating magnetic massive stars
Authors:
Stanley P. Owocki,
Asif ud-Doula,
Jon O. Sundqvist,
Veronique Petit,
David H. Cohen,
Richard H. D. Townsend
Abstract:
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 spa…
▽ More
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.
△ Less
Submitted 28 July, 2016;
originally announced July 2016.
-
The Spectral Temperature of Optically Thick Outflows with Application to Light Echo Spectra from $η$~Carinae's Giant Eruption
Authors:
Stanley P. Owocki,
Nir J. Shaviv
Abstract:
The detection by Rest et al. (2012) of light echoes from $η$ Carinae has provided important new observational constraints on the nature of its 1840's era giant eruption. Spectra of the echoes suggest a relatively cool spectral temperature of about 5500K, lower than the lower limit of about 7000K suggested in the optically thick wind outflow analysis of Davidson (1987). This has lead to a debate ab…
▽ More
The detection by Rest et al. (2012) of light echoes from $η$ Carinae has provided important new observational constraints on the nature of its 1840's era giant eruption. Spectra of the echoes suggest a relatively cool spectral temperature of about 5500K, lower than the lower limit of about 7000K suggested in the optically thick wind outflow analysis of Davidson (1987). This has lead to a debate about the viability of this steady wind model relative to alternative, explosive scenarios. Here we present an updated analysis of the wind outflow model using newer low-temperature opacity tabulations and accounting for the stronger mass loss implied by the $>$10 Msun mass now inferred for the Homunculus. A major conclusion is that, because of the sharp drop in opacity due to free electron recombination for $T<$6500K, a low temperature of about 5000K is compatible with, and indeed expected from, a wind with the extreme mass loss inferred for the eruption. Within a spherical gray model in radiative equilibrium, we derive spectral energy distributions for various assumptions for the opacity variation of the wind, providing a basis for comparisons with observed light echo spectra. The scaling results here are also potentially relevant for other highly optically thick outflows, including those from classical novae, giant eruptions of LBVs and SN Type IIn precursors. A broader issue therefore remains whether the complex, variable features observed from such eruptions are better understood in terms of a steady or explosive paradigm, or perhaps a balance of these idealizations.
△ Less
Submitted 12 June, 2016;
originally announced June 2016.
-
Accretion, Outflows, and Winds of Magnetized Stars
Authors:
M. M. Romanova,
S. P. Owocki
Abstract:
Many types of stars have strong magnetic fields that can dynamically influence the flow of circumstellar matter. In stars with accretion disks, the stellar magnetic field can truncate the inner disk and determine the paths that matter can take to flow onto the star. These paths are different in stars with different magnetospheres and periods of rotation. External field lines of the magnetosphere m…
▽ More
Many types of stars have strong magnetic fields that can dynamically influence the flow of circumstellar matter. In stars with accretion disks, the stellar magnetic field can truncate the inner disk and determine the paths that matter can take to flow onto the star. These paths are different in stars with different magnetospheres and periods of rotation. External field lines of the magnetosphere may inflate and produce favorable conditions for outflows from the disk-magnetosphere boundary. Outflows can be particularly strong in the propeller regime, wherein a star rotates more rapidly than the inner disk. Outflows may also form at the disk-magnetosphere boundary of slowly rotating stars, if the magnetosphere is compressed by the accreting matter. In isolated, strongly magnetized stars, the magnetic field can influence formation and/or propagation of stellar wind outflows. Winds from low-mass, solar-type stars may be either thermally or magnetically driven, while winds from massive, luminous O and B type stars are radiatively driven. In all of these cases, the magnetic field influences matter flow from the stars and determines many observational properties. In this chapter we review recent studies of accretion, outflows, and winds of magnetized stars with a focus on three main topics: (1) accretion onto magnetized stars; (2) outflows from the disk-magnetosphere boundary; and (3) winds from isolated massive magnetized stars. We show results obtained from global magnetohydrodynamic simulations and, in a number of cases compare global simulations with observations.
△ Less
Submitted 16 May, 2016;
originally announced May 2016.
-
Modelling the Central Constant Emission X-ray component of eta Carinae
Authors:
Christopher M. P. Russell,
Michael F. Corcoran,
Kenji Hamaguchi,
Thomas I. Madura,
Stanley P. Owocki,
D. John Hillier
Abstract:
The X-ray emission of $η$ Carinae shows multiple features at various spatial and temporal scales. The central constant emission (CCE) component is centred on the binary and arises from spatial scales much smaller than the bipolar Homunculus nebula, but likely larger than the central wind--wind collision region between the stars as it does not vary over the $\sim$2-3 month X-ray minimum when it can…
▽ More
The X-ray emission of $η$ Carinae shows multiple features at various spatial and temporal scales. The central constant emission (CCE) component is centred on the binary and arises from spatial scales much smaller than the bipolar Homunculus nebula, but likely larger than the central wind--wind collision region between the stars as it does not vary over the $\sim$2-3 month X-ray minimum when it can be observed. Using large-scale 3D smoothed particle hydrodynamics (SPH) simulations, we model both the colliding-wind region between the stars, and the region where the secondary wind collides with primary wind ejected from the previous periastron passage. The simulations extend out to one hundred semimajor axes and make two limiting assumptions (strong coupling and no coupling) about the influence of the primary radiation field on the secondary wind. We perform 3D radiative transfer calculations on the SPH output to synthesize the X-ray emission, with the aim of reproducing the CCE spectrum. For the preferred primary mass-loss rate $\dot{M}_A\approx8.5\times10^{-4}$ M$_\odot$ yr$^{-1}$, the model spectra well reproduce the observation as the strong- and no-coupling spectra bound the CCE observation for longitude of periastron $ω\approx252^\circ$, and bound/converge on the observation for $ω\approx90^\circ$. This suggests that $η$ Carinae has moderate coupling between the primary radiation and secondary wind, that both the region between the stars and the comoving collision on the backside of the secondary generate the CCE, and that the CCE cannot place constraints on the binary's line of sight. We also discuss comparisons with common X-ray fitting parameters.
△ Less
Submitted 4 March, 2016;
originally announced March 2016.
-
Line-driven ablation of circumstellar disks: I. Optically thin decretion disks of classical Oe/Be stars
Authors:
N. D. Kee,
S. P. Owocki,
J. O. Sundqvist
Abstract:
The extreme luminosities of hot, massive stars drive strong stellar winds through UV line-scattering. For OB stars with an orbiting circumstellar disk, we explore the effect of such line-scattering in ablating disk material, initially focusing on the marginally optically thin decretion disks of classical Oe and Be stars. For this we apply a multi-dimensional radiation-hydrodynamics code, assuming…
▽ More
The extreme luminosities of hot, massive stars drive strong stellar winds through UV line-scattering. For OB stars with an orbiting circumstellar disk, we explore the effect of such line-scattering in ablating disk material, initially focusing on the marginally optically thin decretion disks of classical Oe and Be stars. For this we apply a multi-dimensional radiation-hydrodynamics code, assuming optically thin ray tracing for the stellar continuum and a multi-ray Sobolev treatment of the line transfer. This accounts for desaturation of line-absorption by Keplerian shear in the disk, and associated driving by non-radial photons. Results show dense, intermediate-speed surface ablation, consistent with the strong, blue-shifted absorption seen in UV wind lines of Be shell stars. The asymptotic ablation rate is typically an order-unity factor times the stellar wind mass loss rate, leading to disk destruction times of order months to years for Be disks, consistent with observations. The much stronger radiative forces of O stars reduce this time to order days, making sustaining a disk difficult, and so providing a natural explanation for the rarity of Galactic Oe stars. Additionally, the weakened line-driving at lower metallicity implies both a reduction in the winds that help spin-down stars from near-critical rotation, and in the ablation of decretion disks, thus providing a natural explanation for the higher fraction of Classical Be stars, and the presence of Oe stars, in the Magellanic Clouds. We conclude with a discussion of future extensions to study line-driven ablation of denser, optically thick, accretion disks around pre-main-sequence massive stars.
△ Less
Submitted 25 February, 2016;
originally announced February 2016.
-
The MiMeS Survey of Magnetism in Massive Stars: Introduction and overview
Authors:
G. A. Wade,
C. Neiner,
E. Alecian,
J. H. Grunhut,
V. Petit,
B. de Batz,
D. A. Bohlender,
D. H. Cohen,
H. F. Henrichs,
O. Kochukhov,
J. D. Landstreet,
N. Manset,
F. Martins,
S. Mathis,
M. E. Oksala,
S. P. Owocki,
Th. Rivinius,
M. E. Shultz,
J. O. Sundqvist,
R. H. D. Townsend,
A. ud-Doula,
J. -C. Bouret,
J. Braithwaite,
M. Briquet,
A. C. Carciofi
, et al. (25 additional authors not shown)
Abstract:
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 s…
▽ More
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élescope 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.
△ Less
Submitted 26 November, 2015;
originally announced November 2015.
-
Hydrodynamic and radiative transfer modeling of X-ray emission from colliding WR winds: WR 140 & the Galactic center
Authors:
Christopher M. P. Russell,
Michael F. Corcoran,
Jorge Cuadra,
Stanley P. Owocki,
Q. Daniel Wang,
Kenji Hamaguchi,
Yasuharu Sugawara,
Andrew M. T. Pollock,
Timothy R. Kallman
Abstract:
Colliding Wolf-Rayet (WR) winds produce thermal X-ray emission widely observed by X-ray telescopes. In wide WR+O binaries, such as WR 140, the X-ray flux is tied to the orbital phase, and is a direct probe of the winds' properties. In the Galactic center, $\sim$30 WRs orbit the super massive black hole (SMBH) within $\sim$10", leading to a smorgasbord of wind-wind collisions. To model the X-ray em…
▽ More
Colliding Wolf-Rayet (WR) winds produce thermal X-ray emission widely observed by X-ray telescopes. In wide WR+O binaries, such as WR 140, the X-ray flux is tied to the orbital phase, and is a direct probe of the winds' properties. In the Galactic center, $\sim$30 WRs orbit the super massive black hole (SMBH) within $\sim$10", leading to a smorgasbord of wind-wind collisions. To model the X-ray emission of WR 140 and the Galactic center, we perform 3D hydrodynamic simulations to trace the complex gaseous flows, and then carry out 3D radiative transfer calculations to compute the variable X-ray spectra. The model WR 140 RXTE light curve matches the data well for all phases except the X-ray minimum associated with periastron, while the model spectra agree with the RXTE hardness ratio and the shape of the Suzaku observations throughout the orbit. The Galactic center model of the Chandra flux and spectral shape match well in the region r$<$3", but the model flux falls off too rapidly beyond this radius.
△ Less
Submitted 3 November, 2015;
originally announced November 2015.
-
Effect of scattering on the transonic solution topology and intrinsic variability of line-driven stellar winds
Authors:
Jon O. Sundqvist,
Stanley P. Owocki
Abstract:
For line-driven winds from hot, luminous OB stars, we examine the subtle but important role of diffuse, scattered radiation in determining both the topology of steady-state solutions and intrinsic variability in the transonic wind base. We use a smooth source function formalism to obtain nonlocal, integral expressions for the direct and diffuse components of the line-force that account for deviati…
▽ More
For line-driven winds from hot, luminous OB stars, we examine the subtle but important role of diffuse, scattered radiation in determining both the topology of steady-state solutions and intrinsic variability in the transonic wind base. We use a smooth source function formalism to obtain nonlocal, integral expressions for the direct and diffuse components of the line-force that account for deviations from the usual localized, Sobolev forms. As the scattering source function is reduced, we find the solution topology in the transonic region transitions from X-type, with a unique wind solution, to a nodal type, characterized by a degenerate family of solutions.
Specifically, in the idealized case of an optically thin source function and a uniformly bright stellar disk, the unique X-type solution proves to be a stable attractor to which time-dependent numerical radiation-hydrodynamical simulations relax. But in models where the scattering strength is only modestly reduced, the topology instead turns nodal, with the associated solution degeneracy now manifest by intrinsic structure and variability that persist down to the photospheric wind base. This highlights the potentially crucial role of diffuse radiation for the dynamics and variability of line-driven winds, and seriously challenges the use of Sobolev theory in the transonic wind region. Since such Sobolev-based models are commonly used in broad applications like stellar evolution and feedback, this prompts development of new wind models, not only for further quantifying the intrinsic variability found here, but also for computing new theoretical predictions of global properties like velocity laws and mass-loss rates.
△ Less
Submitted 12 August, 2015;
originally announced August 2015.
-
X-ray emission from the giant magnetosphere of the magnetic O-type star NGC 1624-2
Authors:
V. Petit,
D. H. Cohen,
G. A. Wade,
Y. Nazé,
S. P. Owocki,
J. O. Sundqvist,
A. ud-Doula,
A. Fullerton,
M. Leutenegger,
M. Gagné
Abstract:
We observed NGC 1624-2, the O-type star with the largest known magnetic field Bp~20 kG), in X-rays with the ACIS-S camera onboard the Chandra X-ray Observatory. Our two observations were obtained at the minimum and maximum of the periodic Halpha emission cycle, corresponding to the rotational phases where the magnetic field is the closest to equator-on and pole-on, respectively. With these observa…
▽ More
We observed NGC 1624-2, the O-type star with the largest known magnetic field Bp~20 kG), in X-rays with the ACIS-S camera onboard the Chandra X-ray Observatory. Our two observations were obtained at the minimum and maximum of the periodic Halpha emission cycle, corresponding to the rotational phases where the magnetic field is the closest to equator-on and pole-on, respectively. With these observations, we aim to characterise the star's magnetosphere via the X-ray emission produced by magnetically confined wind shocks. Our main findings are:
(i) The observed spectrum of NGC 1624-2 is hard, similar to the magnetic O-type star Theta 1 Ori C, with only a few photons detected below 0.8 keV. The emergent X-ray flux is 30% lower at the Halpha minimum phase.
(ii) Our modelling indicated that this seemingly hard spectrum is in fact a consequence of relatively soft intrinsic emission, similar to other magnetic Of?p stars, combined with a large amount of local absorption (~1-3 x 10^22 cm^-2). This combination is necessary to reproduce both the prominent Mg and Si spectral features, and the lack of flux at low energies. NGC 1624-2 is intrinsically luminous in X-rays (log LX emission ~ 33.4) but 70-95% of the X-ray emission produced by magnetically confined wind shocks is absorbed before it escapes the magnetosphere (log LX ISM corrected ~ 32.5).
(iii) The high X-ray luminosity, its variation with stellar rotation, and its large attenuation are all consistent with a large dynamical magnetosphere with magnetically confined wind shocks.
△ Less
Submitted 30 July, 2015;
originally announced July 2015.
-
The X-ray Lightcurve of the Supermassive star eta Carinae, 1996--2014
Authors:
M. F. Corcoran,
K. Hamaguchi,
J. K. Liburd,
D. Morris,
T. R. Gull,
T. I. Madura,
M. Teodoro,
A. F. J. Moffat,
N. D. Richardson,
C. M. P. Russell,
A. M. T. Pollock,
S. P. Owocki
Abstract:
Eta Carinae is the nearest example of a supermassive, superluminous, unstable star. Mass loss from the system is critical in shaping its circumstellar medium and in determining its ultimate fate. Eta Car currently loses mass via a dense, slow stellar wind and possesses one of the largest mass loss rates known. It is prone to episodes of extreme mass ejection via eruptions from some as-yet unspecif…
▽ More
Eta Carinae is the nearest example of a supermassive, superluminous, unstable star. Mass loss from the system is critical in shaping its circumstellar medium and in determining its ultimate fate. Eta Car currently loses mass via a dense, slow stellar wind and possesses one of the largest mass loss rates known. It is prone to episodes of extreme mass ejection via eruptions from some as-yet unspecified cause; the best examples of this are the large-scale eruptions which occurred in 19th century. Eta Car is a colliding wind binary in which strong variations in X-ray emission and in other wavebands are driven by the violent collision of the wind of eta Car-A and the fast, less dense wind of an otherwise hidden companion star. X-ray variations are the simplest diagnostic we have to study the wind-wind collision and allow us to measure the state of the stellar mass loss from both stars. We present the X-ray lightcurve over the last 20 years from ROSAT observations and monitoring with the Rossi X-ray Timing Explorer and the X-ray Telescope on the Swift satellite. We compare and contrast the behavior of the X-ray emission from the system over that timespan, including surprising variations during the 2014 X-ray minimum.
△ Less
Submitted 28 July, 2015;
originally announced July 2015.
-
A Coordinated X-ray and Optical Campaign on the Nearest Massive Eclipsing Binary, Delta Ori Aa: I. Overview of the X-ray Spectrum
Authors:
M. F. Corcoran,
J. S. Nichols,
H. Pablo,
T. Shenar,
A. M. T. Pollock,
W. L. Waldron,
A. F. J. Moffat,
N. D. Richardson,
C. M. P. Russell,
K. Hamaguchi,
D. P. Huenemoerder,
L. Oskinova,
W. -R. Hamann,
Y. Naze,
R. Ignace,
N. R. Evans,
J. R. Lomax,
J. L. Hoffman,
K. Gayley,
S. P. Owocki,
M. Leutenegger,
T. R. Gull,
K. T. Hole,
J. Lauer,
R. C. Iping
Abstract:
We present an overview of four phase-constrained Chandra HETGS X-ray observations of Delta Ori A. Delta Ori A is actually a triple system which includes the nearest massive eclipsing spectroscopic binary, Delta Ori Aa, the only such object which can be observed with little phase-smearing with the Chandra gratings. Since the fainter star, Delta Ori Aa2, has a much lower X-ray luminosity than the br…
▽ More
We present an overview of four phase-constrained Chandra HETGS X-ray observations of Delta Ori A. Delta Ori A is actually a triple system which includes the nearest massive eclipsing spectroscopic binary, Delta Ori Aa, the only such object which can be observed with little phase-smearing with the Chandra gratings. Since the fainter star, Delta Ori Aa2, has a much lower X-ray luminosity than the brighter primary, Delta Ori A provides a unique system with which to test the spatial distribution of the X-ray emitting gas around Delta Ori Aa1 via occultation by the photosphere of and wind cavity around the X-ray dark secondary. Here we discuss the X-ray spectrum and X-ray line profiles for the combined observation, having an exposure time of nearly 500 ksec and covering nearly the entire binary orbit. Companion papers discuss the X-ray variability seen in the Chandra spectra, present new space-based photometry and ground-based radial velocities simultaneous with the X-ray data to better constrain the system parameters, and model the effects of X-rays on the optical and UV spectrum. We find that the X-ray emission is dominated by embedded wind shock emission from star Aa1, with little contribution from the tertiary star Ab or the shocked gas produced by the collision of the wind of Aa1 against the surface of Aa2. We find a similar temperature distribution to previous X-ray spectrum analyses. We also show that the line half-widths are about $0.3-0.5\times$ the terminal velocity of the wind of star Aa1. We find a strong anti-correlation between line widths and the line excitation energy, which suggests that longer-wavelength, lower-temperature lines form farther out in the wind. Our analysis also indicates that the ratio of the intensities of the strong and weak lines of \ion{Fe}{17} and \ion{Ne}{10} are inconsistent with model predictions, which may be an effect of resonance scattering
△ Less
Submitted 22 July, 2015; v1 submitted 17 July, 2015;
originally announced July 2015.
-
Revisiting the Rigidly Rotating Magnetosphere model for $σ$ Ori E - II. Magnetic Doppler imaging, arbitrary field RRM, and light variability
Authors:
M. E. Oksala,
O. Kochukhov,
J. Krticka,
R. H. D. Townsend,
G. A. Wade,
M. Prvak,
Z. Mikulasek,
J. Silvester,
S. P. Owocki
Abstract:
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 d…
▽ More
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.
△ Less
Submitted 18 May, 2015;
originally announced May 2015.
-
Measuring the shock-heating rate in the winds of O stars using X-ray line spectra
Authors:
David H. Cohen,
Zequn Li,
Kenneth G. Gayley,
Stanley P. Owocki,
Jon O. Sundqvist,
Veronique Petit,
Maurice A. Leutenegger
Abstract:
We present a new method for using measured X-ray emission line fluxes from O stars to determine the shock-heating rate due to instabilities in their radiation-driven winds. The high densities of these winds means that their embedded shocks quickly cool by local radiative emission, while cooling by expansion should be negligible. Ignoring for simplicity any non-radiative mixing or conductive coolin…
▽ More
We present a new method for using measured X-ray emission line fluxes from O stars to determine the shock-heating rate due to instabilities in their radiation-driven winds. The high densities of these winds means that their embedded shocks quickly cool by local radiative emission, while cooling by expansion should be negligible. Ignoring for simplicity any non-radiative mixing or conductive cooling, the method presented here exploits the idea that the cooling post-shock plasma systematically passes through the temperature characteristic of distinct emission lines in the X-ray spectrum. In this way, the observed flux distribution among these X-ray lines can be used to construct the cumulative probability distribution of shock strengths that a typical wind parcel encounters as it advects through the wind. We apply this new method (Gayley 2014) to Chandra grating spectra from five O stars with X-ray emission indicative of embedded wind shocks in effectively single massive stars. Correcting for wind absorption of the X-ray line emission is a crucial component of our analysis, and we use wind optical depth values derived from X-ray line-profile fitting (Cohen et al. 2014) in order to make that correction. The shock-heating rate results we derive for all the stars are quite similar: the average wind mass element passes through roughly one shock that heats it to at least $10^6$ K as it advects through the wind, and the cumulative distribution of shock strengths is a strongly decreasing function of temperature, consistent with a negative power-law of index $n \approx 3$, implying a marginal distribution of shock strengths that scales as $T^{-4}$, and with hints of an even steeper decline or cut-off above $10^7$ K.
△ Less
Submitted 2 September, 2014;
originally announced September 2014.
-
Mass loss from inhomogeneous hot star winds III. An effective-opacity formalism for line radiative transfer in accelerating, clumped two-component media, and first results on theory and diagnostics
Authors:
J. O. Sundqvist,
J. Puls,
S. P. Owocki
Abstract:
[Abridged] We develop and benchmark a fast and easy-to-use effective-opacity formalism for line and continuum radiative transfer in an accelerating two-component clumpy medium. The formalism bridges the limits of optically thin and thick clumps, and is here used to i) design a simple vorosity-modified Sobolev with exact integration (vmSEI) method for analyzing UV wind resonance lines in hot, massi…
▽ More
[Abridged] We develop and benchmark a fast and easy-to-use effective-opacity formalism for line and continuum radiative transfer in an accelerating two-component clumpy medium. The formalism bridges the limits of optically thin and thick clumps, and is here used to i) design a simple vorosity-modified Sobolev with exact integration (vmSEI) method for analyzing UV wind resonance lines in hot, massive stars, and ii) derive simple correction factors to the line force driving the outflows of such stars. We show that (for a given ionization factor) UV resonance doublets may be used to analytically predict the upward corrections in empirically inferred mass-loss rates associated with porosity in velocity space (a.k.a. velocity-porosity, or vorosity), but that severe solution degeneracies exist. For an inter-clump density set to 1 % of the mean density, we for O and B supergiants derive upward empirical mass-loss corrections of typically factors of either ~5 or ~50, depending on which of the two applicable solutions is chosen. Overall, our results indicate this solution dichotomy severely limits the use of UV resonance lines as direct mass-loss indicators of clumped hot stellar winds. We next apply the effective-opacity formalism to the standard CAK theory of line-driven winds. By analytic and numerical hydrodynamics calculations, we show that in cases where vorosity is important at the critical point setting the mass-loss rate, the reduced line-force leads to a lower theoretical mass loss, by a factor scaling with the normalized velocity filling factor fvel. On the other hand, if vorosity is important only above this critical point, the predicted mass loss is not affected, but the wind terminal speed is reduced. This shows that porosity in velocity space can have a significant impact not only on the diagnostics, but also on the dynamics and theory of radiatively driven winds.
△ Less
Submitted 30 May, 2014;
originally announced May 2014.
-
3D Hydrodynamic & Radiative Transfer Models of X-ray Emission from Colliding Wind Binaries
Authors:
Christopher M. P. Russell,
Atsuo T. Okazaki,
Stanley P. Owocki,
Michael F. Corcoran,
Kenji Hamaguchi,
Yasuharu Sugawara
Abstract:
Colliding wind binaries (CWBs) are unique laboratories for X-ray astrophysics. The massive stars in these systems possess powerful stellar winds with speeds up to $\sim$3000 km s$^{-1}$, and their collision leads to hot plasma (up to $\sim10^8$K) that emit thermal X-rays (up to $\sim$10 keV). Many X-ray telescopes have observed CWBs, including Suzaku, and our work aims to model these X-ray observa…
▽ More
Colliding wind binaries (CWBs) are unique laboratories for X-ray astrophysics. The massive stars in these systems possess powerful stellar winds with speeds up to $\sim$3000 km s$^{-1}$, and their collision leads to hot plasma (up to $\sim10^8$K) that emit thermal X-rays (up to $\sim$10 keV). Many X-ray telescopes have observed CWBs, including Suzaku, and our work aims to model these X-ray observations. We use 3D smoothed particle hydrodynamics (SPH) to model the wind-wind interaction, and then perform 3D radiative transfer to compute the emergent X-ray flux, which is folded through X-ray telescopes' response functions to compare directly with observations. In these proceedings, we present our models of Suzaku observations of the multi-year-period, highly eccentric systems $η$ Carinae and WR 140. The models reproduce the observations well away from periastron passage, but only $η$ Carinae's X-ray spectrum is reproduced at periastron; the WR 140 model produces too much flux during this more complicated phase.
△ Less
Submitted 19 May, 2014;
originally announced May 2014.
-
Instabilities in the Envelopes and Winds of Very Massive Stars
Authors:
Stanley P. Owocki
Abstract:
The high luminosity of Very Massive Stars (VMS) means that radiative forces play an important, dynamical role both in the structure and stability of their stellar envelope, and in driving strong stellar-wind mass loss. Focusing on the interplay of radiative flux and opacity, with emphasis on key distinctions between continuum vs. line opacity, this chapter reviews instabilities in the envelopes an…
▽ More
The high luminosity of Very Massive Stars (VMS) means that radiative forces play an important, dynamical role both in the structure and stability of their stellar envelope, and in driving strong stellar-wind mass loss. Focusing on the interplay of radiative flux and opacity, with emphasis on key distinctions between continuum vs. line opacity, this chapter reviews instabilities in the envelopes and winds of VMS. Specifically, we discuss how: 1) the iron opacity bump can induce an extensive inflation of the stellar envelope; 2) the density dependence of mean opacity leads to strange mode instabilities in the outer envelope; 3) desaturation of line-opacity by acceleration of near-surface layers initiates and sustains a line-driven stellar wind outflow; 4) an associated line-deshadowing instability leads to extensive small-scale structure in the outer regions of such line-driven winds; 5) a star with super-Eddington luminosity can develop extensive atmospheric structure from photon bubble instabilities, or from stagnation of flow that exceeds the "photon tiring" limit; 6) the associated porosity leads to a reduction in opacity that can regulate the extreme mass loss of such continuum-driven winds. Two overall themes are the potential links of such instabilities to Luminous Blue Variable (LBV) stars, and the potential role of radiation forces in establishing the upper mass limit of VMS.
△ Less
Submitted 26 March, 2014;
originally announced March 2014.
-
Measuring mass-loss rates and constraining shock physics using X-ray line profiles of O stars from the Chandra archive
Authors:
David H. Cohen,
Emma E. Wollman,
Maurice A. Leutenegger,
Jon O. Sundqvist,
Alex W. Fullerton,
Janos Zsargo,
Stanley P. Owocki
Abstract:
We quantitatively investigate the extent of wind absorption signatures in the X-ray grating spectra of all non-magnetic, effectively single O stars in the Chandra archive via line profile fitting. Under the usual assumption of a spherically symmetric wind with embedded shocks, we confirm previous claims that some objects show little or no wind absorption. However, many other objects do show asymme…
▽ More
We quantitatively investigate the extent of wind absorption signatures in the X-ray grating spectra of all non-magnetic, effectively single O stars in the Chandra archive via line profile fitting. Under the usual assumption of a spherically symmetric wind with embedded shocks, we confirm previous claims that some objects show little or no wind absorption. However, many other objects do show asymmetric and blue shifted line profiles, indicative of wind absorption. For these stars, we are able to derive wind mass-loss rates from the ensemble of line profiles, and find values lower by an average factor of 3 than those predicted by current theoretical models, and consistent with H-alpha if clumping factors of f_cl ~ 20 are assumed. The same profile fitting indicates an onset radius of X-rays typically at r ~ 1.5 R_star, and terminal velocities for the X-ray emitting wind component that are consistent with that of the bulk wind. We explore the likelihood that the stars in the sample that do not show significant wind absorption signatures in their line profiles have at least some X-ray emission that arises from colliding wind shocks with a close binary companion. The one clear exception is zeta Oph, a weak-wind star that appears to simply have a very low mass-loss rate. We also reanalyse the results from the canonical O supergiant zeta Pup, using a solar-metallicity wind opacity model and find Mdot = 1.8 \times 10^{-6} M_sun/yr, consistent with recent multi-wavelength determinations.
△ Less
Submitted 30 January, 2014;
originally announced January 2014.
-
Constraints on decreases in Eta Carinae's mass loss from 3D hydrodynamic simulations of its binary colliding winds
Authors:
Thomas I. Madura,
Theodore R. Gull,
Atsuo T. Okazaki,
Christopher M. P. Russell,
Stanley P. Owocki,
Jose H. Groh,
Michael F. Corcoran,
Kenji Hamaguchi,
Mairan Teodoro
Abstract:
Recent work suggests that the mass-loss rate of the primary star (Eta A) in the massive colliding wind binary Eta Carinae dropped by a factor of 2-3 between 1999 and 2010. We present results from large- (r=1545au) and small- (r=155au) domain, 3D smoothed particle hydrodynamic (SPH) simulations of Eta Car's colliding winds for 3 Eta A mass-loss rates (2.4, 4.8, and 8.5 x 10^-4 M_sun/yr), investigat…
▽ More
Recent work suggests that the mass-loss rate of the primary star (Eta A) in the massive colliding wind binary Eta Carinae dropped by a factor of 2-3 between 1999 and 2010. We present results from large- (r=1545au) and small- (r=155au) domain, 3D smoothed particle hydrodynamic (SPH) simulations of Eta Car's colliding winds for 3 Eta A mass-loss rates (2.4, 4.8, and 8.5 x 10^-4 M_sun/yr), investigating the effects on the dynamics of the binary wind-wind collision (WWC). These simulations include orbital motion, optically thin radiative cooling, and radiative forces. We find that Eta A's mass-loss rate greatly affects the time-dependent hydrodynamics at all spatial scales investigated. The simulations also show that the post-shock wind of the companion star (Eta B) switches from the adiabatic to the radiative-cooling regime during periastron passage. The SPH simulations together with 1D radiative transfer models of Eta A's spectra reveal that a factor of 2 or more drop in Eta A's mass-loss rate should lead to substantial changes in numerous multiwavelength observables. Recent observations are not fully consistent with the model predictions, indicating that any drop in Eta A's mass-loss rate was likely by a factor < 2 and occurred after 2004. We speculate that most of the recent observed changes in Eta Car are due to a small increase in the WWC opening angle that produces significant effects because our line-of-sight to the system lies close to the dense walls of the WWC zone. A modest decrease in Eta A's mass-loss rate may be responsible, but changes in the wind/stellar parameters of Eta B cannot yet be fully ruled out. We suggest observations during Eta Car's next periastron in 2014 to further test for decreases in Eta A's mass-loss rate. If Eta A's mass-loss rate is declining and continues to do so, the 2014 X-ray minimum should be even shorter than that of 2009.
△ Less
Submitted 1 October, 2013;
originally announced October 2013.
-
Constraints on porosity and mass loss in O-star winds from modeling of X-ray emission line profile shapes
Authors:
Maurice A. Leutenegger,
David H. Cohen,
Jon O. Sundqvist,
Stanley P. Owocki
Abstract:
We fit X-ray emission line profiles in high resolution XMM-Newton and Chandra grating spectra of the early O supergiant Zeta Pup with models that include the effects of porosity in the stellar wind. We explore the effects of porosity due to both spherical and flattened clumps. We find that porosity models with flattened clumps oriented parallel to the photosphere provide poor fits to observed line…
▽ More
We fit X-ray emission line profiles in high resolution XMM-Newton and Chandra grating spectra of the early O supergiant Zeta Pup with models that include the effects of porosity in the stellar wind. We explore the effects of porosity due to both spherical and flattened clumps. We find that porosity models with flattened clumps oriented parallel to the photosphere provide poor fits to observed line shapes. However, porosity models with isotropic clumps can provide acceptable fits to observed line shapes, but only if the porosity effect is moderate. We quantify the degeneracy between porosity effects from isotropic clumps and the mass-loss rate inferred from the X-ray line shapes, and we show that only modest increases in the mass-loss rate (<~ 40%) are allowed if moderate porosity effects (h_infinity <~ R_*) are assumed to be important. Large porosity lengths, and thus strong porosity effects, are ruled out regardless of assumptions about clump shape. Thus, X-ray mass-loss rate estimates are relatively insensitive to both optically thin and optically thick clumping. This supports the use of X-ray spectroscopy as a mass-loss rate calibration for bright, nearby O stars.
△ Less
Submitted 23 May, 2013;
originally announced May 2013.
-
On magnetic inhibition of photospheric macro-turbulence generated in the iron-bump opacity zone of O-stars
Authors:
J. O. Sundqvist,
V. Petit,
S. P. Owocki,
G. A. Wade,
J. Puls,
the MiMeS Collaboration
Abstract:
Massive, hot OB-stars show clear evidence of strong macroscopic broadening (in addition to rotation) in their photospheric spectral lines. This paper examines the occurrence of such "macro-turbulence" in slowly rotating O-stars with strong, organised surface magnetic fields. Focusing on the CIV 5811A line, we find evidence for significant macro-turbulent broadening in all stars except NGC1624-2, w…
▽ More
Massive, hot OB-stars show clear evidence of strong macroscopic broadening (in addition to rotation) in their photospheric spectral lines. This paper examines the occurrence of such "macro-turbulence" in slowly rotating O-stars with strong, organised surface magnetic fields. Focusing on the CIV 5811A line, we find evidence for significant macro-turbulent broadening in all stars except NGC1624-2, which also has (by far) the strongest magnetic field. Instead, the very sharp CIV lines in NGC1624-2 are dominated by magnetic Zeeman broadening, from which we estimate a dipolar field of approximately 20 kG. By contrast, magnetic broadening is negligible in the other stars (due to their weaker field strengths, on order 1 kG), and their CIV profiles are typically very broad and similar to corresponding lines observed in non-magnetic O-stars. Quantifying this by an isotropic, Gaussian macro-turbulence, we derive vmac = 2.2 (+- 0.9/2.2) km/s for NGC1624-2, and vmac = 20-65 km/s for the rest of the magnetic sample. We use these observational results to test the hypothesis that the field can stabilise the atmosphere and suppress the generation of macro-turbulence down to stellar layers where the magnetic pressure PB and the gas pressure Pg are comparable. Using a simple grey atmosphere to estimate the temperature T0 at which PB = Pg, we find that T0 > Teff for all investigated magnetic stars, but that T0 reaches the ~160000 K layers associated with the iron opacity-bump in hot stars only for NGC1624-2. This is consistent with the view that the responsible physical mechanism for photospheric O-star macro-turbulence may be stellar gravity-mode oscillations excited by sub-surface convection zones, and suggests that a sufficiently strong magnetic field can suppress such iron-bump generated convection and associated pulsational excitation.
△ Less
Submitted 23 May, 2013;
originally announced May 2013.
-
Thin-shell mixing in radiative wind-shocks and the Lx-Lbol scaling of O-star X-rays
Authors:
Stanley P. Owocki,
Jon O. Sundqvist,
David H. Cohen,
Kenneth G. Gayley
Abstract:
X-ray satellites since Einstein have empirically established that the X-ray luminosity from single O-stars scales linearly with bolometric luminosity, Lx ~ 10^{-7} Lbol. But straightforward forms of the most favored model, in which X-rays arise from instability-generated shocks embedded in the stellar wind, predict a steeper scaling, either with mass loss rate Lx ~ Mdot ~ Lbol^{1.7} if the shocks…
▽ More
X-ray satellites since Einstein have empirically established that the X-ray luminosity from single O-stars scales linearly with bolometric luminosity, Lx ~ 10^{-7} Lbol. But straightforward forms of the most favored model, in which X-rays arise from instability-generated shocks embedded in the stellar wind, predict a steeper scaling, either with mass loss rate Lx ~ Mdot ~ Lbol^{1.7} if the shocks are radiative, or with Lx ~ Mdot^{2} ~ Lbol^{3.4} if they are adiabatic. This paper presents a generalized formalism that bridges these radiative vs. adiabatic limits in terms of the ratio of the shock cooling length to the local radius. Noting that the thin-shell instability of radiative shocks should lead to extensive mixing of hot and cool material, we propose that the associated softening and weakening of the X-ray emission can be parametrized as scaling with the cooling length ratio raised to a power m$, the "mixing exponent". For physically reasonable values m ~= 0.4, this leads to an X-ray luminosity Lx ~ Mdot^{0.6} ~ Lbol that matches the empirical scaling. To fit observed X-ray line profiles, we find such radiative-shock-mixing models require the number of shocks to drop sharply above the initial shock onset radius. This in turn implies that the X-ray luminosity should saturate and even decrease for optically thick winds with very high mass-loss rates. In the opposite limit of adiabatic shocks in low-density winds (e.g., from B-stars), the X-ray luminosity should drop steeply with Mdot^2. Future numerical simulation studies will be needed to test the general thin-shell mixing ansatz for X-ray emission.
△ Less
Submitted 18 December, 2012;
originally announced December 2012.
-
A Magnetic Confinement vs. Rotation Classification of Massive-Star Magnetospheres
Authors:
V. Petit,
S. P. Owocki,
G. A. Wade,
D. H. Cohen,
J. O. Sundqvist,
M. Gagné,
J. Maíz Apellániz,
M. E. Oksala,
D. A. Bohlender,
Th. Rivinius,
H. F. Henrichs,
E. Alecian,
R. H. D. Townsend,
A. ud-Doula,
the MiMeS Collaboration
Abstract:
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én 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 compi…
▽ More
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én 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]
△ Less
Submitted 1 November, 2012;
originally announced November 2012.