Search NASASearch

Engineering topics

Owocki, Stanley P.

Publications and source records attributed to Owocki, Stanley P..

At least 19 records

Constraints on Porosity and Mass Loss in O-star Winds from the 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(sub infinity) less than approximately R(sub ∗)) 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

porosity

Atomic Physics of Shocked Plasma in Winds of Massive Stars

High resolution diffraction grating spectra of X-ray emission from massive stars obtained with Chandra and XMM-Newton have revolutionized our understanding of their powerful, radiation-driven winds. Emission line shapes and line ratios provide diagnostics on a number of key wind parameters. Modeling of resolved emission line velocity profiles allows us to derive independent constraints on stellar mass-loss rates, leading to downward revisions of a factor of a few from previous measurements. Line ratios in He-like ions strongly constrain the spatial distribution of Xray emitting plasma, confirming the expectations of radiation hydrodynamic simulations that X-ray emission begins moderately close to the stellar surface and extends throughout the wind. Some outstanding questions remain, including the possibility of large optical depths in resonance lines, which is hinted at by differences in line shapes of resonance and intercombination lines from the same ion. Resonance scattering leads to nontrivial radiative transfer effects, and modeling it allows us to place constraints on shock size, density, and velocity structure

Physics

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 approximately 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.

Madura, Thomas I.

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. As in previous work, we concentrate on developing bridging laws between the limits of optically thin and thick clumps. 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.

Sundqvist, Jon O.

Constraining the Properties of the Eta Carinae System via 3-D SPH Models of Space-Based Observations: The Absolute Orientation of the Binary Orbit

The extremely massive (> 90 Stellar Mass) and luminous (= 5 x 10(exp 6) Stellar Luminosity) star Eta Carinae, with its spectacular bipolar "Homunculus" nebula, comprises one of the most remarkable and intensely observed stellar systems in the Galaxy. However, many of its underlying physical parameters remain unknown. Multiwavelength variations observed to occur every 5.54 years are interpreted as being due to the collision of a massive wind from the primary star with the fast, less dense wind of a hot companion star in a highly elliptical (e approx. 0.9) orbit. Using three-dimensional (3-D) Smoothed Particle Hydrodynamics (SPH) simulations of the binary wind-wind collision, together with radiative transfer codes, we compute synthetic spectral images of [Fe III] emission line structures and compare them to existing Hubble Space Telescope/Space Telescope Imaging Spectrograph (HST/STIS) observations. We are thus able, for the first time, to tightly constrain the absolute orientation of the binary orbit on the sky. An orbit with an inclination of approx. 40deg, an argument of periapsis omega approx. 255deg, and a projected orbital axis with a position angle of approx. 312deg east of north provides the best fit to the observations, implying that the orbital axis is closely aligned in 3-D space with the Homunculus symmetry axis, and that the companion star orbits clockwise on the sky relative to the primary.

Madura, Thomas I.

Constraining the Properties of the Eta Carinae System via 3-D SPH Models of Space-Based Observations: The Absolute Orientation of the Binary Orbit

The extremely massive (> 90 Solar Mass) and luminous (= 5 x 10(exp 6) Solar Luminosity) star Eta Carinae, with its spectacular bipolar "Homunculus" nebula, comprises one of the most remarkable and intensely observed stellar systems in the galaxy. However, many of its underlying physical parameters remain a mystery. Multiwavelength variations observed to occur every 5.54 years are interpreted as being due to the collision of a massive wind from the primary star with the fast, less dense wind of a hot companion star in a highly elliptical (e approx. 0.9) orbit. Using three-dimensional (3-D) Smoothed Particle Hydrodynamics (SPH) simulations of the binary wind-wind collision in Eta Car, together with radiative transfer codes, we compute synthetic spectral images of [Fe III] emission line structures and compare them to existing Hubble Space Telescope/Space Telescope Imaging Spectrograph (HST/STIS) observations. We are thus able, for the first time, to constrain the absolute orientation of the binary orbit on the sky. An orbit with an inclination of i approx. 40deg, an argument of periapsis omega approx. 255deg, and a projected orbital axis with a position angle of approx. 312deg east of north provides the best fit to the observations, implying that the orbital axis is closely aligned in 3-1) space with the Homunculus symmetry axis, and that the companion star orbits clockwise on the sky relative to the primary.

Madura, Thomas I.

Modeling Broadband X-Ray Absorption of Massive Star Winds

We present a method for computing the net transition of X-rays emitted by shock-heated plasma distributed throughout a partially optically thick stellar wind from a massive star. We find the transmission by an exact integration of the formal solution, assuming the emitting plasma and absorbing plasma are mixed at a constant mass ratio above some minimum radius, below which there is assumed to be no emission. This model is more realistic than either the slab absorption associated with a corona at the base of the wind or the exospheric approximation that assumes all observed X-rays are emitted without attenuation from above the radius of optical depth unity. Our model is implemented in XSPEC as a pre-calculated table that can be coupled to a user-defined table of the wavelength dependent wind opacity. We provide a default wind opacity model that is more representative of real wind opacities than the commonly used neutral ISM tabulation. Preliminary modeling of Chandra grating data indicates that the X-ray hardness trend of OB stars with spectral subtype cars largely be understood as a wind absorption effect.

Leutenegger, Maurice A.

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. As in previous work, we concentrate on developing bridging laws between the limits of optically thin and thick clumps. 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.

generalised

Momentum deposition on Wolf-Rayet winds: Nonisotropic diffusion with effective gray opacity

We derive the velocity and mass-loss rate of a steady state Wolf-Rayet (WR) wind, using a nonisotropic diffusion approximation applied to the transfer between strongly overlapping spectral lines. Following the approach of Friend & Castor (1983), the line list is assumed to approximate a statistically parameterized Poisson distribution in frequency, so that photon transport is controlled by an angle-dependent, effectively gray opacity. We show the nonisotropic diffusion approximation yields good agreement with more accurate numerical treatments of the radiative transfer, while providing analytic insight into wind driving by multiple scattering. We illustrate, in particular, that multiple radiative momentum deposition does not require that potons be repeatedly reflected across substantial distances within the spherical envelope, but indeed is greatest when photons undergo a nearly local diffusion, e.g., through scattering by many lines closely spaced in frequency. Our results reiterate the view that the so-called 'momentum problem' of Wolf-Rayet winds is better characterized as an 'opacity problem' of simply identfying enough lines. One way of increasing the number of thick lines in Wolf-Rayet winds is to transfer opacity from saturated to unsaturated lines, yielding a steeper opacity distribution than that found in OB winds. We discuss the implications of this perspective for extending our approach to W-R wind models that incorporate a more fundamental treatment of the ionization and excitation processes that determine the line opacity. In particular, we argue that developing statistical descriptions of the lines to allow an improved effective opacity for the line ensemble would offer several advantages for deriving such more fundamental W-R wind models.

Gayley, Kenneth G.

The effect of oblateness and gravity darkening on the radiation driving in winds from rapidly rotating B stars

We calculate the radiative driving force for winds around rapidly rotating oblate B stars, and we estimate the impact these forces should have on the production of a wind compressed disk. The effects of limb darkening, gravity darkening, oblateness, and an arbitrary wind velocity field are included in the computation of vector 'oblate finite disk' (OFD) factors, which depend on both radius and colatitude in the wind. The impact of limb darkening alone, with or without rotation, can increase the mass loss by as much as 10% over values computed using the standard uniformly bright spherical finite disk factor. For rapidly rotating stars, limb darkening makes 'sub-stellar' gravity darkening the dominant effect in the radial and latitudinal OFD factors, and lessens the impact of gravity darkening at other visible latitudes (nearer to the oblate limb). Thus, the radial radiative driving is generally stronger over the poles and weaker over the equator, following the gravity darkening at these latitudes. The nonradial radiative driving is considerably smaller in magnitude than the radial component, but is directed both away from the equatorial plane and in a retrograde azimuthal direction, acting to decrease the effective stellar rotation velocity. These forces thus weaken the equatorward wind compression compared to wind models computed with nonrotating finite disk factors.

Cranmer, Steven R.

Theory review: Line-driven instability and other causes of structure and variability in hot-star winds

The winds of the hot, luminous, OB stars are driven by the line-scattering of the star's continuum radiation flux. Several kinds of observational evidence indicate that such winds are highly structured and variable. This paper will review possible theoretical causes of such wind structure. For relatively small-scale, stochastic variability, I review the role of the strong intrinsic instability of the line-driving process itself. for larger scale structure, I describe recent efforts to examine how disturbances from the underlying, rotating star can be translated outward into propagating features in the wind.

Owocki, Stanley P.

Acceleration efficiency in line-driven flows

We reexamine the physics of flow driving by line scattering of a continuum radiation source to determine the degree to which such line scattering can heat as well as accelerate the flow. Within the framework of the Sobolev theory for line transfer, we argue that the finite thermal width of the line scattering profile can lead to a significant 'Doppler heating' via photon frequency redistribution within a Sobolev resonance layer. Quantitative computation of this heating shows, however, that it is largely canceled by a corresponding cooling by the diffuse radiation. The resulting reduction in net Doppler heating or cooling means that the overall effect is only of limited importance in the energy balance of line-driven stellar winds. Through simple scaling relations, we compare the effect to other competing heating or cooling terms, including the ion-drag frictional heating recently discussed by Springmann and Pauldrach. We also provide a physical explanation of the unexpected cooling effect, and comment that its near cancellation of the anticipated heating provides another example of the tendency for ideal Sobolev theory to apply to a higher order than expected.

Gayley, Kenneth G.

Two-dimensional hydrodynamical simulations of wind-compressed disks around rapidly rotating B stars

We use a two-dimensional piecewise parabolic method (PPM) code to simulate numerically the hydrodynamics of a radiation-driven stellar wind from a rapidly rotating Be star. The results generally confirm predictions of the semianalytic 'wind-compressed disk' model recently proposed by Bjorkman and Cassinelli to explain the circumstellar disks inferred observationally to exist around such rapidly rotating stars. However, this numerical simulation is able to incorporate several important effects not accounted for in the simple model, including a dynamical treatment of the outward radiative driving and gas pressure, as well as a rotationally distorted, oblate stellar surface. This enables us to model quantitatively the compressed wind and shock that forms the equatorial disk. The simulation results thus do differ in several important details from the simple method, showing, for example, an inner disk inflow not possible in the heuristic approach of assuming a fixed outward velocity law. There is also no evidence for the predicted detachment of the disk that arises in the fixed outflow picture. The peak equatorward velocity in the dynamical models is furthermore about a factor of 2 smaller than the lytically predicted value of approximately 50% of the stellar equatorial rotation speed. As a result, the dynamical disks are somewhat weaker than predicted, with a wider opening angle, lower disk/pole density ratio, and smaller shock velocity jump. The principal cause of these latter differences appears to be an artificially strong equatorward drift of the subsonic outflow in the original analytic model. Much better agreement with the dynamical results can be obtained, however, from a slightly modified, analytic wind-compression model with a more detailed specification of the fixed wind outflow and a lower boundary set to the sonic radius along a rotationally oblate stellar surface. Hence, despite these detailed differences, the general predicted effect of disk formation by wind compression toward the equator is substantially confirmed.

Owocki, Stanley P.

Fluctuations at the blue edge of saturated wind lines in IUE spectra of O-type stars

We examine basic issues involved in synthesizing resonance-line profiles from 1-D, dynamical models of highly structured hot-star winds. Although these models exhibit extensive variations in density as well as velocity, the density scale length is still typically much greater than the Sobolev length. The line transfer is thus treated using a Sobolev approach, as generalized by Rybicki & Hummer (1978) to take proper account of the multiple Sobolev resonances arising from the nonmonotonic velocity field. The resulting reduced-Lambda-matrix equation describing nonlocal coupling of the source function is solved by iteration, and line profiles and then derived from formal solution integration using this source function. The more appropriate methods that instead use either a stationary or a structured, local source function yield qualitatively similar line-profiles, but are found to violate photon conservation by 10 percent or more. The full results suggest that such models may indeed be able to reproduce naturally some of the qualitative properties long noted in observed UV line profiles, such as discrete absorption components in unsaturated lines, or the blue-edge variability in saturated lines. However, these particular models do not yet produce the black absorption troughs commonly observed in saturated lines, and it seems that this and other important discrepancies (e.g., in acceleration time scale of absorption components) may require development of more complete models that include rotation and other 2-D and/or 3-D effects.

Owocki, Stanley P.

Radiation driven winds of hot stars: Some remarks on stationary models and spectrum synthesis in time-dependent simulations

The basic concept of some stationary models and the way in which these models are used for and in how far they influence the determination of stellar and wind parameters of massive hot stars are discussed. However, as follows directly from the assumption of stationarity, these models are inherently incapable of describing a number of additional observational features, which immediately show that the nonstationary aspects of the wind must be of significant importance. A number of these features are discussed: the soft X-ray emission; the discrete absorption components; the temporal variability of the wind lines, where the blue edges are varying most significantly while the red emission part remains relatively constant; the electron scattering wings of recombination lines in Wolf-Rayet stars; the black troughs in saturated P-Cygni profiles; the variability of optical lines; and the nonthermal radio emission observed in massive stars. The basic picture resulting from the present dynamical models is discussed.

Puls, Joachim

X ray emission from dynamical shock models in hot-star winds

The principal aim of this project was to determine whether x ray emission from instability-generated shocks in dynamical models of highly unstable hot-star winds could explain the x ray flux spectrum observed from such hot stars by Einstein and other x ray satellites. Our initial efforts focused on extending the earlier isothermal simulations of wind instabilities to include an explicit treatment of the energy balance between shock heating and simplified radiative cooling. It was found, however, that direct resolution of cooling regions behind shocks is often impractical, and thus additional, indirect methods for determining this shock x ray emission were also developed. The results indicate that the reverse shocks that dominate simple 1-D instability models typically have too little material undergoing a strong shock to produce the observed x ray emission. Other models with more strongly driven variability from the wind base sometimes show high-speed collisions between relatively dense clumps, and in these instances the computed x ray flux spectrum matches the observed spectrum quite well. This suggests that collisions between relatively large scale wind streams of different speeds may be more suited to producing the observed x rays than the reverse shocks arising from small-scale instabilities.

Owocki, Stanley P.

The effect of viscosity on steady transonic flow with a nodal solution topology

The effect of viscosity on a steady, transonic flow for which the inviscid limit has a nodal solution topology near the critical point is investigated. For the accelerating case, viscous solutions tend to repel each other, so that a very delicate choice of initial conditions is required to prevent them from diverging. Only the two critical solutions extend to arbitrarily large distances into both the subsonic and supersonic flows. For the decelerating case, the solutions tend to attract, and so an entire two-parameter family of solutions now extends over large distances. The general effect of viscosity on the solution degeneracy of a nodal topology is thus to reduce or limit it for the accelerating case and to enhance it for the decelerating case. The astrophysical implications of these findings are addressed.

Owocki, Stanley P.

Instabilities of line-driven stellar winds. V - Effect of an optically thick continuum

Earlier analyses of the linear instability of line-driven stellar winds are extended to the case, relevant to Wolf-Rayet stars, in which the continuum remains optically thick well above the sonic point. It is found that an optically thick flow driven by pure scattering lines is stabilized by the drag effect of the diffuse, scattered radiation. However, even a relatively small photon destruction probability can cause a flow with continuum optical thickness much greater than 1 to remain unstable, with a given growth rate. The implications of these results for the variability characteristics of winds from Wolf-Rayet stars are briefly discussed.

Owocki, Stanley P.