-
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.
-
Winds and magnetospheres of stars and planets: similarities and differences
Authors:
Stan Owocki
Abstract:
Both stars and planets can lose mass through an expansive wind outflow, often constrained or channeled by magnetic fields that form a surrounding magnetosphere. The very strong winds of massive stars are understood to be driven by line-scattering of the star's radiative momentum, while in the Sun and even lower-mass stars a much weaker mass loss arises from the thermal expansion of a mechanically…
▽ More
Both stars and planets can lose mass through an expansive wind outflow, often constrained or channeled by magnetic fields that form a surrounding magnetosphere. The very strong winds of massive stars are understood to be driven by line-scattering of the star's radiative momentum, while in the Sun and even lower-mass stars a much weaker mass loss arises from the thermal expansion of a mechanically heated corona. In exoplanets around such low-mass stars, the radiative heating and wind interaction can lead to thermal expansion or mechanical ablation of their atmospheres. Stellar magnetospheres result from the internal trapping of the wind outflow, while planetary magnetospheres are typically shaped by the external impact from the star's wind. But in both cases the stressing can drive magnetic reconnection that results in observable signatures such as X-ray flares and radio outbursts. This review will aim to give an overview of the underlying physics of these processes with emphasis on their similarities and distinctions for stars vs. planets.
△ Less
Submitted 31 October, 2022;
originally announced October 2022.
-
Magnetically confined wind shock
Authors:
Asif ud-Doula,
Stan Owocki
Abstract:
Many stars across all classes possess strong enough magnetic fields to influence dynamical flow of material off the stellar surface. For the case of massive stars (O and B types), about 10\% of them harbour strong, globally ordered (mostly dipolar) magnetic fields. The trapping and channeling of their stellar winds in closed magnetic loops leads to {\it magnetically confined wind shocks} (MCWS), w…
▽ More
Many stars across all classes possess strong enough magnetic fields to influence dynamical flow of material off the stellar surface. For the case of massive stars (O and B types), about 10\% of them harbour strong, globally ordered (mostly dipolar) magnetic fields. The trapping and channeling of their stellar winds in closed magnetic loops leads to {\it magnetically confined wind shocks} (MCWS), with pre-shock flow speeds that are some fraction of the wind terminal speed that can be a few thousand km s$^{-1}$. These shocks generate hot plasma, a source of X-rays. In the last decade, several developments took place, notably the determination of the hot plasma properties for a large sample of objects using \xmm\ and \ch, as well as fully self-consistent MHD modelling and the identification of shock retreat effects in weak winds. In addition, these objects are often sources of H$α$ emission which is controlled by either sufficiently high mass loss rate or centrifugal breakout. Here we review the theoretical aspects of such magnetic massive star wind dynamics.
△ Less
Submitted 18 September, 2022;
originally announced September 2022.
-
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.
-
Getting started: How a supersonic stellar wind is initiated from a hydrostatic surface
Authors:
Stan Owocki
Abstract:
Most of a star's mass is bound in a hydrostatic equilibrium in which pressure balances gravity. But if at some near-surface layer additional outward forces overcome gravity, this can transition to a supersonic, outflowing wind, with the sonic point, where the outward force cancels gravity, marking the division between hydrostatic atmosphere and wind outflow. This talk will review general issues wi…
▽ More
Most of a star's mass is bound in a hydrostatic equilibrium in which pressure balances gravity. But if at some near-surface layer additional outward forces overcome gravity, this can transition to a supersonic, outflowing wind, with the sonic point, where the outward force cancels gravity, marking the division between hydrostatic atmosphere and wind outflow. This talk will review general issues with such transonic initiation of a stellar wind outflow, and how this helps set the wind mass loss rate. The main discussion contrasts the flow initiation in four prominent classes of steady-state winds: (1) the pressure-driven coronal wind of the sun and other cool stars; (2) line-driven winds from OB stars; (3) a two-stage initiation model for the much denser winds from Wolf-Rayet (WR) stars; and (4) the slow "overflow" mass loss from highly evolved giant stars. A follow on discussion briefly reviews eruptive mass loss, with particular focus on the giant eruption of eta Carinae.
△ 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 Magnetic Early B-type Stars IV: Breakout or Leakage? H$α$ emission as a diagnostic of plasma transport in centrifugal magnetospheres
Authors:
M. E. Shultz,
S. Owocki,
Th. Rivinius,
G. A. Wade,
C. Neiner,
E. Alecian,
O. Kochukhov,
D. Bohlender,
A. ud-Doula,
J. D. Landstreet,
J. Sikora,
A. David-Uraz,
V. Petit,
P. Cerrahoğlu,
R. Fine,
G. Henson,
the MiMeS,
BinaMIcS Collaborations
Abstract:
Rapidly rotating early-type stars with strong magnetic fields frequently show H$α$ emission originating in Centrifugal Magnetospheres (CMs), circumstellar structures in which centrifugal support due to magnetically enforced corotation of the magnetically confined plasma enables it to accumulate to high densities. It is not currently known whether the CM plasma escapes via Centrifugal Breakout (CB)…
▽ More
Rapidly rotating early-type stars with strong magnetic fields frequently show H$α$ emission originating in Centrifugal Magnetospheres (CMs), circumstellar structures in which centrifugal support due to magnetically enforced corotation of the magnetically confined plasma enables it to accumulate to high densities. It is not currently known whether the CM plasma escapes via Centrifugal Breakout (CB), or by an unidentified leakage mechanism. We have conducted the first comprehensive examination of the H$α$ emission properties of all stars currently known to display CM-pattern emission. We find that the onset of emission is dependent primarily on the area of the CM, which can be predicted simply by the value $B_{\rm K}$ of the magnetic field at the Kepler corotation radius $R_{\rm K}$. Emission strength is strongly sensitive to both CM area and $B_{\rm K}$. Emission onset and strength are {\em not} dependent on effective temperature, luminosity, or mass-loss rate. These results all favour a CB scenario, however the lack of intrinsic variability in any CM diagnostics indicates that CB must be an essentially continuous process, i.e.\ it effectively acts as a leakage mechanism. We also show that the emission profile shapes are approximately scale-invariant, i.e.\ they are broadly similar across a wide range of emission strengths and stellar parameters. While the radius of maximum emission correlates closely as expected to $R_{\rm K}$, it is always larger, contradicting models that predict that emission should peak at $R_{\rm K}$.
△ 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.
-
PolStar -- An Explorer-Class FUV Spectropolarimetry Mission to Map the Environments of Massive Stars
Authors:
Paul Scowen,
Richard Ignace,
Coralie Neiner,
Gregg Wade,
Matt Beasley,
Jon Bjorkman,
Jean-Claude Bouret,
Roberto Casini,
Tanausu del Pino Alemán,
Samantha Edgington,
Ken Gayley,
Ed Guinan,
Jennifer Hoffman,
Ian Howarth,
Tony Hull,
Rafael Manso Sainz,
Yael Naze,
Alison Nordt,
Stan Owocki,
Steve Petrinec,
Raman Prinja,
Hugues Sana,
Matt Shultz,
William Sparks,
Nicole St-Louis
, et al. (4 additional authors not shown)
Abstract:
PolStar is an Explorer-class far ultraviolet (FUV) spectropolarimetry mission designed to target massive stars and their environments. PolStar will take advantage of resonance lines only available in the FUV to measure for the first time the magnetic and wind environment around massive stars to constrain models of rotation and mass loss.
PolStar is an Explorer-class far ultraviolet (FUV) spectropolarimetry mission designed to target massive stars and their environments. PolStar will take advantage of resonance lines only available in the FUV to measure for the first time the magnetic and wind environment around massive stars to constrain models of rotation and mass loss.
△ Less
Submitted 15 November, 2019;
originally announced November 2019.
-
The Magnetic Early B-type Stars III: A main sequence magnetic, rotational, and magnetospheric biography
Authors:
M. E. Shultz,
G. A. Wade,
Th. Rivinius,
E. Alecian,
C. Neiner,
V. Petit,
S. Owocki,
A. ud-Doula,
O. Kochukhov,
D. Bohlender,
Z. Keszthelyi,
the MiMeS,
BinaMIcS Collaborations
Abstract:
Magnetic confinement of stellar winds leads to the formation of magnetospheres, which can be sculpted into Centrifugal Magnetospheres (CMs) by rotational support of the corotating plasma. The conditions required for the CMs of magnetic early B-type stars to yield detectable emission in H$α$ -- the principal diagnostic of these structures -- are poorly constrained. A key reason is that no detailed…
▽ More
Magnetic confinement of stellar winds leads to the formation of magnetospheres, which can be sculpted into Centrifugal Magnetospheres (CMs) by rotational support of the corotating plasma. The conditions required for the CMs of magnetic early B-type stars to yield detectable emission in H$α$ -- the principal diagnostic of these structures -- are poorly constrained. A key reason is that no detailed study of the magnetic and rotational evolution of this population has yet been performed. Using newly determined rotational periods, modern magnetic measurements, and atmospheric parameters determined via spectroscopic modelling, we have derived fundamental parameters, dipolar oblique rotator models, and magnetospheric parameters for 56 early B-type stars. Comparison to magnetic A- and O-type stars shows that the range of surface magnetic field strength is essentially constant with stellar mass, but that the unsigned surface magnetic flux increases with mass. Both the surface magnetic dipole strength and the total magnetic flux decrease with stellar age, with the rate of flux decay apparently increasing with stellar mass. We find tentative evidence that multipolar magnetic fields may decay more rapidly than dipoles. Rotational periods increase with stellar age, as expected for a magnetic braking scenario. Without exception, all stars with H$α$ emission originating in a CM are 1) rapid rotators, 2) strongly magnetic, and 3) young, with the latter property consistent with the observation that magnetic fields and rotation both decrease over time.
△ Less
Submitted 5 September, 2019;
originally announced September 2019.
-
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.
-
Line-driven ablation of circumstellar discs: III. Accounting for and analyzing the effects of continuum optical depth
Authors:
Nathaniel Kee,
Stan Owocki,
Rolf Kuiper
Abstract:
In studying the formation of massive stars, it is essential to consider the strong radiative feedback on the stars' natal environments from their high luminosities ($10^4 \sim 10^6 L_\odot$). Given that massive stars contract to main-sequence-like radii before accretion finishes, one form this feedback takes is UV line-acceleration, resulting in outflows much like those expected from main-sequence…
▽ More
In studying the formation of massive stars, it is essential to consider the strong radiative feedback on the stars' natal environments from their high luminosities ($10^4 \sim 10^6 L_\odot$). Given that massive stars contract to main-sequence-like radii before accretion finishes, one form this feedback takes is UV line-acceleration, resulting in outflows much like those expected from main-sequence massive stars. As shown by the prior papers in this series, in addition to driving stellar winds, such line forces also ablate the surface layers off of circumstellar discs within a few stellar radii of the stellar photosphere. This removal of material from an accretion disc in turn results in a decreased accretion rate onto the forming star. Quantifying this, however, requires accounting for the continuum optical depth of the disc along the non-radial rays required for the three-dimensional line-acceleration prescription used in this paper series. We introduce the "thin disc approximation", allowing these continuum optical depths arising from an optically thick but geometrically thin disc to be dynamically treated in the context of radiation-hydrodynamics simulations. Using this approximation in full dynamical simulations, we show that such continuum optical depth effects only reduce the disc ablation by 30 percent or less relative to previous simulations that ignored continuum absorption.
△ Less
Submitted 22 June, 2018;
originally announced June 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.
-
Disruption of circumstellar discs by large-scale stellar magnetic fields
Authors:
Asif ud-Doula,
Stanley Owocki,
Nathaniel Kee
Abstract:
Spectropolarimetric surveys reveal that 8-10\% of OBA stars harbor large-scale magnetic fields, but thus far no such fields have been detected in any classical Be stars. Motivated by this, we present here MHD simulations for how a pre-existing Keplerian disc -- like that inferred to form from decretion of material from rapidly rotating Be stars -- can be disrupted by a rotation-aligned stellar dip…
▽ More
Spectropolarimetric surveys reveal that 8-10\% of OBA stars harbor large-scale magnetic fields, but thus far no such fields have been detected in any classical Be stars. Motivated by this, we present here MHD simulations for how a pre-existing Keplerian disc -- like that inferred to form from decretion of material from rapidly rotating Be stars -- can be disrupted by a rotation-aligned stellar dipole field. For characteristic stellar and disc parameters of a near-critically rotating B2e star, we find that a polar surface field strength of just 10 G can significantly disrupt the disc, while a field of 100 G, near the observational upper limit inferred for most Be stars, completely destroys the disc over just a few days. Our parameter study shows that the efficacy of this magnetic disruption of a disc scales with the characteristic plasma beta {(defined as the ratio between thermal and magnetic pressure)} in the disc, but is surprisingly insensitive to other variations, e.g. in stellar rotation speed, or the mass loss rate of the star's radiatively driven wind. The disc disruption seen here for even a modest field strength suggests that the presumed formation of such Be discs by decretion of material from the star would likely be strongly inhibited by such fields; this provides an attractive explanation for why no large-scale fields are detected from such Be stars.
△ Less
Submitted 8 May, 2018;
originally announced May 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.
-
Rotationally modulated photometric variations in B supergiants?
Authors:
Alexandre David-Uraz,
Gregg Wade,
Anthony Moffat,
Stan Owocki,
Véronique Petit,
the BRITE team
Abstract:
In this contribution, we present BRITE observations of the early-B supergiants $ε$ Ori and $κ$ Ori. We perform a preliminary analysis of the data acquired over the first two Orion observing runs. We evaluate whether they are compatible with co-rotating bright spots and discuss the challenges of such an approach.
In this contribution, we present BRITE observations of the early-B supergiants $ε$ Ori and $κ$ Ori. We perform a preliminary analysis of the data acquired over the first two Orion observing runs. We evaluate whether they are compatible with co-rotating bright spots and discuss the challenges of such an approach.
△ Less
Submitted 22 November, 2017;
originally announced November 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.
-
Line-driven ablation of circumstellar discs: II. Analyzing the role of multiple resonances
Authors:
Nathaniel Kee,
Stan Owocki,
Rolf Kuiper
Abstract:
We extend our previous study of radiative ablation of circumstellar disks by line-scattering of the star's radiation, accounting now for the effect of multiple line resonances off the stellar limb. For an analytic, three-dimensional model of the velocity structure of an equatorial Keplerian disk bounded at higher latitudes by a radially accelerating stellar wind outflow, we use root-finding method…
▽ More
We extend our previous study of radiative ablation of circumstellar disks by line-scattering of the star's radiation, accounting now for the effect of multiple line resonances off the stellar limb. For an analytic, three-dimensional model of the velocity structure of an equatorial Keplerian disk bounded at higher latitudes by a radially accelerating stellar wind outflow, we use root-finding methods to identify multiple resonances from a near-disk circumstellar location along starward rays both on and off the stellar core. Compared to our previous study that accounted only for the effect of on-core resonances in reducing the radiative driving through the scattering of radiation away from a near-disk circumstellar location, including off-limb resonances leads to additional radiative driving from scattering toward this location. Instead of the up-to-50% reduction in line-acceleration previously inferred from multiple resonance effects, we now find a more modest 15-20% net reduction.
△ Less
Submitted 9 November, 2017; v1 submitted 25 October, 2017;
originally announced October 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.
-
Eta Carinae's Thermal X-ray Tail Measured with XMM-Newton and NuSTAR
Authors:
Kenji Hamaguchi,
Michael F. Corcoran,
Theodore R. Gull,
Hiromitsu Takahashi,
Brian W. Grefenstette,
Takayuki Yuasa,
Martin Stuhlinger,
Christopher M. P. Russell,
Anthony F. J. Moffat,
Neetika Sharma,
Thomas I. Madura,
Noel D. Richardson,
Jose Groh,
Julian M. Pittard,
Stanley Owocki
Abstract:
The evolved, massive highly eccentric binary system, eta Carinae, underwent a periastron passage in the summer of 2014. We obtained two coordinated X-ray observations with XMM-Newton and NuSTAR during the elevated X-ray flux state and just before the X-ray minimum flux state around this passage. These NuSTAR observations clearly detected X-ray emission associated with eta Car extending up to ~50 k…
▽ More
The evolved, massive highly eccentric binary system, eta Carinae, underwent a periastron passage in the summer of 2014. We obtained two coordinated X-ray observations with XMM-Newton and NuSTAR during the elevated X-ray flux state and just before the X-ray minimum flux state around this passage. These NuSTAR observations clearly detected X-ray emission associated with eta Car extending up to ~50 keV for the first time. The NuSTAR spectrum above 10 keV can be fit with the bremsstrahlung tail from a kT ~6 keV plasma. This temperature is Delta kT ~2 keV higher than those measured from the iron K emission line complex, if the shocked gas is in collisional ionization equilibrium. This result may suggest that the companion star's pre-shock wind velocity is underestimated. The NuSTAR observation near the X-ray minimum state showed a gradual decline in the X-ray emission by 40% at energies above 5 keV in a day, the largest rate of change of the X-ray flux yet observed in individual eta Car observations. The column density to the hardest emission component, NH ~1e24 cm-2, marked one of the highest values ever observed for eta Car, strongly suggesting increased obscuration of the wind-wind colliding X-ray emission by the thick primary stellar wind prior to superior conjunction. Neither observation detected the power-law component in the extremely hard band that INTEGRAL and Suzaku observed prior to 2011. If the non-detection by NuSTAR is caused by absorption, the power-law source must be small and located very near the WWC apex. Alternatively, it may be that the power-law source is not related to either eta Car or the GeV gamma-ray source.
△ Less
Submitted 2 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.