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.
We have updated our publicly available dust radiative transfer code (Hochunk3d) to include new emission processes and various three-dimensional (3D) geometries appropriate for forming stars. The 3D ...geometries include warps and spirals in disks, accretion hotspots on the central star, fractal clumping density enhancements, and misaligned inner disks. Additional axisymmetric (2D) features include gaps in disks and envelopes, "puffed-up inner rims" in disks, multiple bipolar cavity walls, and iteration of disk vertical structure assuming hydrostatic equilibrium (HSEQ). We include the option for simple power-law envelope geometry, which, combined with fractal clumping and bipolar cavities, can be used to model evolved stars as well as protostars. We include non-thermal emission from polycyclic aromatic hydrocarbons (PAHs) and very small grains, and external illumination from the interstellar radiation field. The grid structure was modified to allow multiple dust species in each cell; based on this, a simple prescription is implemented to model dust stratification. We describe these features in detail, and show example calculations of each. Some of the more interesting results include the following: (1) outflow cavities may be more clumpy than infalling envelopes. (2) PAH emission in high-mass stars may be a better indicator of evolutionary stage than the broadband spectral energy distribution slope; and related to this, (3) externally illuminated clumps and high-mass stars in optically thin clouds can masquerade as young stellar objects. (4) Our HSEQ models suggest that dust settling is likely ubiquitous in T Tauri disks, in agreement with previous observations.
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.
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, ..., and exponent, n. We found that the n values range mainly between 1.5 and 3.5, and ... varies between 10 super( -12) and 10 super( -10) g cm super( -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 super( -12) and 10 super( -9) M yr super( -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 ... 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 time-scale for disc growth, during which the disc is being actively fed by mass injection episodes, is shorter than the time-scale 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. (ProQuest: ... denotes formulae/symbols omitted.)
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.)