No Arabic abstract
In this work we explore the effect of binarity in the decretion disc of Be stars, in order to explain their variability. To this aim, we performed smoothed particle hydrodynamics (SPH) simulations on Be binary systems, following the matter ejected isotropically from the equator of the Be star towards the base of an isothermal decretion disc. We let the system evolve for time long enough to be considered at steady state, and focus on the effect of viscosity for coplanar prograde binary orbits. The disc structure is found to be locked to the orbital phase, exhibiting also a dependence on the azimuthal angle. Additionally, we present the first results from detailed non-local thermodynamic equilibrium (non-LTE) radiative transfer calculations of the disc structure computed with the SPH code. This is achieved by the use of the three-dimensional (3D) Monte Carlo code HDUST, which can produce predictions with respect to a series of observables.
We apply the viscous decretion disc (VDD) model to interpret the infrared disc continuum emission of 80 Be stars observed in different epochs. In this way, we determined 169 specific disc structures, namely their density scale, $rho_0$, and exponent, $n$. We found that the $n$ values range mainly between $1.5$ and $3.5$, and $rho_0$ varies between $10^{-12}$ and $10^{-10},mathrm{g,cm^{-3}}$, with a peak close to the lower value. Our large sample also allowed us to firmly establish that the discs around early-type stars are denser than in late-type stars. Additionally, we estimated the disc mass decretion rates and found that they range between $10^{-12}$ and $10^{-9},mathrm{M_{odot},yr^{-1}}$. These values are compatible with recent stellar evolution models of fast-rotating stars. One of the main findings of this work is a correlation between the $rho_0$ and $n$ values. In order to find out whether these relations can be traced back to the evolution of discs or have some other origin, we used the VDD model to calculate temporal sequences under different assumptions for the time profile of the disc mass injection. The results support the hypothesis that the observed distribution of disc properties is due to a common evolutionary path. In particular, our results suggest that the timescale for disc growth, during which the disc is being actively fed by mass injection episodes, is shorter than the timescale for disc dissipation, when the disc is no longer fed by the star and dissipates as a result of the viscous diffusion of the disc material.
The first results of radiative transfer calculations on decretion discs of binary Be stars are presented. A smoothed particle hydrodynamics code computes the structure of Be discs in coplanar circular binary systems for a range of orbital and disc parameters. The resulting disc configuration consists of two spiral arms, and can be given as input into a Monte Carlo code, which calculates the radiative transfer along the line of sight for various observational coordinates. Making use of the property of steady disc structure in coplanar circular binaries, observables are computed as functions of the orbital phase. Orbital-phase series of line profiles are given for selected parameter sets under various viewing angles, to allow comparison with observations. Flat-topped profiles with and without superimposed multiple structures are reproduced, showing, for example, that triple-peaked profiles do not have to be necessarily associated with warped discs and misaligned binaries. It is demonstrated that binary tidal effects give rise to phase-locked variability of the violet-to-red (V/R) ratio of hydrogen emission lines. The V/R ratio exhibits two maxima per cycle; in certain cases those maxima are equal, leading to a clear new V/R cycle every half orbital period. This study opens a way in identifying binaries and in constraining the parameters of binary systems that exhibit phase-locked variations induced by tidal interaction with a companion star.
Evolutionary models of fast-rotating stars show that the stellar rotational velocity may approach the critical speed. Critically rotating stars cannot spin up more, therefore they lose their excess angular momentum through an equatorial outflowing disk. The radial extension of such disks is unknown, partly because we lack information about the radial variations of the viscosity. We study the magnetorotational instability, which is considered to be the origin of anomalous viscosity in outflowing disks. We used analytic calculations to study the stability of outflowing disks submerged in the magnetic field. The magnetorotational instability develops close to the star if the plasma parameter is large enough. At large radii the instability disappears in the region where the disk orbital velocity is roughly equal to the sound speed. The magnetorotational instability is a plausible source of anomalous viscosity in outflowing disks. This is also true in the region where the disk radial velocity approaches the sound speed. The disk sonic radius can therefore be roughly considered as an effective outer disk radius, although disk material may escape from the star to the insterstellar medium. The radial profile of the angular momentum-loss rate already flattens there, consequently, the disk mass-loss rate can be calculated with the sonic radius as the effective disk outer radius. We discuss a possible observation determination of the outer disk radius by using Be and Be/X-ray binaries.
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 report optical spectroscopic observations of the Be/gamma-ray binaries LSI+61303, MWC 148 and MWC 656. The peak separation and equivalent widths of prominent emission lines (H-alpha, H-beta, H-gamma, HeI, and FeII) are measured. We estimated the circumstellar disc size, compared it with separation between the components, and discussed the disc truncation. We find that in LSI+61303 the compact object comes into contact with the outer parts of the circumstellar disc at periastron, in MWC 148 the compact object goes deeply into the disc during the periastron passage, and in MWC 656 the black hole is accreting from the outer parts of the circumstellar disc along the entire orbit. The interstellar extinction was estimated using interstellar lines. The rotation of the mass donors appears to be similar to the rotation of the mass donors in Be/X-ray binaries. We suggest that X-ray/optical periodicity of about 1 day deserves to be searched for.