No Arabic abstract
The observed emission lines of Be stars originate from a circumstellar Keplerian disk that are generally well explained by the Viscous Decretion Disk model. In an earlier work we performed the modeling of the full light curve of the bright Be star $omega$ CMa (Ghoreyshi et al. 2018) with the 1-D time-dependent hydrodynamics code SINGLEBE and the Monte Carlo radiative-transfer code HDUST. We used the V -band light curve that probes the inner disk through four disk formation and dissipation cycles. This new study compares predictions of the same set of model parameters with time-resolved photometry from the near UV through the mid-infrared, comprehensive series of optical spectra, and optical broad-band polarimetry, that overall represent a larger volume of the disk. Qualitatively, the models reproduce the trends in the observed data due to the growth and decay of the disk. However, quantitative differences exist, e.g., an overprediction of the flux increasing with wavelength, too slow decreases in Balmer emission-line strength that are too slow during disk dissipation, and the discrepancy between the range of polarimetric data and the model. We find that a larger value of the viscosity parameter alone, or a truncated disk by a companion star, reduces these discrepancies by increasing the dissipation rate in the outer regions of the disk.
We present interferometric observations of the Be star Zeta Tau obtained using the MIRC beam combiner at the CHARA Array. We resolved the disk during four epochs in 2007-2009. We fit the data with a geometric model to characterize the circumstellar disk as a skewed elliptical Gaussian and the central Be star as a uniform disk. The visibilities reveal a nearly edge-on disk with a FWHM major axis of ~ 1.8 mas in the H-band. The non-zero closure phases indicate an asymmetry within the disk. Interestingly, when combining our results with previously published interferometric observations of Zeta Tau, we find a correlation between the position angle of the disk and the spectroscopic V/R ratio, suggesting that the tilt of the disk is precessing. This work is part of a multi-year monitoring campaign to investigate the development and outward motion of asymmetric structures in the disks of Be stars.
We use a time-dependent hydrodynamic code and a non-LTE Monte Carlo code to model disk dissipation for the Be star 66 Ophiuchi. We compiled 63 years of observations from 1957 to 2020 to encompass the complete history of the growth and subsequent dissipation of the stars disk. Our models are constrained by new and archival photometry, spectroscopy and polarization observations, allowing us to model the disk dissipation event. Using Markov chain Monte Carlo methods, we find 66 Oph is consistent with standard B2Ve stellar properties. We computed a grid of 61568 Be star disk models to constrain the density profile of the disk before dissipation using observations of the H$alpha$ line profile and SED. We find at the onset of dissipation the disk has a base density of $2.5times10^{-11} rm{g cm^{-3}}$ with a radial power-law index of $n=2.6$. Our models indicate that after 21 years of disk dissipation 66 Ophs outer disk remained present and bright in the radio. We find an isothermal disk with constant viscosity with an $alpha = 0.4$ and an outer disk radius of $sim$115 stellar radii best reproduces the rate of 66 Ophs disk dissipation. We determined the interstellar polarization in the direction of the star in the V-band is $p=0.63 pm 0.02%$ with a polarization position angle of $theta_{IS}approx85.7 pm 0.7^circ$. Using the Stokes QU diagram, we find the intrinsic polarization position angle of 66 Ophs disk is $theta_{int}approx98 pm 3^circ$.
The viscous decretion disk (VDD) model is able to explain most of the currently observable properties of the circumstellar disks of Be stars. However, more stringent tests, focusing on reproducing multitechnique observations of individual targets via physical modeling, are needed to study the predictions of the VDD model under specific circumstances. In the case of nearby, bright Be star $beta$ CMi, these circumstances are a very stable low-density disk and a late-type (B8Ve) central star. The aim is to test the VDD model thoroughly, exploiting the full diagnostic potential of individual types of observations, in particular, to constrain the poorly known structure of the outer disk if possible, and to test truncation effects caused by a possible binary companion using radio observations. We use the Monte Carlo radiative transfer code HDUST to produce model observables, which we compare with a very large set of multitechnique and multiwavelength observations that include ultraviolet and optical spectra, photometry covering the interval between optical and radio wavelengths, optical polarimetry, and optical and near-IR (spectro)interferometry. Due to the absence of large scale variability, data from different epochs can be combined into a single dataset. A parametric VDD model with radial density exponent of $n$ = 3.5, which is the canonical value for isothermal flaring disks, is found to explain observables typically formed in the inner disk, while observables originating in the more extended parts favor a shallower, $n$ = 3.0, density falloff. Modeling of radio observations allowed for the first determination of the physical extent of a Be disk (35$^{+10}_{-5}$ stellar radii), which might be caused by a binary companion. Finally, polarization data allowed for an indirect measurement of the rotation rate of the star, which was found to be $W gtrsim 0.98$, i.e., very close to critical.
The Cepheus B (CepB) molecular cloud and a portion of the nearby CepOB3b OB association, one of the most active regions of star formation within 1 kpc, have been observed with the IRAC detector on board the Spitzer Space Telescope. The goals are to study protoplanetary disk evolution and processes of sequential triggered star formation in the region. Out of ~400 pre-main sequence (PMS) stars selected with an earlier Chandra X-ray Observatory observation, 95% are identified with mid-infrared sources and most of these are classified as diskless or disk-bearing stars. The discovery of the additional >200 IR-excess low-mass members gives a combined Chandra+Spitzer PMS sample complete down to 0.5 Mo outside of the cloud, and somewhat above 1 Mo in the cloud. Analyses of the nearly disk-unbiased combined Chandra+Spitzer selected stellar sample give several results. Our major finding is a spatio-temporal gradient of young stars from the hot molecular core towards the primary ionizing O star HD 217086. This strongly supports the radiation driven implosion (RDI) model of triggered star formation in the region. The empirical estimate for the shock velocity of 1 km/s is very similar to theoretical models of RDI in shocked molecular clouds...ABRIDGED... Other results include: 1. agreement of the disk fractions, their mass dependency, and fractions of transition disks with other clusters; 2. confirmation of the youthfulness of the embedded CepB cluster; 3. confirmation of the effect of suppression of time-integrated X-ray emission in disk-bearing versus diskless systems.
A global disk oscillation implemented in the viscous decretion disk (VDD) model has been used to reproduce most of the observed properties of the well known Be star $zeta$ Tau. 48 Librae shares several similarities with $zeta$ Tau -- they are both early-type Be stars, they display shell characteristics in their spectra, and they exhibit cyclic $V/R$ variations -- but has some marked differences as well, such as a much denser and more extended disk, a much longer $V/R$ cycle, and the absence of the so-called triple-peak features. We aim to reproduce the photometric, polarimetric, and spectroscopic observables of 48 Librae with a self-consistent model, and to test the global oscillation scenario for this target. Our calculations are carried out with the three-dimensional NLTE radiative transfer code HDUST. We employ a rotationally deformed, gravity-darkened central star, surrounded by a disk whose unperturbed state is given by the VDD model. A two-dimensional global oscillation code is then used to calculate the disk perturbation, and superimpose it on the unperturbed disk. A very good, self-consistent fit to the time-averaged properties of the disk is obtained with the VDD. The calculated perturbation has a period $P = 12$ yr, which agrees with the observed period, and the behaviour of the $V/R$ cycle is well reproduced by the perturbed model. The perturbed model improves the fit to the photometric data and reproduces some features of the observed spectroscopic data. Some suggestions to improve the synthesized spectroscopy in a future work are given.