Source author record

S. P. Owocki

S. P. Owocki appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

28works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2025arXiv

The effect of collisional cooling of energetic electrons on radio emission from the centrifugal magnetospheres of magnetic hot stars

This paper extends our previous study of gyro-emission by energetic electrons in the magnetospheres of rapidly rotating, magnetic massive stars, through a quantitative analysis of the role of Coulomb collisions with thermal electrons from stellar wind material trapped within the centrifugal magnetosphere (CM). For a dipolar field with aligned magnetic and rotational axes, we show that both gyro-cooling along magnetic loops and Coulomb cooling in the CM layer have nearly the same dependence on the magnitude and radial variation of magnetic field, implying that their ratio is a global parameter that is largely independent of the field. Analytic analysis shows that, for electrons introduced near the CM layer around a magnetic loop apex, collisional cooling is more important for electrons with high pitch angle, while more field-aligned electrons cool by gyro-emission near their mirror point close to the loop base. Numerical models that assume a gyrotropic initial deposition with a gaussian distribution in both radius and loop co-latitude show the residual gyro-emission is generally strongest near the loop base, with highly relativistic electrons suffering much lower collisional losses than lower-energy electrons that are only mildly relativistic. Finally, we briefly discuss the potential applicability of this formalism to magnetic ultracool dwarfs, for which VLBI observations indicate incoherent radio emission to be concentrated around the magnetic equator, in contrast to our predictions here for magnetic hot stars. We suggest that this difference could be attributed to either a lower ambient density of thermal electrons, or more highly relativistic non-thermal electrons, both of which would reduce the relative importance of the collisional cooling explored here.

preprint2022arXiv

Discovery and origin of the radio emission from the multiple stellar system KQVel

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.

preprint2022arXiv

Electron Scattering Emission in the Light Curves of Stars with Centrifugal Magnetospheres

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.

preprint2022arXiv

MOBSTER -- VI. The crucial influence of rotation on the radio magnetospheres of hot stars

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.

preprint2020arXiv

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

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.

preprint2016arXiv

Accretion, Outflows, and Winds of Magnetized Stars

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.

preprint2016arXiv

Line-driven ablation of circumstellar disks: I. Optically thin decretion disks of classical Oe/Be stars

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.

preprint2015arXiv

Revisiting the Rigidly Rotating Magnetosphere model for $σ$ Ori E - II. Magnetic Doppler imaging, arbitrary field RRM, and light variability

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

preprint2015arXiv

The X-ray Lightcurve of the Supermassive star eta Carinae, 1996--2014

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.

preprint2014arXiv

Mass loss from inhomogeneous hot star winds III. An effective-opacity formalism for line radiative transfer in accelerating, clumped two-component media, and first results on theory and diagnostics

[Abridged] We develop and benchmark a fast and easy-to-use effective-opacity formalism for line and continuum radiative transfer in an accelerating two-component clumpy medium. The formalism bridges the limits of optically thin and thick clumps, and is here used to i) design a simple vorosity-modified Sobolev with exact integration (vmSEI) method for analyzing UV wind resonance lines in hot, massive stars, and ii) derive simple correction factors to the line force driving the outflows of such stars. We show that (for a given ionization factor) UV resonance doublets may be used to analytically predict the upward corrections in empirically inferred mass-loss rates associated with porosity in velocity space (a.k.a. velocity-porosity, or vorosity), but that severe solution degeneracies exist. For an inter-clump density set to 1 % of the mean density, we for O and B supergiants derive upward empirical mass-loss corrections of typically factors of either ~5 or ~50, depending on which of the two applicable solutions is chosen. Overall, our results indicate this solution dichotomy severely limits the use of UV resonance lines as direct mass-loss indicators of clumped hot stellar winds. We next apply the effective-opacity formalism to the standard CAK theory of line-driven winds. By analytic and numerical hydrodynamics calculations, we show that in cases where vorosity is important at the critical point setting the mass-loss rate, the reduced line-force leads to a lower theoretical mass loss, by a factor scaling with the normalized velocity filling factor fvel. On the other hand, if vorosity is important only above this critical point, the predicted mass loss is not affected, but the wind terminal speed is reduced. This shows that porosity in velocity space can have a significant impact not only on the diagnostics, but also on the dynamics and theory of radiatively driven winds.

preprint2013arXiv

On magnetic inhibition of photospheric macro-turbulence generated in the iron-bump opacity zone of O-stars

Massive, hot OB-stars show clear evidence of strong macroscopic broadening (in addition to rotation) in their photospheric spectral lines. This paper examines the occurrence of such "macro-turbulence" in slowly rotating O-stars with strong, organised surface magnetic fields. Focusing on the CIV 5811A line, we find evidence for significant macro-turbulent broadening in all stars except NGC1624-2, which also has (by far) the strongest magnetic field. Instead, the very sharp CIV lines in NGC1624-2 are dominated by magnetic Zeeman broadening, from which we estimate a dipolar field of approximately 20 kG. By contrast, magnetic broadening is negligible in the other stars (due to their weaker field strengths, on order 1 kG), and their CIV profiles are typically very broad and similar to corresponding lines observed in non-magnetic O-stars. Quantifying this by an isotropic, Gaussian macro-turbulence, we derive vmac = 2.2 (+- 0.9/2.2) km/s for NGC1624-2, and vmac = 20-65 km/s for the rest of the magnetic sample. We use these observational results to test the hypothesis that the field can stabilise the atmosphere and suppress the generation of macro-turbulence down to stellar layers where the magnetic pressure PB and the gas pressure Pg are comparable. Using a simple grey atmosphere to estimate the temperature T0 at which PB = Pg, we find that T0 > Teff for all investigated magnetic stars, but that T0 reaches the ~160000 K layers associated with the iron opacity-bump in hot stars only for NGC1624-2. This is consistent with the view that the responsible physical mechanism for photospheric O-star macro-turbulence may be stellar gravity-mode oscillations excited by sub-surface convection zones, and suggests that a sufficiently strong magnetic field can suppress such iron-bump generated convection and associated pulsational excitation.

preprint2012arXiv

A Magnetic Confinement vs. Rotation Classification of Massive-Star Magnetospheres

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

preprint2012arXiv

Diagnosing small- and large-scale structure in the winds of hot, massive OB-star

It is observationally as well as theoretically well established that the winds of hot, massive OB-stars are highly structured on a broad range of spatial scales. This paper first discusses consequences of the small-scale structures associated with the strong instability inherent to the line-driving of these winds. We demonstrate the importance of a proper treatment of such wind clumping to obtain reliable estimates of mass-loss rates, and also show that instability simulations that are perturbed at the lower boundary indeed display significant clumping quite close to the wind base, in general agreement with observations. But a growing subset of massive stars has also been found to possess strong surface magnetic fields, which may channel the star's outflow and induce also large-scale wind structures and cyclic behavior of spectral diagnostics. The paper concludes by showing that multi-dimensional, magneto-hydrodynamical wind simulations, together with detailed radiative-transfer modeling, can reproduce remarkably well the periodic Balmer line emission observed in slowly rotating magnetic O stars like HD191612.

preprint2012arXiv

Modeling high-energy light curves of the PSR B1259-63/LS 2883 binary based on 3-D SPH simulations

Temporal changes of X-ray to very-high-energy gamma-ray emissions from the pulsar-Be star binary PSR B1259-63/LS 2883 are studied based on 3-D SPH simulations of pulsar wind interaction with Be-disk and wind. We focus on the periastron passage of the binary and calculate the variation of the synchrotron and inverse-Compton emissions using the simulated shock geometry and pressure distribution of the pulsar wind. The characteristic double-peaked X-ray light curve from observations is reproduced by our simulation under a dense Be disk condition (base density ~10^{-9} g cm^{-3}). We interpret the pre- and post-periastron peaks as being due to a significant increase in the conversion efficiency from pulsar spin down power to the shock-accelerated particle energy at orbital phases when the pulsar crosses the disk before periastron passage, and when the pulsar wind creates a cavity in the disk gas after periastron passage, respectively. On the contrary, in the model TeV light curve, which also shows a double peak feature, the first peak appears around the periastron phase. The possible effects of cooling processes on the TeV light curve are briefly discussed.

preprint2011arXiv

Magnetospheres of massive stars across the EM spectrum

Magnetic massive stars -- which are being discovered with increasing frequency -- represent a new category of wind-shaping mechanism for O and B stars. Magnetic channeling of these stars' radiation-driven winds, the Magnetically Confined Wind Shock paradigm, leads to the formation of a shock-heated magnetosphere, which can radiate X-rays, modify UV resonance lines, and create disks of Halpha emitting material. The dynamical properties of these magnetospheres are well understood from a theoretical point of view as an interplay between the magnetic wind confinement and rotation. However, the manifestations of magnetospheres across the spectrum may be more complex and various than first anticipated. On the other hand, recent advances in modeling these magnetospheres provide a key to better understand massive star winds in general. We will summarize the coordinated observational, theoretical, and modeling efforts from the Magnetism in Massive Star Project, addressing key outstanding questions regarding magnetosphere manifestations across the spectral domain.

preprint2011arXiv

Mass and angular momentum loss via decretion disks

We examine the nature and role of mass loss via an equatorial decretion disk in massive stars with near-critical rotation induced by evolution of the stellar interior. In contrast to the usual stellar wind mass loss set by exterior driving from the stellar luminosity, such decretion-disk mass loss stems from the angular momentum loss needed to keep the star near and below critical rotation, given the interior evolution and decline in the star's moment of inertia. Because the specific angular momentum in a Keplerian disk increases with the square root of the radius, the decretion mass loss associated with a required level of angular momentum loss depends crucially on the outer radius for viscous coupling of the disk, and can be significantly less than the spherical mass loss the spherical, wind-like mass loss commonly assumed in evolutionary calculations. We discuss the physical processes that affect the outer disk radius, including thermal disk outflow, and ablation of the disk material via a line-driven wind induced by the star's radiation. We present parameterized scaling laws for taking account of decretion-disk mass loss in stellar evolution codes, including how these are affected by metallicity, or by presence within a close binary and/or a dense cluster. Effects similar to those discussed here should also be present in accretion disks during star formation, and may play an important role in shaping the distribution of rotation speeds on the ZAMS.

preprint2011arXiv

Mass loss from inhomogeneous hot star winds II. Constraints from a combined optical/UV study

Mass-loss rates currently in use for hot, massive stars have recently been seriously questioned, mainly because of the effects of wind clumping. We investigate the impact of clumping on diagnostic ultraviolet resonance and optical recombination lines. Optically thick clumps, a non-void interclump medium, and a non-monotonic velocity field are all accounted for in a single model. We used 2D and 3D stochastic and radiation-hydrodynamic (RH) wind models, constructed by assembling 1D snapshots in radially independent slices. To compute synthetic spectra, we developed and used detailed radiative transfer codes for both recombination lines (solving the "formal integral") and resonance lines (using a Monte-Carlo approach). In addition, we propose an analytic method to model these lines in clumpy winds, which does not rely on optically thin clumping. Results: Synthetic spectra calculated directly from current RH wind models of the line-driven instability are unable to in parallel reproduce strategic optical and ultraviolet lines for the Galactic O-supergiant LCep. Using our stochastic wind models, we obtain consistent fits essentially by increasing the clumping in the inner wind. A mass-loss rate is derived that is approximately two times lower than predicted by the line-driven wind theory, but much higher than the corresponding rate derived from spectra when assuming optically thin clumps. Our analytic formulation for line formation is used to demonstrate the potential impact of optically thick clumping in weak-winded stars and to confirm recent results that resonance doublets may be used as tracers of wind structure and optically thick clumping. (Abridged)

preprint2011arXiv

Revisiting the Rigidly Rotating Magnetosphere model for sigma Ori E. I. Observations and Data Analysis

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

preprint2011arXiv

Stellar envelope inflation near the Eddington limit. Implications for the radii of Wolf-Rayet stars and luminous blue variables

(shortened) It has been proposed that the envelopes of luminous stars may be subject to substantial radius inflation. The inflation effect has been discussed in relation to the radius problem of WR stars, but has yet failed to explain the large observed radii of Galactic WR stars. We wish to obtain a physical perspective of the inflation effect, and study the consequences for the radii of WR stars, and LBVs. For WR stars the observed radii are up to an order of magnitude larger than predicted by theory, whilst S Doradus-type LBVs are subject to humongous radius variations, which remain as yet ill-explained. We use a dual approach to investigate the envelope inflation, based on numerical models for stars near the Eddington limit, and a new analytic formalism to describe the effect. An additional new aspect is that we take the effect of density inhomogeneities (clumping) within the outer stellar envelopes into account. Due to the effect of clumping we are able to bring the observed WR radii in agreement with theory. Based on our new formalism, we find that the radial inflation is a function of a dimensionless parameter W, which largely depends on the topology of the Fe-opacity peak, i.e., on material properties. For W>1, we discover an instability limit, for which the stellar envelope becomes gravitationally unbound, i.e. there no longer exists a static solution. Within this framework we are also able to explain the S Doradus-type instabilities for LBVs like AG Car, with a possible triggering due to changes in stellar rotation. The stellar effective temperatures in the upper HR diagram are potentially strongly affected by the inflation effect. This may have particularly strong effects on the evolved massive LBV and WR stars just prior to their final collapse, as the progenitors of SNe Ibc, SNe II, and long GRBs.

preprint2010arXiv

Discovery of Rotational Braking in the Magnetic Helium-Strong Star Sigma Orionis E

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

preprint2009arXiv

Gamma-ray variability from wind clumping in HMXBs with jets

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

preprint2006arXiv

Centrifugal Breakout of Magnetically Confined Line-Driven Stellar Winds

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

preprint2005arXiv

The Rigidly Rotating Magnetosphere of Sigma Ori E

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

preprint2004arXiv

A Rigidly Rotating Magnetosphere Model for Circumstellar Emission from Magnetic OB Stars

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

preprint2004arXiv

An XMM-Newton observation of the massive binary HD159176

We report the analysis of an XMM-Newton observation of the close binary HD159176 (O7V + O7V). The observed L_X/L_bol ratio reveals an X-ray luminosity exceeding by a factor ~ 7 the expected value for X-ray emission from single O-stars, therefore suggesting a wind-wind interaction scenario. EPIC and RGS spectra are fitted consistently with a two temperature mekal optically thin thermal plasma model, with temperatures ranging from ~ 2 to 6 10^6 K. At first sight, these rather low temperatures are consistent with the expectations for a close binary system where the winds collide well before reaching their terminal velocities. We also investigate the variability of the X-ray light curve of HD159176 on various short time scales. No significant variability is found and we conclude that if hydrodynamical instabilities exist in the wind interaction region of HD159176, they are not sufficient to produce an observable signature in the X-ray emission. Hydrodynamic simulations using wind parameters from the literature reveal some puzzling discrepancies. The most striking one concerns the predicted X-ray luminosity which is one or more orders of magnitude larger than the observed one. A significant reduction of the mass loss rate of the components compared to the values quoted in the literature alleviates the discrepancy but is not sufficient to fully account for the observed luminosity. Because hydrodynamical models are best for the adiabatic case whereas the colliding winds in HD159176 are most likely highly radiative, a totally new approach has been envisaged... (see paper for complete abstract)

preprint2004arXiv

Kinematic model inversions of hot star recurrent DAC data - tests against dynamical CIR models

The Discrete Absorption Components (DACs) commonly observed in the ultraviolet lines of hot stars have previously been modelled by dynamical simulations of Corotating Interaction Regions (CIRs) in their line-driven stellar winds. Here we apply the kinematic DAC inversion method of Brown et al. to the hydrodynamical CIR models and test the reliability of the results obtained. We conclude that the inversion method is able to recover valuable information on the velocity structure of the mean wind and to trace movement of velocity plateaux in the hydrodynamical data, though the recovered density profile of the stream is correct only very near to the stellar surface.

preprint2003arXiv

On multicomponent effects in stellar winds of stars at extremely low metallicity

We calculate multicomponent line-driven wind models of stars at extremely low metallicity suitable for massive first generation stars. For most of the models we find that the multicomponent wind nature is not important for either wind dynamics or for wind temperature stratification. However, for stars with the lowest metallicities we find that multicomponent effects influence the wind structure. These effects range from pure heating to possible fallback of the nonabsorbing wind component. We present a simple formula for the calculation of metallicity for which the multicomponent effects become important. We show that the importance of the multicomponent nature of winds of low metallicity stars is characterised not only by the low density of driving ions, but also by lower mass-loss rate.