Massive star winds are important contributors to the energy, momentum and chemical enrichment of the interstellar medium. Strong, organized and predominantly dipolar magnetic fields have been firmly detected in a small subset of massive O-type stars. Magnetic massive stars are known to exhibit phase-locked variability of numerous observable quantities that is hypothesized to arise due to the presence of an obliquely rotating magnetosphere formed via the magnetic confinement of their strong outflowing winds. Analyzing the observed modulations of magnetic O-type stars is thus a key step towards the better understanding of the physical processes that occur within their magnetospheres. The dynamical processes that lead to the formation of a magnetosphere are formally solved utilizing complex MHD simulations. Recently, an Analytic Dynamical Magnetosphere (ADM) model has been developed that can quickly be employed to compute the time-averaged density, temperature and velocity gradients within a dynamical magnetosphere. Here, we exploit the ADM model to compute photometric and polarimetric observables of magnetic Of?p stars, to test geometric models inferred from magnetometry. We showcase important results on the prototypical Of?p-type star HD 191612, that lead to a better characterization of massive star wind and magnetic properties.

We present a detailed visible and near-IR spectro-interferometric analysis of the Be-shell star $\omicron$ Aquarii from quasi-contemporaneous CHARA/VEGA and VLTI/AMBER observations. We measured the stellar radius of $\omicron$ Aquarii as 4.0 $\pm$ 0.3 $\mathrm{R_{\odot}}$. We constrained the disk geometry and kinematics using a kinematic model and a MCMC fitting procedure. The disk sizes in H$\alpha$ and Br$\gamma$ were found to be similar, at $\sim$10-12 $\mathrm{D_{\star}}$, which is uncommon since most results for Be stars show a larger extension in H$\alpha$ than in Br$\gamma$. We found that the inclination angle $i$ derived from H$\alpha$ is significantly lower ($\sim$15 deg) than the one derived from Br$\gamma$. The disk kinematics were found to be near to the Keplerian rotation in Br$\gamma$, but not in H$\alpha$. After analyzing all our data using a grid of HDUST models (BeAtlas), we found a common physical description for the disk in both lines: $\Sigma_{0}$ = 0.12 g cm\textsuperscript{-2} and $m$ = 3.0. The stellar rotational rate was found to be very close ($\sim$96\%) to the critical value. Our analysis of multi-epoch H$\alpha$ profiles and imaging polarimetry indicates that the disk has been stable for at least 20 years. Compared to Br$\gamma$, the data in H$\alpha$ shows a substantially different picture that cannot fully be understood using the current physical models of Be star disks. $\omicron$ Aquarii presents a stable disk, but the measured $m$ is lower than the standard value in the VDD model for steady-state. Such long-term stability can be understood in terms of the high rotational rate for this star, the rate being a main source for the mass injection in the disk. Our results on the stellar rotation and disk stability are consistent with results in the literature showing that late-type Be stars are more likely to be fast rotators and have stable disks.

We model the circumstellar disk of $\delta$ Sco using the 3-dimensional Monte Carlo radiative transfer code HDUST in order to quantify the large scale changes in the disk through the years 2000 to 2018, and to see if these changes can be attributed to the secondary star affecting the disk throughout its orbit. We determine our best-fitting models through matching simulated observations to actual H$\rm \alpha$ spectroscopy and V-band photometric observations. Our modelling results confirm previous findings that the disk of $\delta$ Sco was forming early in the century. We also find a period of disk dissipation when the companion is at apastron, as well as a significant growth of the disk between 2009 and 2011, prior to the periastron of 2011. Due to the steady-state nature of the disk after 2011, it is difficult to say whether the variations seen are due to the effect of the close passage of the binary companion.

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 IR 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 ($\lambda \approx$ 10\,cm) for 2 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, 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.

The tidal interaction of a Be star with a binary companion forms two spiral arms that cause orbital modulation of the Be disc structure. The aim of this work is to identify observables in which this modulation is apparent. The structure of a Be disc in a coplanar circular binary system is computed with a smoothed-particle hydrodynamics code, and a radiation transfer code calculates the spectral energy distribution. Line depolarisation was confirmed, with polarisation profiles nearly reverse to emission-line profiles. The continuum flux maximizes for pole-on discs, but photometric variability maximizes for edge-on discs. The linear polarisation exhibits one or two maxima per orbital cycle. While polarisation variability in visible passbands is important only at low inclinations, infrared bands may demonstrate high orbital variability even at large inclinations. More evident is the modulation in the polarisation angle (PA) for low inclinations. The latter can be used to track azimuthal asymmetries for pole-on discs, where the spectroscopic variability in the violet-to-red (V/R) emission-component ratio disappears. PA reversals coincide with phases where V/R=1, tracking lines of sight directed towards regions where the approaching and receding arms overlap. Continuum flux and polarisation are mostly in phase for neighbouring wavelength regions. It is suggested that studies of non-symmetric discs distorted by tidal forces from a secondary star may be used to study disc variabilities of other origins.

In early-type Be stars, groups of nonradial pulsation (NRP) modes with numerically related frequencies may be instrumental for the release of excess angular momentum through mass-ejection events. Difference and sum/harmonic frequencies often form additional groups. The goal of this study is to find out whether a similar frequency pattern occurs in the cooler third-magnitude B7-8\,IIIe shell star $\nu$ Pup. Time-series analyses are performed of space photometry with BRITE-Constellation (2015, 2016/17, and 2017/18), SMEI (2003--011), and Hipparcos (1989-1993). Two IUE SWP and 27 optical echelle spectra spanning 20 years were retrieved from various archives. The optical spectra exhibit no anomalies or well-defined variabilities. A magnetic field was not detected. All three photometry satellites recorded variability near 0.656 c/d which is resolved into three features separated by ~0.0021 c/d. First harmonics form a second frequency group, also spaced by ~0.0021 c/d. The frequency spacing is very nearly but not exactly equidistant. Variability near 0.0021 c/d was not detected. The long-term frequency stability could be used to derive meaningful constraints on the properties of a putative companion star. The IUE spectra do not reveal the presence of a hot subluminous secondary. $\nu$\,Pup is another Be star exhibiting an NRP variability pattern with long-term constancy and underlining the importance of combination frequencies and frequency groups. The star is a good target for efforts to identify an effectively single Be star.

We analyzed V-band photometry of the Be star {\omega} CMa, obtained during the last four decades, during which the star went through four complete cycles of disc formation and dissipation. The data were simulated by hydrodynamic models based on a time-dependent implementation of the viscous decretion disc (VDD) paradigm, in which a disc around a fast-spinning Be star is formed by material ejected by the star and driven to progressively larger orbits by means of viscous torques. Our simulations offer a good description of the photometric variability during phases of disc formation and dissipation, which suggests that the VDD model adequately describes the structural evolution of the disc. Furthermore, our analysis allowed us to determine the viscosity parameter {\alpha}, as well as the net mass and angular momentum (AM) loss rates. We find that {\alpha} is variable, ranging from 0.1 to 1.0, not only from cycle to cycle but also within a given cycle. Additionally, build-up phases usually have larger values of {\alpha} than the dissipation phases. Furthermore, during dissipation the outward AM flux is not necessarily zero, meaning that {\omega} CMa does not experience a true quiescence but, instead, switches between a high to a low AM loss rate during which the disc quickly assumes an overall lower density but never zero. We confront the average AM loss rate with predictions from stellar evolution models for fast-rotating stars, and find that our measurements are smaller by more than one order of magnitude.

High-dispersion spectroscopic monitoring of HESS J0632+057 has been carried out over four orbital cycles in order to search for orbital modulation, covering the entire orbital phase. We have measured radial velocity of H$\alpha$ emission line with the method introduced by Shafter et al. (1986), which has been successfully applied to some Be stars. The velocity is seen to increase much earlier than expected for the orbital period of 315 days, and much more steeply than expected at around "apastron". The period of the H$\alpha$ modulation is found to be as $308^{+26}_{-23}$ days. We have also analyzed Swift/XRT data from 2009 to 2015 to study the orbital modulation, selecting the data with good statistics ($\geq$ 30 counts). With additional two-year data to the previous works, the orbital period has been updated to $313^{+11}_{-8}$ days, which is consistent with the previous X-ray periods and the spectroscopic one. The past XMM-Newton and Chandra observations prefer the period of 313 days. With the new period, assuming that H$\alpha$ velocities accurately trace the motion of the Be star, we have derived a new set of the orbital parameters. In the new orbit, which is less eccentric ($e \sim 0.6$), two outbursts occur after apastron, and just after periastron. Besides, the column density in bright phase ($4.7^{+0.9}_{-08}\times10^{21}\;\mathrm{cm^{-2}}$) is higher than in faint phase ($2.2\pm0.5\times10^{21}\;\mathrm{cm^{-2}}$). These facts suggest that outbursts occur when the compact object passes nearby/through the Be disk. The mass function implies that mass of the compact object is less than 2.5 $\mathrm{M_{sun}}$ assuming that the mass of the Be star is 13.2--18.2 $\mathrm{M_{sun}}$ (Aragona at al. 2010) unless the inclination is extremely small. The photon index indicates that the spectra becomes softer when the system is bright. These suggest that the compact object is a pulsar.

Be stars are main-sequence massive stars with emission features in their spectrum, which originates in circumstellar gaseous discs. Even though the viscous decretion disc (VDD) model can satisfactorily explain most observations, two important physical ingredients, namely the magnitude of the viscosity ($\alpha$) and the disk mass injection rate, remain poorly constrained. The light curves of Be stars that undergo events of disc formation and dissipation offer an opportunity to constrain these quantities. A pipeline was developed to model these events that uses a grid of synthetic light curves, computed from coupled hydrodynamic and radiative transfer calculations. A sample of 54 Be stars from the OGLE survey of the Small Magellanic Cloud (SMC) was selected for this study. Because of the way our sample was selected (bright stars with clear disc events), it likely represents the densest discs in the SMC. Like their siblings in the Galaxy, the mass of the disc in the SMC increases with the stellar mass. The typical mass and angular momentum loss rates associated with the disk events are of the order of $\sim$$10^{-10}\, M_\odot\,\mathrm{yr^{-1}}$ and $\sim$$5\times 10^{36}\, \mathrm{g\, cm^{2}\, s^{-2}}$, respectively. The values of $\alpha$ found in this work are typically of a few tenths, consistent with recent results in the literature and with the ones found in dwarf novae, but larger than current theory predicts. Considering the sample as a whole, the viscosity parameter is roughly two times larger at build-up ($\left\langle\alpha_\mathrm{bu}\right\rangle = 0.63$) than at dissipation ($\left\langle\alpha_\mathrm{d}\right\rangle = 0.26$). Further work is necessary to verify whether this trend is real or a result of some of the model assumptions.

In order to study the growth and evolution of circumstellar disks around classical Be stars, we analyze optical time-series photometry from the KELT survey with simultaneous infrared and visible spectroscopy from the APOGEE survey and BeSS database for a sample of 160 Galactic classical Be stars. The systems studied here show variability including transitions from a diskless to a disk-possessing state (and vice versa), and persistent disks that vary in strength, being replenished at either regularly or irregularly occurring intervals. We detect disk-building events (outbursts) in the light curves of 28\% of our sample. Outbursts are more commonly observed in early- (57\%), compared to mid- (27\%) and late-type (8\%) systems. A given system may show anywhere between 0 -- 40 individual outbursts in its light curve, with amplitudes ranging up to $\sim$0.5 mag and event durations between $\sim$2 -- 1000 days. We study how both the photometry and spectroscopy change together during active episodes of disk growth or dissipation, revealing details about the evolution of the circumstellar environment. We demonstrate that photometric activity is linked to changes in the inner disk, and show that, at least in some cases, the disk growth process is asymmetrical. Observational evidence of Be star disks both growing and clearing from the inside out is presented. The duration of disk buildup and dissipation phases are measured for 70 outbursts, and we find that the average outburst takes about twice as long to dissipate as it does to build up in optical photometry. Our analysis hints that dissipation of the inner disk occurs relatively slowly for late-type Be stars.

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.

Observations of 25 Ori much expand the picture derived of other early-type Be stars with BRITE and SMEI. Two instead of one difference frequencies rule the variability: (a) The lower one, 0.0129 c/d, is the frequency of events with full amplitudes of 100-200 mmag which may signal mass loss possibly driven by the higher one, 0.1777 c/d. (b) Much of the entire power spectrum is a tightly woven network of combination frequencies: (i) Below 0.25\,c/d, numerous frequencies are difference frequencies. (ii) Many frequencies above 2.5 c/d can be represented as sum frequencies and in a few cases as harmonics. (iii) Many frequencies between 1.1 and 1.75\,c/d can be portrayed as parents of combination frequencies. The number and fraction of combination frequencies increases steeply with decreasing amplitude and and accuracy of the frequency matching.

The BRITE Constellation of nanosatellites obtained mmag photometry of 28 Cygni for 11 months in 2014-2016. Observations with the Solar Mass Ejection Imager in 2003-2010 and 118 H$\alpha$ line profiles were added. For decades, 28 Cyg has exhibited four large-amplitude frequencies: two closely spaced frequencies of spectroscopically confirmed $g$ modes near 1.5 c/d, one slightly lower exophotospheric (Stefl) frequency, and at 0.05 c/d the difference frequency between the two g modes. This top-level framework is indistinguishable from eta Cen (Paper I), which is also very similar in spectral type, rotation rate, and viewing angle. The Stefl frequency is the only one that does not seem to be affected by the difference frequency. The amplitude of the latter undergoes large variations; around maximum the amount of near-circumstellar matter is increased, and the amplitude of the Stefl frequency grows by some factor. During such brightenings dozens of transient spikes appear in the frequency spectrum, concentrated in three groups. Only eleven frequencies were common to all years of BRITE observations. Be stars seem to be controlled by several coupled clocks, most of which are not very regular on timescales of weeks to months but function for decades. The combination of g modes to the low difference frequency and/or the atmospheric response to it appears significantly nonlinear. Like in eta Cen, the difference-frequency variability seems the main responsible for the modulation of the star-to-disc mass transfer in 28 Cyg. A hierarchical set of difference frequencies may reach the longest timescales known of the Be phenomenon.

We investigated the properties of $80$ Be star disks, using the time dependent viscous decretion disk (VDD) model to interpret their mid-infrared continuum fluxes. Besides estimating the typical parameter ranges for such disks, this approach allowed us to relate the disk structure to its evolutionary status. In this contribution, we present a brief summary of the recently published results, and discuss some properties and time-scales of the VDD model solution.

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.

We use a smoothed particle hydrodynamics (SPH) code to examine the effects of misaligned binary companions on Be star discs. We systematically vary the degree of misalignment between the disc and the binary orbit, as well as the disc viscosity and orbital period to study their effects on the density in the inner and outer parts of the disc. We find that varying the degree of misalignment, the viscosity, and the orbital period affects both the truncation radius and the density structure of the outer disc, while the inner disc remains mostly unaffected. We also investigate the tilting of the disc in the innermost part of the disc and find the tilt increases with radius until reaching a maximum around 5 stellar radii. The direction of the line of nodes, with respect to the equator of the central star, is found to be offset compared to the orbital line of nodes, and to vary periodically in time, with a period of half a orbital phase. We also compare the scale height of our discs with the analytical scale height of an isothermal disc, which increases with radius as $r^{1.5}$. We find that this formula reproduces the scale height well for both aligned and misaligned systems but underestimates the scale height in regions of the disc where density enhancements develop.

High spatial resolution is the key for the understanding various astrophysical phenomena. But even with the future E-ELT, single dish instruments are limited to a spatial resolution of about 4 mas in the visible. For the closest objects within our Galaxy most of the stellar photosphere remains smaller than 1 mas. With the success of long baseline interferometry these limitations were soom overcome. Today low and high resolution interferometric instruments on the VLTI and CHARA offer an immense range of astrophysical studies. Combining more telescopes and moving to visible wavelengths broadens the science cases even more. With the idea of developing strong science cases for a future visible interferometer, we organized a science group around the following topics: pre-main sequence and main sequence stars, fundamental parameters, asteroseismology and classical pulsating stars, evolved stars, massive stars, active galactic nuclei (AGNs) and imaging techniques. A meeting was organized on the 15th and 16th of January, 2015 in Nice with the support of the Action Specific in Haute Resolution Angulaire (ASHRA), the Programme National en Physique Stellaire (PNPS), the Lagrange Laboratory and the Observatoire de la Cote d'Azur, in order to present these cases and to discuss them further for future visible interferometers. This White Paper presents the outcome of the exchanges. This book is dedicated to the memory of our colleague Olivier Chesneau who passed away at the age of 41.

We revisit the viscous decretion disk (VDD) model of the classical Be star $\beta$ CMi as presented by \citet{klement} using an updated version of the radiative transfer code {\ttfamily HDUST}. A software bug was causing the mean intensities to be slightly underestimated in the equatorial region of the disk, with small but detectable effects on the disk temperature and emergent spectrum. The new model fixes an unexplained feature of the original model, which was able to reproduce the observations only when considering a dual density behavior: a steep density fall-off in the very inner parts of the disk followed by a shallower density profile. The new model is able to reproduce all the observables reasonably well using a single power-law for the density profile throughout the whole disk, as predicted by the VDD model. All the other original conclusions, most importantly the reported truncation of the disk at the distance of 35 stellar equatorial radii from the central star, remain unchanged.

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)

Context. The mechanism of disk formation around fast-rotating Be stars is not well understood. In particular, it is not clear which mechanisms operate, in addition to fast rotation, to produce the observed variable ejection of matter. The star Achernar is a privileged laboratory to probe these additional mechanisms because it is close, presents B-Be phase variations on timescales ranging from 6 yr to 15 yr, a companion star was discovered around it, and probably presents a polar wind or jet. Aims. Despite all these previous studies, the disk around Achernar was never directly imaged. Therefore we seek to produce an image of the photosphere and close environment of the star. Methods. We used infrared long-baseline interferometry with the PIONIER/VLTI instrument to produce reconstructed images of the photosphere and close environment of the star over four years of observations. To study the disk formation, we compared the observations and reconstructed images to previously computed models of both the stellar photosphere alone (normal B phase) and the star presenting a circumstellar disk (Be phase). Results. The observations taken in 2011 and 2012, during the quiescent phase of Achernar, do not exhibit a disk at the detection limit of the instrument. In 2014, on the other hand, a disk was already formed and our reconstructed image reveals an extended H-band continuum excess flux. Our results from interferometric imaging are also supported by several H-alpha line profiles showing that Achernar started an emission-line phase sometime in the beginning of 2013. The analysis of our reconstructed images shows that the 2014 near-IR flux extends to 1.7 - 2.3 equatorial radii. Our model-independent size estimation of the H-band continuum contribution is compatible with the presence of a circumstellar disk, which is in good agreement with predictions from Be-disk models.

We analyze the V-band photometry data of the Be star omega CMa, observed during the last four decades. The data is fitted by hydrodynamic models based on the viscous decretion disk (VDD) theory, in which a disk around a fast-spinning Be star is formed by material ejected by the central star and driven to progressively wider orbits by means of viscous torques. For the first time, we apply the model for both the disk build up and the dissipation phases. Our simulations offer a good description of the photometric variability in both phases, which suggests that the VDD model adequately describes the disk structural evolution. Furthermore, our analysis allowed us to determine the viscosity parameter (alpha) of the gas, as well as the net mass loss rate. We find that alpha is variable, ranging from 0.1 to 1.0, and that build up phases have larger values of alpha than the dissipation phases. Additionally, we find that, contrary to what is generally assumed, even during quiescence the outward mass flux is never zero, suggesting that the star alternates between a high mass-loss phase (outburst) and a low mass-loss phase (quiescence).

The late-type Be star $\beta$ CMi is remarkably stable compared to other Be stars that have been studied. This has led to a realistic model of the outflowing Be disk by Klement et al. These results showed that the disk is likely truncated at a finite radius from the star, which Klement et al.~suggest is evidence for an unseen binary companion in orbit. Here we report on an analysis of the Ritter Observatory spectroscopic archive of $\beta$ CMi to search for evidence of the elusive companion. We detect periodic Doppler shifts in the wings of the H$\alpha$ line with a period of 170 d and an amplitude of 2.25 km s$^{-1}$, consistent with a low-mass binary companion ($M\approx 0.42 M_\odot$). We then compared the small changes in the violet-to-red peak height changes ($V/R$) with the orbital motion. We find weak evidence that it does follow the orbital motion, as suggested by recent Be binary models by Panoglou et al. Our results, which are similar to those for several other Be stars, suggest that $\beta$ CMi may be a product of binary evolution where Roche lobe overflow has spun up the current Be star, likely leaving a hot subdwarf or white dwarf in orbit around the star. Unfortunately, no direct sign of this companion star is found in the very limited archive of {\it International Ultraviolet Explorer} spectra.

New high-cadence observations with BRITE covering many months confirm that coupled pairs of nonradial pulsation modes are widespread among early-type Be stars. With the difference frequency between the parental variations they may form a roughly sinusoidal variability or the amplitude may cyclicly vary. A first - amplified - beat pattern is also found. In all three cases the amplitudes of difference frequencies can exceed the amplitude sum of the base frequencies, and modulations of the star-to-circumstellar-disk mass-transfer rate may be associated with these slow variations. This suggests more strongly than any earlier observations that significant dissipation of pulsational energy in the atmosphere may be a cause of mass ejections from Be stars. A unifying interpretative concept is presented.

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.

Classical Be stars are introduced as object class and their particular potential for space based photometry is highlighted. A brief summary of the various types of variability observed in Be stars makes clear that an interpretation of every single frequency as a pulsation mode falls short, instead there are as well purely circumstellar variations and those that originate in the immediate stellar to circumstellar interaction region. In particular the latter offer great potential, as they are linked to one of the few remaining great riddles of Be stars, namely how they feed their disks.

Space photometry has revealed rich pulsation spectra in classical Be stars. Often frequency pairs can be found that differ by nearly the frequency of a third genuine stellar variability. The lowest currently known of these so-called difference (or $\Delta$) frequencies reach 0.03 c/d, which is the upper limit of outburst frequencies. The potential of such slow variations for pulsation-assisted mass loss is discussed.

We analyze 45 spectropolarimetric observations of the eclipsing, interacting binary star V356 Sgr, obtained over a period of 21 years, to characterize the geometry of the system's circumstellar material. After removing interstellar polarization from these data, we find the system exhibits a large intrinsic polarization signature arising from electron scattering. In addition, the lack of repeatable eclipses in the polarization phase curves indicates the presence of a substantial pool of scatterers not occulted by either star. We suggest that these scatterers form either a circumbinary disk coplanar with the gainer's accretion disk or an elongated structure perpendicular to the orbital plane of V356 Sgr, possibly formed by bipolar outflows. We also observe small-scale, cycle-to-cycle variations in the magnitude of intrinsic polarization at individual phases, which we interpret as evidence of variability in the amount of scattering material present within and around the system. This may indicate a mass transfer or mass loss rate that varies on the time-scale of the system's orbital period. Finally, we compare the basic polarimetric properties of V356 Sgr with those of the well studied beta Lyr system; the significant differences observed between the two systems suggests diversity in the basic circumstellar geometry of Roche-lobe overflow systems

Context: [abbreviated] 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: [abbreviated] 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 \v{S}tefl 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 4 c/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.

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.

Be stars are surrounded by outflowing circumstellar matter structured in the form of decretion discs. They are often members of binary systems, where it is expected that the decretion disc interacts both radiatively and gravitationally with the companion. In this work we study how various orbital (period, mass ratio and eccentricity) and disc (viscosity) parameters affect the disc structure in coplanar systems. We simulate such binaries with the use of a smoothed particle hydrodynamics code. The main effects of the secondary on the disc are its truncation and the accumulation of material inwards of truncation. We find two cases with respect to the effects of eccentricity: (i) In circular or nearly circular prograde orbits, the disc maintains a rotating, constant in shape, configuration, which is locked to the orbital phase. The disc is smaller in size, more elongated and more massive for low viscosity parameter, small orbital separation and/or high mass ratio. (ii) Highly eccentric orbits are more complex, with the disc structure and total mass strongly dependent on the orbital phase and the distance to the secondary. We also study the effects of binarity in the disc continuum emission. Since the infrared and radio SED are sensitive to the disc size and density slope, the truncation and matter accumulation result in considerable modifications in the emergent spectrum. We conclude that binarity can serve as an explanation for the variability exhibited in observations of Be stars, and that our model can be used to detect invisible companions.

In the last decade, the viscous decretion disk model has emerged as the new paradigm for Be star disks. In this contribution, we propose a simple analytical model to estimate the continuum infrared excess arising from these circumstellar disks, in the light of the currently accepted scenario. We demonstrate that the disk can be satisfactorily described by a two component system: an inner optically thick region, which we call the pseudo-photosphere, and a diffuse outer part. In particular, a direct connexion between the disk brightness profile and the thermal structure is derived, and then confronted to realistic numerical simulations. This result quantifies how the non-isothermality of the disk ultimately affects both infrared measured fluxes and visibilities.

Spectroscopic observations of Achernar over the past decades, have shown the photospheric line width, as measured by the rotational parameter $v \sin i$, to vary in correlation with the emission activity. Here we present new observations, covering the most recent activity phase, and further archival data collected from the archives. The $v \sin i$ variation is confirmed. On the basis of the available data it cannot be decided with certainty whether the increased line width precedes the emission activity, i.e. is a signature of the ejection mechanism, or postdates is, which would make it a signature of re-accretion of some of the disk-material. However, the observed evidence leans towards the re-accretion hypothesis. Two further stars showing the effect of variable line width in correlation with emission activity, namely 66 Oph and $\pi$ Aqr, are presented as well.

Empirical evidence for the involvement of nonradial pulsations (NRP's) in the mass loss from Be stars ranges from (i) a singular case (\object{$\mu$ Cen}) of repetitive mass ejections triggered by multi-mode beating to (ii) several photometric reports about enormous numbers of pulsation modes popping up during outbursts and on to (iii) effective single-mode pulsators. The BRITE Constellation of nanosatellites was used to obtain mmag photometry of the Be stars $\eta$ and \object{$\mu$ Cen}. In the low-inclination star \object{$\mu$ Cen}, light pollution by variable amounts of near-stellar matter prevented any new insights into the variability and other properties of the central star. In the equator-on star \object{$\eta$ Cen}, BRITE photometry and {\sc Heros} echelle spectroscopy from the 1990s reveal an intricate clockwork of star-disk interactions. The mass transfer is modulated with the frequency difference of two NRP modes and an amplitude three times as large as the amplitude sum of the two NRP modes. This process feeds a high-amplitude circumstellar activity running with the incoherent and slightly lower so-called \v{S}tefl frequency. The mass loss-modulation cycles are tightly coupled to variations in the value of the \v{S}tefl frequency and in its amplitude, albeit with strongly drifting phase differences. The observations are well described by the decomposition of the mass loss into a pulsation-related engine in the star and a viscosity-dominated engine in the circumstellar disk. Arguments are developed that large-scale gas-circulation flows occur at the interface. The propagation rates of these eddies manifest themselves as \v{S}tefl frequencies. Bursts in power spectra during mass-loss events can be understood as the noise inherent to these gas flows.

The MiMeS project is a large-scale, high resolution, sensitive spectropolarimetric investigation of the magnetic properties of O and early B type stars. Initiated in 2008 and completed in 2013, the project was supported by 3 Large Program allocations, as well as various programs initiated by independent PIs and archival resources. Ultimately, over 4800 circularly polarized spectra of 560 O and B stars were collected with the instruments ESPaDOnS at the Canada-France-Hawaii Telescope, Narval at the T\'elescope Bernard Lyot, and HARPSpol at the European Southern Observatory La Silla 3.6m telescope, making MiMeS by far the largest systematic investigation of massive star magnetism ever undertaken. In this paper, the first in a series reporting the general results of the survey, we introduce the scientific motivation and goals, describe the sample of targets, review the instrumentation and observational techniques used, explain the exposure time calculation designed to provide sensitivity to surface dipole fields above approximately 100 G, discuss the polarimetric performance, stability and uncertainty of the instrumentation, and summarize the previous and forthcoming publications.

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^{\circ}$. 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\%$, 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^{-12}$ to $10^{-9}\,{\rm M_{\odot}\,yr^{-1}}$, which is at least two orders of magnitude smaller than the previous outflowing disc model predictions.

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.

We present spatially and spectrally resolved Br-gamma emission around the planet-hosting, transitional Herbig Ae/Be star HD 100546. Aiming to gain insight into the physical origin of the line in possible relation to accretion processes, we carried out Br-gamma spectro-interferometry using AMBER/VLTI from three different baselines achieving spatial and spectral resolutions of 2-4 mas and 12000. The Br-gamma visibility is larger than that of the continuum for all baselines. Differential phases reveal a shift between the photocentre of the Br-gamma line -displaced 0.6 mas (0.06 au at 100 pc) NE from the star- and that of the K-band continuum emission -displaced 0.3 mas NE from the star. The photocentres of the redshifted and blueshifted components of the Br-gamma line are located NW and SE from the photocentre of the peak line emission, respectively. Moreover, the photocentre of the fastest velocity bins within the spectral line tends to be closer to that of the peak emission than the photocentre of the slowest velocity bins. Our results are consistent with a Br-gamma emitting region inside the dust inner rim (<0.25 au) and extending very close to the central star, with a Keplerian, disc-like structure rotating counter-clockwise, and most probably flared (25 deg). Even though the main contribution to the Br-gamma line does not come from gas magnetically channelled on to the star, accretion on to HD 100546 could be magnetospheric, implying a mass accretion rate of a few 10^-7 Msun/yr. This value indicates that the observed gas has to be replenished on time-scales of a few months to years, perhaps by planet-induced flows from the outer to the inner disc as has been reported for similar systems.

We have used the radiative transfer code HDUST to analyze and interpret the long-term photometric behavior of the Be star omega CMa, considering four complete cycles of disk formation and dissipation. This is the first time in which a full lightcurve of a Be star was investigated and modeled including both disk build-up and dissipation phases. Based on the quite good fit of the observed data we were able to derive the history of stellar mass decretion rates (including long- and short-term changes) during the disk formation and dissipation phases in all four cycles.

We report on monitoring observations of the TeV gamma-ray binary HESS J0632+057, which were carried out to constrain the interaction between the Be circumstellar disk and the compact object of unknown nature, and provide for the first time high-dispersion (R > 50000) optical spectra in the second half of the orbital cycle, from apastron through periastron. The Halpha, Hbeta, and Hgamma line profiles are found to exhibit remarkable short-term variability for ~1 month after the apastron (phase 0.6--0.7), whereas they show little variation near the periastron. These emission lines show "S-shaped" variations with timescale of ~150 days, which is about twice that reported previously. In contrast to the Balmer lines, no profile variability is seen in any FeII emission line. We estimate the radii of emitting regions of the Halpha, Hbeta, Hgamma, and FeII emission lines to be ~30, 11, 7, and 2 stellar radii (R_*), respectively. The amplitudes of the line profile variations in different lines indicate that the interaction with the compact object affects the Be disk down to, at least, the radius of 7 R_* after the apastron. This fact, together with little profile variability near the periastron, rules out the tidal force as the major cause of disk variability. Although this leaves the pulsar wind as the most likely candidate mechanism for disk variations, understanding the details of the interaction, particularly the mechanism for causing a large-scale disk disturbance after the apastron, remains an open question.

Context. Various Be stars exhibit intensity variations of the violet and red emission peaks in their HI lines observed in emission. This so-called $V/R$ phenomenon is usually explained by the precession of a one-armed spiral density perturbation in the circumstellar disk. That global-disk oscillation scenario was confirmed, both observationally and theoretically, in the previous series of two papers analyzing the Be shell star {\zeta} Tauri. The vertically averaged (2D) global-disk oscillation model used at the time was able to reproduce the $V/R$ variations observed in H{\alpha}, as well as the spatially resolved interferometric data from AMBER/VLTI. Unfortunately, that model failed to reproduce the $V/R$ phase of Br15 and the amplitude of the polarization variation, suggesting that the inner disk structure predicted by the model was incorrect. Aims. The first aim of the present paper is to quantify the temporal variations of the shell-line characteristics of {\zeta} Tauri. The second aim is to better understand the physics underlying the $V/R$ phenomenon by modeling the shell-line variations together with the $V/R$ and polarimetric variations. The third aim is to test a new 2.5D disk oscillation model, which solves the set of equations that describe the 3D perturbed disk structure but keeps only the equatorial (i.e., 2D) component of the solution. This approximation was adopted to allow comparisons with the previous 2D model, and as a first step toward a future 3D model. Results. The new 2.5D formalism improves the agreement with the observed $V/R$ variations of H{\alpha} and Br15, under the proviso that a large value of the viscosity parameter, {\alpha} = 0.8, be adopted. Nonetheless, it remains challenging for the models to reproduce consistently the amplitude and the average level of the polarization data, whatever formalism is adopted.

Be stars exhibit variability for a great number of observables. Putting the pieces together of the disk dynamics is not an easy task and requires arduous modeling before achieving a good fit of the observational data. In order to guide the modeling process and make it more efficient, it is very instructive to investigate reference dynamical cases. This paper focuses on continuum polarimetric quantities and is the second of a series that aims to demonstrate the capacity of deriving the dynamical history and fundamental parameters of a classical Be star through the follow-up of various observables. After a detailed study of the different opacities at play in the formation of polarized spectra, we investigate predictions of polarimetric observables in the continuum for different dynamical scenarios. Our models are based on a coupling of a hydrodynamic viscous decretion simulations in a disk and a 3-D non-LTE radiative transfer code. Introducing the polarization color diagram (PCD), we show that certain combinations of polarimetric observables exhibit features that are characteristic of a mass loss history. This diagram also enables to estimate fundamental parameters such as the inclination angle, the disk density scale and the $\alpha$ viscous diffusion parameter. We present the PCD as a powerful diagnosis tool to track the dynamical phases of a Be star such as disk building-up, dissipation, periodic and episodic outbursts. Finally we confront our models with observations of 4 Be stars that exhibited long-term polarimetric activity.

Recent observational and theoretical studies of classical Be stars have established the utility of polarization color diagrams (PCD) in helping to constrain the time-dependent mass decretion rates of these systems. We expand on our pilot observational study of this phenomenon, and report the detailed analysis of a long-term (1989-2004) spectropolarimetric survey of 9 additional classical Be stars, including systems exhibiting evidence of partial disk-loss/disk-growth episodes as well as systems exhibiting long-term stable disks. After carefully characterizing and removing the interstellar polarization along the line of sight to each of these targets, we analyze their intrinsic polarization behavior. We find that many steady-state Be disks pause at the top of the PCD, as predicted by theory. We also observe sharp declines in the Balmer jump polarization for later spectral type, near edge-on steady-state disks, again as recently predicted by theory, likely caused when the base density of the disk is very high, and the outer region of the edge-on disk starts to self absorb a significant number of Balmer jump photons. The intrinsic $V$-band polarization and polarization position angle of $\gamma$ Cas exhibits variations that seem to phase with the orbital period of a known one-armed density structure in this disk, similar to the theoretical predictions of Halonen & Jones. We also observe stochastic jumps in the intrinsic polarization across the Balmer jump of several known Be+sdO systems, and speculate that the thermal inflation of part of the outer region of these disks could be responsible for producing this observational phenomenon. Finally, we estimate the base densities of this sample of stars to be between $\approx 8\times 10^{-11}$ to $\approx 4 \times 10^{-12}\,\rm g cm^{-3}$ during quasi steady state periods given there maximum observed polarization.

In the past decade, a consensus has emerged regarding the nature of classical Be stars: They are very rapidly rotating main sequence B stars, which, through a still unknown, but increasingly constrained process, form an outwardly diffusing gaseous, dust-free Keplerian disk. In this work, first the definition of Be stars is contrasted to similar classes, and common observables obtained for Be stars are introduced and the respective formation mechanisms explained. We then review the current state of knowledge concerning the central stars as non-radially pulsating objects and non-magnetic stars, as far as it concerns large scale, i.e., mostly dipolar, global fields. Localized, weak magnetic fields remain possible, but are as of yet unproven. The Be phenomenon, linked with one or more mass ejection processes, acts on top of a rotation rate of about 75% of critical or above. The properties of the process can be well constrained, leaving only few options, most importantly, but not exclusively, non-radial pulsation and small scale magnetic fields. Of these, it is well possible that all are realized: In different stars, different processes may be acting. Once the material has been lifted into Keplerian orbit, memory of the details of the ejection process is lost, and the material is governed by viscosity. The disks are fairly well understood in the theoretical framework of the viscous decretion disk model. This is not only true for the disk structure, but as well for its variability, both cyclic and secular. Be binaries are reviewed under the aspect of the various types of interactions a companion can have with the circumstellar disk. Finally, extragalactic Be stars, at lower metallicities, seem more common and more rapidly rotating.

The main theoretical problem for the formation of a Keplerian disk around Be stars is how to supply angular momentum from the star to the disk, even more so since Be stars probably rotate somewhat sub-critically. For instance, nonradial pulsation may transport angular momentum to the stellar surface until (part of) this excess supports the disk formation/replenishment. The nearby Be star Achernar is presently building a new disk and offers an excellent opportunity to observe this process from relatively close-up. Spectra from various sources and epochs are scrutinized to identify the salient stellar parameters characterizing the disk life cycle as defined by H\alpha emission. Variable strength of the non-radial pulsation is confirmed, but does not affect the further results. For the first time it is demonstrated that the photospheric line width does vary in a Be star, by as much as \Delta v sin i \lesssim 35kms^{-1}. However, contrary to assumptions in which a photospheric spin-up accumulates during the diskless phase and then is released into the disk as it is fed, the apparent photospheric spin-up is positively correlated with the appearance of H\alpha line emission: The photospheric line widths and circumstellar emission increase together, and the apparent stellar rotation declines to the value at quiescence after the H\alpha line emission becomes undetectable.

Context. The circumstellar disks ejected by many rapidly rotating B stars (so-called Be stars) offer the rare opportunity of studying the structure and dynamics of gaseous disks at high spectral as well as angular resolution. Aims. This paper explores a newly identified effect in spectro-interferometric phase that can be used for probing the inner regions of gaseous edge-on disks on a scale of a few stellar radii. Methods. The origin of this effect (dubbed central quasi-emission phase signature, CQE-PS) lies in the velocity-dependent line absorption of photospheric radiation by the circumstellar disk. At high spectral and marginal interferometric resolution, photocenter displacements between star and isovelocity regions in the Keplerian disk reveal themselves through small interferometric phase shifts. To investigate the diagnostic potential of this effect, a series of models are presented, based on detailed radiative transfer calculations in a viscous decretion disk. Results. Amplitude and detailed shape of the CQE-PS depend sensitively on disk density and size and on the radial distribution of the material with characteristic shapes in differential phase diagrams. In addition, useful lower limits to the angular size of the central stars can be derived even when the system is almost unresolved. Conclusions. The full power of this diagnostic tool can be expected if it can be applied to observations over a full life-cycle of a disk from first ejection through final dispersal, over a full cycle of disk oscillations, or over a full orbital period in a binary system.

Some massive stars possess strong magnetic fields that confine plasma in the circumstellar environment. These \textit{magnetospheres} have been studied spectroscopically, photometrically and, more recently, interferometrically. Here we report on the first firm detection of a magnetosphere in continuum linear polarization, as a result of monitoring of $\sigma$\,Ori\,E at the Pico dos Dias Observatory. {The non-zero intrinsic polarization indicates an asymmetric structure, whose minor elongation axis is oriented $150\fdg0$ east of the celestial north.} A modulation of the polarization was observed, with a period of half of the rotation period, which supports the theoretical prediction of the presence of two diametrally opposed, co-rotating blobs of gas. A phase lag of -0.085 was detected between the polarization minimum and the primary minimum of the light curve, suggestive of a complex shape of the plasma clouds. We present a preliminary analysis of the data with the Rigidly Rotating Magnetosphere model, which could not reproduce simultaneously the photometric and polarimetric data. A toy model comprising two spherical co-rotating blobs {joined by a thin disk} proved more successful in reproducing the polarization modulation. {With this model we were able to determine that the total scattering mass of the thin disk is similar to the mass of the blobs ($2M_{\rm b}/M_{\rm d}=1.2$) and that the blobs are rotating counterclockwise on the plane of the sky.} This result shows that polarimetry can provide a diagnostic of the geometry of clouds, which will serve as an important constraint for {improving} the Rigidly Rotating Magnetosphere model.

We describe the results of the world-wide observing campaign of the highly eccentric Be binary system delta Scorpii 2011 periastron passage which involved professional and amateur astronomers. Our spectroscopic observations provided a precise measurement of the system orbital period at 10.8092+/- 0.0005 years. Fitting of the He II 4686A line radial velocity curve determined the periastron passage time on 2011 July 3, UT 9:20 with a 0.9--day uncertainty. Both these results are in a very good agreement with recent findings from interferometry. We also derived new evolutionary masses of the binary components (13 and 8.2 Msun) and a new distance of 136 pc from the Sun, consistent with the HIPPARCOS parallax. The radial velocity and profile variations observed in the H_alpha line near the 2011 periastron reflected the interaction of the secondary component and the circumstellar disk around the primary component. Using these data, we estimated a disk radius of 150 Rsun. Our analysis of the radial velocity variations measured during the periastron passage time in 2000 and 2011 along with those measured during the 20th century, the high eccentricity of the system, and the presence of a bow shock-like structure around it suggest that delta Sco might be a runaway triple system. The third component should be external to the known binary and move on an elliptical orbit that is tilted by at least 40 degree with respect to the binary orbital plane for such a system to be stable and responsible for the observed long-term radial velocity variations.

The Be star \delta\ Scorpii is an interesting binary system, whose primary companion created a circumstellar disk after the periastron passage of the secondary in 2000, being since then classified as Be. This work presents the results of a long-term monitoring of this star in broad-band imaging polarimetry. The observational data collected since 2006 in the Pico dos Dias Observatory (Brazil) show a variable polarization that seems to correlate with the photometric light curve. From this data we see that the disk density varied since 2006; furthermore, the data suggests that there was some disturbance of the disk during the last periastron passage in July, 2011.

We present prelimiary results from a coordinated spectroscopic campaign in 2011, centered on the \delta Sco periastron passage in July. Data have mostly been obtained with the FEROS/2.2 m at La Silla and ESPaDOnS/CFHT at Mauna Kea echelle instruments. Main results include the absence of tidally induced disturbance to the main \beta Cephei pulsation mode and the absence of tidally triggered mass-ejection at time of periastron proper. The observed (as far as yet analyzed) variations are compatible with the picture of a disk that is disturbed on its outer radius, with the disturbance propagating inwards after the periastron.

We analyse a feature that is detectable at the differential phases of high-resolution spectro-interferometry of Be stars. The origin of this feature (dubbed CQE-PS, Central Quasi Emission Phase Signature) lies in the differential absorption of photospheric radiation by the circumstellar disk, which is spectroscopically detected as a CQE line profile in shell stars. This phenomenon has great diagnostic potential for Be star disks, revealing properties of these disks on the scale of a few stellar radii.

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 build-up 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 3-D non-LTE radiative transfer code. We show that lightcurves for different wavelengths are specific of a mass loss history, inclination angle and $\alpha$ viscosity parameter. The diagnostic potential of those lightcurves is also discussed.

While the presence of discs around classical Be stars is well established, their origin is still uncertain. To understand what processes result in the creation of these discs and how angular momentum is transported within them, their physical properties must be constrained. This requires comparing high spatial and spectral resolution data with detailed radiative transfer modelling. We present a high spectral resolution, R~80,000, sub milli-arcsecond precision, spectroastrometric study of the circumstellar disc around the Be star beta CMi. The data are confronted with three-dimensional, NLTE radiative transfer calculations to directly constrain the properties of the disc. Furthermore, we compare the data to disc models featuring two velocity laws; Keperian, the prediction of the viscous disc model, and angular momentum conserving rotation. It is shown that the observations of beta CMi can only be reproduced using Keplerian rotation. The agreement between the model and the observed SED, polarisation and spectroastrometric signature of beta CMi confirms that the discs around Be stars are well modelled as viscous decretion discs.

Spectroscopic, polarimetric, and high spectral resolution interferometric data covering the period 1995-2011 are analyzed to document the transition into a new phase of circumstellar disk activity in the classical Be-shell star 48 Lib. The long-term disk evolution is described using the V/R ratio of the violet and red emission components of H alpha and Br gamma, radial velocities and profiles of He I and optical metal shell lines, multi-band BVRI polarimetry, broad-band and high-resolution interferometric visibilities and phases. Spectroscopic signatures of disk asymmetries in 48 Lib vanished in the late nineties but recovered some time between 2004 and 2007. Variations in the radial velocity and line profile of conventional shell lines correlate with the V/R behavior. They are shared by narrow absorption cores superimposed on otherwise seemingly photospheric He I lines, which may form in high-density gas at the inner disk close to the photosphere. Large radial velocity variations continued also during the V/R-quiet years. The broad H-band interferometry gives a disk diameter of (1.72+-0.2) mas (equivalent to 15 stellar radii), position angle of the disk (50+-9) deg and a relatively low disk flattening of 1.66+-0.3. Within the errors the same disk position angle is derived from polarimetric observations and from photocenter shifts across Br gamma. The high-resolution interferometric visibility and phase profiles show a double or even multiple-component structure. A preliminary estimate based on the size of the Br gamma emitting region indicates a large diameter for the disk (tens of stellar radii).

SOUTH POL will be a survey of the Southern sky in optical polarized light. It will use a newly designed polarimetric module at an 80cm Robotic Telescope. Telescope and polarimeter will be installed at CTIO, Chile, in late 2012. The initial goal is to cover the sky south of declination -15{\deg} in two years of observing time, aiming at a polarimetric accuracy \lesssim 0.1% down to V=15, with a camera covering a field of about 2.0 square degrees. SOUTH POL will impact areas such as Cosmology, Extragalactic Astronomy, Interstellar Medium of the Galaxy and Magellanic Clouds, Star Formation, Stellar Envelopes, Stellar explosions and Solar System, among others.

Be stars possess gaseous circumstellar decretion disks, which are well described using standard $\alpha$-disk theory. The Be star 28 CMa recently underwent a long outburst followed by a long period of quiescence, during which the disk dissipated. Here we present the first time-dependent models of the dissipation of a viscous decretion disk. By modeling the rate of decline of the V-band excess, we determine that the viscosity parameter $\alpha=1.0\pm0.2$, corresponding to a mass injection rate $\dot{M}=(3.5\pm 1.3) \times 10^{-8}\ M_\sun\,\mathrm{yr}^{-1}$. Such a large value of $\alpha$ suggests that the origin of the turbulent viscosity is an instability in the disk whose growth is limited by shock dissipation. The mass injection rate is more than an order of magnitude larger than the wind mass loss rate inferred from UV observations, implying that the mass injection mechanism most likely is not the stellar wind, but some other mechanism.

We present the basic concept of a new astronomical instrument: SPARC4 - Simultaneous Polarimeter and Rapid Camera in 4 bands. SPARC4 combines in one instrument: (i) photometric and polarimetric modes; (ii) sub-second time-resolution in photometric mode and excellent time-resolution in polarimetric mode; (iii) simultaneous imaging in four broad-bands for both modes. This combination will make SPARC4 a unique facility for ground-based optical observatories. Presently, the project is in its conceptual design phase.

High spectral-resolution Br-gamma VLTI/AMBER observations of classical Be stars reveal complex visibility, phase, and closure phase profiles, which were not seen in previous lower-resolution data. They present new challenges for the modeling of Be stars and may require additional ingredients over and above the conventional central star + circumstellar disk paradigm.

We present mid infrared (Mid-IR) spectra of the Compton-thick Seyfert 2 galaxy NGC\,3281, obtained with the Thermal-Region Camera Spectrograph (T-ReCS) at the Gemini South telescope. The spectra present a very deep silicate absorption at 9.7\,$\mu$m, and [S{\sc\,iv]}\,10.5\,$\mu$m and [Ne{\sc\,ii]}\,12.7\,$\mu$m ionic lines, but no evidence of PAH emission. We find that the nuclear optical extinction is in the range 24 $\leq$ A$_{V}$ $\leq$ 83\,mag. A temperature T = 300\,K was found for the black-body dust continuum component of the unresolved 65\,pc nucleus and at 130\,pc SE, while the region at 130\,pc reveals a colder temperature (200\,K). We describe the nuclear spectrum of NGC\,3281 using a clumpy torus model that suggests that the nucleus of this galaxy hosts a dusty toroidal structure. According to this model, the ratio between the inner and outer radius of the torus in NGC\,3281 is $R_0/R_d$ = 20, with {\bf 14} clouds in the equatorial radius with optical depth of $\tau_{V}$ = 40\,mag. We would be looking in the direction of the torus equatorial radius ($i$ = {\bf 60$^{\circ}$}), which has outer radius of R$_{0}\,\sim$ 11\,pc. The column density is N$_{H}\approx$\,{\bf 1.2}\,$\times\,10^{24}\,cm^{-2}$ and iron K$\alpha$ equivalent width ($\approx$ 0.5 - 1.2\,keV) are used to check the torus geometry. Our findings indicate that the X-ray absorbing column density, which classifies NGC\,3281 as a Compton-thick source, may also be responsible for the absorption at 9.7\,$\mu$m providing strong evidence that the silicate dust responsible for this absorption can be located in the AGN torus.

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 system's stellar decretion rate.

In Wisniewski et al. 2010, paper I, we analyzed 15 years of spectroscopic and spectropolarimetric data from the Ritter and Pine Bluff Observatories of 2 Be stars, 60 Cygni and {\pi} Aquarii, when a transition from Be to B star occurred. Here we anaylize the intrinsic polarization, where we observe loop-like structures caused by the rise and fall of the polarization Balmer Jump and continuum V-band polarization being mismatched temporally with polarimetric outbursts. We also see polarization angle deviations from the mean, reported in paper I, which may be indicative of warps in the disk, blobs injected at an inclined orbit, or spiral density waves. We show our ongoing efforts to model time dependent behavior of the disk to constrain the phenomena, using 3D Monte Carlo radiative transfer codes.

We present a novel theoretical tool to analyze the dynamical behaviour of a Be disk fed by non-constant decretion rates. It is mainly based on the computer code HDUST, a fully three-dimensional radiative transfer code that has been successfully applied to study several Be systems so far, and the SINGLEBE code that solves the 1D viscous diffusion problem. We have computed models of the temporal evolution of different types of Be star disks for different dynamical scenarios. By showing the behaviour of a large number of observables (interferometry, polarization, photometry and spectral line profiles), we show how it is possible to infer from observations some key dynamical parameters of the disk.

B[e] supergiants are luminous, massive post-main sequence stars exhibiting non-spherical winds, forbidden lines, and hot dust in a disc-like structure. The physical properties of their rich and complex circumstellar environment (CSE) are not well understood, partly because these CSE cannot be easily resolved at the large distances found for B[e] supergiants (typically $\ga 1$~kpc). From mid-IR spectro-interferometric observations obtained with VLTI/MIDI we seek to resolve and study the CSE of the Galactic B[e] supergiant CPD-57\degr\,2874. For a physical interpretation of the observables (visibilities and spectrum) we use our ray-tracing radiative transfer code (FRACS), which is optimised for thermal spectro-interferometric observations. Thanks to the short computing time required by FRACS ($<10$~s per monochromatic model), best-fit parameters and uncertainties for several physical quantities of CPD-57\degr\,2874 were obtained, such as inner dust radius, relative flux contribution of the central source and of the dusty CSE, dust temperature profile, and disc inclination. The analysis of VLTI/MIDI data with FRACS allowed one of the first direct determinations of physical parameters of the dusty CSE of a B[e] supergiant based on interferometric data and using a full model-fitting approach. In a larger context, the study of B[e] supergiants is important for a deeper understanding of the complex structure and evolution of hot, massive stars.

Circumstellar discs of Be stars are thought to be formed from material ejected from a fast-spinning central star. This material possesses large amounts of angular momentum and settles in a quasi-Keplerian orbit around the star. This simple description outlines the basic issues that a successful disc theory must address: 1) What is the mechanism responsible for the mass ejection? 2) What is the final configuration of the material? 3) How the disc grows? With the very high angular resolution that can be achieved with modern interferometers operating in the optical and infrared we can now resolve the photosphere and immediate vicinity of nearby Be stars. Those observations are able to provide very stringent tests for our ideas about the physical processes operating in those objects. This paper discusses the basic hydrodynamics of viscous decretion discs around Be stars. The model predictions are quantitatively compared to observations, demonstrating that the viscous decretion scenario is currently the most viable theory to explain the discs around Be stars.

FS CMa type stars are a recently described group of objects with the B[e] phenomenon that exhibit strong emission-line spectra and strong IR excesses. In this paper we report the first attempt for a detailed modeling of IRAS 00470+6429, for which we have the best set of observations. Our modeling is based on two key assumptions: the star has a main-sequence luminosity for its spectral type (B2) and the circumstellar envelope is bimodal, composed of a slowly outflowing disk-like wind and a fast polar wind. Both outflows are assumed to be purely radial. We adopt a novel approach to describe the dust formation site in the wind that employs timescale arguments for grain condensation and a self-consistent solution for the dust destruction surface. With the above assumptions we were able to reproduce satisfactorily many observational properties of IRAS 00470+6429, including the H line profiles and the overall shape of the spectral energy distribution. Our adopted recipe for dust formation proved successful in reproducing the correct amount of dust formed in the circumstellar envelope. Possible shortcomings of our model, as well as suggestions for future improvements, are discussed.

Context. The star HD 87643, exhibiting the "B[e] phenomenon", has one of the most extreme infrared excesses for this object class. It harbours a large amount of both hot and cold dust, and is surrounded by an extended reflection nebula. Aims. One of our major goals was to investigate the presence of a companion in HD87643. In addition, the presence of close dusty material was tested through a combination of multi-wavelength high spatial 5Aresolution observations. Methods. We observed HD 87643 with high spatial resolution techniques, using the near-IR AMBER/VLTI interferometer with baselines ranging from 60 m to 130 m and the mid-IR MIDI/VLTI interferometer with baselines ranging from 25 m to 65 m. These observations are complemented by NACO/VLT adaptive-optics-corrected images in the K and L-bands, ESO-2.2m optical Wide-Field Imager large-scale images in the B, V and R-bands, Results. We report the direct detection of a companion to HD 87643 by means of image synthesis using the AMBER/VLTI instrument. The presence of the companion is confirmed by the MIDI and NACO data, although with a lower confidence. The companion is separated by ~ 34 mas with a roughly north-south orientation. The period must be large (several tens of years) and hence the orbital parameters are not determined yet. Binarity with high eccentricity might be the key to interpreting the extreme characteristics of this system, namely a dusty circumstellar envelope around the primary, a compact dust nebulosity around the binary system and a complex extended nebula witnessing past violent ejections.

Emission lines formed in decretion disks of Be stars often undergo long-term cyclic variations, especially in the violet-to-red (V/R) ratio of their primary components. From observations of the bright Be-shell star zeta Tau, the possibly broadest and longest data set illustrating the prototype of this behaviour was compiled from our own and archival observations. It comprises optical and infrared spectra, broad-band polarimetry, and interferometric observations. From 3 V/R cycles between 1997 and 2008, a mean cycle length in H alpha of 1400-1430 days was derived. After each minimum in V/R, the shell absorption weakens and splits into two components, leading to 3 emission peaks. This phase makes the strongest contribution to the variability in cycle length. V/R curves of different lines are shifted in phase. Lines formed on average closer to the central star are ahead of the others. The shell absorption lines fall into 2 categories differing in line width, ionization/excitation potential, and variability of the equivalent width. The interferometry has resolved the continuum and the line emission in Br gamma and He I 2.06. The phasing of the Br gamma emission shows that the photocenter of the line-emitting region lies within the plane of the disk but is offset from the continuum source. The plane of the disk is constant throughout the observed V/R cycles. The self-consistent, one-armed, disk-oscillation model is developed in Paper II.

FS CMa type stars are a group of Galactic objects with the B[e] phenomenon. They exhibit strong emission-line spectra and infrared excesses, which are most likely due to recently formed circumstellar dust. The group content and identification criteria were described in the first two papers of the series. In this paper we report our spectroscopic and photometric observations of the optical counterpart of IRAS 00470+6429 obtained in 2003--2008. The optical spectrum is dominated by emission lines, most of which have P Cyg type profiles. We detected significant brightness variations, which may include a regular component, and variable spectral line profiles in both shape and position. The presence of a weak Li {\sc I} 6708 \AA line in the spectrum suggests that the object is most likely a binary system with a B2--B3 spectral type primary companion of a luminosity $\log$ L/L$\odot$ = 3.9$\pm$0.3 and a late-type secondary companion. We estimate a distance toward the object to be 2.0$\pm$0.3 kpc from the Sun.

Aims. In this paper we model, in a self-consistent way, polarimetric, photometric, spectrophotometric and interferometric observations of the classical Be star $\zeta$ Tauri. Our primary goal is to conduct a critical quantitative test of the global oscillation scenario. Methods. We have carried out detailed three-dimensional, NLTE radiative transfer calculations using the radiative transfer code HDUST. For the input for the code we have used the most up-to-date research on Be stars to include a physically realistic description for the central star and the circumstellar disc. We adopt a rotationally deformed, gravity darkened central star, surrounded by a disc whose unperturbed state is given by a steady-state viscous decretion disc model. We further assume that disc is in vertical hydrostatic equilibrium. Results. By adopting a viscous decretion disc model for $\zeta$ Tauri and a rigorous solution of the radiative transfer, we have obtained a very good fit of the time-average properties of the disc. This provides strong theoretical evidence that the viscous decretion disc model is the mechanism responsible for disc formation. With the global oscillation model we have successfully fitted spatially resolved VLTI/AMBER observations and the temporal V/R variations of the H$\alpha$ and Br$\gamma$ lines. This result convincingly demonstrates that the oscillation pattern in the disc is a one-armed spiral. Possible model shortcomings, as well as suggestions for future improvements, are also discussed.

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 3-D NLTE Monte Carlo simulations, that viscous keplerian disks can be highly non-isothermal. In this paper we build upon our previous work and solve the full problem of the steady-state non-isothermal 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 non-isothermal disk models must be used for a detailed modeling of Be star disks.

The recent interferometric study of Achernar, leading to the conclusion that its geometrical oblateness cannot be explained in the Roche approximation, has stirred substantial interest in the community, in view of its potential impact in many fields of stellar astrophysics. It is the purpose of this paper to reinterpret the interferometric observations with a fast rotating, gravity darkened central star surrounded by a small equatorial disk, whose presence is consistent with contemporaneous spectroscopic data. We find that we can only fit the available data assuming a critically rotating central star. We identified two different disk models that simultaneously fit the spectroscopic, polarimetric, and interferometric observational constraints: a tenuous disk in hydrostatic equilibrium (i.e., with small scaleheight) and a smaller, scaleheight enhanced disk. We believe that these relatively small disks correspond to the transition region between the photosphere and the circumstellar environment, and that they are probably perturbed by some photospheric mechanism. The study of this interface between photosphere and circumstellar disk for near-critical rotators is crucial to our understanding of the Be phenomenon, and the mass and angular momentum loss of stars in general. This work shows that it is nowadays possible to directly study this transition region from simultaneous multi-technique observations.

We present the results of a high accuracy ($\sigma \approx 0.005%$) polarization monitoring of the Be Star Achernar that was carried out between July 7th and November 5th, 2006. Our results indicate that, after a near quiescent phase from 1998 to 2002, Achernar is presently in an active phase and has built a circumstellar disk. We detect variations both in the polarization level and position angle in timescales as short as one hour and as long as several weeks. Detailed modeling of the observed polarization strongly suggests that the short-term variations originate from discrete mass ejection events which produce transient inhomogeneities in the inner disk. Long-term variations, on the other hand, can be explained by the formation of an inner ring following one or several mass ejection events.

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$ Tau's 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.

We present optical $WBVR$ and infrared $JHKL$ photometric observations of the Be binary system $\delta$ Sco, obtained in 2000--2005, mid-infrared (10 and $18 \mu$m) photometry and optical ($\lambda\lambda$ 3200--10500 \AA) spectropolarimetry obtained in 2001. Our optical photometry confirms the results of much more frequent visual monitoring of $\delta$ Sco. In 2005, we detected a significant decrease in the object's brightness, both in optical and near-infrared brightness, which is associated with a continuous rise in the hydrogen line strenghts. We discuss possible causes for this phenomenon, which is difficult to explain in view of current models of Be star disks. The 2001 spectral energy distribution and polarization are succesfully modeled with a three-dimensional non-LTE Monte Carlo code which produces a self-consistent determination of the hydrogen level populations, electron temperature, and gas density for hot star disks. Our disk model is hydrostatically supported in the vertical direction and radially controlled by viscosity. Such a disk model has, essentially, only two free parameters, viz., the equatorial mass loss rate and the disk outer radius. We find that the primary companion is surrounded by a small (7 $R_\star$), geometrically-thin disk, which is highly non-isothermal and fully ionized. Our model requires an average equatorial mass loss rate of $1.5\times 10^{-9} M_{\sun}$ yr$^{-1}$.

We present self-consistent solutions for the disk structure of classical Be stars. Our disk model is hydrostatically supported in the vertical direction and the radial structure is governed by viscosity ($\alpha$-disks). We perform three-dimensional non-LTE Monte Carlo simulations to calculate simultaneously both the equilibrium temperature and Hydrogen level populations and to solve self-consistently for the density structure of the disk. We discuss the general properties of the solution for the disk structure and test our model against observations of $\delta$ Scorpii. Our results confirm that a viscous decretion disk model is consistent with these observations.

We discuss the nature of the circumstellar envelopes around the B[e] supergiants (B[e]SG) in the Magellanic Clouds (MC). Contrary to those in the Galaxy, the MC B[e]SG have a well defined luminosity and can be considered members of a well defined class. We discuss spectroscopy and optical broadband polarimetry and spectropolarimetry data. These data show for the first time detailed changes in the polarization across several spectral features. We show that the envelopes of the B[e]SG are generally variable. Broadband polarimetry data show that the envelopes are definitely non-spherically symmetric and large non-axisymmetric ejections may occur. In addition to that, spectropolarimetry is coming of age as a tool to study the B[e]SG envelope structure.

We present a 3-D NLTE 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 similiar to the Monte Carlo transition probability method developed by Lucy. Here we present a simplification of his method that avoids the use of the macro atom 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. We also find that the disks are fully ionized, as expected, but there is an ionization minimum in the vicinity of the temperature minimum. 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.

We study the effects of an annular gap induced by an embedded proto-planet on disk scattered light images and the infrared spectral energy distribution. We find that the outer edge of a gap is brighter in the scattered light images than a similar location in a gap-free disk. The stellar radiation that would have been scattered by material within in the gap is instead scattered by the disk wall at the outer edge of the gap, producing a bright ring surrounding the dark gap in the images. Given sufficient resolution, such gaps can be detected by the presence of this bright ring in scattered light images. A gap in a disk also changes the shape of the SED. Radiation that would have been absorbed by material in the gap is instead reprocessed by the outer gap wall. This leads to a decrease in the SED at wavelengths corresponding to the temperature at the radius of the missing gap material and a corresponding flux increase at longer wavelengths corresponding to the temperature of the outer wall. We note that, unlike an inner hole in the disk, the presence of an annular gap does not change the bolometric IR flux; it simply redistributes the radiation, previously produced by material within the gap, to longer wavelengths. This implies that the changes in the SED generally will be smaller for gaps than holes. Although it will be difficult on the basis of the SED alone to distinguish between the presence of a gap and other physical effects, the level of changes can be sufficiently large to be measurable with current instruments (e.g., Spitzer).

We investigate the linear polarization in the light of extrasolar planetary systems that may arise as a result of an occultation of the star by a transiting planet. Such an occultation breaks any spherical symmetry over the projected stellar disk and thus results in a non-vanishing linear polarization. This polarization will furthermore vary as the occultation progresses. We present both analytical and numerical results for the occultation of G-K-M-T dwarf stars by planets with sizes ranging from the one of Earth to two times the size of Jupiter. We find that the occultation polarization may result in an observable signal and provide additional means to characterize various parameters of the system. A particularly interesting result is that, for the later spectral types (i.e., smaller stellar radii), this polarization signature may be observable even for Earth-like planets. This suggests polarization as a possible tool to detect such planets. Departs from symmetry around mid-transit in the time dependence of the polarization signature may provide an estimate of the orbital eccentricity.

We present a model for the disk of the classical Be star $\zeta$ Tauri. The model consists of a Keplerian rotating disk with a power-law surface density and a vertical density distribution that follows from the balance between the thermal gas pressure and the z-component of the stellar gravitation. The opening angle of such a disk is not a fixed value but increases with the distance to the star (flared disk). We use a Monte Carlo code that solves simultaneously the thermal equilibrium, the statistical equilibrium and the hydrostatic equilibrium equations. The result is a self-consistent solution for the electron temperature, level populations and gas density distribution. It follows from this model that the only free parameter is the disk density. All other parameters (e.g. shape of the disk surface) are self-consistently calculated given the disk density. In this paper we present a successful fit for $\zeta$ Tauri data using this model.

We study the effects of dust grain size on the spectral energy distribution (SED) of spherical circumstellar envelopes. Based on the self-similarity relations of dusty SEDs derived by Ivezic & Elitzur (1997), we expect an approximate invariance of the IR SED for models with different grain sizes. Approximate invariance follows from the fact that differently sized grains have similar optical properties at long wavelengths where the dust reprocesses the starlight. In this paper, we discuss what are the physical requirements on the model parameters to maintain the approximate invariance of the IR SED. Single grain size models are studied for a wide range of grain sizes in three optical depth regimes. In this study, we find limits for the cases where the IR SED is and is not capable of conveying information about grain sizes, and to what extent it does so. We find that approximate invariance occurs for a much larger range of grain sizes than previously believed, and, when approximate invariance holds, the SED is controlled mainly by one parameter, the reprocessing optical depth, a quantity that measures the fraction of starlight that is absorbed by the dust grains. Models with a grain size distribution are studied as well. For these models, we find that, in many instances, the concept of approximate invariance may be extended from the IR SED to all wavelengths. This means that, for a wide range of optical depths, models with different grain size distributions will produce very similar SEDs and, hence, the reprocessing optical depth is the only quantity that can be unambiguously obtained from the SED. The observational consequences of this result are discussed in detail.

We have obtained linear polarization measurements of the Large Magellanic Cloud B[e] supergiant S 111 using optical imaging polarimetry. The intrinsic polarization found is consistent with the presence of an axisymmetric circumstellar envelope. We have additionally estimated the electron density for S 111 using data from the literature and revisited the correlation between polarization and envelope parameters of the B[e] supergiant stars using more recent IR calibration color data. The data suggest that the polarization can be indeed explained by electron scattering. We have used Monte Carlo codes to model the continuum polarization of the Magellanic B[e] supergiants. The results indicate that the electron density distribution in their envelopes is closer to a homogeneous distribution rather than an r$^{-2}$ dependence. At the same time, the data are best fitted by a spherical distribution with density contrast than a cylindrical distribution. The data and the model results support the idea of the presence of an equatorial disk and of the two-component wind model for the envelopes of the B[e] supergiants. Spectropolarimetry would help further our knowledge of these envelopes.