No Arabic abstract
The first results from a near-contemporaneous optical and infrared spectroscopic observing program designed to probe the detailed density structure of classical Be circumstellar disks are presented. We report the discovery of asymmetrical infrared emission lines of He I, O I, Fe II, and the Brackett, Paschen, and Pfund series lines of H I which exhibit an opposite V/R orientation (V $>$ R) to that observed for the optical Balmer H$alpha$ line (V $<$ R) in the classical Be star $zeta$ Tau. We interpret these data as evidence that the density wave which characterizes $zeta$ Taus disk has a significantly different average azimuthal morphology in the inner disk region as compared to the outer disk region. A follow-up multi-wavelength observational campaign to trace the temporal evolution of these line profile morphologies, along with detailed theoretical modeling, is suggested to test this hypothesis.
The extreme luminosities of hot, massive stars drive strong stellar winds through UV line-scattering. For OB stars with an orbiting circumstellar disk, we explore the effect of such line-scattering in ablating disk material, initially focusing on the marginally optically thin decretion disks of classical Oe and Be stars. For this we apply a multi-dimensional radiation-hydrodynamics code, assuming optically thin ray tracing for the stellar continuum and a multi-ray Sobolev treatment of the line transfer. This accounts for desaturation of line-absorption by Keplerian shear in the disk, and associated driving by non-radial photons. Results show dense, intermediate-speed surface ablation, consistent with the strong, blue-shifted absorption seen in UV wind lines of Be shell stars. The asymptotic ablation rate is typically an order-unity factor times the stellar wind mass loss rate, leading to disk destruction times of order months to years for Be disks, consistent with observations. The much stronger radiative forces of O stars reduce this time to order days, making sustaining a disk difficult, and so providing a natural explanation for the rarity of Galactic Oe stars. Additionally, the weakened line-driving at lower metallicity implies both a reduction in the winds that help spin-down stars from near-critical rotation, and in the ablation of decretion disks, thus providing a natural explanation for the higher fraction of Classical Be stars, and the presence of Oe stars, in the Magellanic Clouds. We conclude with a discussion of future extensions to study line-driven ablation of denser, optically thick, accretion disks around pre-main-sequence massive stars.
We completed a search for circumstellar disks around Herbig Be stars using the NRAO Very Large Array (VLA) and the IRAM Plateau de Bure (PdB) interferometers. We present our new VLA and PdBI data for the three objects MWC 297, Z CMa, and LKHa 215. We constructed the SED from near-IR to centimeter wavelengths by adding our millimeter and centimeter data to the available data at other wavelengths, mainly Spitzer images. The entire SED was fitted using a disk+envelope model. In addition, we compiled all the disk millimeter observations in the literature and completed a statistical analysis of all the data. We show that the disk mass is usually only a small percentage (less than 10%) of the mass of the entire envelope in HBe stars. For the disks, there are large source-to-source variations. Two disks in our sample, R Mon and Z CMa, have similar sizes and masses to those found in T Tauri and Herbig Ae stars. The disks around MWC 1080 and MWC 297 are, however, smaller (rout<100 AU). We did not detect the disks towards MWC 137 and LkHa 215 at millimeter wavelengths, which limits the mass and the size of the possible circumstellar disks. A comparison between our data and previous results for T Tauri and Herbig Ae stars indicates that although massive disks (0.1 Msun) are found in young objects (10^4 yr), the masses of the disks around Herbig Be stars are usually 5-10 times lower than those around lower mass stars. We propose that disk photoevaporation is responsible for this behavior. In Herbig Be stars, the UV radiation disperses the gas in the outer disk on a timescale of a few 10^5 yr. Once the outer part of the disk has vanished, the entire gaseous disk is photoevaporated on a very short timescale (10^5 yr) and only a small, dusty disk consisting of large grains remains.
We analyze the intrinsic polarization of two classical Be stars in the process of losing their circumstellar disks via a Be to normal B star transition originally reported by Wisniewski et al. During each of five polarimetric outbursts which interrupt these disk-loss events, we find that the ratio of the polarization across the Balmer jump (BJ+/BJ-) versus the V-band polarization traces a distinct loop structure as a function of time. Since the polarization change across the Balmer jump is a tracer of the innermost disk density whereas the V-band polarization is a tracer of the total scattering mass of the disk, we suggest such correlated loop structures in Balmer jump-V band polarization diagrams (BJV diagrams) provide a unique diagnostic of the radial distribution of mass within Be disks. We use the 3-D Monte Carlo radiation transfer code HDUST to reproduce the observed clockwise loops simply by turning on/off the mass decretion from the disk. We speculate that counter-clockwise loop structures we observe in BJV diagrams might be caused by the mass decretion rate changing between subsequent on/off sequences. Applying this new diagnostic to a larger sample of Be disk systems will provide insight into the time-dependent nature of each systems stellar decretion rate.
Context. The structure of the inner parts of Be star disks (20 stellar radii) is well explained by the viscous decretion disk (VDD) model, which is able to reproduce the observable properties of most of the objects studied so far. The outer parts, on the ther hand, are not observationally well-explored, as they are observable only at radio wavelengths. A steepening of the spectral slope somewhere between infrared and radio wavelengths was reported for several Be stars that were previously detected in the radio, but a convincing physical explanation for this trend has not yet been provided. Aims. We test the VDD model predictions for the extended parts of a sample of six Be disks that have been observed in the radio to address the question of whether the observed turndown in the spectral energy distribution (SED) can be explained in the framework of the VDD model, including recent theoretical development for truncated Be disks in binary systems. Methods. We combine new multi-wavelength radio observations from the Karl. G. Jansky Very Large Array (JVLA) and Atacama Pathfinder Experiment (APEX) with previously published radio data and archival SED measurements at ultraviolet, visual, and infrared wavelengths. The density structure of the disks, including their outer parts, is constrained by radiative transfer modeling of the observed spectrum using VDD model predictions. In the VDD model we include the presumed effects of possible tidal influence from faint binary companions. Results. For 5 out of 6 studied stars, the observed SED shows strong signs of SED turndown between far-IR and radio wavelengths. A VDD model that extends to large distances closely reproduces the observed SEDs up to far IR wavelengths, but fails to reproduce the radio SED. ... (abstract continues but did not fit here)
We report on a high angular resolution survey of circumstellar disks around 24 northern sky Be stars. The K-band continuum survey was made using the CHARA Array long baseline interferometer (baselines of 30 to 331 m). The interferometric visibilities were corrected for the flux contribution of stellar companions in those cases where the Be star is a member of a known binary or multiple system. For those targets with good uv coverage, we used a four parameter Gaussian elliptical disk model to fit the visibilities and to determine the axial ratio, position angle, K-band photospheric flux contribution, and angular diameter of the disk major axis. For the other targets with relatively limited uv coverage, we constrained the axial ratio, inclination angle, and or disk position angle where necessary in order to resolve the degeneracy between possible model solutions. We also made fits of the ultraviolet and infrared spectral energy distributions to estimate the stellar angular diameter and infrared flux excess of each target. The mean ratio of the disk diameter (measured in K-band emission) to stellar diameter (from SED modeling) is 4.4 among the 14 cases where we reliably resolved the disk emission, a value which is generally lower than the disk size ratio measured in the higher opacity Halpha emission line. We estimated the equatorial rotational velocity from the projected rotational velocity and disk inclination for 12 stars, and most of these stars rotate close to or at the critical rotational velocity.