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.

PolStar is an Explorer-class far ultraviolet (FUV) spectropolarimetry mission designed to target massive stars and their environments. PolStar will take advantage of resonance lines only available in the FUV to measure for the first time the magnetic and wind environment around massive stars to constrain models of rotation and mass loss.

We use a 3D Monte Carlo radiative transfer code to study the polarimetric and photometric variability from stationary corotating interaction regions (CIR) in the wind of massive stars. Our CIRs are approximated by Archimedean spirals of higher (or lower) density formed in a spherical wind originating from the star and we also made allowance for a bright gaussian spot at the base of the CIR. Comparing results from our code to previous analytical calculations in the optically thin case, we find differences which we attribute mainly to a better estimation of the total unpolarized flux reaching the observer. In the optically thick case, the differences with the analytical calculations are much larger, as multiple scattering introduces extra complexities including occultation effects. The addition of a gaussian spot does not alter the shape of the polarization curve significantly but does create a small excess in polarization. On the other hand, the effect can be larger on the light curve and can become dominant over the resulting CIR, depending on the spot parameters and density of the wind.

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.

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.

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.

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)

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

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.

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 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.

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.

We have updated our publicly available dust radiative transfer code (HOCHUNK3D) to include new emission processes and various 3-D geometries appropriate for forming stars. The 3-D geometries include warps and spirals in disks, accretion hotspots on the central star, fractal clumping density enhancements, and misaligned inner disks. Additional axisymmetric (2-D) features include gaps in disks and envelopes, "puffed-up inner rims" in disks, multiple bipolar cavity walls, and iteration of disk vertical structure assuming hydrostatic equilibrium. We include the option for simple power-law envelope geometry, which combined with fractal clumping, and bipolar cavities, can be used to model evolved stars as well as protostars. We include non-thermal emission from PAHs and very small grains, and external illumination from the interstellar radiation field. The grid structure was modified to allow multiple dust species in each cell; based on this, a simple prescription is implemented to model dust stratification. We describe these features in detail, and show example calculations of each. Some of the more interesting results include the following: 1) Outflow cavities may be more clumpy than infalling envelopes. 2) PAH emission in high-mass stars may be a better indicator of evolutionary stage than the broadband SED slope; and related to this, 3) externally illuminated clumps and high-mass stars in optically thin clouds can masquerade as YSOs. 4) Our hydrostatic equilibrium models suggest that dust settling is likely ubiquitous in T Tauri disks, in agreement with previous observations.

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 report the detection of a unique CO2 ice band toward the deeply embedded, low-mass protostar HOPS-68. Our spectrum, obtained with the Infrared Spectrograph onboard the Spitzer Space Telescope, reveals a 15.2 micron CO2 ice bending mode profile that cannot modeled with the same ice structure typically found toward other protostars. We develop a modified CO2 ice profile decomposition, including the addition of new high-quality laboratory spectra of pure, crystalline CO2 ice. Using this model, we find that 87-92% of the CO2 is sequestered as spherical, CO2-rich mantles, while typical interstellar ices show evidence of irregularly-shaped, hydrogen-rich mantles. We propose that (1) the nearly complete absence of unprocessed ices along the line-of-sight is due to the flattened envelope structure of HOPS-68, which lacks cold absorbing material in its outer envelope, and possesses an extreme concentration of material within its inner (10 AU) envelope region and (2) an energetic event led to the evaporation of inner envelope ices, followed by cooling and re-condensation, explaining the sequestration of spherical, CO2 ice mantles in a hydrogen-poor mixture. The mechanism responsible for the sublimation could be either a transient accretion event or shocks in the interaction region between the protostellar outflow and envelope. The proposed scenario is consistent with the rarity of the observed CO2 ice profile, the formation of nearly pure CO2 ice, and the production of spherical ice mantles. HOPS-68 may therefore provide a unique window into the protostellar feedback process, as outflows and heating shape the physical and chemical structure of protostellar envelopes and molecular clouds.

In this paper we analyse a set of 33 optical spectra, which were acquired with the ARCES echelle spectrograph (R = 38,000) on the 3.5-m telescope at the Apache Point Observatory. We examine the H{\alpha} profile in each of these observations in order to determine the geometry and velocity structure of the previously discovered bipolar jet, which originates from the secondary star of HD 44179 located at the centre of the Red Rectangle nebula. Using a 3D geometric model we are able to determine the orbital coverage during which the jet occults the primary star. During the occultation, part of the H{\alpha} line profile appears in absorption. The velocity structure of the jet was determined by modelling the absorption line profile using the Sobolev approximation for each orbital phase during which we have observations. The results indicate the presence of a wide angle jet, likely responsible for observed biconical structure of the outer nebula. Furthermore, we were able to determine a likely velocity structure and rule out several others. We find that the jet is comprised of low-density, high-velocity, central region and a higher-density, lower-velocity, conical shell.

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.

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.

In this paper we examine the profiles of the complex Na I D-lines in the Red Rectangle. The spectra were acquired with the ARCES \'{e}chelle spectrograph R = 38,000 on the 3.5-m telescope at the Apache Point Observatory. Additional spectra taken with STIS were acquired from the Hubble Legacy Archive (HLA) and were used to independently confirm the spatial origin of the spectral features. The profile of a single D-line consists of double-peaked emission, red-shifted absorption and blue-shifted absorption. We find that the double-peaked emission originates from the bipolar outflow, the red-shifted absorption feature is due to the photospheric line, and the blue-shifted absorption arises from the bipolar outflow as seen against the photosphere of the luminous post-AGB component in HD 44179. In order to better understand the Na I D-line profile, we examined the periodically variable asymmetric photospheric absorption lines. The asymmetric lines are interpreted as a signature of slow self-accretion following enhanced mass-loss around periastron. An empirical model was constructed to remove the photospheric component from the Na I D-line profile in order to study the nebular emission and absorption of sodium along the line-of-sight to the primary. This paper also discusses the different origins of the single-peaked emission, the double-peaked emission and the blue-shifted and red-shifted absorption components.

We present the Spitzer Space Telescope Infrared Spectrograph spectrum of the Orion A protostar HOPS-68. The mid-infrared spectrum reveals crystalline substructure at 11.1, 16.1, 18.8, 23.6, 27.9, and 33.6 microns superimposed on the broad 9.7 and 18 micron amorphous silicate features; the substructure is well matched by the presence of the olivine end-member forsterite. Crystalline silicates are often observed as infrared emission features around the circumstellar disks of Herbig Ae/Be stars and T Tauri stars. However, this is the first unambiguous detection of crystalline silicate absorption in a cold, infalling, protostellar envelope. We estimate the crystalline mass fraction along the line-of-sight by first assuming that the crystalline silicates are located in a cold absorbing screen and secondly by utilizing radiative transfer models. The resulting crystalline mass fractions of 0.14 and 0.17, respectively, are significantly greater than the upper limit found in the interstellar medium (< 0.02-0.05). We propose that the amorphous silicates were annealed within the hot inner disk and/or envelope regions and subsequently transported outward into the envelope by entrainment in a protostellar outflow

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.

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.

(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.

Disk fragmentation resulting from the gravitational instability has been proposed as an efficient mechanism for forming giant planets. We use the planet Fomalhaut b, the triple-planetary system HR 8799, and the potential protoplanet associated with HL Tau to test the viability of this mechanism. We choose the above systems since they harbor planets with masses and orbital characteristics favored by the fragmentation mechanism. We do not claim that these planets must have formed as the result of fragmentation, rather the reverse: if planets can form from disk fragmentation, then these systems are consistent with what we should expect to see. We use the orbital characteristics of these recently discovered planets, along with a new technique to more accurately determine the disk cooling times, to place both lower and upper limits on the disk surface density--and thus mass--required to form these objects by disk fragmentation. Our cooling times are over an order of magnitude shorter than those of Rafikov (2005),which makes disk fragmentation more feasible for these objects. We find that the required mass interior to the planet's orbital radius is ~0.1 Msun for Fomalhaut b, the protoplanet orbiting HL Tau, and the outermost planet of HR 8799. The two inner planets of HR 8799 probably could not have formed in situ by disk fragmentation.

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.

Microlensing surveys have proven to be tremendously fruitful in providing valuable data products for many fields of astrophysics, from eclipse lightcurves for substellar candidates to limb darkening in stellar atmospheres. We report on a program of modeling observables from microlensing of circumstellar envelopes, particularly those of red giant stars that are the most likely to show finite source effects. We will summarize work for how polarization light curves can be used to infer envelope properties and will describe recent modeling of the time dependent spectral energy distributions (SEDs) for microlensing of dusty winds. One of the most exciting developments is the possibility of measuring variable polarization from microlensing in a suitable source using the RINGO polarimeter at La Palma. Also quite interesting is the possibility of probing a dusty wind using IRAC data for a suitable source in the event that Spitzer has a ``warm'' cycle.

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 medium resolution spectropolarimetry and long term photo-polarimetry of two massive post-red supergiants, IRC +10420 and HD 179821. The data provide new information on their circumstellar material as well as their evolution. In IRC +10420, the polarization of the Halpha line is different to that of the continuum, which indicates that the electron-scattering region is not spherically symmetric. The observed long term changes in the polarimetry can be associated with an axi-symmetric structure, along the short axis of the extended reflection nebulosity. Long term photometry reveals that the star increased in temperature until the mid-nineties, after which the photospheric flux in the optical levelled off. As the photometric changes are mostly probed in the red, they do not trace high stellar temperatures sensitively. And so, it is not obvious whether the star has halted its increase in temperature or not. For HD 179821 we find no polarization effects across any absorption or emission lines, but observe very large polarization changes of order 5% over 15 years. Unexpectedly, during the same period, the optical photometry displayed modest variability at the 0.2 magnitude level. Several explanations for this puzzling fact are discussed. Most of which, involving asymmetries in the circumstellar material, seem to fail as there is no evidence for the presence of hot, dusty material close to the star. Alternatively, the variations can be explained by the presence of a non-radially pulsating photosphere. Changes in the photometry hint at an increase in temperature corresponding to a change through two spectral subclasses over the past ten years.

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.

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.

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 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 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.

The high gradient of magnification across the source during both small impact parameter events and caustic crossings offers a unique opportunity for determining the surface brightness profile of the source. Furthermore, models indicate that these events can also provide an appreciable polarisation signal -- arising from differential magnification across the otherwise symmetric source. In a previous paper, polarimetric variations from point lensing of a circumstellar envelope were considered, as would be suitable for an extended envelope around a red giant. In this work we examine the polarisation in the context of caustic crossing events, the scenario which represents the most easily accessible situation for actually observing a polarisation signal in Galactic microlensing.

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).

2-D and 3-D radiation transfer models of forming stars generally produce bluer 1-10 micron colors than 1-D models of the same evolutionary state and envelope mass. Therefore, 1-D models of the shortwave radiation will generally estimate a lower envelope mass and later evolutionary state than multidimensional models. 1-D models are probably reasonable for very young sources, or longwave analysis (wavelengths > 100 microns). In our 3-D models of high-mass stars in clumpy molecular clouds, we find no correlation between the depth of the 10 micron silicate feature and the longwave (> 100 micron) SED (which sets the envelope mass), even when the average optical extinction of the envelope is >100 magnitudes. This is in agreement with the observations of Faison et al. (1998) of several UltraCompact HII (UCHII) regions, suggesting that many of these sources are more evolved than embedded protostars. We have calculated a large grid of 2-D models and find substantial overlap between different evolutionary states in the mid-IR color-color diagrams. We have developed a model fitter to work in conjunction with the grid to analyze large datasets. This grid and fitter will be expanded and tested in 2005 and released to the public in 2006.

We present Monte Carlo scattered light models of a warped disk that reproduce the observed photopolarimetric variability of the classical T Tauri star AA Tau. For a system inclination of 75 deg and using an analytic description for a warped inner disk, we find that the shape and amplitude of the photopolarimetric variability are reproduced with a warp which occults the star, located at 0.07 au, amplitude 0.016 au, extending over radial and azimuthal ranges 0.0084 au and 145 deg. We also show a time sequence of high spatial resolution scattered light images, showing a dark shadow cast by the warp sweeping round the disk. Using a modified smooth particle hydrodynamics code, we find that a stellar dipole magnetic field of strength 5.2 kG, inclined at 30 deg to the stellar rotation axis can reproduce the required disk warping to explain AA Tau's photopolarimetric variability.

We present observations of four T Tauri stars using long baseline infrared interferometry from the Palomar Testbed Interferometer. The target sources, T Tau N, SU Aur, RY Tau and DR Tau, are all known to be surrounded by dusty circumstellar disks. The observations directly trace the inner regions (< 1 AU) of the disk and can be used to constrain the physical properties of this material. For three of the sources observed, the infrared emission is clearly resolved. We first use geometric models to characterize the emission region size, which ranges from 0.04 to 0.3 AU in radius. We then use Monte Carlo radiation transfer models of accretion disks to jointly model the spectral energy distribution and the interferometric observations with disk models including accretion and scattering. With these models, we are able to reproduce the data set with extended emission arising from structures larger than 10 milliarcseconds contributing less than 6% of the K band emission, consistent with there being little or no envelope remaining for these Class II sources (d log(lambda*F_lambda)/d log(lambda) ~ -2 to 0 in the infrared). The radiation transfer models have inner radii for the dust similar to the geometric models; however, for RY Tau emission from gas within the inner dust radius contributes significantly to the model flux and visibility at infrared wavelengths. The main conclusion of our modeling is that emission from inner gas disks (between the magnetic truncation radius and the dust destruction radius) can be a significant component in the inner disk flux for sources with large inner dust radii.

We model how the mid-infrared colors of Young Stellar Objects (YSOs) vary with stellar temperature. The spectral energy distribution (SED) of each object has contributions from thermal emission of circumstellar dust, from direct stellar photospheric emission, and from scattered stellar emission. We first isolate the effects of stellar contributions (direct+scattered) to the SED using homologous "Class I" models, where the shape of the thermal contribution to the SED remains constant. Stellar light contributes relatively more to the MIR emission of YSOs with lower temperature stars (Tstar ~4000 K) because the emission peak wavelength of the star is closer to that of the thermal radiation. Next we consider realistic Class 0, I, and II source models and find that the other dominant effect of varying stellar temperature on YSO SEDs is that of the inner disk wall: In high-Tstar models, the dust destruction radius is much further out with a consequently larger inner disk wall that contributes relatively more to the 2-10 um flux. This effect partially offsets that of the stellar contribution leading to varying behaviors of the 2-10 um flux. Our modeled MIR protostellar colors have implications for interpretations of Spitzer IRAC observations of star formation regions: It is commonly assumed that the slope of the SED at 1-10 um is directly related to evolutionary state. We show that inclination effects, aperture size, scattered light, and stellar temperature cause a broad spread in the colors of a source at a single evolutionary state. Color magnitude diagrams can help sort out these effects by separating sources with different Tstar based on their different brightness (for sources at the same distance).

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 present synthetic spectra for circumstellar disks that are heated by radiation from a central brown dwarf. Under the assumption of vertical hydrostatic equilibrium, our models yield scaleheights for brown dwarf disks in excess of three times those derived for classical T Tauri (CTTS) disks. If the near-IR excess emission observed from brown dwarfs is indeed due to circumstellar disks, then the large scaleheights we find could have a significant impact on the optical and near-IR detectability of such systems. Our radiation transfer calculations show that such highly flared disks around brown dwarfs will result in a large fraction of obscured sources due to extinction of direct starlight by the disk over a wide range of sightlines. The obscured fraction for a 'typical' CTTS is less than 20%. We show that the obscured fraction for brown dwarfs may be double that for CTTS, but this depends on stellar and disk mass. We also comment on possible confusion in identifying brown dwarfs via color-magnitude diagrams: edge-on CTTS display similar colors and magnitudes as a face-on brown dwarf plus disk systems.

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 present model spectral energy distributions, colors, polarization, and images for an evolutionary sequence of a low-mass protostar from the early collapse stage (Class 0) to the remnant disk stage (Class III). We find a substantial overlap in colors and SEDs between protostars embedded in envelopes (Class 0-I) and T Tauri disks (Class II), especially at mid-IR wavelengths. Edge-on Class I-II sources show double-peaked spectral energy distributions, with a short-wavelength hump due to scattered light and the long-wavelength hump due to thermal emission. These are the bluest sources in mid-IR color-color diagrams. Since Class 0 and I sources are diffuse, the size of the aperture over which fluxes are integrated has a substantial effect on the computed colors, with larger aperture results showing significantly bluer colors. This causes overlap in color-color diagrams between all evolutionary states, especially in the mid-IR. However the near-IR polarization of the Class 0 sources is much higher than the Class I-II sources, providing a means to separate these evolutionary states. We varied the grain properties in the circumstellar envelope, allowing for larger grains in the disk midplane and smaller in the envelope. We find that grain growth in disks of Class I sources can be detected at wavelengths greater than 100 $\mu$m. Our image calculations predict that the diffuse emission from edge-on Class I and II sources should be detectable in the mid-IR with the Space Infrared Telescope Facility (SIRTF) in nearby star forming regions (out to several hundred parsecs).

We present 2-D radiation transfer models of Class I Protostars and show the effect of including more realistic geometries on the resulting spectral energy distributions and images. We begin with a rotationally flattened infalling envelope as our comparison model, and add a flared disk and bipolar cavity. The disk affects the spectral energy distribution most strongly at edge-on inclinations, causing a broad dip at about 10 um (independent of the silicate feature) due to high extinction and low scattering albedo in this wavelength region. The bipolar cavities allow more direct stellar+disk radiation to emerge into polar directions, and more scattering radiation to emerge into all directions. The wavelength-integrated flux, often interpreted as luminosity, varies with viewing angle, with pole-on viewing angles seeing 2-4 times as much flux as edge-on, depending on geometry. Thus, observational estimates of luminosity should take into account the inclination of a source. The envelopes with cavities are significantly bluer in near-IR and mid-IR color-color plots than those without cavities. Using 1-D models to interpret Class I sources with bipolar cavities would lead to an underestimate of envelope mass and an overestimate of the implied evolutionary state. We compute images at near-, mid-, and far-IR wavelengths. We find that the mid-IR colors and images are sensitive to scattering albedo, and that the flared disk shadows the midplane on large size scales at all wavelengths plotted. Finally, our models produce polarization spectra which can be used to diagnose dust properties, such as albedo variations due to grain growth. Our results of polarization across the 3.1 um ice feature agree well with observations for ice mantles covering 5% of the radius of the grains.

[Abridged] We present ISO far-infrared (IR) observations of 21 hard X-ray selected AGN from the HEAO-1 A2 sample. We compare the far-IR to X-ray spectral energy distributions (SEDs) of this sample with various radio and optically selected AGN samples. The hard-X-ray selected sample shows a wider range of optical/UV shapes extending to redder near-IR colors. The bluer objects are Seyfert 1s, while the redder AGN are mostly intermediate or type 2 Seyferts. This is consistent with a modified unification model in which the amount of obscuring material increases with viewing angle and may be clumpy. Such a scenario, already suggested by differing optical/near-IR spectroscopic and X-ray AGN classifications, allows for different amounts of obscuration of the continuum emission in different wavebands and of the broad emission line region which results in a mixture of behaviors for AGN with similar optical emission line classifications. The resulting limits on the column density of obscuring material through which we are viewing the redder AGN are 100 times lower than for the standard optically thick torus models. The resulting decrease in optical depth of the obscuring material allows the AGN to heat more dust at larger radial distances. We show that an AGN-heated, flared, dusty disk with mass 10^9 solar and size of few hundred pc is able to generate optical-far-IR SEDs which reproduce the wide range of SEDs present in our sample with no need for an additional starburst component to generate the long-wavelength, cooler part of the IR continuum.

The unusual spectral energy distribution (SED) of the classical T Tauri star GM Aurigae provides evidence for the presence of an inner disc hole extending to several au. Using a combination of hydrodynamical simulations and Monte Carlo radiative transport, we investigate whether the observed SED is consistent with the inner hole being created and maintained by an orbiting planet. We show that an ~ 2 M_Jupiter planet, orbiting at 2.5 au in a disc with mass 0.047 M_sun and radius 300 au, provides a good match both to the SED and to CO observations which constrain the velocity field in the disc. A range of planet masses is allowed by current data, but could in principle be distinguished with further observations between 3 and ~ 20 microns. Future high precision astrometric instruments should also be able to detect the motion of the central star due to an orbiting Jupiter mass planet. We argue that the small number of T Tauri stars with SEDs resembling that of GM Aur is broadly consistent with the expected statistics of embedded migrating planets.

We report the discovery in NIR with SofI at the NTT of a resolved circumstellar dust disk around a 2MASS source at the periphery of the rho Ophiuchi dark cloud. We present follow-up observations in J, H, and Ks-band obtained with ISAAC at the VLT, under 0.4"-seeing conditions, which unveil a dark dust lane oriented East-West between two characteristic northern and southern reflection nebulae. This new circumstellar dust disk has a radius of 2.15" (300 AU at 140 pc), and a width of 1.2" (170 AU at 140 pc). Thanks to its location at the periphery of the dense cores, it suffers small foreground visual extinction (Av=2.1\pm2.6 mag). Although this disk is seen close to edge-on, the two reflection nebulae display very different colors. We introduce a new NIR data visualization called ``Pixel NIR Color Mapping'', which allows to visualize directly the NIR colors of the nebula pixels. Thanks to this method we identify a ridge, 0.3" (40 AU at 140 pc) to the north of the dark lane and parallel to it, which displays a NIR color excess. This ridge corresponds to an unusual increase of brightness from J to Ks, which is also visible in the NTT observation obtained 130 days before the VLT one. We also find that the northern nebula shows ~3 mag more extinction than the southern nebula. We compute axisymmetric disk models to reproduce the VLT scattered light images and the spectral energy distribution from optical to NIR. Our best model, with a disk inclination i=86\pm1 deg, correctly reproduces the extension of the southern reflection nebula, but it is not able to reproduce either the observed NIR color excess in the northern nebula or the extinction difference between the two reflection nebulae. We discuss the possible origin of the peculiar asymmetrical NIR color properties of this object.

We investigate the observational signatures of a straightforward evolutionary scenario for protoplanetary disks in which the disk mass of small (50 micron) particles decreases homologously with time, but the disk structure and stellar parameters do not change. Our goal is to identify optimal infrared spectral indicators of the existence of disks, their structure, and mass evolution that may be tested with the upcoming SIRTF mission. We present simulated spectral energy distributions and colors over a wide range of masses. The SED is most sensitive to disk mass in the far-IR and longer wavelengths, which is already known from millimeter and radio observations. As the disk mass decreases, the excess emission of the disk over the stellar photosphere diminishes more rapidly at the longest rather than at short wavelengths. At near-infrared wavelengths, the disk remains optically thick to stellar radiation over a wide range of disk mass, resulting in a slower decline in the SED in this spectral regime. Therefore, near-IR excesses (K-L) provide a robust means of detecting disks in star clusters down to 1E-7 solar masses, while the far-IR excess probes the disk mass. Reducing the disk mass results in a clear progression in color-color diagrams with low mass disks displaying the bluest colors. We interpret color-color diagrams for Taurus-Auriga sources in the context of decreasing disk mass.

We present spectral energy distribution (SED) models for the edge-on classical T Tauri star HH30 IRS that indicate dust grains have grown to larger than 50 microns within its circumstellar disk. The disk geometry and inclination are known from previous modeling of multiwavelength Hubble Space Telescope images and we use the SED to constrain the dust size distribution. Model spectra are shown for different circumstellar dust models: a standard ISM mixture and larger grain models. As compared to ISM grains, the larger dust grain models have a shallower wavelength dependent opacity. Models with the larger dust grains provide a good match to the currently available data, but mid and far-IR observations are required to more tightly constrain the dust size distribution. The accretion luminosity in our models is L_acc<0.2 L_star corresponding to an accretion rate of 4E-9M_sun/yr. Dust size distributions that are simple power-law extensions (i.e., no exponential cutoff) yield acceptable fits to the optical/near-IR but too much emission at mm wavelengths and require larger disk masses. Such a simple size distribution would not be expected in an environment such as the disk of HH30 IRS, particularly over such a large range in grain sizes. However, its ability to adequately characterize the grain populations may be determined from more complete observational sampling of the SED in the mid to far-IR.

We describe a general radiative equilibrium and temperature correction procedure for use in Monte Carlo radiation transfer codes with sources of temperature-independent opacity, such as astrophysical dust. The technique utilizes the fact that Monte Carlo simulations track individual photon packets, so we may easily determine where their energy is absorbed. When a packet is absorbed, it heats a particular cell within the envelope, raising its temperature. To enforce radiative equilibrium, the absorbed packet is immediately re-emitted. To correct the cell temperature, the frequency of the re-emitted packet is chosen so that it corrects the temperature of the spectrum previously emitted by the cell. The re-emitted packet then continues being scattered, absorbed, and re-emitted until it finally escapes from the envelope. As the simulation runs, the envelope heats up, and the emergent spectral energy distribution (SED) relaxes to its equilibrium value, without iteration. This implies that the equilibrium temperature calculation requires no more computation time than the SED calculation of an equivalent pure scattering model with fixed temperature. In addition to avoiding iteration, our method conserves energy exactly, because all injected photon packets eventually escape. Furthermore, individual packets transport energy across the entire system because they are never destroyed. This long-range communication, coupled with the lack of iteration, implies that our method does not suffer the convergence problems commonly associated with lambda-iteration. To verify our temperature correction procedure, we compare our results to standard benchmark tests, and finally we present the results of simulations for two-dimensional axisymmetric density structures.

Observations of the B0.2V star tau Scorpii have revealed unusual stellar wind characteristics: red-shifted absorption in the far-ultraviolet O VI resonance doublet up to +250 km/s, and extremely hard X-ray emission implying gas at temperatures in excess of 10^7 K. We describe a phenomenological model to explain these properties. We assume the wind of tau Sco consists of two components: ambient gas in which denser clumps are embedded. The clumps are optically thick in the UV resonance lines primarily responsible for accelerating the ambient wind. The reduced acceleration causes the clumps to slow and even infall, all the while being confined by the ram pressure of the outflowing ambient wind. We calculate detailed trajectories of the clumps in the ambient stellar wind, accounting for a line radiation driving force and the momentum deposited by the ambient wind in the form of drag. We show these clumps will fall back towards the star with velocities of several hundred km/sec for a broad range of initial conditions. The infalling clumps produce X-ray emitting plasmas with temperatures in excess of (1-6)x10^7 K in bow shocks at their leading edge. The infalling material explains the peculiar red-shifted absorption wings seen in the O VI doublet. The required mass loss in clumps is 3% - 30% ofthe total mass loss rate. The model developed here can be generally applied to line-driven outflows with clumps or density irregularities. (Abstract Abridged)

In this paper we provide a possible link between the structure of the bipolar nebula surrounding SN1987A and the properties of its progenitor star. A Wind Blwon Bubble (WBB) scenario is emplyed, in which a fast, tenuous wind from a Blue Supergiant expands into a slow, dense wind, expelled during an earlier Red Supergiant phase. The bipolar shapre develops due to a pole-to-equator density contrast in the slow wind (ie, the slow wind forms a slow torus). We use the Wind Compressed Disk (WCD) model of Bjorkman & Cassinelli (1992) to determine the shape of the slow torus. In the WCD scenario, the shape of the torus is determined by the rotation of the progenitor star. We then use a self-similar semi-analytical method for wind blown bubble evolution to determine the shape of the resulting bipolar nebula. We find that the union of the wind-compressed-disk and bipolar-wind-blown- bubble models allows us to recover the salient properties of SN1987A's circumstellar nebula. In particular, the size, speed and density of SN1987A's inner ring are easily reproduced in our calculations. An exploration of parameter space shows the the red supergiant progenitor must be been rotating at > 0.3 of its breakup speed. We conclude that the progenitor was most likely spun up by a merger with a binary companion. Using a simple model for the binary merger we find that the companion is likely to have had a mass > 0.5 M_sun.