Abstract We examine the evolution of the disk surrounding the Be star in the highly eccentric 10.8 yr binary system δ Scorpii over its three most recent periastron passages. V -band and B − V ...photometry, along with H α spectroscopy, are combined with a new set of extensive multiband polarimetry data to produce a detailed comparison of the disk's physical conditions during the time periods surrounding each closest approach of the secondary star. We use the three-dimensional Monte Carlo radiative transfer code HDUST and a smoothed particle hydrodynamics code to support our observations with models of disk evolution, discussing the behavior of the H α and He i 6678 Å lines, V -band magnitude, and polarization degree. We compare the characteristics of the disk immediately before each periastron passage to create a baseline for the unperturbed disk. We find that the extent of the H α emitting region increased between each periastron passage, and that transient asymmetries in the disk become more pronounced with each successive encounter. Asymmetries of the H α and He i 6678 Å lines in 2011 indicate that perturbations propagate inward through the disk near periastron. When the disk’s direction of orbit is opposite to that of the secondary, the parameters used in our models do not produce spiral density enhancements in the H α emitting region because the tidal interaction time is short due to the relative velocities of the disk particles with the secondary. The effects of the secondary star on the disk are short-lived and the disk shows independent evolution between each periastron event.
ABSTRACT
We use a three-dimensional smoothed particle hydrodynamics code to simulate growth and dissipation of Be star discs in systems where the binary orbit is misaligned with respect to the spin ...axis of the primary star. We investigate six different scenarios of varying orbital period and misalignment angle, feeding the disc at a constant rate for 100 orbital periods, and then letting the disc dissipate for 100 orbital periods. During the disc growth phase, we find that the binary companion tilts the disc away from its initial plane at the equator of the primary star before settling to a constant orientation after 40–50 orbital periods. While the mass-injection into the disc is ongoing, the tilting of the disc can cause material to reaccrete on to the primary star prematurely. Once disc dissipation begins, usually the disc precesses about the binary companion’s orbital axis with precession periods ranging from 20 to 50 orbital periods. In special cases, we detect phenomena of disc tearing, as well as Kozai–Lidov oscillations of the disc. These oscillations reach a maximum eccentricity of about 0.6, and a minimum inclination of about 20○ with respect to the binary’s orbit. We also find the disc material to have highly eccentric orbits beyond the transition radius, where the disc changes from being dominated by viscous forces, to heavily controlled by the companion star, in contrast to its nearly circular motion inwards of the transition radius. Finally, we offer predictions to how these changes will affect Be star observables.
We present a three-dimensional non-LTE Monte Carlo radiative transfer code that we use to study the temperature and ionization structure of Keplerian disks around classical Be stars. The method we ...employ is largely similar to the Monte Carlo transition probability method developed by Lucy. Here we present a simplification of his method that avoids the use of the macroatom concept. Our investigations of the temperature structure of Be star disks show that the disk temperature behavior is a hybrid between the behavior of young stellar object (YSO) disks and hot star winds. The optically thick inner parts of Be star disks have temperatures that are similar to YSO disks, while the optically thin outer parts are like stellar winds. Thus, the temperature at the disk midplane initially drops, reaching a minimum at 3-5 stellar radii, after which it rises back to the optically thin radiative equilibrium temperature at large distances. On the other hand, the optically thin upper layers of the disk are approximately isothermal--a behavior that is analogous to the hot upper layers of YSO disks. Interestingly, unlike the case of YSO disks, we find that disk flaring has little effect on the temperature structure of Be star disks. We also find that the disks are fully ionized, as expected, but that there is an ionization minimum in the vicinity of the temperature minimum. The deficit of photoionization at this location makes it the most likely site for the low ionization state lines (e.g., Fe II) that produce the shell features observed in Be stars. Finally, we find that despite the complex temperature structure, the infrared excess is well approximated by an equivalent isothermal disk model whose temperature is about 60% of the stellar temperature. This is largely because at long wavelengths, the effective photosphere of the disk is located in its isothermal regions.
Rapid rotation is a fundamental characteristic of classical Be stars and a crucial property allowing for the formation of their circumstellar disks. Past evolution in a mass and angular momentum ...transferring binary system offers a plausible solution to how Be stars attained their fast rotation. Although the subdwarf remnants of mass donors in such systems should exist in abundance, only a few have been confirmed due to tight observational constraints. An indirect method of detecting otherwise hidden companions is offered by their effect on the outer parts of Be star disks, which are expected to be disrupted or truncated. In the context of the infrared and radio continuum excess radiation originating in the disk, the disk truncation can be revealed by a turndown in the spectral energy distribution due to reduced radio flux levels. In this work, we search for signs of spectral turndown in a sample of 57 classical Be stars with radio data, which include new data for 23 stars and the longest-wavelength detections so far (λ 10 cm) for two stars. We confidently detect the turndown for all 26 stars with sufficient data coverage (20 of which are not known to have close binary companions). For the remaining 31 stars, the data are inconclusive as to whether the turndown is present or not. The analysis suggests that many if not all Be stars have close companions influencing their outer disks. If confirmed to be subdwarf companions, the mass transfer spin-up scenario might explain the existence of the vast majority of classical Be stars.
Be stars possess gaseous circumstellar disks that modify in many ways the spectrum of the central B star. Furthermore, they exhibit variability at several timescales and for a large number of ...observables. Putting the pieces together of this dynamical behavior is not an easy task and requires a detailed understanding of the physical processes that control the temporal evolution of the observables. There is an increasing body of evidence that suggests that Be disks are well described by standard alpha -disk theory. This paper is the first of a series that aims at studying the possibility of inferring several disk and stellar parameters through the follow-up of various observables. Here we study the temporal evolution of the disk density for different dynamical scenarios, including the disk buildup as a result of a long and steady mass injection from the star, the disk dissipation that occurs after mass injection is turned off, as well as scenarios in which active periods are followed by periods of quiescence. For those scenarios, we investigate the temporal evolution of continuum photometric observables using a three-dimensional non-LTE radiative transfer code. We show that light curves for different wavelengths are specific of a mass loss history, inclination angle, and alpha viscosity parameter. The diagnostic potential of those light curves is also discussed.
Context. Classical Be stars have been established as pulsating stars. Space-based photometric monitoring missions contributed significantly to that result. However, whether Be stars are just rapidly ...rotating SPB or beta Cep stars, or whether they have to be understood differently, remains debated in the view of their highly complex power spectra. Aims. Kepler data of three known Be stars are re-visited to establish their pulsational nature and assess the properties of additional, non-pulsational variations. The three program stars turned out to be one inactive Be star, one active, continuously outbursting Be star, and one Be star transiting from a non-outbursting into an outbursting phase, thus forming an excellent sample to distill properties of Be stars in the various phases of their life-cycle. Methods. The Kepler data was first cleaned from any long-term variability with Lomb-Scargle based pre-whitening. Then a Lomb-Scargle analysis of the remaining short-term variations was compared to a wavelet analysis of the cleaned data. This offers a new view on the variability, as it enables us to see the temporal evolution of the variability and phase relations between supposed beating phenomena, which are typically not visualized in a Lomb-Scargle analysis. Results. The short-term photometric variability of Be stars must be disentangled into a stellar and a circumstellar part. The stellar part is on the whole not different from what is seen in non-Be stars. However, some of the observed phenomena might be to be due to resonant mode coupling, a mechanism not typically considered for B-type stars. Short-term circumstellar variability comes in the form of either a group of relatively well-defined, short-lived frequencies during outbursts, which are called Stefl frequencies, and broad bumps in the power spectra, indicating aperiodic variability on a time scale similar to typical low-order g-mode pulsation frequencies, rather than true periodicity. Conclusions. From a stellar pulsation perspective, Be stars are rapidly rotating SPB stars, that is they pulsate in low order g-modes, even if the rapid rotation can project the observed frequencies into the traditional high-order p-mode regime above about 4c/d. However, when a circumstellar disk is present, Be star power spectra are complicated by both cyclic, or periodic, and aperiodic circumstellar phenomena, possibly even dominating the power spectrum.
ABSTRACT
We build on our previous work involving smoothed particle hydrodynamic simulations of Be stars, by using the model that exhibited disc tearing as input into the three-dimensional Monte Carlo ...radiative transfer code hdust to predict observables from a variety of viewing angles throughout the disc tearing process. We run one simulation at the start of each orbital period from 20 to 72 orbital periods, which covers two complete disc tearing events. The resulting trends in observables are found to be dependent on the relative position of the observer and the tearing disc. The H $\rm \alpha$ equivalent width, V magnitude, and polarization can all increase or decrease in any combination depending on the viewpoint of the observer. The H $\rm \alpha$ line profile also displays changes in strength and peak separation throughout the tearing process. We show how the outer disc of the torn system can have a large effect on the H $\rm \alpha$ line profile, and also contributes to a wavelength-dependent polarization position angle, resulting in a similar sawtooth shape to the polarization percentage. Finally, we compare our predictions to Pleione (28 Tau) where evidence has suggested that a disc tearing event has occurred in the past. We find that our tearing disc model can broadly match the trends seen in Pleione’s observables, as well as produce the two-component H $\rm \alpha$ lines observed in Pleione. This is the strongest evidence, thus far, of Pleione’s disc having indeed experienced a tearing event.
We discuss the basic hydrodynamics that determines the density structure of the disks around hot stars. Observational evidence supports the idea that these disks are Keplerian (rotationally ...supported) gaseous disks. A popular scenario in the literature, which naturally leads to the formation of Keplerian disks, is the viscous decretion model. According to this scenario, the disks are hydrostatically supported in the vertical direction, while the radial structure is governed by the viscous transport. This suggests that the temperature is one primary factor that governs the disk density structure. In a previous study we demonstrated, using three-dimensional non-LTE Monte Carlo simulations, that viscous Keplerian disks can be highly nonisothermal. In this paper we build on our previous work and solve the full problem of the steady state nonisothermal viscous diffusion and vertical hydrostatic equilibrium. We find that the self-consistent solution departs significantly from the analytic isothermal density, with potentially large effects on the emergent spectrum. This implies that nonisothermal disk models must be used for a detailed modeling of Be star disks.
Abstract
We acquired H
α
spectroscopic observations from 2005 to 2019 showing Pleione has transitioned from a Be phase to a Be-shell phase during this period. Using the radiative transfer code
hdust
..., we created a grid of ∼100,000 disk models for Pleione. We successfully reproduced the observed transition with a disk model that varies in inclination while maintaining an equatorial density of
ρ
0
(
r
)
=
3
×
10
−
11
(
r
/
R
eq
)
−
2.7
g
cm
−
3
, and an H
α
-emitting region extending to 15
R
eq
. We use a precessing disk model to extrapolate the changing disk inclination over 120 yr and follow the variability in archival observations. The best-fit disk model precesses over a line-of-sight inclination between ∼25° and ∼144° with a precessional period of ∼80.5 yr. Our precessing models match some of the observed variability but fail to reproduce all of the historical data available. Therefore, we propose an ad hoc model based on our precessing disk model inspired by recent smoothed particle hydrodynamics simulations of similar systems, where the disk tears due to the tidal influence of a companion star. In this model, a single disk is slowly tilted to an angle of 30° from the stellar equator over 34 yr. Then, the disk is torn by the companion’s tidal torque, with the outer region separating from the innermost disk. The small inner disk returns to the stellar equator as mass injection remains constant. The outer disk precesses for ∼15 yr before gradually dissipating. The process repeats every 34 yr and reproduces all trends in Pleione’s variability.
We investigate the continuum emission of viscous decretion discs around Be stars in this paper. The results obtained from non-LTE (local thermodynamic equilibrium) radiative transfer models show two ...regimes in the disc surface brightness profile: an inner optically thick region, which behaves as a pseudo-photosphere with a wavelength-dependent size, and an optically thin tenuous outer part, which contributes with about a third of the total flux. The isophotal shape of the surface brightness is well described by elliptical contours with an axial ratio b/a = cos...i for inclinations i < 75... Based on these properties, a semi-analytical model was developed to describe the continuum emission of gaseous discs. It provides fluxes and spectral slopes at the infrared within an accuracy of 10 and 5 per cent, respectively, when compared to the numerical results. The model indicates that the infrared spectral slope is mainly determined by both the density radial slope and the disc flaring exponent, being practically independent of disc inclination and base density. As a first application, the density structure of 15 Be stars was investigated, based on the infrared flux excess, and the results compared to previous determinations in the literature. Our results indicate that the decretion rates are in the range of 10...-10... M...yr..., which is at least two orders of magnitude smaller than the previous outflowing disc model predictions. (ProQuest: ... denotes formulae/symbols omitted.)