Search

We examine the evolution of the disk surrounding the Be star in the highly eccentric binary system $\delta$ Scorpii over its three most recent periastron passages. $V$-band and $B-V$ photometry, along with H$\alpha$ spectroscopy are combined with a new set of extensive multi-band polarimetry data to produce a detailed comparison of the disk's physical conditions during the time periods surrounding each closest approach of the secondary star. We use the three-dimensional Monte Carlo radiative transfer code \textsc{HDUST} and smoothed particle hydrodynamics (\textsc{SPH}) code to support our observations with models of disk evolution, discussing the behaviour of the H$\alpha$ and He\,\textsc{i} 6678 lines, $V$-band magnitude, and polarization degree. We compare the characteristics of the disk immediately before each periastron passage to create a baseline for the unperturbed disk. We find that the extent of the H$\alpha$ emitting region increased between each periastron passage, and that transient asymmetries in the disk become more pronounced with each successive encounter. Asymmetries of the H$\alpha$ and He\,\textsc{i} 6678 lines in 2011 indicate that perturbations propagate inward through the disk near periastron. When the disk's direction of orbit is opposite to that of the secondary, the parameters used in our models do not produce spiral density enhancements in the H$\alpha$ emitting region because the tidal interaction time is short due to the relative velocities of the disk particles with the secondary. The effects of the secondary star on the disk are short-lived and the disk shows independent evolution between each periastron event., Comment: 20 pages, 13 figures, accepted for publication in ApJ

Altough both components of the bright eclipsing binary $\alpha$ Dra having been resolved using long baseline interferometry and the secondary component shown to contribute some 15\% of the total flux, a spectroscopic detection of the companion star was so far unsuccessful. To achieve our goals, we use a combined data set from interferometry with the Navy Precision Optical Interferometer (NPOI), photometry with the TESS space observatory, and high-resolution spectroscopy with the HERMES fibre-fed spectrograph at the La Palma observatory. We use the method of spectral disentangling to search for the contribution of a companion star in the observed composite HERMES spectra, to separate the spectral contributions of both components, and to determine orbital elements of the $\alpha$ Dra system. TESS light curves are analysed in an iterative fashion with spectroscopic inference of stellar atmospheric parameters to determine fundamental stellar properties and their uncertainties. Finally, NPOI interferometric measurements are used for determination of the orbital parameters of the system and angular diameters of both binary components. We report the first firm spectroscopic detection of the secondary component in $\alpha$ Dra and deliver disentangled spectra of both binary components. The inferred near-core mixing properties of both components do not support a dependence of the convective core overshooting on the stellar mass. At the same time, the $\alpha$ Dra system provides extra support to hypothesise that the mass discrepancy in eclipsing spectroscopic double-lined binaries is associated with inferior atmospheric modelling of intermediate- and high-mass stars, and less so with the predictive power of stellar structure and evolution models as to the amount of near-core mixing and mass of the convective core. (abridged abstract), Comment: 23 pages, 15 figures, 8 tables, accepted for publication in Astronomy & Astrophysics

We present the results of a multiplicity survey for a magnitude-limited sample of 31 classical Be stars conducted with the Navy Precision Optical Interferometer and the Mark III Stellar Interferometer. The interferometric observations were used to detect companions in ten previously known binary systems. For two of these sources (66 Oph and $\beta$ Cep) new orbital solutions were obtained, while for a third source ($\upsilon$ Sgr) our observations provide the first direct, visual detection of the hot sdO companion to the Be primary star. Combining our interferometric observations with an extensive literature search, we conclude that an additional four sources (o Cas, 15 Mon, $\beta$ Lyr, and $\beta$ Cep) also contain wider binary components that are physical companions to the narrow binaries, thus forming hierarchical multiple systems. Among the sources not previously confirmed as spectroscopic or visual binaries, BK Cam was resolved on a number of nights within a close physical proximity of another star with relative motion possibly suggesting a physical binary. Combining our interferometric observations with an extensive literature search, we provide a detailed listing of companions known around each star in the sample, and discuss the multiplicity frequency in the sample. We also discuss the prospects for future multiplicity studies of classical Be stars by long baseline optical interferometry., Comment: 58 pages, 14 figures, 8 tables

We demonstrate that the angle between star's rotation axis and the observer's line-of-sight, usually called the inclination angle, can be reliably determined for Be stars via H$\alpha$ emission-line profile fitting. We test our method on a sample of 11 Be stars with available long-baseline interferometric data from the Navy Precision Optical Interferometer~(NPOI). We fit the H$\alpha$ emission line profile of each star to obtain a spectroscopic inclination angle $i_{\rm H\alpha}$. We then obtain an independent inclination angle estimate, $i_{\rm V^2}$, by fitting the observed interferometric visibilities with model visibilities based on a purely geometric representation of the light distribution on the sky. The sample differences, $\Delta i \equiv i_{\rm H\alpha} - i_{\rm V^2}$, are normally distributed with a mean of zero and a standard deviation of $6.7$ degrees, and the linear correlation coefficient between $i_{\rm H\alpha}$ and $i_{\rm V^2}$ is $r=0.93$. As Be stars comprise upwards of one fifth of all main-sequence B-type stars, this H$\alpha$ line profile fitting technique has the potential to provide an efficient method for detecting correlated stellar spin axes in young open clusters. Furthermore, if the orientation of the Be star circumstellar disk on the plane of the sky can be constrained by polarization measurements, it is possible to determine the full 3D stellar rotation vector of each Be star., Comment: 19 pages, 23 figures, accepted for publication in The Astrophysical Journal

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., Comment: 16 pages, 12 figures, 4 tables, accepted for publishing in ApJ

We utilize a multi-step modelling process to produce synthetic interferometric and spectroscopic observables, which are then compared to their observed counterparts. Our extensive set of interferometric observations of the Be star 48 Per, totaling 291 data points, were obtained at the Navy Precision Optical Interferometer from 2006 Nov 07 to 2006 Nov 23. Our models were further constrained by comparison with contemporaneous H$\alpha$ line spectroscopy obtained at the John S. Hall Telescope at Lowell Observatory on 2006 Nov 1. Theoretical spectral energy distributions, SEDs, for 48~Per were confirmed by comparison with observations over a wavelength regime of 0.4 to 60 microns from Touhami et al. (2010) and Vieira et al. (2017). Our best-fitting combined model from H$\alpha$ spectroscopy, H$\alpha$ interferometry and SED fitting has a power law density fall off, $n$, of 2.3 and an initial density at the stellar surface of $\rho_0 = 1.0 \times 10^{-11} \rm{g\,cm}^{-3}$ with a inclination constrained by H$\alpha$ spectroscopy and interferometry of 45$^{\circ} \pm 5^{\circ}$. The position angle for the system, measured east from north, is 121$^{\circ}$ $\pm$ 1$^{\circ}$. Our best-fit model shows that the disk emission originates in a moderately large disk with a radius of 25 R$_{*}$ which is consistent with a disk mass of approximately 5 $\times$ 10$^{24}$ g or 3 $\times$ 10$^{-10}$ M$_{*}$. Finally, we compare our results with previous studies of 48~Per by Quirrenbach et al. (1997) with agreement but find, with our much larger data set, that our disk size contradicts the findings of Delaa et al. (2011)., Comment: 26 pages, 9 figures

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.

Omicron Aquarii is late-type, Be shell star with a stable and nearly symmetric H$\alpha$ emission line. We combine H$\alpha$ interferometric observations obtained with the Navy Precision Optical Interferometer (NPOI) covering 2007 through 2014 with H$\alpha$ spectroscopic observations over the same period and a 2008 observation of the system's near-infrared spectral energy distribution to constrain the properties of $o$~Aqr's circumstellar disk. All observations are consistent with a circumstellar disk seen at an inclination of $75\pm\,3^{\circ}$ with a position angle on the sky of $110\pm\,8^{\circ}$ measured E from N. From the best-fit disk density model, we find that 90\% of the H$\alpha$ emission arises from within $9.5$ stellar radii, and the mass associated with this H$\alpha$ disk is $\sim 1.8\times10^{-10}$ of the stellar mass and the associated angular momentum, assuming Keplerian rotation for the disk, is $\sim 1.6\times10^{-8}$ of the total stellar angular momentum. The occurrence of a central quasi-emission (CQE) feature in Mg\,{\sc ii} $\lambda\,4481$ is also predicted by this best-fit disk model and the computed profile compares successfully with observations from 1999. To obtain consistency between the H$\alpha$ line profile modelling and the other constraints, it was necessary in the profile fitting to weight the line core (emission peaks and central depression) more heavily than the line wings, which were not well reproduced by our models. This may reflect the limitation of assuming a single power-law for the disk's equatorial density variation. The best-fit disk density model for $o$~Aqr predicts that H$\alpha$ is near its maximum strength as a function of disk density, and hence the H$\alpha$ equivalent width and line profile change only weakly in response to large (factor of $\sim 5$) changes in the disk density., Comment: 16 pages; 16 figures

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., Comment: 19 pages (35 including online material), 17 figures, 2 online figures, 2 online tables with data

We acquired Johnson BV photometry of the binary Be disk system delta~Scorpii during its 2009, 2010, 2011, and 2012 observing seasons and used it to probe the innermost regions of the disk. We found that several disk building events have occurred during this time, resulting in an overall brightening in the V-band and reddening of the system. In addition to these long-term trends, we found cyclical variability in each observing season on timescales between 60 and 100 days. We were able to reproduce the changes in the magnitude and colour of delta Sco using our theoretical models and found that variable mass-loss rates in the range 2.5-7.0x10^{-9} M_{sun}/yr over ~35 days can reproduce the observed increase in brightness., Comment: 15 pages, emulateapj

Binary post-AGB stars are interesting laboratories to study both the evolution of binaries as well as the structure of circumstellar disks. A multiwavelength high angular resolution study of the prototypical object 89 Herculis is performed with the aim of identifying and locating the different emission components seen in the SED. A large interferometric data set, collected over the past decade and covering optical and near-IR wavelengths, is analyzed with simple geometric models. Combining the interferometric constraints with the photometry and the optical spectra, we reassess the energy budget of the post-AGB star and its circumstellar environment. We report the first (direct) detection of a large (35-40%) optical circumstellar flux contribution and spatially resolve its emission region. Given this large amount of reprocessed and/or redistributed optical light, the fitted size of the emission region is rather compact and fits with(in) the inner rim of the circumbinary dust disk. This rim dominates our K band data through thermal emission and is rather compact, emitting significantly already at a radius of twice the orbital separation. We interpret the circumstellar optical flux as due to a scattering process, with the scatterers located in the extremely puffed-up inner rim of the disk and possibly also in a bipolar outflow seen pole-on. A non-LTE gaseous origin in an inner disk cannot be excluded but is considered highly unlikely. This direct detection of a significant amount of circumbinary light at optical wavelengths poses several significant questions regarding our understanding of both post-AGB binaries and the physics in their circumbinary disks. Although the identification of the source of emission/scattering remains inconclusive without further study on this and similar objects, the implications are manifold., Comment: Accepted for publication in A&A, 16 pages, 15 figures

Focusing on B-emission stars, we investigated a set of H$\alpha$ equivalent widths calculated from observed spectra acquired over a period of about 4 years from 2003 to 2007. During this time, changes in equivalent width for our program stars were monitored. We have found a simple statistical method to quantify these changes in our observations. This statistical test, commonly called the F ratio, involves calculating the ratio of the external and internal error. We show that the application of this technique can be used to place bounds on the degree of variability of Be stars. This observational tool provides a quantitative way to find Be stars at particular stages of variability requiring relatively little observational data., Comment: 8 figures accepted by The Astronomical Journal

In anticipation of the possible collision between a circumstellar disk and the secondary star in the highly eccentric binary system Delta Scorpii, high angular resolution interferometric observations have been acquired aimed at revising the binary parameters. The Navy Prototype Optical Interferometer (NPOI) was used to spatially resolve the binary components in 2000 and over a period between 2005 and 2010. The interferometric observations are used to obtain the angular separations and orientations of the two stellar components at all epochs for which data has been obtained, including 2005 and 2006, for which based on previous studies there was some uncertainty as to if the signature of binarity can be clearly detected. The results of this study represent the most complete and accurate coverage of the binary orbit of this system to date and allow for the revised timing of the upcoming periastron passage that will occur in 2011 to be obtained., Comment: Accepted to ApJ Letters, 10 pages, 3 figures, 1 table

Interacting binaries typically have separations in the milli-arcsecond regime and hence it has been challenging to resolve them at any wavelength. However, recent advances in optical interferometry have improved our ability to discern the components in these systems and have now enabled the direct determination of physical parameters. We used the Navy Prototype Optical Interferometer to produce for the first time images resolving all three components in the well-known Algol triple system. Specifically, we have separated the tertiary component from the binary and simultaneously resolved the eclipsing binary pair, which represents the nearest and brightest eclipsing binary in the sky. We present revised orbital elements for the triple system, and we have rectified the 180-degree ambiguity in the position angle of Algol C. Our directly determined magnitude differences and masses for this triple star system are consistent with earlier light curve modeling results., Comment: Accepted to Astrophysical Journal Letters, 15 pages, 2 eps figures.

High spatial resolution observations of the H\alpha-emitting wind structure associated with the Luminous Blue Variable star P Cygni were obtained with the Navy Prototype Optical Interferometer (NPOI). These observations represent the most comprehensive interferometric data set on P Cyg to date. We demonstrate how the apparent size of the H\alpha-emitting region of the wind structure of P Cyg compares between the 2005, 2007 and 2008 observing seasons and how this relates to the H\alpha line spectroscopy. Using the data sets from 2005, 2007 and 2008 observing seasons, we fit a circularly symmetric Gaussian model to the interferometric signature from the H\alpha-emitting wind structure of P Cyg. Based on our results we conclude that the radial extent of the H\alpha-emitting wind structure around P Cyg is stable at the 10% level. We also show how the radial distribution of the H\alpha flux from the wind structure deviates from a Gaussian shape, whereas a two-component Gaussian model is sufficient to fully describe the H\alpha-emitting region around P Cyg., Comment: 27 pages, 6 figues, accepted for publication in AJ

Recovering images from optical interferometric observations is one of the major challenges in the field. Unlike the case of observations at radio wavelengths, in the optical the atmospheric turbulence changes the phases on a very short time scale, which results in corrupted phase measurements. In order to overcome these limitations, several groups developed image reconstruction techniques based only on squared visibility and closure phase information, which are unaffected by atmospheric turbulence. We present the results of two techniques used by our group, which employed coherently integrated data from the Navy Prototype Optical Interferometer. Based on these techniques we were able to recover complex visibilities for several sources and image them using standard radio imaging software. We describe these techniques, the corrections applied to the data, present the images of a few sources, and discuss the implications of these results., Comment: 8 pages, 2008 SPIE Astronomical Telescopes and Instrumentation, Eds. M. Schoeller, W. C. Danchi, F. Delplancke

We present a numerical model describing a circularly symmetric gaseous disk around the Be star chi Ophiuchi. The model is constrained by long-baseline interferometric observations that are sensitive to the H-alpha Balmer line emission from the disk. For the first time our interferometric observations spatially resolve the inner region of the circumstellar disk around chi Ophiuchi and we use these results to place a constraint on the physical extent of the H-alpha-emitting region. We demonstrate how this in turn results in very specific constraints on the parameters that describe the variation of the gas density as a function of radial distance from the central star., Comment: Accepted for publication in ApJ

In this paper we use coherently integrated visibilities (see separate paper in these proceedings, Jorgensen et al. 2008) to measure the properties of binary stars. We use only the phase of the complex visibility and not the amplitude. The reason for this is that amplitudes suffer from the calibration effect (the same for coherent and incoherent averages) and thus effectively provide lower accuracy measurements. We demonstrate that the baseline phase alone can be used to measure the separation, orientation and brightness ratio of a binary star, as a function of wavelength., Comment: 2008 SPIE Astronomical Telescopes and Instrumentation

We have computed theoretical models of circumstellar disks for the classical Be stars $\kappa$ Dra, $\beta$ Psc, and $\upsilon$ Cyg. Models were constructed using a non-LTE radiative transfer code developed by \citet{sig07} which incorporates a number of improvements over previous treatments of the disk thermal structure, including a realistic chemical composition. Our models are constrained by direct comparison with long baseline optical interferometric observations of the H$\alpha$ emitting regions and by contemporaneous H$\alpha$ line profiles. Detailed comparisons of our predictions with H$\alpha$ interferometry and spectroscopy place very tight constraints on the density distributions for these circumstellar disks., Comment: 10 figures,28 pages, accepted by ApJ

We present infrared photometric observations of the Be binary system delta Scorpii obtained in 2006. The J,H and K magnitudes are the same within the errors compared to observations taken 10 months earlier. We derive the infrared excess from the observation and compare this to the color excess predicted by a radiative equilibrium model of the primary star and its circumstellar disk. We use a non-LTE computational code to model the gaseous envelope concentrated in the star's equatorial plane and calculate the expected spectral energy distribution and Halpha emission profile of the star with its circumstellar disk. Using the observed infrared excess of delta Sco, as well as Halpha spectroscopy bracketing the IR observations in time, we place constraints on the radial density distribution in the circumstellar disk. Because the disk exhibits variability in its density distribution, this work will be helpful in understanding its dynamics., Comment: 12 pages, 14 figures, to be published in PASP May 2008

We present the results of an experiment to image the interacting binary star beta Lyrae with data from the Navy Prototype Optical Interferometer (NPOI), using a differential phase technique to correct for the effects of the instrument and atmosphere on the interferometer phases. We take advantage of the fact that the visual primary of beta Lyrae and the visibility calibrator we used are both nearly unresolved and nearly centrally symmetric, and consequently have interferometric phases near zero. We used this property to detect and correct for the effects of the instrument and atmosphere on the phases of beta Lyrae and to obtain differential phases in the channel containing the Halpha emission line. Combining the Halpha-channel phases with information about the line strength, we recovered complex visibilities and imaged the Halpha emission using standard radio interferometry methods. We find that the results from our differential phase technique are consistent with those obtained from a more-standard analysis using squared visibilities (V^2's). Our images show the position of the Halpha emitting regions relative to the continuum photocenter as a function of orbital phase and indicate that the major axis of the orbit is oriented along p.a.=248.8+/-1.7 deg. The orbit is smaller than previously predicted, a discrepancy that can be alleviated if we assume that the system is at a larger distance from us, or that the contribution of the stellar continuum to the Halpha channel is larger than estimated. Finally, we also detected a differential phase signal in the channels containing HeI emission lines at 587.6 and 706.5nm, with orbital behavior different from that of the Halpha, indicating that it originates from a different part of this interacting system., Comment: Submitted to ApJ

With the Navy Prototype Optical Interferometer (NPOI), the binary system Theta 1 Orionis C, the most massive member of the Trapezium, was spatially resolved over a time period extending from February 2006 to March 2007. The data show significant orbital motion over the 14 months, and, after combining the NPOI data with previous measurements of the system from the literature, the observations span 10 years of the orbit. Our results indicate that the secondary did not experience an unusually close periastron passage this year, in contradiction to the prediction of a recently published, highly eccentric ~11 year orbit. Future observations of this source will be required to improve the orbital solution. Possible implications of the results in terms of system distance are discussed, although a main conclusion of this work is that a definitive orbit solution will require more time to obtain sufficient phase coverage, and that the interaction effects expected at periastron did not occur in 2007., Comment: 13 pages, 2 tables, 3 figures. Accepted for publication in ApJL

The secondary photometric standard star #79 for the FS Aur field (Henden & Honeycutt 1997) designated as [HH97] FS Aur-79 (GSC 1874 399) is a short period (0.2508 days) eclipsing binary whose light curve is a combination of the $\beta$ Lyr and BY Dra type variables. High signal-to-noise multi-color photometry were obtained using the USNO 1-m telescope. These light curves show asymmetry at quadrature phases (O'Connell effect), which can be modeled with the presence of star spots. A low resolution spectrum obtained with the 3.5-m WIYN telescope at orbital phase 0.76 is consistent with a spectral type of dK7e and dM3e. A radial velocity curve for the primary star was constructed using twenty-four high resolution spectra from the 9.2 m HET. Spectra show H-alpha and H-beta in emission confirming chromospheric activity and possibly the presence of circumstellar material. Binary star models that simultaneously fit the U, B, V, R and RV curves are those with a primary star of mass 0.59+-0.02 Msun, temperature 4100+-25 K, mean radius of 0.67 Rsun, just filling its Roche lobe and a secondary star of mass 0.31+-0.09 Msun, temperature 3425+-25 K, mean radius of 0.48 Rsun, just within its Roche lobe. An inclination angle of 83+-2 degrees with a center of mass separation of 1.62 Rsun is also derived. Star spots, expected for a rotation period of less than a day, had to be included in the modeling to fit the O'Connell effect.

Interferometric observations of two well-known Be stars, gamma Cas and phi Per, were collected and analyzed to determine the spatial characteristics of their circumstellar regions. The observations were obtained using the Navy Prototype Optical Interferometer equipped with custom-made narrowband filters. The filters isolate the H-alpha emission line from the nearby continuum radiation, which results in an increased contrast between the interferometric signature due to the H-alpha-emitting circumstellar region and the central star. Because the narrowband filters do not significantly attenuate the continuum radiation at wavelengths 50 nm or more away from the line, the interferometric signal in the H-alpha channel is calibrated with respect to the continuum channels. The observations used in this study represent the highest spatial resolution measurements of the H-alpha-emitting regions of Be stars obtained to date. These observations allow us to demonstrate for the first time that the intensity distribution in the circumstellar region of a Be star cannot be represented by uniform disk or ring-like structures, whereas a Gaussian intensity distribution appears to be fully consistent with our observations., Comment: 23 pages, 14 figures, accepted for publication in AJ