Search

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, $\eta_{\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 $\sigma$ Ori E., Comment: 5 pages, 3 figures, accepted by ApJL

This paper describes ZEUS-MP, a multi-physics, massively parallel, message- passing implementation of the ZEUS code. ZEUS-MP differs significantly from the ZEUS-2D code, the ZEUS-3D code, and an early "version 1" of ZEUS-MP distributed publicly in 1999. ZEUS-MP offers an MHD algorithm better suited for multidimensional flows than the ZEUS-2D module by virtue of modifications to the Method of Characteristics scheme first suggested by Hawley and Stone (1995), and is shown to compare quite favorably to the TVD scheme described by Ryu et. al (1998). ZEUS-MP is the first publicly-available ZEUS code to allow the advection of multiple chemical (or nuclear) species. Radiation hydrodynamic simulations are enabled via an implicit flux-limited radiation diffusion (FLD) module. The hydrodynamic, MHD, and FLD modules may be used in one, two, or three space dimensions. Self gravity may be included either through the assumption of a GM/r potential or a solution of Poisson's equation using one of three linear solver packages (conjugate-gradient, multigrid, and FFT) provided for that purpose. Point-mass potentials are also supported. Because ZEUS-MP is designed for simulations on parallel computing platforms, considerable attention is paid to the parallel performance characteristics of each module. Strong-scaling tests involving pure hydrodynamics (with and without self-gravity), MHD, and RHD are performed in which large problems (256^3 zones) are distributed among as many as 1024 processors of an IBM SP3. Parallel efficiency is a strong function of the amount of communication required between processors in a given algorithm, but all modules are shown to scale well on up to 1024 processors for the chosen fixed problem size., Comment: Accepted for publication in the ApJ Supplement. 42 pages with 29 inlined figures; uses emulateapj.sty. Discussions in sections 2 - 4 improved per referee comments; several figures modified to illustrate grid resolution. ZEUS-MP source code and documentation available from the Laboratory for Computational Astrophysics at http://lca.ucsd.edu/codes/currentcodes/zeusmp2/

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., Comment: 52 pages, 14 figures (1 color), 6 tables. To appear in the Astrophysical Journal, 1 August 2005, v628, issue 2. New version corrects e-mail address, figure and table formatting problems

This talk summarizes results from recent MHD simulations of the role of a dipole magnetic field in inducing large-scale structure in the line-driven stellar winds of hot, luminous stars. Unlike previous fixed-field analyses, the MHD simulations here take full account of the dynamical competition between the field and the flow. 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', $\eta_\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, $\eta_\ast \le 1$, the field is fully opened by wind outflow, but nonetheless, for confinement as small as $\eta_\ast=1/10$ it can have significant back-influence in enhancing the density and reducing the flow speed near the magnetic equator. For stronger confinement, $\eta_\ast > 1$, the magnetic field remains closed over limited range of latitude and height above the equatorial surface, but eventually is opened into nearly radial configuration at large radii. Within the closed loops, the flow is channeled toward loop tops into shock collisions that are strong enough to produce hard X-rays. Within the open field region, the equatorial channeling leads to oblique shocks that are again strong enough to produce X-rays and also lead to a thin, dense, slowly outflowing ``disk'' at the magnetic equator., Comment: To be published in the proceedings of International Conference on magnetic fields in O, B and A stars, ASP Conference Series, Vol. 216, 1003, L.A. Balona, H. Henrichs & T. Medupe, eds.; for alternate upload go to: http://wonka.physics.ncsu.edu/~auddoul/publications/udDoula.ps

We summarize recent 2D MHD simulations of line-driven stellar winds from rotating hot-stars with a dipole magnetic field aligned to the star's rotation axis. For moderate to strong fields, much wind outflow is initially along closed magnetic loops that nearly corotate as a solid body with the underlying star, thus providing a torque that results in an effective angular momentum spin-up of the outflowing material. But instead of forming the ``magnetically torqued disk'' (MTD) postulated in previous phenemenological analyses, the dynamical simulations here show that material trapped near the tops of such closed loops tends either to fall back or break out, depending on whether it is below or above the Keplerian corotation radius. Overall the results raise serious questions about whether magnetic torquing of a wind outflow could naturally result in a Keplerian circumstellar disk. However, for very strong fields, it does still seem possible to form a %``magnetically confined, centrifugally supported, rigid-body disk'', centrifugally supported, ``magnetically rigid disk'' (MRD), in which the field not only forces material to maintain a rigid-body rotation, but for some extended period also holds it down against the outward centrifugal force at the loop tops. We argue that such rigid-body disks seem ill-suited to explain the disk emission from Be stars, but could provide a quite attractive paradigm for circumstellar emission from the magnetically strong Bp and Ap stars., Comment: To be published in the proceedings of International Conference on magnetic fields in O, B and A stars, ASP Conference Series, Vol. 216, 1003, L.A. Balona, H. Henrichs & T. Medupe, eds.; for alternate upload go to: http://www.bartol.udel.edu/~owocki/preprints/Owocki-SA.pdf

We carry out an extended analytic study of how the tilt and faster-than-radial expansion from a magnetic field affect the mass flux and flow speed of a line-driven stellar wind. A key motivation is to reconcile results of numerical MHD simulations with previous analyses that had predicted non-spherical expansion would lead to a strong speed enhancement. By including finite-disk correction effects, a dynamically more consistent form for the non-spherical expansion, and a moderate value of the line-driving power index $\alpha$, we infer more modest speed enhancements that are in good quantitative agreement with MHD simulations, and also are more consistent with observational results. Our analysis also explains simulation results that show the latitudinal variation of the surface mass flux scales with the square of the cosine of the local tilt angle between the magnetic field and the radial direction. Finally, we present a perturbation analysis of the effects of a finite gas pressure on the wind mass loss rate and flow speed in both spherical and magnetic wind models, showing that these scale with the ratio of the sound speed to surface escape speed, $a/v_{esc}$, and are typically 10-20% compared to an idealized, zero-gas-pressure model., Comment: Accepted for publication in ApJ, for the full version of the paper go to: http://www.bartol.udel.edu/~owocki/preprints/btiltdiv-mdotvinf.pdf

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., Comment: 4 pages, 2 figures, to appear in Review of Scientific Instruments, March 2003

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', $\eta_{\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 $\eta_{\ast} \le 1$, the field is fully opened by the wind outflow, but nonetheless for confinements as small as $\eta_{\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 $\eta_{\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., Comment: 39 pages, 9 figures, for better quality figures visit: http://www.bartol.udel.edu/~owocki

High-resolution X-ray spectra of high-mass stars and low-mass T-Tauri stars obtained during the first year of the Chandra mission are providing important clues about the mechanisms which produce X-rays on very young stars. For zeta Puppis (O4 If) and zeta Ori (O9.5 I), the broad, blue-shifted line profiles, line ratios, and derived temperature distribution suggest that the X-rays are produced throughout the wind via instability-driven wind shocks. For some less luminous OB stars, like theta^1 Ori C (O7 V) and tau Sco (B0 V), the line profiles are symmetric and narrower. The presence of time-variable emission and very high-temperature lines in theta^1 Ori C and tau Sco suggest that magnetically confined wind shocks may be at work. The grating spectrum of the classical T-Tauri star TW Hya is remarkable because the forbidden-line emission of He-like Ne IX and O VII is very weak, implying that the X-ray emitting region is very dense, n = 6E+12 cgs, or that the X-rays are produced very close to the ultraviolet hotspot at the base of an accretion funnel. ACIS light curves and spectra of flares and low-mass and high-mass young stellar objects in Orion and rho Ophiuchus further suggest that extreme magnetic activity is a general property of many very young stars., Comment: 12 pages, 5 figures, to appear in "The X-ray Universe at Sharp Focus", eds. S. Vrtilek, E. Schlegel and L. Kuhi, ASP Conference Series References updated 2001 Sept 17

In this paper, we present results from the spectropolarimetric follow-up of photometrically selected candidate magnetic B stars from the MOBSTER project. Out of four observed targets, one (HD 38170) is found to host a detectable surface magnetic field, with a maximum longitudinal field measurement of 105±14 G. This star is chemically peculiar and classified as a^2 CVn variable. Its detection validates the use of TESS to perform a photometric selection of magnetic candidates. Furthermore, upper limits on the strength of a putative dipolar magnetic field are derived for the remaining three stars, and we report the discovery of a previously unknown spectroscopic binary system, HD 25709. Finally, we use our non-detections as case studies to further inform the criteria to be used for the selection of a larger sample of stars to be followed up using high-resolution spectropolarimetry.

We conducted an observational campaign towards one of the most massive and luminous colliding wind binaries in the Galaxy, HD 93129A, close to its periastron passage in 2018. During this time the source was predicted to be in its maximum of high-energy emission. Here we present our data analysis from the X-ray satellites Chandra and NuSTAR and the γ-ray satellite AGILE. High-energy emission coincident with HD 93129A was detected in the X-ray band up to 18 keV, whereas in the γ-ray band only upper limits were obtained. We interpret the derived fluxes using a non-thermal radiative model for the wind-collision region. We establish a conservative upper limit for the fraction of the wind kinetic power that is converted into relativistic electron acceleration, f(NT,e) < 0.02. In addition, we set a lower limit for the magnetic field in the wind-collision region as BWCR > 0.3 G. We also argue a putative interpretation of the emission from which we estimate f(NT,e) ≈ 0.006 and B(WCR) ≈ 0.5 G. We conclude that multi-wavelength, dedicated observing campaigns during carefully selected epochs are a powerful tool for characterising the relativistic particle content and magnetic field intensity in colliding wind binaries.