THE FIRST DETERMINATION OF THE VISCOSITY PARAMETER IN THE CIRCUMSTELLAR DISK OF A Be STAR

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 +/- 0.2, corresponding to a mass injection rate (M) over dot = (3.5 +/- 1.3) x 10(-8) M-circle dot yr(-1). Such a large value of a 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.

New activity in the large circumstellar disk of the Be-shell star 48 Librae

Aims. Spectroscopic, polarimetric, and high spectral resolution interferometric data covering the period 1995-2011 are analyzed to document the transition into a new phase of circumstellar disk activity in the classical Be-shell star 48 Lib. The objective is to use this broad data set to additionally test disk oscillations as the basic underlying dynamical process. Methods. The long-term disk evolution is described using the V/R ratio of the violet and red emission components of H alpha and Br gamma, radial velocities and profiles of He I and optical metal shell lines, as well as multi-band BVRI polarimetry. Single-epoch broad-band and high-resolution interferometric visibilities and phases are discussed with respect to a classical disk model and the given baseline orientations. Results. Spectroscopic signatures of disk asymmetries in 48 Lib vanished in the late nineties but recovered some time between 2004 and 2007, as shown by a new large-amplitude and long-duration V/R cycle. Variations in the radial velocity and line profile of conventional shell lines correlate with the V/R behavior. They are shared by narrow absorption cores superimposed on otherwise seemingly photospheric He I lines, which may form in high-density gas at the inner disk close to the photosphere. Large radial velocity variations continued also during the V/R-quiet years, suggesting that V/R may not always be a good indicator of global density waves in the disk. The comparison of the polarization after the recovery of the V/R activity shows a slight increase, while the polarization angle has been constant for more than 20 years, placing tight limits on any 3-D precession or warping of the disk. The broad H-band interferometry gives a disk diameter of (1.72 +/- 0.2) mas (equivalent to 15 stellar radii), position angle of the disk (50 +/- 9)degrees and a relatively low disk flattening of 1.66 +/- 0.3. Within the errors the same disk position angle is derived from polarimetric observations and from photocenter shifts across Br gamma. The high-resolution interferometric visibility and phase profiles show a double or even multiple-component structure. A preliminary estimate based on the size of the Br gamma emitting region indicates a large diameter for the disk (tens of stellar radii). Overall, no serious contradiction between the observations and the disk-oscillation model could be construed.

Circumstellar emission in Be/X-ray binaries of the Magellanic Clouds and the Milky Way

Aims. We study the optical and near-infrared colour excesses produced by circumstellar emission in a sample of Be/X-ray binaries. Our main goals are exploring whether previously published relations, valid for isolated Be stars, are applicable to Be/X-ray binaries and computing the distance to these systems after correcting for the effects of the circumstellar contamination. Methods. Simultaneous UBVRI photometry and spectra in the 3500−7000 Å spectral range were obtained for 11 optical counterparts to Be/X-ray binaries in the LMC, 5 in the SMC and 12 in the Milky Way. As a measure of the amount of circumstellar emission we used the Hα equivalent width corrected for photospheric absorption. Results. We find a linear relationship between the strength of the Hα emission line and the component of E(B − V) originating from the circumstellar disk. This relationship is valid for stars with emission lines weaker than EW ≈ −15   Å. Beyond this point, the circumstellar contribution to E(B − V) saturates at a value ≈0.17   mag. A similar relationship is found for the (V − I) near infrared colour excess, albeit with a steeper slope and saturation level. The circumstellar excess in (B − V) is found to be about five times higher for Be/X-ray binaries than for isolated Be stars with the same equivalent width EW(Hα), implying significant differences in the physical properties of their circumstellar envelopes. The distance to Be/X-ray binaries (with non-shell Be star companions) can only be correctly estimated by taking into account the excess emission in the V band produced by free-free and free-bound transitions in the circumstellar envelope. We provide a simple method to determine the distances that includes this effect.

New Insights from Aperiodic Variability of Young Stars

Nearly all young stars are variable, with the variability traditionally divided into two classes: periodic variables and aperiodic or "irregular" variables. Periodic variables have been studied extensively, typically using periodograms, while aperiodic variables have received much less attention due to a lack of standard statistical tools. However, aperiodic variability can serve as a powerful probe of young star accretion physics and inner circumstellar disk structure. For my dissertation, I analyzed data from a large-scale, long-term survey of the nearby North America Nebula complex, using Palomar Transient Factory photometric time series collected on a nightly or every few night cadence over several years. This survey is the most thorough exploration of variability in a sample of thousands of young stars over time baselines of days to years, revealing a rich array of lightcurve shapes, amplitudes, and timescales.

I have constrained the timescale distribution of all young variables, periodic and aperiodic, on timescales from less than a day to ~100 days. I have shown that the distribution of timescales for aperiodic variables peaks at a few days, with relatively few (~15%) sources dominated by variability on tens of days or longer. My constraints on aperiodic timescale distributions are based on two new tools, magnitude- vs. time-difference (Δm-Δt) plots and peak-finding plots, for describing aperiodic lightcurves; this thesis provides simulations of their performance and presents recommendations on how to apply them to aperiodic signals in other time series data sets. In addition, I have measured the error introduced into colors or SEDs from combining photometry of variable sources taken at different epochs. These are the first quantitative results to be presented on the distributions in amplitude and time scale for young aperiodic variables, particularly those varying on timescales of weeks to months.

Pre-Main-Sequence Accretion in the Low Metallicity Galactic Star-forming Region Sh 2-284

We present optical spectra of pre-main-sequence (PMS) candidates around the Ha region taken with the Southern African Large Telescope in the low metallicity (Z) Galactic region Sh 2-284, which includes the open cluster Dolidze 25 with an atypical low metallicity of Z similar to 1/5 Z(circle dot). It has been suggested on the basis of both theory and observations that PMS mass-accretion rates, (M) over dot(acc), are a function of Z. We present the first sample of spectroscopic estimates of mass-accretion rates for PMS stars in any low-Z star-forming region. Our data set was enlarged with literature data of H alpha emission in intermediate-resolution R-band spectroscopy. Our total sample includes 24 objects spanning a mass range between 1 and 2 M-circle dot and with a median age of approximately 3.5 Myr. The vast majority (21 out of 24) show evidence for a circumstellar disk on the basis of Two Micron All Sky Survey and Spitzer infrared photometry. We find (M) over dot(acc) in the 1-2 M-circle dot interval to depend quasi-quadratically on stellarmass, with (M) over dot(acc) proportional to M-*(2.4 +/- 0.35), and inversely with stellar age, with (M) over dot(acc) proportional to t(*)(-0.7 +/- 0.4). Furthermore, we compare our spectroscopic (M) over dot(acc) measurements with solar Z Galactic PMS stars in the same mass range, but, surprisingly find no evidence for a systematic change in (M) over dot(acc) with Z. We show that literature accretion-rate studies are influenced by detection limits, and we suggest that (M) over dot(acc) may be controlled by factors other than Z(*), M-*, and age.

Non-LTE Monte Carlo radiative transfer. II. Nonisothermal solutions for viscous Keplerian disks

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 three-dimensional non-LTE Monte Carlo simulations, that viscous Keplerian disks can be highly nonisothermal. In this paper we build on our previous work and solve the full problem of the steady state nonisothermal 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 nonisothermal disk models must be used for a detailed modeling of Be star disks.

TOWARD UNDERSTANDING THE B[e] PHENOMENON. IV. MODELING OF IRAS 00470+6429

FS CMa type stars are a recently described group of objects with the B[e] phenomenon which exhibits 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 (CS) envelope is bimodal, composed of a slowly outflowing disklike 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 satisfactorily reproduce many observational properties of IRAS 00470+6429, including the Hi 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 CS envelope. Possible shortcomings of our model, as well as suggestions for future improvements, are discussed.

DYNAMICAL EVOLUTION OF VISCOUS DISKS AROUND Be STARS. I. PHOTOMETRY

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 buildup 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 three-dimensional non-LTE radiative transfer code. We show that light curves for different wavelengths are specific of a mass loss history, inclination angle, and alpha viscosity parameter. The diagnostic potential of those light curves is also discussed.

Observational effects of stellar rotation on Be stars' disks

Be stars are known to be fast rotators. At high rotation rates a profound modification of the radiation field reaching the circumstellar environment is expected. The origin of this modification is the decrease of the effective gravity on stellar surface leading to the stellar geometrical flattening and the gravity darkening effect predicted by Von Zeipel. Making use of the radiative transfer code HDUST we discuss the consequences of such stellar rotation on the structure of Be star disks based on the Viscous Decretion Disk model. Observational predictions are also made, as SED, IR-excess and Hydrogen line profiles. The modified illumination of the circumstellar disk generates significant changes in these quantities. Ascertaining these changes is useful to set some of the fundamental parameters of the Be system and to unveil the role of stellar rotation over the stellar evolution.

Super-orbital variability of LS I+61 degrees 303 at TeV energies

Context. The gamma-ray binary LS I +61º303 is a well-established source from centimeter radio up to very high energy (VHE; E > 100 GeV). The broadband emission shows a periodicity of ∼26.5 days, coincident with the orbital period. A longer (super-orbital) period of 1667 ± 8 days was proposed from radio variability and confirmed using optical and high-energy (HE; E ¿ 100 MeV) gamma-ray observations. In this paper, we report on a four-year campaign performed by MAGIC together with archival data concentrating on a search for a long-timescale signature in the VHE emission from LS I +61º303. Aims. We focus on the search for super-orbital modulation of the VHE emission, similar to that observed at other energies, and on the search for correlations between TeV emission and an optical determination of the extension of the circumstellar disk. Methods. A four-year campaign has been carried out using the MAGIC telescopes. The source was observed during the orbital phases when the periodic VHE outbursts have occurred (φ = 0.55 – 0.75, one orbit = 26.496 days). Additionally, we included archival MAGIC observations and data published by the VERITAS collaboration in these studies. For the correlation studies, LS I +61◦303 has also been observed during the orbital phases where sporadic VHE emission had been detected in the past (φ = 0.75 – 1.0). These MAGIC observations were simultaneous with optical spectroscopy from the LIVERPOOL telescope. Results. The TeV flux of the periodical outburst in orbital phases φ = 0.5 – 0.75 was found to show yearly variability consistent with the long-term modulation of ∼4.5 years found in the radio band. This modulation of the TeV flux can be well described by a sine function with a best-fit period of 1610±58 days. The complete data, including archival observations, span two super-orbital periods. There is no evidence for a correlation between the TeV emission and the mass-loss rate of the Be star, but this may be affected by the strong, short-timescale (as short as intra-day) variation displayed by the Hα fluxes.

CSI 2264: Simultaneous Optical and Infrared Light Curves of Young Disk-bearing Stars in NGC 2264 with CoRoT and Spitzer—Evidence for Multiple Origins of Variability

We present the Coordinated Synoptic Investigation of NGC 2264, a continuous 30 day multi-wavelength photometric monitoring campaign on more than 1000 young cluster members using 16 telescopes. The unprecedented combination of multi-wavelength, high-precision, high-cadence, and long-duration data opens a new window into the time domain behavior of young stellar objects. Here we provide an overview of the observations, focusing on results from Spitzer and CoRoT. The highlight of this work is detailed analysis of 162 classical T Tauri stars for which we can probe optical and mid-infrared flux variations to 1% amplitudes and sub-hour timescales. We present a morphological variability census and then use metrics of periodicity, stochasticity, and symmetry to statistically separate the light curves into seven distinct classes, which we suggest represent different physical processes and geometric effects. We provide distributions of the characteristic timescales and amplitudes and assess the fractional representation within each class. The largest category (>20%) are optical "dippers" with discrete fading events lasting ~1-5 days. The degree of correlation between the optical and infrared light curves is positive but weak; notably, the independently assigned optical and infrared morphology classes tend to be different for the same object. Assessment of flux variation behavior with respect to (circum)stellar properties reveals correlations of variability parameters with Hα emission and with effective temperature. Overall, our results point to multiple origins of young star variability, including circumstellar obscuration events, hot spots on the star and/or disk, accretion bursts, and rapid structural changes in the inner disk. Based on data from the Spitzer and CoRoT missions. The CoRoT space mission was developed and is operated by the French space agency CNES, with participation of ESA's RSSD and Science Programmes, Austria, Belgium, Brazil, Germany, and Spain.

The relationship between gamma Cassiopeiae's X-ray emission and its circumstellar environment

gamma Cas is the prototypical classical Be star and is recently 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 four XMM-Newton observations. The observational techniques included long baseline optical interferometry (LBOI) from two instruments at CHARA, photometry carried out by an automated photometric telescope and H alpha observations. Because gamma Cas is also known to be in a binary, we measured radial velocities from the H alpha line and redetermined its period as 203.55 +/- 0.20 days and its eccentricity as near zero. The LBOI observations suggest that the star's decretion disk was axisymmetric in 2010, has an system inclination angle near 45 degrees, and a larger radius than previously reported. In addition, the Be star began an "outburst" at the beginning of our campaign, made visible by a brightening and reddening of the disk during our campaign and beyond. Our analyses of the new high resolution spectra disclosed many attributes also found from spectra obtained in 2001 (Chandra) and 2004 (XMM-Newton). As well as a dominant hot (approximate to 14 keV) thermal component, the familiar attributes included: (i) a fluorescent feature of Fe K even 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 two absorption columns are required to fit the continuum. While the first one maintained its historical average of 1 x 10(21) cm(-2), the second was very large and doubled to 7.4 x 10(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.

The Gemini planet-finding campaign: the frequency of giant planets aroud debris disk stars

We have completed a high-contrast direct imaging survey for giant planets around 57 debris disk stars as part of the Gemini NICI Planet-Finding Campaign. We achieved median H-band contrasts of 12.4 mag at 0.''5 and 14.1 mag at 1'' separation. Follow-up observations of the 66 candidates with projected separation <500 AU show that all of them are background objects. To establish statistical constraints on the underlying giant planet population based on our imaging data, we have developed a new Bayesian formalism that incorporates (1) non-detections, (2) single-epoch candidates, (3) astrometric and (4) photometric information, and (5) the possibility of multiple planets per star to constrain the planet population. Our formalism allows us to include in our analysis the previously known β Pictoris and the HR 8799 planets. Our results show at 95% confidence that <13% of debris disk stars have a ≥5 M Jup planet beyond 80 AU, and <21% of debris disk stars have a ≥3 M Jup planet outside of 40 AU, based on hot-start evolutionary models. We model the population of directly imaged planets as d 2 N/dMdavpropm α a β, where m is planet mass and a is orbital semi-major axis (with a maximum value of a max). We find that β < –0.8 and/or α > 1.7. Likewise, we find that β < –0.8 and/or a max < 200 AU. For the case where the planet frequency rises sharply with mass (α > 1.7), this occurs because all the planets detected to date have masses above 5 M Jup, but planets of lower mass could easily have been detected by our search. If we ignore the β Pic and HR 8799 planets (should they belong to a rare and distinct group), we find that <20% of debris disk stars have a ≥3 M Jup planet beyond 10 AU, and β < –0.8 and/or α < –1.5. Likewise, β < –0.8 and/or a max < 125 AU. Our Bayesian constraints are not strong enough to reveal any dependence of the planet frequency on stellar host mass. Studies of transition disks have suggested that about 20% of stars are undergoing planet formation; our non-detections at large separations show that planets with orbital separation >40 AU and planet masses >3 M Jup do not carve the central holes in these disks.

The circumstellar environment and evolutionary status of the supergiant B[e] star Wd1-9

Context. Historically, supergiant (sg)B[e] stars have been difficult to include in theoretical schemes for the evolution of massive OB stars. Aims. The location of Wd1-9 within the coeval starburst cluster Westerlund 1 means that it may be placed into a proper evolutionary context and we therefore aim to utilise a comprehensive multiwavelength dataset to determine its physical properties and consequently its relation to other sgB[e] stars and the global population of massive evolved stars within Wd1. Methods. Multi-epoch R- and I-band VLT/UVES and VLT/FORS2 spectra are used to constrain the properties of the circumstellar gas, while an ISO-SWS spectrum covering 2.45−45μm is used to investigate the distribution, geometry and composition of the dust via a semi-analytic irradiated disk model. Radio emission enables a long term mass-loss history to be determined, while X-ray observations reveal the physical nature of high energy processes within the system. Results. Wd1-9 exhibits the rich optical emission line spectrum that is characteristic of sgB[e] stars. Likewise its mid-IR spectrum resembles those of the LMC sgB[e] stars R66 and 126, revealing the presence of equatorially concentrated silicate dust, with a mass of ~10−4M⊙. Extreme historical and ongoing mass loss (≳ 10−4M⊙yr−1) is inferred from the radio observations. The X-ray properties of Wd1-9 imply the presence of high temperature plasma within the system and are directly comparable to a number of confirmed short-period colliding wind binaries within Wd1. Conclusions. The most complete explanation for the observational properties of Wd1-9 is that it is a massive interacting binary currently undergoing, or recently exited from, rapid Roche-lobe overflow, supporting the hypothesis that binarity mediates the formation of (a subset of) sgB[e] stars. The mass loss rate of Wd1-9 is consistent with such an assertion, while viable progenitor and descendent systems are present within Wd1 and comparable sgB[e] binaries have been identified in the Galaxy. Moreover, the rarity of sgB[e] stars - only two examples are identified from a census of ~ 68 young massive Galactic clusters and associations containing ~ 600 post-Main Sequence stars - is explicable given the rapidity (~ 104yr) expected for this phase of massive binary evolution.

Note on porous rotating disk flow

We revisit the classical Karman rotating disk problem. A series analysis is used to derive estimates of boundary conditions at the surface. Using these estimates, computed thermal and flow fields for large mass transfer through the disk are readily obtained using a shooting method. The relevance of the problem to practical flows is discussed briefly.