Be stars have generally been characterized by the emission lines in their spectra, and especially the time variability of those spectroscopic features. They are known to also exhibit photometric variability at multiple timescales, but have not been broadly compared and analyzed by that behavior. We have taken advantage of the advent of wide-field, long-baseline, and high-cadence photometric surveys that search for transiting exoplanets to perform a comprehensive analysis of brightness variations among a large number of known Be stars. The photometric data comes from the KELT transit survey, with a typical cadence of 30 minutes, baseline of up to ten years, photometric precision of about 1%, and coverage of about 60% of the sky. We analyze KELT light curves of 610 known Be stars in both the Northern and Southern hemispheres in an effort to study their variability. Consistent with other studies of Be star variability, we find most of the stars to be photometrically variable. We derive lower limits on the fraction of stars in our sample that exhibit features consistent with non-radial pulsations (25%), outbursts (36%), and long term trends in the circumstellar disk (37%), and show how these are correlated with spectral sub-type. Other types of variability, such as those owing to binarity, are also explored. Simultaneous spectroscopy for some of these systems from the Be Star Spectral Database (BeSS) allow us to better understand the physical causes for the observed variability, especially in cases of outbursts and changes in the disk.

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.

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

Supernova (SN) 2015bh (or SNhunt275) was discovered in NGC 2770 on 2015 February with an absolute magnitude of Mr ~ -13.4 mag, and was initially classified as a SN impostor. Here we present the photometric and spectroscopic evolution of SN 2015bh from discovery to late phases (~ 1 yr after). In addition, we inspect archival images of the host galaxy up to ~ 21 yr before discovery, finding a burst ~ 1 yr before discovery, and further signatures of stellar instability until late 2014. Later on, the luminosity of the transient slowly increases, and a broad light curve peak is reached after about three months. We propose that the transient discovered in early 2015 could be a core-collapse SN explosion. The pre-SN luminosity variability history, the long-lasting rise and faintness first light curve peak suggests that the progenitor was a very massive, unstable and blue star, which exploded as a faint SN because of severe fallback of material. Later on, the object experiences a sudden brightening of 3 mag, which results from the interaction of the SN ejecta with circumstellar material formed through repeated past mass-loss events. Spectroscopic signatures of interaction are however visible at all epochs. A similar chain of events was previously proposed for the similar interacting SN 2009ip.

We present the results of 938 speckle measures of double stars and suspected double stars drawn mainly from the Hipparcos Catalogue, as well as 208 observations where no companion was noted. One hundred fourteen pairs have been resolved for the first time. The data were obtained during four observing runs in 2014 using the Differential Speckle Survey Instrument (DSSI) at Lowell Observatory's Discovery Channel Telescope. The measurement precision obtained when comparing to ephemeris positions of binaries with very well-known orbits is generally less than 2 mas in separation and 0.5 degrees in position angle. Differential photometry is found to have internal precision of approximately 0.1 magnitudes and to be in very good agreement with Hipparcos measures in cases where the comparison is most relevant. We also estimate the detection limit in the cases where no companion was found. Visual orbital elements are derived for 6 systems.

The eruptive variable V838 Monocerotis gained notoriety in 2002 when it brightened nine magnitudes in a series of three outbursts and then rapidly evolved into an extremely cool supergiant. We present optical, near-IR, and mid-IR spectroscopic and photometric observations of V838 Monocerotis obtained between 2008 and 2012 at the Apache Point Observatory 3.5m, NASA IRTF 3m, and Gemini South 8m telescopes. We contemporaneously analyze the optical & IR spectroscopic properties of V838 Monocerotis to arrive at a revised spectral type L3 supergiant and effective temperature Teff~2000--2200 K. Because there are no existing optical observational data for L supergiants in the optical, we speculate that V838 Monocerotis may represent the prototype for L supergiants in this wavelength regime. We find a low level of Halpha emission present in the system, consistent with interaction between V838 Monocerotis and its B3V binary; however, we cannot rule out a stellar collision as the genesis event, which could result in the observed Halpha activity. Based upon a two-component blackbody fit to all wavelengths of our data, we conclude that, as of 2009, a shell of ejecta surrounded V838 Monocerotis at a radius of R=263+/-10 AU with a temperature of T=285+/-2 K. This result is consistent with IR interferometric observations from the same era and predictions from the Lynch et al. model of the expanding system, which provides a simple framework for understanding this complicated system.

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.

The work is aimed at a study of the circumstellar disk of the bright classical binary Be star {\pi} Aqr. We analysed variations of a double-peaked profile of the H{\alpha} emission line in the spectrum of {\pi} Aqr that was observed in many phases during ~40 orbital cycles in 2004--2013. We applied the Discrete Fourier Transform (DFT) method to search for periodicity in the peak intensity (V/R) ratio. Doppler tomography was used to study the structure of the disk around the primary. The dominant frequency in the power spectrum of the H{\alpha} V/R ratio is 0.011873 day^-1 that correspond to a period of 84.2(2) days and is in agreement with the earlier determined orbital period of the system, Porb=84.1 days. The V/R ratio shows a sinusoidal variation phase-locked with the orbital period. Doppler maps of all our spectra show a non-uniform structure of the disk around the primary: a ring with the inner and outer radii at Vin~ 450 km/s and Vout~ 200km/s, respectively, along with an extended stable region (spot) at V_x ~ 225 km/s and V_y~100 km/s. The disk radius of ~ 65 Rsun = 0.33 AU was estimated assuming Keplerian motion of a particle on a circular orbit at the disk outer edge.

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.

\gamma Cas is the prototypical classical Be star and is best known for its variable hard X-ray emission. To elucidate the reasons for this emission, we mounted a multiwavelength campaign in 2010 centered around 4 XMM observations. The observational techniques included long baseline optical interferometry (LBOI), monitoring by an Automated Photometric Telescope and Halpha observations. Because gamma Cas is also known to be in a binary, we measured Halpha radial velocities and redetermined its period as 203.55+/-0.2 days and an eccentricity near zero. The LBOI observations suggest that the star's decretion disk was axisymmetric in 2010, has an inclination angle near 45^o, and a larger radius than previously reported. The Be star began an "outburst" at the beginning of our campaign, made visible by a disk brightening and reddening during our campaign. Our analyses of the new high resolution spectra disclosed many attributes found from spectra obtained in 2001 (Chandra) and 2004 (XMM). As well as a dominant hot 14 keV thermal component, these familiar ones included: (i) a fluorescent feature of Fe K stronger than observed at previous times, (ii) strong lines of N VII and Ne XI lines indicative of overabundances, and (iii) a subsolar Fe abundance from K-shell lines but a solar abundance from L-shell ions. We also found that 2 absorption columns are required to fit the continuum. While the first one maintained its historical average of 1X10^21 cm^-2, the second was very large and doubled to 7.4X10^23 cm^-2 during our X-ray observations. Although we found no clear relation between this column density and orbital phase, it correlates well with the disk brightening and reddening both in the 2010 and earlier observations. Thus, the inference from this study is that much (perhaps all?) of the X-ray emission from this source originates behind matter ejected by gamma Cas into our line of sight.

We present near-IR (J,H,Ks) photometry for 27 of the 28 candidate Herbig Ae/Be stars in the Small and Large Magellanic Clouds identified via the EROS1 and EROS2 surveys as well as near-contemporaneous optical (H-alpha) spectroscopy for 21 of these 28 candidates. Our observations extend previous efforts to determine the evolutionary status of these objects. We compare the IR brightness and colors of a subset of our sample with archival ground-based IR data and find evidence of statistically significant photometric differences for ELHC 5, 7, 12, 18, and 21 in one or more filter. In all cases, these near-IR photometric variations exhibit a grey color as compared to earlier epoch data. The ~1 magnitude IR brightening and minimal change in the H-alpha emission strength we observe in ELHC 7 is consistent with previous claims that it is a UX Ori type HAe/Be star, which is occasionally obscurred by dust clouds. We also detect a ~1 magnitude IR brightening of ELHC 12, but find little evidence of a similar large-scale change in its H-alpha line strength, suggesting that its behavior could also be caused by a UX Ori-like event. The ~0.5 magnitude IR variability we observe for ELHC 21, which also exhibited little evidence of a change in its H-alpha emission strength, could conceivably be caused by a major recent enhancement in the density of the inner disk region of a classical Be star. We also report the first near-IR photometry for two ESHC stars and the first H-alpha spectroscopy for one ELHC and five ESHC stars. Although H-alpha emission is detected in all of these new observations, they do not exhibit a strong near-IR excess. It is therefore possible that many of these objects may be classical Be stars rather than Herbig Ae/Be stars.

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.

(Abridged) Classical Be stars occasionally transition from having a gaseous circumstellar disk (''Be phase'') to a state in which all observational evidence for the presence of these disks disappears (''normal B-star phase''). We present one of the most comprehensive spectropolarimetric views to date of such a transition for two Be stars, pi Aquarii and 60 Cygni. 60 Cyg's disk loss episode was characterized by a monotonic decrease in emission strength over a time-scale of 1000 days, consistent with the viscous time-scale of the disk, assuming alpha is 0.14. pi Aqr's disk loss was episodic in nature and occurred over a time-scale of 2440 days. An observed time lag between the behavior of the polarization and H-alpha in both stars indicates the disk clearing proceeded in an ''inside-out'' manner. We determine the position angle of the intrinsic polarization to be 166.7 +/- 0.1 degrees for pi Aqr and 107.7 +/- 0.4 degrees for 60 Cyg, and model the observed polarization during the quiescent diskless phase of each star to determine the interstellar polarization along the line of sight. Minor outbursts observed during the quiescent phase of each star shared similar lifetimes as those previously reported for mu Cen, suggesting that the outbursts represent the injection and subsequent viscous dissipation of individual blobs of material into the inner circumstellar environments of these stars. We also observe deviations from the mean intrinsic polarization position angle during polarization outbursts in each star, indicating deviations from axisymmetry. We propose that these deviations might be indicative of the injection (and subsequent circularization) of new blobs into the inner disk, either in the plane of the bulk of the disk material or in a slightly inclined (non-coplanar) orbit.

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.

Science opportunities and recommendations concerning optical/infrared polarimetry for the upcoming decade in the fields of stars and stellar evolution. Community-based White Paper to Astro2010 in response to the call for such papers.

We report high spatial resolution 11.2 and 18.1 micron imaging of the eruptive variable V838 Monocerotis, obtained with Gemini Observatory's Michelle in 2007 March. The 2007 flux density of the unresolved stellar core is roughly 2 times brighter than that observed in 2004. We interpret these data as evidence that V838 Mon has experienced a new circumstellar dust creation event. We also report a gap of spatially extended thermal emission over radial distances of 1860-93000 AU from the central source, which suggests that no prior significant circumstellar dust production events have occurred within the past 900-1500 years.

We report high spatial resolution 11.2 and 18.1 micron imaging of V838 Monocerotis obtained with Gemini Observatory's Michelle instrument in 2007 March. Strong emission is observed from the unresolved stellar core of V838 Mon in our Gemini imagery, and is confirmed by Spitzer MIPS 24 micron imaging obtained in 2007 April. The 2007 flux density of the unresolved mid-infrared emission component is 2 times brighter than that observed in 2004. No clear change in the net amount of 24 micron extended emission is observed between the 2004 and 2007 epoch Spitzer imagery. We interpret these data as evidence that V838 Mon has experienced a new circumstellar dust creation event. We suggest that this newly created dust has condensed from the expanding ejecta produced from V838 Mon's 2002 outburst events, and is most likely clumpy. We speculate that one (or more) of these clumps might have passed through the line-of-sight in late 2006, producing the brief multi-wavelength photometric event reported by Bond (2006) and Munari et al (2007). We detect no evidence of extended emission above a level of 1 mJy at 11.2 microns and 7 mJy at 18.1 microns over radial distances of 1860 - 93000 AU (0.3-15.0 arcsec) from the central source. Using the simple assumption that ejecta material expands at a constant velocity of 300-500 km/s, this gap of thermal emission suggests that no significant prior circumstellar dust production events have occurred within the past 900-1500 years.

We present new optical spectroscopic and archival Spitzer IRAC photometric observations of a B-type star in the SMC cluster NGC 346, NGC 346:KWBBe 200. We detect numerous Fe II, [O I], and [Fe II] lines, as well as strong P-Cygni profile H I emission lines in its optical spectrum. The star's near-IR color and optical to IR SED clearly indicate the presence of an infrared excess, consistent with the presence of gas and warm, T ~800 K, circumstellar dust. Based on a crude estimate of the star's luminosity and the observed spectroscopic line profile morphologies, we find that the star is likely to be a B-type supergiant. We suggest that NGC 346:KWBBe 200 is a newly discovered B[e] supergiant star, and represents the fifth such object to be identified in the SMC.

(Abridged version) We present the first detailed imaging polarization observations of six SMC and six LMC clusters, known to have large populations of B-type stars which exhibit excess H-alpha emission, to constrain the evolutionary status of these stars and hence better establish links between the onset of disk formation in classical Be stars and cluster age and/or metallicity. The wavelength dependence of our intrinsic polarization data provides a diagnostic of the dominant and any secondary polarigenic agents present, enabling us to discriminate pure gas disk systems, i.e. classical Be stars, from composite gas plus dust disk systems, i.e. Herbig Ae/Be or B[e] stars. Our intrinsic polarization results, along with available near-IR color information, strongly supports the suggestion of Wisniewski et al. that classical Be stars are present in clusters of age 5-8 Myr, and contradict assertions that the Be phenomenon only develops in the second half of a B star's main sequence lifetime, i.e. no earlier than 10 Myr. Comparing the polarimetric properties of our dataset to a similar survey of Galactic classical Be stars, we find that the prevalence of polarimetric Balmer jump signatures decreases with metallicity. We speculate that these results might indicate that either it is more difficult to form large disk systems in low metallicity environments, or that the average disk temperature is higher in these low metallicity environments. We have characterized the polarimetric signatures of all candidate Be stars in our data sample and find ~25% are unlikely to arise from true classical Be star-disk systems.

We present U,B,V,R,I imaging polarimetry of NGC 2100 and its surrounding environment, which comprise a part of the LMC 2 supershell. The morphology of the observed position angle distribution provides a tracer of the projected magnetic field in this environment. Our polarization maps detail regions exhibiting similarly aligned polarization position angles, as well as more complex position angle patterns. We observe regions of coherent fields on spatial scales of 42 x 24 pc to 104 x 83 pc, and infer projected field strengths of $\sim14-30\mu$G. We propose that the superposition of global outflows from the LMC 2 environment, as well as outflows created within NGC 2100, produce the unique field geometry in the region.

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 the first K'-band, long-baseline interferometric observations of the northern Be stars gamma Cas, phi Per, zeta Tau, and kappa Dra. The measurements were made with multiple telescope pairs of the CHARA Array interferometer, and in every case the observations indicate that the circumstellar disks of the targets are resolved. We fit the interferometric visibilities with predictions from a simple disk model that assumes an isothermal gas in Keplerian rotation. We derive fits of the four model parameters (disk base density, radial density exponent, disk normal inclination, and position angle) for each of the targets. The resulting densities are in broad agreement with prior studies of the IR excess flux and the resulting orientations generally agree with those from interferometric H-alpha and continuum polarimetric observations. We find that the angular size of the K' disk emission is smaller than that determined for the H-alpha emission, and we argue that the difference is the result of a larger H-alpha opacity and the relatively larger neutral hydrogen fraction with increasing disk radius. All the targets are known binaries with faint companions, and we find that companions appear to influence the interferometric visibilities in the cases of phi Per and kappa Dra. We also present contemporaneous observations of the H-alpha, H-gamma, and Br-gamma emission lines. Synthetic model profiles of these lines that are based on the same disk inclination and radial density exponent as derived from the CHARA Array observations match the observed emission line strength if the disk base density is reduced by approximately 1.7 dex.

We present B, V, R, and H alpha photometry of 8 clusters in the Small Magellanic Cloud, 5 in the Large Magellanic Cloud, and 3 Galactic clusters, and use 2 color diagrams (2-CDs) to identify candidate Be star populations in these clusters. We find evidence that the Be phenomenon is enhanced in low metallicity environments, based on the observed fractional early-type candidate Be star content of clusters of age 10-25 Myr. Numerous candidate Be stars of spectral types B0 to B5 were identified in clusters of age 5-8 Myr, challenging the suggestion of Fabregat & Torrejon (2000) that classical Be stars should only be found in clusters at least 10 Myr old. These results suggest that a significant number of B-type stars must emerge onto the zero-age-main-sequence as rapid rotators. We also detect an enhancement in the fractional content of early-type candidate Be stars in clusters of age 10-25 Myr, suggesting that the Be phenomenon does become more prevalent with evolutionary age. We briefly discuss the mechanisms which might contribute to such an evolutionary effect. A discussion of the limitations of utilizing the 2-CD technique to investigate the role evolutionary age and/or metallicity play in the development of the Be phenomenon is offered, and we provide evidence that other B-type objects of very different nature, such as candidate Herbig Ae/Be stars may contaminate the claimed detections of ``Be stars'' via 2-CDs.

We present the results of long-term spectropolarimetric and spectroscopic monitoring of MWC 314, a candidate Luminous Blue Variable star. We detect the first evidence of H$\alpha$ variability in MWC 314, and find no apparent periodicity in this emission. The total R-band polarization is observed to vary between 2.21% and 3.00% at a position angle consistently around $\sim0^{\circ}$, indicating the presence of a time-variable intrinsic polarization component, hence an asymmetrical circumstellar envelope. We find suggestive evidence that MWC 314's intrinsic polarization exhibits a wavelength-independent magnitude varying between 0.09% and 0.58% at a wavelength-independent position angle covering all four quadrants of the Stokes Q-U plane. Electron scattering off of density clumps in MWC 314's wind is considered as the probable mechanism responsible for these variations.

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 present the results of optical and near-IR spectroscopic and near-IR photometric observations of the emission-line stars Hen3-298 and Hen3-303. Strong emission in the H-alpha line is found in both objects. The presence of Fe II and [O I] emission lines in the spectrum of Hen3-298 indicates that it is a B[e] star. The double-peaked CO line profiles, found in the infrared spectrum of Hen3-298, along with the optical line profiles suggest that the star is surrounded by a rotating circumstellar disk. Both objects also show infrared excesses, similar to those of B[e] stars. The radial velocities of the absorption and emission lines as well as a high reddening level suggest that the objects are located in the Norma spiral arm at a distance of 3-4.5 kpc. We estimated a luminosity of log (L/L_sun) ~ 5.1 and a spectral type of no earlier than B3 for Hen3-298. Hen3-303 seems to be a less luminous B-type object (log (L/L_sun) ~ 4.3), located in the same spiral arm.

We discuss three different observational diagnostics related to disks around hot stars: absorption line determinations of rotational velocities of Be stars; polarization diagnostics of circumstellar disks; and X-ray line diagnostics of one specific magnetized hot star, $\theta^1$ Ori C. Some common themes that emerge from these studies include (a) the benefits of having a specific physical model as a framework for interpreting diagnostic data; (b) the importance of combining several different types of observational diagnostics of the same objects; and (c) that while there is often the need to reinterpret traditional diagnostics in light of new theoretical advances, there are many new and powerful diagnostics that are, or will soon be, available for the study of disks around hot stars.

We present imaging polarimetry observations of the eruptive variable V838 Monocerotis and its neighboring field obtained in 2002 October. The polarization of field stars confirms the previously determined interstellar polarization along the line of sight to V838 Mon. While V838 Mon showed intrinsic polarization shortly after its second outburst on 2002 February 8, all subsequent observations only showed a quiescent interstellar polarization component. We find V838 Mon once again showed significant intrinsic polarization in 2002 October, suggesting the presence of an asymmetrical geometry of scattering material close to the star. Furthermore, an observed 90 degree position angle flip in the intrinsic polarization from 2002 February to 2002 October suggests that the distribution of nearby circumstellar material has experienced significant changes. We discuss the opacity changes in the evolving circumstellar cloud around V838 Mon that may explain these observations.

The spectroscopic and spectropolarimetric variability of the peculiar variable V838 Monocerotis during the brighter phases of its multiple outbursts in 2002 is presented. Significant line profile variability of H$\alpha$ and Si II 6347.10\AA & 6371.36\AA occurred in spectra obtained between 2002 February 5 and 2002 March 14, and a unique secondary absorption component was observed near the end of this time period. Our observations also suggest that multiple shifts in ionization states occurred during the outbursts. Spectropolarimetric observations reveal that V838 Mon exhibited both intrinsic and interstellar polarization components during the initial stages of the second outburst, indicating the presence of an asymmetric geometry; however, the intrinsic component had significantly declined by February 14. We determine the interstellar polarization to be $P_{max} = 2.746 \pm 0.011 %$, $\lambda_{max} = 5790 \pm 37\AA$, $PA = 153.43 \pm 0.12 ^{\circ} $, and we find the integrated intrinsic V band polarization on February 5 to be $P = 0.983 \pm 0.012 %$ at a position angle of $127.0 \pm 0.5^{\circ}$. The implications of these observations for the nature of V838 Monocerotis, its distance, and its ejecta are discussed.

We present the results of recent multicolor photometric and high-resolution spectroscopic observations of the bright Be star Pi Aquarii. Observational data collected from the literature were used to study the star's variations over the last four decades. The star is identified with the IR sources F22227+0107 in the IRAS Faint Point Source catalog and MSX5_G066.0066-44.7392 in the MSX catalog. The variations in near-IR brightness of Pi Aqr are found to be among the largest reported for Be stars. Since 1996, the star has shown only weak signs of circumstellar emission, which has allowed us to refine the fundamental stellar parameters: A_V=0.15 mag., T_eff=24000K, log g=3.9, and M_V=-2.95 mag. A weak emission component of the H-alpha line has been detected during the recent quasi-normal star phase. From analysis of the H-alpha line profiles, we find anti-phased radial velocity variations of the emission component and the photospheric absorption, with a period of 84.1 days and semi-amplitudes of 101.4 and 16.7 km/s, respectively. This result suggests that Pi Aqr may be a binary system consisting of stars with masses of M_1 sin^{3}i = 12.4 M_sun, M_2 sin^{3}i = 2.0 M_sun. We also estimate the orbital inclination angle to be between 50 and 75 degrees. We suggest that the photometric, spectroscopic, and polarimetric variations observed during the second half of the 20th century may be due to variable mass transfer between the binary components.

Ground based imaging, imaging polarimetry, and recent Hubble Space Telescope WFPC2 and NICMOS images of protostars have revealed very complex scattered light patterns that cannot be entirely explained by 2-D radiation transfer models. We present here for the first time radiation transfer models of T Tau and IRAS 04016+2610 that are fully 3-D, with the aim of investigating the effects on image morphology of multiple illuminating sources and infalling envelopes that have been shaped by multiple outflows. For T Tau we have constructed scattered light models where the illumination of the surrounding envelope is by a binary with each source surrounded by its own small circumstellar disk or envelope. We find that the asymmetries in the WFPC2 image of T Tau can be reproduced if the disks in the binary system are misaligned, consistent with a recently discovered bipolar outflow believed to originate from the secondary. For IRAS 04016+2610 we find that the observed scattered light pattern can be reproduced by scattering in an envelope with cavities carved by two sets of bipolar outflows, suggestive of an embedded binary system.

The bright star delta Sco has been considered a typical B0-type object for many years. Spectra of the star published prior to 1990 showed no evidence of emission, but only of short-term line profile variations attributed to nonradial pulsations. Speckle interferometric observations show that delta Sco is a binary system with a highly-eccentric orbit and a period of 10.6 years. Weak emission in the H-alpha line was detected in its spectrum for the first time during a periastron passage in 1990. Shortly before the next periastron passage in the summer of 2000, the binary entered a strong H-alpha emission and enhanced mass loss phase. We monitored the spectroscopic development of the Be outburst from July 2000 through March 2001. In this paper we present results from our spectroscopy, refine elements of the binary orbit, and discuss possible mechanisms for the mass loss.

We report on new data, including four vacuum ultraviolet spectropolarimetric observations by the Wisconsin Ultraviolet Photo-Polarimeter Experiment ("WUPPE") on the Astro-1 (one observation) and Astro-2 (three observations) shuttle missions, and 15 new visible-wavelength observations obtained by the HPOL CCD spectropolarimeter at the Pine Bluff Observatory of the University of Wisconsin. This includes three HPOL observations made within 12 hours of each of the three Astro-2 WUPPE observations, giving essentially simultaneous observations extending from 1500 to 10500 Angstroms. An analysis of these data yields estimates on the properties of wind "clumps" when they are detectable by the polarization of their scattered light. We find that the clumps must be detected near the base of the wind, r/R* = 1.3 - 2.5. At this time the clump density is at least 10^13 cm^-3, 20 times the mean wind, and the temperature is roughly 10,000 K, about 20% cooler than the mean wind. The clumps leading to the largest observed polarization must have electron optical depth of 0.1 - 1 and an angular extent of 0.1 - 1 ster, and account for at most 2% of the wind mass loss. The H and HeI emission lines from the wind are unpolarized, but their P Cygni absorption is enhanced by a factor of 4 in the scattered light. We speculate on the possible relationship of the clumps detected polarimetrically to those seen as Discrete Absorption Features (DACs) and in interferometry.

We have monitored the polarized radiation of the Type IIb SN 1993J in M81 over a period of 41 days, starting from 7 days after the explosion on 1993 March 27.5 (UT). Our data show clear evidence that the intrinsic continuum polarization of SN 1993J evolved from being essentially negligible on April 3-4, to a peak value of ~ 1% in late April 1993, and started to decline by the middle of May. The polarized flux spectrum in late April strongly resembled spectra of Type Ib supernovae, with prominent He I lines but redshifted ~ 3380 km/s relative to the total flux spectrum. These data are consistent with models of Hoflich; they suggest that the polarization was most likely produced by either an asymmetric helium core configuration of material and/or flux, or scattering from an asymmetric circumstellar distribution of dusty material. A combination of electron and dust scattering, as well as a clumpy or stratified distribution of the emitting gas, are possible as the polarization mechanism of the continuum and emission lines. The latter interpretation is supported by the fact that 1-2 months after the explosion, the observed rotations of polarization position angle across prominent line features remain even after correction for effects of interstellar polarization. This indicates that emission lines of He I, Fe II, [O I], and H are all intrinsically polarized at position angles different from that of the continuum, with the non-Balmer lines generally being most highly polarized. If the supernova had an oblate geometry, our data are consistent with a small viewing angle (i.e., more or less equator-on), although the degree of asphericity that gave rise to the polarization at early times is probably smaller (minor to major axis ratio > 0.7) than has been previously suggested.