BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Stanley P. Owocki

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

ResearcherUnclaimed source record
astro-ph.SRastro-phastro-ph.HEastro-ph.IM

Catalog footprint

What is connected

38works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

38 published item(s)

preprint2016arXiv

An `Analytic Dynamical Magnetosphere' formalism for X-ray and optical emission from slowly rotating magnetic massive stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Modelling the Central Constant Emission X-ray component of eta Carinae

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Spectral Temperature of Optically Thick Outflows with Application to Light Echo Spectra from $η$~Carinae's Giant Eruption

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Effect of scattering on the transonic solution topology and intrinsic variability of line-driven stellar winds

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hydrodynamic and radiative transfer modeling of X-ray emission from colliding WR winds: WR 140 & the Galactic center

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

3D Hydrodynamic & Radiative Transfer Models of X-ray Emission from Colliding Wind Binaries

Colliding wind binaries (CWBs) are unique laboratories for X-ray astrophysics. The massive stars in these systems possess powerful stellar winds with speeds up to $\sim$3000 km s$^{-1}$, and their collision leads to hot plasma (up to $\sim10^8$K) that emit thermal X-rays (up to $\sim$10 keV). Many X-ray telescopes have observed CWBs, including Suzaku, and our work aims to model these X-ray observations. We use 3D smoothed particle hydrodynamics (SPH) to model the wind-wind interaction, and then perform 3D radiative transfer to compute the emergent X-ray flux, which is folded through X-ray telescopes' response functions to compare directly with observations. In these proceedings, we present our models of Suzaku observations of the multi-year-period, highly eccentric systems $η$ Carinae and WR 140. The models reproduce the observations well away from periastron passage, but only $η$ Carinae's X-ray spectrum is reproduced at periastron; the WR 140 model produces too much flux during this more complicated phase.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Instabilities in the Envelopes and Winds of Very Massive Stars

The high luminosity of Very Massive Stars (VMS) means that radiative forces play an important, dynamical role both in the structure and stability of their stellar envelope, and in driving strong stellar-wind mass loss. Focusing on the interplay of radiative flux and opacity, with emphasis on key distinctions between continuum vs. line opacity, this chapter reviews instabilities in the envelopes and winds of VMS. Specifically, we discuss how: 1) the iron opacity bump can induce an extensive inflation of the stellar envelope; 2) the density dependence of mean opacity leads to strange mode instabilities in the outer envelope; 3) desaturation of line-opacity by acceleration of near-surface layers initiates and sustains a line-driven stellar wind outflow; 4) an associated line-deshadowing instability leads to extensive small-scale structure in the outer regions of such line-driven winds; 5) a star with super-Eddington luminosity can develop extensive atmospheric structure from photon bubble instabilities, or from stagnation of flow that exceeds the "photon tiring" limit; 6) the associated porosity leads to a reduction in opacity that can regulate the extreme mass loss of such continuum-driven winds. Two overall themes are the potential links of such instabilities to Luminous Blue Variable (LBV) stars, and the potential role of radiation forces in establishing the upper mass limit of VMS.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Measuring mass-loss rates and constraining shock physics using X-ray line profiles of O stars from the Chandra archive

We quantitatively investigate the extent of wind absorption signatures in the X-ray grating spectra of all non-magnetic, effectively single O stars in the Chandra archive via line profile fitting. Under the usual assumption of a spherically symmetric wind with embedded shocks, we confirm previous claims that some objects show little or no wind absorption. However, many other objects do show asymmetric and blue shifted line profiles, indicative of wind absorption. For these stars, we are able to derive wind mass-loss rates from the ensemble of line profiles, and find values lower by an average factor of 3 than those predicted by current theoretical models, and consistent with H-alpha if clumping factors of f_cl ~ 20 are assumed. The same profile fitting indicates an onset radius of X-rays typically at r ~ 1.5 R_star, and terminal velocities for the X-ray emitting wind component that are consistent with that of the bulk wind. We explore the likelihood that the stars in the sample that do not show significant wind absorption signatures in their line profiles have at least some X-ray emission that arises from colliding wind shocks with a close binary companion. The one clear exception is zeta Oph, a weak-wind star that appears to simply have a very low mass-loss rate. We also reanalyse the results from the canonical O supergiant zeta Pup, using a solar-metallicity wind opacity model and find Mdot = 1.8 \times 10^{-6} M_sun/yr, consistent with recent multi-wavelength determinations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Measuring the shock-heating rate in the winds of O stars using X-ray line spectra

We present a new method for using measured X-ray emission line fluxes from O stars to determine the shock-heating rate due to instabilities in their radiation-driven winds. The high densities of these winds means that their embedded shocks quickly cool by local radiative emission, while cooling by expansion should be negligible. Ignoring for simplicity any non-radiative mixing or conductive cooling, the method presented here exploits the idea that the cooling post-shock plasma systematically passes through the temperature characteristic of distinct emission lines in the X-ray spectrum. In this way, the observed flux distribution among these X-ray lines can be used to construct the cumulative probability distribution of shock strengths that a typical wind parcel encounters as it advects through the wind. We apply this new method (Gayley 2014) to Chandra grating spectra from five O stars with X-ray emission indicative of embedded wind shocks in effectively single massive stars. Correcting for wind absorption of the X-ray line emission is a crucial component of our analysis, and we use wind optical depth values derived from X-ray line-profile fitting (Cohen et al. 2014) in order to make that correction. The shock-heating rate results we derive for all the stars are quite similar: the average wind mass element passes through roughly one shock that heats it to at least $10^6$ K as it advects through the wind, and the cumulative distribution of shock strengths is a strongly decreasing function of temperature, consistent with a negative power-law of index $n \approx 3$, implying a marginal distribution of shock strengths that scales as $T^{-4}$, and with hints of an even steeper decline or cut-off above $10^7$ K.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Constraints on decreases in Eta Carinae's mass loss from 3D hydrodynamic simulations of its binary colliding winds

Recent work suggests that the mass-loss rate of the primary star (Eta A) in the massive colliding wind binary Eta Carinae dropped by a factor of 2-3 between 1999 and 2010. We present results from large- (r=1545au) and small- (r=155au) domain, 3D smoothed particle hydrodynamic (SPH) simulations of Eta Car's colliding winds for 3 Eta A mass-loss rates (2.4, 4.8, and 8.5 x 10^-4 M_sun/yr), investigating the effects on the dynamics of the binary wind-wind collision (WWC). These simulations include orbital motion, optically thin radiative cooling, and radiative forces. We find that Eta A's mass-loss rate greatly affects the time-dependent hydrodynamics at all spatial scales investigated. The simulations also show that the post-shock wind of the companion star (Eta B) switches from the adiabatic to the radiative-cooling regime during periastron passage. The SPH simulations together with 1D radiative transfer models of Eta A's spectra reveal that a factor of 2 or more drop in Eta A's mass-loss rate should lead to substantial changes in numerous multiwavelength observables. Recent observations are not fully consistent with the model predictions, indicating that any drop in Eta A's mass-loss rate was likely by a factor < 2 and occurred after 2004. We speculate that most of the recent observed changes in Eta Car are due to a small increase in the WWC opening angle that produces significant effects because our line-of-sight to the system lies close to the dense walls of the WWC zone. A modest decrease in Eta A's mass-loss rate may be responsible, but changes in the wind/stellar parameters of Eta B cannot yet be fully ruled out. We suggest observations during Eta Car's next periastron in 2014 to further test for decreases in Eta A's mass-loss rate. If Eta A's mass-loss rate is declining and continues to do so, the 2014 X-ray minimum should be even shorter than that of 2009.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Constraints on porosity and mass loss in O-star winds from modeling of X-ray emission line profile shapes

We fit X-ray emission line profiles in high resolution XMM-Newton and Chandra grating spectra of the early O supergiant Zeta Pup with models that include the effects of porosity in the stellar wind. We explore the effects of porosity due to both spherical and flattened clumps. We find that porosity models with flattened clumps oriented parallel to the photosphere provide poor fits to observed line shapes. However, porosity models with isotropic clumps can provide acceptable fits to observed line shapes, but only if the porosity effect is moderate. We quantify the degeneracy between porosity effects from isotropic clumps and the mass-loss rate inferred from the X-ray line shapes, and we show that only modest increases in the mass-loss rate (<~ 40%) are allowed if moderate porosity effects (h_infinity <~ R_*) are assumed to be important. Large porosity lengths, and thus strong porosity effects, are ruled out regardless of assumptions about clump shape. Thus, X-ray mass-loss rate estimates are relatively insensitive to both optically thin and optically thick clumping. This supports the use of X-ray spectroscopy as a mass-loss rate calibration for bright, nearby O stars.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A dynamical magnetosphere model for periodic Halpha emission from the slowly rotating magnetic O star HD191612

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An X-ray Survey of Colliding Wind Binaries

We have compiled a list of 36 O+O and 89 Wolf-Rayet binary candidates in the Milky Way and Magellanic clouds detected with the Chandra, XMM-Newton and ROSAT satellites to probe the connection between their X-ray properties and their system characteristics. Of the WR binaries with published parameters, all but two have kT > 0.9 keV. The most X-ray luminous WR binaries are typically very long period systems. The WR binaries show a nearly four-order of magnitude spread in X-ray luminosity, even among among systems with very similar WR primaries. Among the O+O binaries, short-period systems generally have soft X-ray spectra and longer period systems show harder X-ray spectra, again with a large spread in LX/Lbol.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Clumping in the inner winds of hot, massive stars from hydrodynamical line-driven instability simulations

We investigate the effects of stellar limb-darkening and photospheric perturbations for the onset of wind structure arising from the strong, intrinsic line-deshadowing instability (LDI) of a line-driven stellar wind. A linear perturbation analysis shows that including limb-darkening reduces the stabilizing effect of the diffuse radiation, leading to a net instability growth rate even at the wind base. Numerical radiation-hydrodynamics simulations of the non-linear evolution of this instability then show that, in comparison with previous models assuming a uniformly bright star without base perturbations, wind structure now develops much closer ($r \la 1.1 R_\star$) to the photosphere. This is in much better agreement with observations of O-type stars, which typically indicate the presence of strong clumping quite near the wind base.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

First 3D MHD simulation of a massive-star magnetosphere with application to Hαemission from θ^1 Ori C

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Thin-shell mixing in radiative wind-shocks and the Lx-Lbol scaling of O-star X-rays

X-ray satellites since Einstein have empirically established that the X-ray luminosity from single O-stars scales linearly with bolometric luminosity, Lx ~ 10^{-7} Lbol. But straightforward forms of the most favored model, in which X-rays arise from instability-generated shocks embedded in the stellar wind, predict a steeper scaling, either with mass loss rate Lx ~ Mdot ~ Lbol^{1.7} if the shocks are radiative, or with Lx ~ Mdot^{2} ~ Lbol^{3.4} if they are adiabatic. This paper presents a generalized formalism that bridges these radiative vs. adiabatic limits in terms of the ratio of the shock cooling length to the local radius. Noting that the thin-shell instability of radiative shocks should lead to extensive mixing of hot and cool material, we propose that the associated softening and weakening of the X-ray emission can be parametrized as scaling with the cooling length ratio raised to a power m$, the "mixing exponent". For physically reasonable values m ~= 0.4, this leads to an X-ray luminosity Lx ~ Mdot^{0.6} ~ Lbol that matches the empirical scaling. To fit observed X-ray line profiles, we find such radiative-shock-mixing models require the number of shocks to drop sharply above the initial shock onset radius. This in turn implies that the X-ray luminosity should saturate and even decrease for optically thick winds with very high mass-loss rates. In the opposite limit of adiabatic shocks in low-density winds (e.g., from B-stars), the X-ray luminosity should drop steeply with Mdot^2. Future numerical simulation studies will be needed to test the general thin-shell mixing ansatz for X-ray emission.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A generalised porosity formalism for isotropic and anisotropic effective opacity and its effects on X-ray line attenuation in clumped O star winds

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Chandra X-ray spectroscopy of the very early O supergiant HD 93129A: constraints on wind shocks and the mass-loss rate

We present analysis of both the resolved X-ray emission line profiles and the broadband X-ray spectrum of the O2 If* star HD 93129A, measured with the Chandra HETGS. This star is among the earliest and most massive stars in the Galaxy, and provides a test of the embedded wind shock scenario in a very dense and powerful wind. A major new result is that continuum absorption by the dense wind is the primary cause of the hardness of the observed X-ray spectrum, while intrinsically hard emission from colliding wind shocks contributes less than 10% of the X-ray flux. We find results consistent with the predictions of numerical simulations of the line-driving instability, including line broadening indicating an onset radius of X-ray emission of several tenths Rstar. Helium-like forbidden-to-intercombination line ratios are consistent with this onset radius, and inconsistent with being formed in a wind-collision interface with the star's closest visual companion at a distance of ~100 AU. The broadband X-ray spectrum is fit with a dominant emission temperature of just kT = 0.6 keV along with significant wind absorption. The broadband wind absorption and the line profiles provide two independent measurements of the wind mass-loss rate: Mdot = 5.2_{-1.5}^{+1.8} \times 10^{-6} Msun/yr and Mdot = 6.8_{-2.2}^{+2.8} \times 10^{-6} Msun/yr, respectively. This is the first consistent modeling of the X-ray line profile shapes and broadband X-ray spectral energy distribution in a massive star, and represents a reduction of a factor of 3 to 4 compared to the standard H-alpha mass-loss rate that assumes a smooth wind.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Constraining the Absolute Orientation of Eta Carinae's Binary Orbit: A 3-D Dynamical Model for the Broad [Fe III] Emission

We present a three-dimensional (3-D) dynamical model for the broad [Fe III] emission observed in Eta Carinae using the Hubble Space Telescope/Space Telescope Imaging Spectrograph (HST/STIS). This model is based on full 3-D Smoothed Particle Hydrodynamics (SPH) simulations of Eta Car's binary colliding winds. Radiative transfer codes are used to generate synthetic spectro-images of [Fe III] emission line structures at various observed orbital phases and STIS slit position angles (PAs). Through a parameter study that varies the orbital inclination i, the PA θ that the orbital plane projection of the line-of-sight makes with the apastron side of the semi-major axis, and the PA on the sky of the orbital axis, we are able, for the first time, to tightly constrain the absolute 3-D orientation of the binary orbit. To simultaneously reproduce the blue-shifted emission arcs observed at orbital phase 0.976, STIS slit PA = +38 degrees, and the temporal variations in emission seen at negative slit PAs, the binary needs to have an i \approx 130 to 145 degrees, θ \approx -15 to +30 degrees, and an orbital axis projected on the sky at a PA \approx 302 to 327 degrees east of north. This represents a system with an orbital axis that is closely aligned with the inferred polar axis of the Homunculus nebula, in 3-D. The companion star, Eta B, thus orbits clockwise on the sky and is on the observer's side of the system at apastron. This orientation has important implications for theories for the formation of the Homunculus and helps lay the groundwork for orbital modeling to determine the stellar masses.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hydrodynamic Interaction between the Be Star and the Pulsar in the TeV Binary PSR B1259-63/LS 2883

We study the interaction between the Be star and the pulsar in the TeV binary PSR B1259-63/LS 2883, using 3-D SPH simulations of the tidal and wind interactions in this Be-pulsar system. We first run a simulation without pulsar wind nor Be wind, taking into account only the gravitational effect of the pulsar on the Be disk. In this simulation, the gas particles are ejected at a constant rate from the equatorial surface of the Be star, which is tilted in a direction consistent with multi-waveband observations. We run the simulation until the Be disk is fully developed and starts to repeat a regular tidal interaction with the pulsar. Then, we turn on the pulsar wind and the Be wind. We run two simulations with different wind mass-loss rates for the Be star, one for a B2V type and the other for a significantly earlier spectral type. Although the global shape of the interaction surface between the pulsar wind and the Be wind agrees with the analytical solution, the effect of the pulsar wind on the Be disk is profound. The pulsar wind strips off an outer part of the Be disk, truncating the disk at a radius significantly smaller than the pulsar orbit. Our results, therefore, rule out the idea that the pulsar passes through the Be disk around periastron, which has been assumed in the previous studies. It also turns out that the location of the contact discontinuity can be significantly different between phases when the pulsar wind directly hits the Be disk and those when the pulsar wind collides with the Be wind. It is thus important to adequately take into account the circumstellar environment of the Be star, in order to construct a satisfactory model for this prototypical TeV binary.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Multi-Wavelength Implications of the Companion Star in Eta Carinae

Eta Carinae is considered to be a massive colliding wind binary system with a highly eccentric (e \sim 0.9), 5.54-yr orbit. However, the companion star continues to evade direct detection as the primary dwarfs its emission at most wavelengths. Using three-dimensional (3-D) SPH simulations of Eta Car's colliding winds and radiative transfer codes, we are able to compute synthetic observables across multiple wavebands for comparison to the observations. The models show that the presence of a companion star has a profound influence on the observed HST/STIS UV spectrum and H-alpha line profiles, as well as the ground-based photometric monitoring. Here, we focus on the Bore Hole effect, wherein the fast wind from the hot secondary star carves a cavity in the dense primary wind, allowing increased escape of radiation from the hotter/deeper layers of the primary's extended wind photosphere. The results have important implications for interpretations of Eta Car's observables at multiple wavelengths.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The nature and consequences of clumping in hot, massive star winds

This review describes the evidence for small-scale structure, `clumping', in the radiation line-driven winds of hot, massive stars. In particular, we focus on examining to what extent simulations of the strong instability inherent to line-driving can explain the multitude of observational evidence for wind clumping, as well as on how to properly account for extensive structures in density and velocity when interpreting the various wind diagnostics used to derive mass-loss rates.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

X-ray Modeling of η Carinae and WR140 from SPH Simulations

The colliding wind binary (CWB) systems η Carinae and WR140 provide unique laboratories for X-ray astrophysics. Their wind-wind collisions produce hard X-rays that have been monitored extensively by several X-ray telescopes, including RXTE. To interpret these X-ray light curves and spectra, we apply 3D hydrodynamic simulations of the wind-wind collision using smoothed particle hydrodynamics (SPH), with the recent improvements of radiative cooling and the acceleration of the stellar winds according to a β law. For both systems, the 2-10 keV RXTE light curves are well-reproduced in absolute units for most phases, but the light curve dips associated with the periastron passages are not well matched. In WR140, the dip is too weak, and in η Carinae, the large difference in wind speeds of the two stars leads to a hot, post-periastron bubble that produces excess emission toward the end of the X-ray minimum.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

X-ray Modeling of η Carinae and WR140 from SPH Simulations

The colliding wind binary (CWB) systems η Carinae and WR140 provide unique laboratories for X-ray astrophysics. Their wind-wind collisions produce hard X-rays that have been monitored extensively by several X-ray telescopes, including RXTE. To interpret these RXTE X-ray light curves, we model the wind-wind collision using 3D smoothed particle hydrodynamics (SPH) simulations. Adiabatic simulations that account for the absorption of X-rays from an assumed point source at the apex of the wind-collision shock cone by the distorted winds can closely match the observed 2-10keV RXTE light curves of both η Car and WR140. This point-source model can also explain the early recovery of η Car's X-ray light curve from the 2009.0 minimum by a factor of 2-4 reduction in the mass loss rate of η Car. Our more recent models relax the point-source approximation and account for the spatially extended emission along the wind-wind interaction shock front. For WR140, the computed X-ray light curve again matches the RXTE observations quite well. But for η Car, a hot, post-periastron bubble leads to an emission level that does not match the extended X-ray minimum observed by RXTE. Initial results from incorporating radiative cooling and radiatively-driven wind acceleration via a new anti-gravity approach into the SPH code are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A mass-loss rate determination for zeta Puppis from the quantitative analysis of X-ray emission line profiles

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A proper description of clumping in hot star winds: the key to obtaining reliable mass-loss rates?

Small-scale inhomogeneities, or `clumping', in the winds of hot, massive stars are conventionally included in spectral analyses by assuming optically thin clumps. To reconcile investigations of different diagnostics using this microclumping technique, very low mass-loss rates must be invoked for O stars. Recently it has been suggested that by using the microclumping approximation one may actually drastically underestimate the mass-loss rates. Here we demonstrate this, present a new, improved description of clumpy winds, and show how corresponding models, in a combined UV and optical analysis, can alleviate discrepancies between previously derived rates and those predicted by the line-driven wind theory. Furthermore, we show that the structures obtained in time-dependent, radiation-hydrodynamic simulations of the intrinsic line-driven instability of such winds, which are the basis to our current understanding of clumping, in their present-day form seem unable to provide a fully self-consistent, simultaneous fit to both UV and optical lines. The reasons for this are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Is Eta Carinae a fast rotator, and how much does the companion influence the inner wind structure?

We analyze interferometric measurements of the Luminous Blue Variable Eta Carinae with the goal of constraining the rotational velocity of the primary star and probing the influence of the companion. Using 2-D radiative transfer models of latitude-dependent stellar winds, we find that prolate wind models with a ratio of the rotational velocity (vrot) to the critical velocity (vcrit) of W=0.77-0.92, inclination angle of i=60-90 degrees, and position angle PA=108-142 degrees reproduce simultaneously K-band continuum visibilities from VLTI/VINCI and closure phase measurements from VLTI/AMBER. Interestingly, oblate models with W=0.73-0.90 and i=80-90 degrees produce similar fits to the interferometric data, but require PA=210-230 degrees. Therefore, both prolate and oblate models suggest that the rotation axis of the primary star is not aligned with the Homunculus polar axis. We also compute radiative transfer models of the primary star allowing for the presence of a cavity and dense wind-wind interaction region created by the companion star. We find that the wind-wind interaction has a significant effect on the K-band image mainly via free-free emission from the compressed walls and, for reasonable model parameters, can reproduce the VLTI/VINCI visibilities taken at phase 0.92-0.93. We conclude that the density structure of the primary wind can be sufficiently disturbed by the companion, thus mimicking the effects of fast rotation in the interferometric observables. Therefore, fast rotation may not be the only explanation for the interferometric observations. Intense temporal monitoring and 3-D modeling are needed to resolve these issues.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerical models of collisions between core-collapse supernovae and circumstellar shells

Recent observations of luminous Type IIn supernovae (SNe) provide compelling evidence that massive circumstellar shells surround their progenitors. In this paper we investigate how the properties of such shells influence the SN lightcurve by conducting numerical simulations of the interaction between an expanding SN and a circumstellar shell ejected a few years prior to core collapse. Our parameter study explores how the emergent luminosity depends on a range of circumstellar shell masses, velocities, geometries, and wind mass-loss rates, as well as variations in the SN mass and energy. We find that the shell mass is the most important parameter, in the sense that higher shell masses (or higher ratios of M_shell/M_SN) lead to higher peak luminosities and higher efficiencies in converting shock energy into visual light. Lower mass shells can also cause high peak luminosities if the shell is slow or if the SN ejecta are very fast, but only for a short time. Sustaining a high luminosity for durations of more than 100 days requires massive circumstellar shells of order 10 M_sun or more. This reaffirms previous comparisons between pre-SN shells and shells produced by giant eruptions of luminous blue variables (LBVs), although the physical mechanism responsible for these outbursts remains uncertain. The lightcurve shape and observed shell velocity can help diagnose the approximate size and density of the circumstellar shell, and it may be possible to distinguish between spherical and bipolar shells with multi-wavelength lightcurves. These models are merely illustrative. One can, of course, achieve even higher luminosities and longer duration light curves from interaction by increasing the explosion energy and shell mass beyond values adopted here.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

A Nozzle Analysis of Slow-Acceleration Solutions in One-Dimensional Models of Rotating Hot-Star Winds

We analyze the steady 1D flow equations for a rotating stellar wind based on a ``nozzle'' analogy for terms that constrain the local mass flux. For low rotation, we find the nozzle minimum occurs near the stellar surface, allowing a transition to a standard, CAK-type steep-acceleration solution; but for rotations > 75% of the critical rate, this inner nozzle minimum exceeds the global minimum, implying near-surface supercritical solutions would have an overloaded mass loss rate. In steady, analytic models in which the acceleration is assumed to be monotonically positive, this leads the solution to switch to a slow acceleration mode. However, time-dependent simulations using a numerical hydrodynamics code show that, for rotation rates 75 - 85% of critical, the flow can develop abrupt "kink" transitions from a steep acceleration to a decelerating solution. For rotations above 85% of critical, the hydrodynamic simulations confirm the slow acceleration, with the lower flow speed implying densities 5 - 30 times higher than the polar (or a nonrotating) wind. Still, when gravity darkening and 2D flow effects are accounted for, it seems unlikely that rotationally modified equatorial wind outflows could account for the very large densities inferred for the equatorial regions around B[e] supergiants.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

On the role of continuum-driven eruptions in the evolution of very massive stars and Population III stars

We suggest that the mass lost during the evolution of very massive stars may be dominated by optically thick, continuum-driven outbursts or explosions, instead of by steady line-driven winds. In order for a massive star to become a WR star, it must shed its H envelope, but new estimates of the effects of clumping in winds indicate that line driving is vastly insufficient. We discuss massive stars above roughly 40-50 Msun, for which the best alternative is mass loss during brief eruptions of luminous blue variables (LBVs). Our clearest example of this phenomenon is the 19th century outburst of eta Car, when the star shed 12-20 Msun or more in less than a decade. Other examples are circumstellar nebulae of LBVs, extragalactic eta Car analogs (``supernova impostors''), and massive shells around SNe and GRBs. We do not yet fully understand what triggers LBV outbursts, but they occur nonetheless, and present a fundamental mystery in stellar astrophysics. Since line opacity from metals becomes too saturated, the extreme mass loss probably arises from a continuum-driven wind or a hydrodynamic explosion, both of which are insensitive to metallicity. As such, eruptive mass loss could have played a pivotal role in the evolution and fate of massive metal-poor stars in the early universe. If they occur in these Population III stars, such eruptions would profoundly affect the chemical yield and types of remnants from early SNe and hypernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Chandra HETGS Multi-Phase Spectroscopy of the Young Magnetic O Star theta^1 Orionis C

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

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

A Porosity-Length Formalism for Photon-Tiring-Limited Mass Loss from Stars Above the Eddington Limit

We examine radiatively driven mass loss from stars near and above the Eddington limit (Ledd). We begin by reviewing the instabilities that are expected to form extensive structure near Ledd. We investigate how this "porosity" can reduce the effective coupling between the matter and radiation. Introducing a new "porosity-length'' formalism, we derive a simple scaling for the reduced effective opacity, and use this to derive an associated scaling for the porosity-moderated, continuum-driven mass loss rate from stars that formally exceed Ledd. For a simple super-Eddington model with a single porosity length that is assumed to be on the order of the gravitational scale height, the overall mass loss is similar to that derived in previous porosity work. This is much higher than is typical of line-driven winds, but is still only a few percent of the photon tiring limit--for which the luminosity becomes insufficient to carry the flow out of the gravitational potential. To obtain still stronger mass loss that approaches observationally inferred values near this limit, we introduce a power-law-porosity model in which the associated structure has a broad range of scales. We show that the mass loss rate can be enhanced by a factor that increases with the Eddington parameter Gamma, such that for moderately large Gamma (> 3-4), mass loss rates could approach the photon tiring limit. Together with the ability to drive quite fast outflow speeds, the derived mass loss could explain the large inferred mass loss and flow speeds of giant outbursts in eta Carinae and other LBV stars.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Be-star rotation: how close to critical?

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

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

High-Resolution {\it Chandra} Spectroscopy of tau Scorpii: A Narrow-Line X-ray Spectrum From a Hot Star

Long known to be an unusual early-type star by virtue of its hard and strong X-ray emission, tau Scorpii poses a severe challenge to the standard picture of O-star wind-shock X-ray emission. The Chandra HETGS spectrum now provides significant direct evidence that this B0.2 star does not fit this standard wind-shock framework. The many emission lines detected with the Chandra gratings are significantly narrower than what would be expected from a star with the known wind properties of tau Sco, although they are broader than the corresponding lines seen in late-type coronal sources. While line ratios are consistent with the hot plasma on this star being within a few stellar radii of the photosphere, from at least one He-like complex there is evidence that the X-ray-emitting plasma is located more than a stellar radius above the photosphere. The Chandra spectrum of tau Sco is harder and more variable than those of other hot stars, with the exception of the young magnetized O star theta Ori C. We discuss these new results in the context of wind, coronal, and hybrid wind-magnetic models of hot-star X-ray emission.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

On the Steady Nature of Line-Driven Disk Winds

We perform an analytic investigation of the stability of line-driven disk winds, independent of hydrodynamic simulations. Our motive is to determine whether or not line-driven disk winds can account for the wide/broad UV resonance absorption lines seen in cataclysmic variables (CVs) and quasi-stellar objects (QSOs). In both CVs and QSOs observations generally indicate that the absorption arising in the outflowing winds has a steady velocity structure on time scales exceeding years (for CVs) and decades (for QSOs). However, published results from hydrodynamic simulations of line-driven disk winds are mixed, with some researchers claiming that the models are inherently unsteady, while other models produce steady winds. The analytic investigation presented here shows that if the accretion disk is steady, then the line-driven disk wind emanating from it can also be steady. In particular, we show that a gravitational force initially increasing along the wind streamline, which is characteristic of disk winds, does not imply an unsteady wind. The steady nature of line-driven disk winds is consistent with the 1D streamline disk-wind models of Murray and collaborators and the 2.5D time-dependent models of Pereyra and collaborators. This paper emphasizes the underlying physics behind the steady nature of line-driven disk winds using mathematically simple models that mimic the disk environment

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

X-ray Emission Line Profile Modeling of O stars: Fitting a Spherically-Symmetric Analytic Wind-Shock Model to the Chandra Spectrum of zeta Puppis

X-ray emission line profiles provide the most direct insight into the dynamics and spatial distribution of the hot, X-ray-emitting plasma above the surfaces of OB stars. The O supergiant zeta Puppis shows broad, blueshifted, and asymmetric line profiles, generally consistent with the wind-shock picture. We model the profiles of eight lines in the Chandra HETGS spectrum of this prototypical hot star. The fitted lines indicate that the plasma is distributed throughout the wind starting close to the photosphere, that there is significantly less attenuation of the X-rays by the overlying wind than is generally supposed, and that there is not a strong trend in wind absorption with wavelength.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Dynamical Simulations of Magnetically Channeled Line-Driven Stellar Winds: I. Isothermal, Nonrotating, Radially Driven Flow

We present numerical magnetohydrodynamic (MHD) simulations of the effect of stellar dipole magnetic fields on line-driven wind outflows from hot, luminous stars. Unlike previous fixed-field analyses, the simulations here take full account of the dynamical competition between field and flow, and thus apply to a full range of magnetic field strength, and within both closed and open magnetic topologies. A key result is that the overall degree to which the wind is influenced by the field depends largely on a single, dimensionless, `wind magnetic confinement parameter', $η_{\ast}$ ($ = B_{eq}^2 R_\ast^2/{\dot M} v_\infty$), which characterizes the ratio between magnetic field energy density and kinetic energy density of the wind. For weak confinement $η_{\ast} \le 1$, the field is fully opened by the wind outflow, but nonetheless for confinements as small as $η_{\ast}=1/10$ can have a significant back-influence in enhancing the density and reducing the flow speed near the magnetic equator. For stronger confinement $η_{\ast} > 1$, the magnetic field remains closed over a limited range of latitude and height about the equatorial surface, but eventually is opened into a nearly radial configuration at large radii.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

X ray emission line profile modeling of hot stars

The launch of high-spectral-resolution x-ray telescopes (Chandra, XMM) has provided a host of new spectral line diagnostics for the astrophysics community. In this paper we discuss Doppler-broadened emission line profiles from highly supersonic outflows of massive stars. These outflows, or winds, are driven by radiation pressure and carry a tremendous amount of kinetic energy, which can be converted to x rays by shock-heating even a small fraction of the wind plasma. The unshocked, cold wind is a source of continuum opacity to the x rays generated in the shock-heated portion of the wind. Thus the emergent line profiles are affected by transport through a two-component, moving, optically thick medium. While complicated, the interactions among these physical effects can provide quantitative information about the spatial distribution and velocity of the x-ray-emitting and absorbing plasma in stellar winds. We present quantitative models of both a spherically-symmetric wind and a wind with hot plasma confined in an equatorial disk by a dipole magnetic field.

0 citations0 reviewsNo code linkedOpen access