Theoretical models of the ionization state in protoplanetary disks suggest the existence of large areas with low ionization and weak coupling between the gas and magnetic fields. In this regime hydrodynamical instabilities may become important. In this work we investigate the gas and dust structure and dynamics for a typical T Tauri system under the influence of the vertical shear instability (VSI). We use global 3D radiation hydrodynamics simulations covering all $360^\circ$ of azimuth with embedded particles of 0.1 and 1mm size, evolved for 400 orbits. Stellar irradiation heating is included with opacities for 0.1- to 10-$\mu$m-sized dust. Saturated VSI turbulence produces a stress-to-pressure ratio of $\alpha \simeq 10^{-4}$. The value of $\alpha$ is lowest within 30~au of the star, where thermal relaxation is slower relative to the orbital period and approaches the rate below which VSI is cut off. The rise in $\alpha$ from 20 to 30~au causes a dip in the surface density near 35~au, leading to Rossby wave instability and the generation of a stationary, long-lived vortex spanning about 4~au in radius and 40~au in azimuth. Our results confirm previous findings that mm size grains are strongly vertically mixed by the VSI. The scale height aspect ratio for 1mm grains is determined to be 0.037, much higher than the value $H/r=0.007$ obtained from millimeter-wave observations of the HL~Tau system. The measured aspect ratio is better fit by non-ideal MHD models. In our VSI turbulence model, the mm grains drift radially inwards and many are trapped and concentrated inside the vortex. The turbulence induces a velocity dispersion of $\sim 12$~m/s for the mm grains, indicating that grain-grain collisions could lead to fragmentation.

Outflows driven by large-scale magnetic fields likely play an important role in the evolution and dispersal of protoplanetary disks, and in setting the conditions for planet formation. We extend our 2-D axisymmetric non-ideal MHD model of these outflows by incorporating radiative transfer and simplified thermochemistry, with the twin aims of exploring how heating influences wind launching, and illustrating how such models can be tested through observations of diagnostic spectral lines. Our model disks launch magnetocentrifugal outflows primarily through magnetic tension forces, so the mass-loss rate increases only moderately when thermochemical effects are switched on. For typical field strengths, thermochemical and irradiation heating are more important than magnetic dissipation. We furthermore find that the entrained vertical magnetic flux diffuses out of the disk on secular timescales as a result of non-ideal MHD. Through post-processing line radiative transfer, we demonstrate that spectral line intensities and moment-1 maps of atomic oxygen, the HCN molecule, and other species show potentially observable differences between a model with a magnetically driven outflow and one with a weaker, photoevaporative outflow. In particular, the line shapes and velocity asymmetries in the moment-1 maps could enable the identification of outflows emanating from the disk surface.

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

The Habitable Exoplanet Observatory, or HabEx, has been designed to be the Great Observatory of the 2030s. For the first time in human history, technologies have matured sufficiently to enable an affordable space-based telescope mission capable of discovering and characterizing Earthlike planets orbiting nearby bright sunlike stars in order to search for signs of habitability and biosignatures. Such a mission can also be equipped with instrumentation that will enable broad and exciting general astrophysics and planetary science not possible from current or planned facilities. HabEx is a space telescope with unique imaging and multi-object spectroscopic capabilities at wavelengths ranging from ultraviolet (UV) to near-IR. These capabilities allow for a broad suite of compelling science that cuts across the entire NASA astrophysics portfolio. HabEx has three primary science goals: (1) Seek out nearby worlds and explore their habitability; (2) Map out nearby planetary systems and understand the diversity of the worlds they contain; (3) Enable new explorations of astrophysical systems from our own solar system to external galaxies by extending our reach in the UV through near-IR. This Great Observatory science will be selected through a competed GO program, and will account for about 50% of the HabEx primary mission. The preferred HabEx architecture is a 4m, monolithic, off-axis telescope that is diffraction-limited at 0.4 microns and is in an L2 orbit. HabEx employs two starlight suppression systems: a coronagraph and a starshade, each with their own dedicated instrument.

Young planets with masses approaching Jupiter's have tides strong enough to clear gaps around their orbits in the protostellar disk. Gas flow through the gaps regulates the planets' further growth and governs the disks' evolution. Magnetic forces may drive that flow if the gas is sufficiently ionized to couple to the fields. We compute the ionizing effects of the X-rays from the central young star, using Monte Carlo radiative transfer calculations to find the spectrum of Compton-scattered photons reaching the planet's vicinity. The scattered X-rays ionize the gas at rates similar to or greater than the interstellar cosmic ray rate near planets the mass of Saturn and of Jupiter, located at 5 au and at 10 au, in disks with the interstellar mass fraction of sub-micron dust and with the dust depleted a factor 100. Solving a gas-grain recombination reaction network yields charged particle populations whose ability to carry currents is sufficient to partly couple the magnetic fields to the gas around the planet. Most cases can undergo Hall shear instability, and some can launch magnetocentrifugal winds. However the material on the planet's orbit has diffusivities so large in all the cases we examine, that magneto-rotational turbulence is prevented and the non-ideal terms govern the magnetic field's evolution. Thus the flow of gas in the gaps opened by the young giant planets depends crucially on the finite conductivity.

The increasing number of newly detected exoplanets at short orbital periods raises questions about their formation and migration histories. A particular puzzle that requires explanation arises from one of the key results of the Kepler mission, namely the increase in the planetary occurrence rate with orbital period up to 10 days for F, G, K and M stars. We investigate the conditions for planet formation and migration near the dust sublimation front in protostellar disks around young Sun-like stars. For this analysis we use iterative 2D radiation hydrostatic disk models which include irradiation by the star, and dust sublimation and deposition depending on the local temperature and vapor pressure. We perform a parameter study by varying the magnetized turbulence onset temperature, the accretion stress, the dust mass fraction, and the mass accretion rate. Our models feature a gas-only inner disk, a silicate sublimation front and dust rim starting at around 0.08 au, an ionization transition zone with a corresponding density jump, and a pressure maximum which acts as a pebble trap at around 0.12 au. Migration torque maps show Earth- and super-Earth-mass planets halt in our model disks at orbital periods ranging from 10 to 22 days. Such periods are in good agreement with both the inferred location of the innermost planets in multiplanetary systems, and the break in planet occurrence rates from the Kepler sample at 10 days. In particular, models with small grains depleted produce a trap located at a 10-day orbital period, while models with a higher abundance of small grains present a trap at around a 17-day orbital period. The snow line lies at 1.6 au, near where the occurrence rate of the giant planets peaks. We conclude that the dust sublimation zone is crucial for forming close-in planets, especially when considering tightly packed super-Earth systems.

Hydrodynamical simulations show that a pair of spiral arms can form in the disk around a rapidly-growing young star and that the arms are crucial in transporting angular momentum as the disk accretes material from the surrounding envelope. Here we report the detection of a pair of symmetric spiral structures in a protostellar disk, supporting the formation of spiral arms in the disk around a forming star. The HH 111 VLA 1 source is a young Class I source embedded in a massive infalling protostellar envelope and is actively accreting, driving the prominent HH 111 jet. Previous observations showed a ring of shock emission around the disk's outer edge, indicating accretion of the envelope material onto the disk at a high rate. Now with ALMA observations of thermal emission from dust particles, we detect a pair of spiral arms extending from the inner region to the disk's outer edge, similar to that seen in many simulations. Additionally, the disk is massive, with Toomre's Q parameter near unity in the outer parts where the spiral structures are detected, supporting the notion that envelope accretion is driving the outer disk gravitationally unstable. In our observations, another source, HH 111 VLA 2, is spatially resolved for the first time, showing a disk-like structure with a diameter of ~ 26 au and an orientation nearly orthogonal to that of the HH 111 VLA 1 disk.

The recent discovery of a staggering diversity of planets beyond the Solar System has brought with it a greatly expanded search space for habitable worlds. The Kepler exoplanet survey has revealed that most planets in our interstellar neighborhood are larger than Earth and smaller than Neptune. Collectively termed super-Earths and mini-Neptunes, some of these planets may have the conditions to support liquid water oceans, and thus Earth-like biology, despite differing in many ways from our own planet. In addition to their quantitative abundance, super-Earths are relatively large and are thus more easily detected than true Earth twins. As a result, super-Earths represent a uniquely powerful opportunity to discover and explore a panoply of fascinating and potentially habitable planets in 2020 - 2030 and beyond.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

For the first time in human history, technologies have matured sufficiently to enable a mission capable of discovering and characterizing habitable planets like Earth orbiting sunlike stars other than the Sun. At the same time, such a platform would enable unique science not possible from ground-based facilities. This science is broad and exciting, ranging from new investigations of our own solar system to a full range of astrophysics disciplines. The Habitable Exoplanet Observatory, or HabEx, is one of four studies currently being undertaken by NASA in preparation for the 2020 Astrophysics Decadal Survey. HabEx has been designed to be the Great Observatory of the 2030s, with community involvement through a competed and funded Guest Observer (GO) program. This interim report describes the HabEx baseline concept, which is a space-based 4-meter diameter telescope mission concept with ultraviolet (UV), optical, and near-infrared (near-IR) imaging and spectroscopy capabilities. More information on HabEx can be found at https://www.jpl.nasa.gov/habex

In spite of making a small contribution to total protoplanetary disk mass, dust affects the disk temperature by controlling absorption of starlight. As grains grow from their initial ISM-like size distribution, settling depletes the disk's upper layers of dust and decreases the optical depth, cooling the interior. Here we investigate the effect of collisional growth of dust grains and their dynamics on the thermal and optical profile of the disk, and explore the possibility that cooling induced by grain growth and settling could lead to gravitational instability. We develop a Monte Carlo dust collision model with a weighting technique and allow particles to collisionally evolve through sticking and fragmentation, along with vertical settling and turbulent mixing. We explore two disk models, the MMEN (minimum-mass extrasolar nebula), and a "heavy" disk with higher surface density than the MMEN, and perform simulations for both constant and spatially variable turbulence efficiency profile $\alpha(R,z)$. We then calculate mean wavelength-dependent opacities for the evolving disks and perform radiative transfer to calculate the temperature profile $T(R,z)$. Finally, we calculate the Toomre Q parameter, a measure of the disk's stability against self-gravity, for each disk model after it reaches a steady state dust-size distribution. We find that even weak turbulence can keep sub-micron sized particles stirred in the disk's upper layer, affecting its optical and thermal profiles, and the growth of large particles in the midplane can make a massive disk optically thick at millimeter wavelengths, making it difficult to calculate the surface density of dust available for planet formation in the inner disk. Also, for an initially massive disk, grain settling and growth can produce a drop in the Toomre Q parameter, driving the disk to $Q < 1.4$ and possibly triggering spiral instabilities.

The Planet Formation Imager (PFI, www.planetformationimager.org) is a next-generation infrared interferometer array with the primary goal of imaging the active phases of planet formation in nearby star forming regions. PFI will be sensitive to warm dust emission using mid-infrared capabilities made possible by precise fringe tracking in the near-infrared. An L/M band combiner will be especially sensitive to thermal emission from young exoplanets (and their disks) with a high spectral resolution mode to probe the kinematics of CO and H2O gas. In this paper, we give an overview of the main science goals of PFI, define a baseline PFI architecture that can achieve those goals, point at remaining technical challenges, and suggest activities today that will help make the Planet Formation Imager facility a reality.

Many protostellar disks show central cavities, rings, or spiral arms likely caused by low-mass stellar or planetary companions, yet few such features are conclusively tied to bodies embedded in the disks. We note that even small features on the disk's surface cast shadows, because the starlight grazes the surface. We therefore focus on accurately computing the disk's thickness, which depends on its temperature. We present models with temperatures set by the balance between starlight heating and radiative cooling, that are also in vertical hydrostatic equilibrium. The planet has 20, 100, or 1000~M$_\oplus$, ranging from barely enough to perturb the disk significantly, to clearing a deep tidal gap. The hydrostatic balance strikingly alters the model disk's appearance. The planet-carved gap's outer wall puffs up under starlight heating, throwing a shadow across the disk beyond. The shadow appears in scattered light as a dark ring that could be mistaken for a gap opened by another more distant planet. The surface brightness contrast between outer wall and shadow for the 1000-M$_\oplus$ planet is an order of magnitude greater than a model neglecting the temperature disturbances. The shadow is so deep it largely hides the planet-launched spiral wave's outer arm. Temperature gradients are such that outer low-mass planets undergoing orbital migration will converge within the shadow. Furthermore the temperature perturbations affect the shape, size, and contrast of features at millimeter and centimeter wavelengths. Thus radiative heating and cooling are key to the appearance of protostellar disks with embedded planets.

Dozens of habitable zone, approximately earth-sized exoplanets are known today. An emerging frontier of exoplanet studies is identifying which of these habitable zone, small planets are actually habitable (have all necessary conditions for life) and, of those, which are earth-like. Many parameters and processes influence habitability, ranging from the orbit through detailed composition including volatiles and organics, to the presence of geological activity and plate tectonics. While some properties will soon be directly observable, others cannot be probed by remote sensing for the foreseeable future. Thus, statistical understanding of planetary systems' formation and evolution is a key supplement to the direct measurements of planet properties. Probabilistically assessing parameters we cannot directly measure is essential to reliably assess habitability, to prioritizing habitable-zone planets for follow-up, and for interpreting possible biosignatures.

The discovery of thousands of exoplanets over the last couple of decades has shown that the birth of planets is a very efficient process in nature. Theories invoke a multitude of mechanisms to describe the assembly of planets in the disks around pre-main-sequence stars, but observational constraints have been sparse on account of insufficient sensitivity and resolution. Understanding how planets form and interact with their parental disk is crucial also to illuminate the main characteristics of a large portion of the full population of planets that is inaccessible to current and near-future observations. This White Paper describes some of the main issues for our current understanding of the formation and evolution of planets, and the critical contribution expected in this field by the Next Generation Very Large Array.

HH 212 is a Class 0 protostellar system found to host a "hamburger"-shaped dusty disk with a rotating disk atmosphere and a collimated SiO jet at a distance of ~ 400 pc. Recently, a compact rotating outflow has been detected in SO and SO2 toward the center along the jet axis at ~ 52 au (0.13") resolution. Here we resolve the compact outflow into a small-scale wide-opening rotating outflow shell and a collimated jet, with the observations in the same S-bearing molecules at ~ 16 au (0.04") resolution. The collimated jet is aligned with the SiO jet, tracing the shock interactions in the jet. The wide-opening outflow shell is seen extending out from the inner disk around the SiO jet and has a width of ~ 100 au. It is not only expanding away from the center, but also rotating around the jet axis. The specific angular momentum of the outflow shell is ~ 40 au km/s. Simple modeling of the observed kinematics suggests that the rotating outflow shell can trace either a disk wind or disk material pushed away by an unseen wind from the inner disk or protostar. We also resolve the disk atmosphere in the same S-bearing molecules, confirming the Keplerian rotation there.

We report 885$\mu$m ALMA continuum flux densities for 24 Taurus members spanning the stellar/substellar boundary, with spectral types from M4 to M7.75. Of the 24 systems, 22 are detected at levels ranging from 1.0-55.6 mJy. The two non-detections are transition disks, though other transition disks in the sample are detected. Converting ALMA continuum measurements to masses using standard scaling laws and radiative transfer modeling yields dust mass estimates ranging from $\sim$0.3-20M$_{\oplus}$. The dust mass shows a declining trend with central object mass when combined with results from submillimeter surveys of more massive Taurus members. The substellar disks appear as part of a continuous sequence and not a distinct population. Compared to older Upper Sco members with similar masses across the substellar limit, the Taurus disks are brighter and more massive. Both Taurus and Upper Sco populations are consistent with an approximately linear relationship in $M_{dust}$ to $M_{star}$, although derived power-law slopes depend strongly upon choices of stellar evolutionary model and dust temperature relation. The median disk around early M-stars in Taurus contains a comparable amount of mass in small solids as the average amount of heavy elements in Kepler planetary systems on short-period orbits around M-dwarf stars, with an order of magnitude spread in disk dust mass about the median value. Assuming a gas:dust ratio of 100:1, only a small number of low-mass stars and brown dwarfs have a total disk mass amenable to giant planet formation, consistent with the low frequency of giant planets orbiting M-dwarfs.

Turbulence is the leading candidate for angular momentum transport in protoplanetary disks and therefore influences disk lifetimes and planet formation timescales. However, the turbulent properties of protoplanetary disks are poorly constrained observationally. Recent studies have found turbulent speeds smaller than what fully-developed MRI would produce (Flaherty et al. 2015, 2017). However, existing studies assumed a constant CO/H2 ratio of 0.0001 in locations where CO is not frozen-out or photo-dissociated. Our previous studies of evolving disk chemistry indicate that CO is depleted by incorporation into complex organic molecules well inside the freeze-out radius of CO. We consider the effects of this chemical depletion on measurements of turbulence. Simon et al. (2015) suggested that the ratio of the peak line flux to the flux at line center of the CO J=3-2 transition is a reasonable diagnostic of turbulence, so we focus on that metric, while adding some analysis of the more complex effects on spatial distribution. We simulate the emission lines of CO based on chemical evolution models presented in Yu et al. (2016), and find that the peak-to-trough ratio changes as a function of time as CO is destroyed. Specifically, a CO-depleted disk with high turbulent velocity mimics the peak-to-trough ratios of a non-CO-depleted disk with lower turbulent velocity. We suggest that disk observers and modelers take into account the possibility of CO depletion when using line peak-to-trough ratios to constrain the degree of turbulence in disks. Assuming that CO/H2 = 0.0001 at all disk radii can lead to underestimates of turbulent speeds in the disk by at least 0.2 km/s.

Planets are born in protostellar disks, which are now observed with enough resolution to address questions about internal gas flows. Candidates for driving the flows include magnetic forces, but ionization state estimates suggest much of the gas mass decouples from magnetic fields. Thus, hydrodynamical instabilities could play a major role. We investigate disk dynamics under conditions typical for a T Tauri system, using global 3D radiation hydrodynamics simulations with embedded particles and a resolution of 70 cells per scale height. Stellar irradiation heating is included with realistic dust opacities. The disk starts in joint radiative balance and hydrostatic equilibrium. The vertical shear instability (VSI) develops into turbulence that persists up to at least 1600 inner orbits (143 outer orbits). Turbulent speeds are a few percent of the local sound speed at the midplane, increasing to 20%, or 100 m/s, in the corona. These are consistent with recent upper limits on turbulent speeds from optically thin and thick molecular line observations of TW Hya and HD 163296. The predominantly vertical motions induced by the VSI efficiently lift particles upwards. Grains 0.1 and 1 mm in size achieve scale heights greater than expected in isotropic turbulence. We conclude that while kinematic constraints from molecular line emission do not directly discriminate between magnetic and nonmagnetic disk models, the small dust scale heights measured in HL Tau and HD 163296 favor turbulent magnetic models, which reach lower ratios of the vertical kinetic energy density to the accretion stress.

Recent observations of the protoplanetary disk around the Herbig Be star HD 100546 show two bright features in infrared (H and L' bands) at about 50 AU, with one so far unexplained. We explore the observational signatures of a high mass planet causing shock heating in order to determine if it could be the source of the unexplained infrared feature in HD 100546. More fundamentally, we identify and characterize planetary shocks as an extra, hitherto ignored, source of luminosity in transition disks. The RADMC-3D code is used to perform dust radiative transfer calculations on the hydrodynamical disk models, including volumetric heating. A stronger shock heating rate by a factor 20 would be necessary to qualitatively reproduce the morphology of the second infrared source. Instead, we find that the outer edge of the gap carved by the planet heats up by about 50% relative to the initial reference temperature, which leads to an increase in the scale height. The bulge is illuminated by the central star, producing a lopsided feature in scattered light, as the outer gap edge shows an asymmetry in density and temperature attributable to a secondary spiral arm launched not from the Lindblad resonances but from the 2:1 resonance. We conclude that high-mass planets lead to shocks in disks that may be directly observed, particularly at wavelengths of 10 mm or longer, but that they are more likely to reveal their presence in scattered light by puffing up their outer gap edges and exciting multiple spiral arms.

Stellar rotation is a key in our understanding of both mass-loss and evolution of intermediate and massive stars. It can lead to anisotropic mass-loss in the form of radiative wind or an excretion disk. We wished to spatially resolve the photosphere and gaseous environment of 51 Oph, a peculiar star with a very high vsin(i) of 267km s$^{-1}$ and an evolutionary status that remains unsettled. It has been classified by different authors as a Herbig, a $\beta$ Pic, or a classical Be star. We used the VEGA visible beam combiner installed on the CHARA array that reaches a submilliarcsecond resolution. Observation were centered on the H$\alpha$ emission line. We derived, for the first time, the extension and flattening of 51 Oph photosphere. We found a major axis of $\theta_{{\mathrm{eq}}}$=8.08$\pm$0.70$R_\odot$ and a minor axis of $\theta_{{\mathrm{pol}}}$=5.66$\pm$0.23$R_\odot$ . This high photosphere distortion shows that the star is rotating close to its critical velocity. Finally, using spectro-interferometric measurements in the H$ \alpha$ line, we constrained the circumstellar environment geometry and kinematics and showed that the emission is produced in a 5.2$\pm$2R$_{*}$ disk in Keplerian rotation. From the visible point of view, 51 Oph presents all the features of a classical Be star: near critical-rotation and double-peaked H$\alpha $ line in emission produced in a gaseous disk in Keplerian rotation. However, this does not explain the presence of dust as seen in the mid-infrared and millimeter spectra, and the evolutionary status of 51 Oph remains unsettled.

We report the results of high-angular-resolution observations that search for exozodiacal light in a sample of main sequence stars and sub-giants. Using the "jouvence" of the fiber linked unit for optical recombination (JouFLU) at the center for high angular resolution astronomy (CHARA) telescope array, we have observed a total of 44 stars. Out of the 44 stars, 33 are new stars added to the initial, previously published survey of 42 stars performed at CHARA with the fiber linked unit for optical recombiation (FLUOR). Since the start of the survey extension, we have detected a K-band circumstellar excess for six new stars at the ~ 1\% level or higher, four of which are known or candidate binaries, and two for which the excess could be attributed to exozodiacal dust. We have also performed follow-up observations of 11 of the stars observed in the previously published survey and found generally consistent results. We do however detect a significantly larger excess on three of these follow-up targets: Altair, $\upsilon$ And and $\kappa$ CrB. Interestingly, the last two are known exoplanet host stars. We perform a statistical analysis of the JouFLU and FLUOR samples combined, which yields an overall exozodi detection rate of $21.7^{+5.7}_{-4.1}\%$. We also find that the K-band excess in FGK-type stars correlates with the existence of an outer reservoir of cold ($\lesssim 100\,$K) dust at the $99\%$ confidence level, while the same cannot be said for A-type stars.

We report the results of high-angular-resolution observations that search for exozodiacal light in a sample of main sequence stars and sub-giants. Using the "jouvence" of the fiber linked unit for optical recombination (JouFLU) at the center for high angular resolution astronomy (CHARA) telescope array, we have observed a total of 44 stars. Out of the 44 stars, 33 are new stars added to the initial, previously published survey of 42 stars performed at CHARA with the fiber linked unit for optical recombiation (FLUOR). Since the start of the survey extension, we have detected a K-band circumstellar excess for six new stars at the ~ 1\% level or higher, four of which are known or candidate binaries, and two for which the excess could be attributed to exozodiacal dust. We have also performed follow-up observations of 11 of the stars observed in the previously published survey and found generally consistent results. We do however detect a significantly larger excess on three of these follow-up targets: Altair, $\upsilon$ And and $\kappa$ CrB. Interestingly, the last two are known exoplanet host stars. We perform a statistical analysis of the JouFLU and FLUOR samples combined, which yields an overall exozodi detection rate of $21.7^{+5.7}_{-4.1}\%$. We also find that the K-band excess in FGK-type stars correlates with the existence of an outer reservoir of cold ($\lesssim 100\,$K) dust at the $99\%$ confidence level, while the same cannot be said for A-type stars.

The mechanism of angular momentum transport in protoplanetary disks is fundamental to understand the distributions of gas and dust in the disks. The unprecedented, high spatial resolution ALMA observations taken toward HL Tau and subsequent radiative transfer modeling reveal that a high degree of dust settling is currently achieved at the outer part of the HL Tau disk. Previous observations however suggest a high disk accretion rate onto the central star. This configuration is not necessarily intuitive in the framework of the conventional viscous disk model, since efficient accretion generally requires a high level of turbulence, which can suppress dust settling considerably. We develop a simplified, semi-analytical disk model to examine under what condition these two properties can be realized in a single model. Recent, non-ideal MHD simulations are utilized to realistically model the angular momentum transport both radially via MHD turbulence and vertically via magnetically induced disk winds. We find that the HL Tau disk configuration can be reproduced well when disk winds are properly taken into account. While the resulting disk properties are likely consistent with other observational results, such an ideal situation can be established only if the plasma $\beta$ at the disk midplane is $\beta_0 \simeq 2 \times 10^4$ under the assumption of steady accretion. Equivalently, the vertical magnetic flux at 100 au is about 0.2 mG. More detailed modeling is needed to fully identify the origin of the disk accretion and quantitatively examine plausible mechanisms behind the observed gap structures in the HL Tau disk.

Gas in protostellar disks provides the raw material for giant planet formation and controls the dynamics of the planetesimal-building dust grains. Accurate gas mass measurements help map the observed properties of planet-forming disks onto the formation environments of known exoplanets. Rare isotopologues of carbon monoxide (CO) have been used as gas mass tracers for disks in the Lupus star-forming region, with an assumed interstellar CO/H$_2$ abundance ratio. Unfortunately, observations of T-Tauri disks show that CO abundance is not interstellar---a finding reproduced by models that show CO abundance decreasing both with distance from the star and as a function of time. Here we present radiative transfer simulations that assess the accuracy of CO-based disk mass measurements. We find that the combination of CO chemical depletion in the outer disk and optically thick emission from the inner disk leads observers to underestimate gas mass by more than an order of magnitude if they use the standard assumptions of interstellar CO/H$_2$ ratio and optically thin emission. Furthermore, CO abundance changes on million-year timescales, introducing an age/mass degeneracy into observations. To reach factor of a few accuracy for CO-based disk mass measurements, we suggest that observers and modelers adopt the following strategies: (1) select the low-$J$ transitions; (2) observe multiple CO isotopologues and use either intensity ratios or normalized line profiles to diagnose CO chemical depletion; and (3) use spatially resolved observations to measure the CO abundance distribution.

Many planets orbit within an AU of their stars, raising questions about their origins. Particularly puzzling are the planets found near the silicate sublimation front. We investigate conditions near the front in the protostellar disk around a young intermediate-mass star, using the first global 3-D radiation non-ideal MHD simulations in this context. We treat the starlight heating; the silicate grains sublimation and deposition at the local, time-varying temperature and density; temperature-dependent Ohmic dissipation; and various initial magnetic fields. The results show magnetorotational turbulence around the sublimation front at 0.5 AU. The disk interior to 0.8 AU is turbulent, with velocities exceeding 10% of the sound speed. Beyond 0.8 AU is the dead zone, cooler than 1000 K and with turbulence orders of magnitude weaker. A local pressure maximum just inside the dead zone concentrates solid particles, favoring their growth. Over many orbits, a vortex develops at the dead zones inner edge, increasing the disks thickness locally by around 10%. We synthetically observe the results using Monte Carlo transfer calculations, finding the sublimation front is near-infrared bright. The models with net vertical magnetic fields develop extended, magnetically-supported atmospheres that reprocess extra starlight, raising the near-infrared flux 20%.The vortex throws a nonaxisymmetric shadow on the outer disk. At wavelengths > 2 micron, the flux varies several percent on monthly timescales. The variations are more regular when the vortex is present. The vortex is directly visible as an arc at ultraviolet through near-infrared wavelengths, given sub-AU spatial resolution.

We present interferometric observations of the sigma Orionis triple system using the CHARA Array, NPOI, and VLTI. Using these measurements, we spatially resolve the orbit of the close spectroscopic binary (Aa,Ab) for the first time and present a revised orbit for the wide pair (A,B). Combining the visual orbits with previously published radial velocity measurements and new radial velocities measured at CTIO, we derive dynamical masses for the three massive stars in the system of M_Aa = 16.99 +/- 0.20 Msun, M_Ab = 12.81 +/- 0.18 Msun, and M_B = 11.5 +/- 1.2 Msun. The inner and outer orbits in the triple are not coplanar, with a relative inclination of 120-127 deg. The orbital parallax provides a precise distance of 387.5 +/- 1.3 pc to the system. This is a significant improvement over previous estimates of the distance to the young sigma Orionis cluster.

We perform two-dimensional axially symmetric radiation-hydrodynamic simulations to assess the impact of outflows and radiative force feedback from massive protostars by varying when the protostellar outflow starts, the ratio of ejection to accretion rates, and the strength of the wide angle disk wind component. The star formation efficiency, i.e. the ratio of final stellar mass to initial core mass, is dominated by radiative forces and the ratio of outflow to accretion rates. Increasing this ratio has three effects: First, the protostar grows slower with a lower luminosity at any given time, lowering radiative feedback. Second, bipolar cavities cleared by the outflow are larger, further diminishing radiative feedback on disk and core scales. Third, the higher momentum outflow sweeps up more material from the collapsing envelope, decreasing the protostar's potential mass reservoir via entrainment. The star formation efficiency varies with the ratio of ejection to accretion rates from 50% in the case of very weak outflows to as low as 20% for very strong outflows. At latitudes between the low density bipolar cavity and the high density accretion disk, wide angle disk winds remove some of the gas, which otherwise would be part of the accretion flow onto the disk; varying the strength of these wide angle disk winds, however, alters the final star formation efficiency by only +/-6%. For all cases, the opening angle of the bipolar outflow cavity remains below 20 degree during early protostellar accretion phases, increasing rapidly up to 65 degree at the onset of radiation pressure feedback.

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

We report new CHARA/MIRC interferometric observations of the Cepheid archetype $\delta$ Cep, which aimed at detecting the newly discovered spectroscopic companion. We reached a maximum dynamic range $\Delta H $ = 6.4, 5.8, and 5.2 mag, respectively within the relative distance to the Cepheid $r < 25$ mas, $25 < r < 50$ mas and $50 < r < 100$ mas. Our observations did not show strong evidence of a companion. We have a marginal detection at $3\sigma$ with a flux ratio of 0.21\%, but nothing convincing as we found other possible probable locations. We ruled out the presence of companion with a spectral type earlier than F0V, A1V and B9V, respectively for the previously cited ranges $r$. From our estimated sensitivity limits and the Cepheid light curve, we derived lower-limit magnitudes in the $H$ band for this possible companion to be $H_\mathrm{comp} > 9.15, 8.31$ and 7.77 mag, respectively for $r < 25$ mas, $25 < r < 50$ mas and $50 < r < 100$ mas. We also found that to be consistent with the predicted orbital period, the companion has to be located at a projected separation $< 24$ mas with a spectral type later than a F0V star.

Many stars host planets orbiting within a few astronomical units (AU). The occurrence rate and distributions of masses and orbits vary greatly with the host stars mass. These close planets origins are a mystery that motivates investigating protoplanetary disks central regions. A key factor governing the conditions near the star is the silicate sublimation front, which largely determines where the starlight is absorbed, and which is often called the inner rim. We present the first radiation hydrodynamical modeling of the sublimation front in the disks around the young intermediate-mass stars called Herbig Ae stars. The models are axisymmetric, and include starlight heating, silicate grains sublimating and condensing to equilibrium at the local, time-dependent temperature and density, and accretion stresses parametrizing the results of MHD magneto-rotational turbulence models. The results compare well with radiation hydrostatic solutions, and prove to be dynamically stable. Passing the model disks into Monte Carlo radiative transfer calculations, we show that the models satisfy observational constraints on the inner rims location. A small optically-thin halo of hot dust naturally arises between the inner rim and the star. The inner rim has a substantial radial extent, corresponding to several disk scale heights. While the fronts overall position varies with the stellar luminosity, its radial extent depends on the mass accretion rate. A pressure maximum develops near the location of thermal ionization at temperatures about 1000 K. The pressure maximum is capable of halting solid pebbles radial drift and concentrating them in a zone where temperatures are sufficiently high for annealing to form crystalline silicates.

The gas near the midplanes of planet-forming protostellar disks remains largely unprobed by observations due to the high optical depth of commonly observed molecules such as CO and H$_2$O. However, rotational emission lines from rare molecules may have optical depths near unity in the vertical direction, so that the lines are strong enough to be detected, yet remain transparent enough to trace the disk midplane. Here we present a chemical model of an evolving T-Tauri disk and predict the optical depths of rotational transitions of $^{12}$C$^{16}$O, $^{13}$C$^{16}$O, $^{12}$C$^{17}$O and $^{12}$C$^{18}$O. The MRI-active disk is primarily heated by the central star due to the formation of the dead zone. CO does not freeze out in our modeled region within $70$AU around a sunlike star. However, the abundance of CO decreases because of the formation of complex organic molecules (COM), producing an effect that can be misinterpreted as the "snow line". These results are robust to variations in our assumptions about the evolution of the gas to dust ratio. The optical depths of low-order rotational lines of C$^{17}$O are around unity, making it possible to see into the disk midplane using C$^{17}$O. Combining observations with modeled C$^{17}$O$/$H$_2$ ratios, like those we provide, can yield estimates of protoplanetary disks' gas masses.

Chondritic meteorites provide valuable opportunities to investigate the origins of the solar system. We explore impact jetting as a mechanism of chondrule formation and subsequent pebble accretion as a mechanism of accreting chondrules onto parent bodies of chondrites, and investigate how these two processes can account for the currently available meteoritic data. We find that when the solar nebula is $\le 5$ times more massive than the minimum-mass solar nebula at $a \simeq 2-3$ AU and parent bodies of chondrites are $\le 10^{24}$ g ($\le$ 500 km in radius) in the solar nebula, impact jetting and subsequent pebble accretion can reproduce a number of properties of the meteoritic data. The properties include the present asteroid belt mass, the formation timescale of chondrules, and the magnetic field strength of the nebula derived from chondrules in Semarkona. Since this scenario requires a first generation of planetesimals that trigger impact jetting and serve as parent bodies to accrete chondrules, the upper limit of parent bodies' masses leads to the following implications: primordial asteroids that were originally $\ge 10^{24}$ g in mass were unlikely to contain chondrules, while less massive primordial asteroids likely had a chondrule-rich surface layer. The scenario developed from impact jetting and pebble accretion can therefore provide new insights into the origins of the solar system.

We provide CoRoT and Spitzer light curves, as well as broad-band multi-wavelength photometry and high resolution, multi- and single-epoch spectroscopy for 17 classical T Tauris in NGC 2264 whose CoRoT light curves (LCs) exemplify the "stochastic" LC class as defined in Cody et al. (2014). The most probable physical mechanism to explain the optical variability in this LC class is time-dependent mass accretion onto the stellar photosphere, producing transient hot spots. As evidence in favor of this hypothesis, multi-epoch high resolution spectra for a subset of these stars shows that their veiling levels also vary in time and that this veiling variability is consistent in both amplitude and timescale with the optical LC morphology. Furthermore, the veiling variability is well-correlated with the strength of the HeI 6678A emission line, a feature predicted by models to arise in accretion shocks on or near the stellar photosphere. Stars with accretion burst LC morphology (Stauffer et al. 2014) are also attributed to variable mass accretion. Both the stochastic and accretion burst LCs can be explained by a simple model of randomly occurring flux bursts, with the stochastic LC class having a higher frequency of lower amplitude events. Based on their UV excesses, veiling, and mean Ha equivalent widths, members of the stochastic LC class have only moderate time-averaged mass accretion rates. The most common feature of their Ha profiles is for them to exhibit blue-shifted absorption features, most likely originating in a disk wind. The lack of periodic signatures in the LCs suggests that little of the variability is due to long-lived hot spots rotating into or out of our line of sight; instead, the primary driver of the observed photometric variability is likely to be instabilities in the inner disk that lead to variable mass accretion.

This is a joint summary of the reports from the three Astrophysics Program Analysis Groups (PAGs) in response to the "Planning for the 2020 Decadal Survey" charge given by the Astrophysics Division Director Paul Hertz. This joint executive summary contains points of consensus across all three PAGs. Additional findings specific to the individual PAGs are reported separately in the individual reports. The PAGs concur that all four large mission concepts identified in the white paper as candidates for maturation prior to the 2020 Decadal Survey should be studied in detail. These include the Far-IR Surveyor, the Habitable-Exoplanet Imaging Mission, the UV/Optical/IR Surveyor, and the X-ray Surveyor. This finding is predicated upon assumptions outlined in the white paper and subsequent charge, namely that 1) major development of future large flagship missions under consideration are to follow the implementation phases of JWST and WFIRST; 2) NASA will partner with the European Space Agency on its L3 Gravitational Wave Surveyor; 3) the Inflation Probe be classified as a probe-class mission to be developed according to the 2010 Decadal Survey report. If these key assumptions were to change, this PAG finding would need to be re-evaluated. The PAGs find that there is strong community support for the second phase of this activity - maturation of the four proposed mission concepts via Science and Technology Definition Teams (STDTs). The PAGs find that there is strong consensus that all of the STDTs contain broad and interdisciplinary representation of the science community. Finally, the PAGs find that there is community support for a line of Probe-class missions within the Astrophysics mission portfolio (condensed).

Recent high-resolution, near-infrared images of protoplanetary disks have shown that these disks often present spiral features. Spiral arms are among the structures predicted decades ago by numerical simulations of disk-planet interaction and thus it is tempting to suspect that planetary perturbers are responsible for the observed signatures. However, such interpretation is not free of problems. The spirals are found to have large pitch angles, and in at least one case (HD 100546) the spiral feature appears effectively unpolarized, which implies thermal emission of the order of 1000K (465$\pm$40K at closer inspection). We have recently shown in two-dimensional models that shock dissipation in the supersonic wake of high-mass planets can lead to significant heating if the disk is sufficiently adiabatic. In this paper we extend this analysis to three dimensions in thermodynamically evolving disks. We use the Pencil Code in spherical coordinates for our models, with a prescription for thermal cooling based on the optical depth of the local vertical gas column. We use a 5$M_J$ planet, and show that shocks in the region around the planet where the Lindblad resonances occur heat the gas to substantially higher temperatures than the ambient disk gas at that radius. The gas is accelerated vertically away from the midplane by the shocks to form shock bores, and the gas falling back toward the midplane breaks up into a turbulent surf near the Lindblad resonances. This turbulence, although localized, has high $\alpha$ values, reaching 0.05 in the inner Lindblad resonance, and 0.1 in the outer one. We also find evidence that the disk regions heated up by the planetary shocks eventually becomes superadiabatic, generating convection far from the planet's orbit.

The parallax of pulsation, and its implementations such as the Baade-Wesselink method and the infrared surface bright- ness technique, is an elegant method to determine distances of pulsating stars in a quasi-geometrical way. However, these classical implementations in general only use a subset of the available observational data. Freedman & Madore (2010) suggested a more physical approach in the implementation of the parallax of pulsation in order to treat all available data. We present a global and model-based parallax-of-pulsation method that enables including any type of observational data in a consistent model fit, the SpectroPhoto-Interferometric modeling of Pulsating Stars (SPIPS). We implemented a simple model consisting of a pulsating sphere with a varying effective temperature and a combina- tion of atmospheric model grids to globally fit radial velocities, spectroscopic data, and interferometric angular diameters. We also parametrized (and adjusted) the reddening and the contribution of the circumstellar envelopes in the near-infrared photometric and interferometric measurements. We show the successful application of the method to two stars: delta Cep and eta Aql. The agreement of all data fitted by a single model confirms the validity of the method. Derived parameters are compatible with publish values, but with a higher level of confidence. The SPIPS algorithm combines all the available observables (radial velocimetry, interferometry, and photometry) to estimate the physical parameters of the star (ratio distance/ p-factor, Teff, presence of infrared excess, color excess, etc). The statistical precision is improved (compared to other methods) thanks to the large number of data taken into account, the accuracy is improved by using consistent physical modeling and the reliability of the derived parameters is strengthened thanks to the redundancy in the data.

MWC 314 is a bright candidate luminous blue variable that resides in a fairly close binary system, with an orbital period of 60.753$\pm$0.003 d. We observed MWC 314 with a combination of optical spectroscopy, broad-band ground- and space-based photometry, as well as with long baseline, near-infrared interferometry. We have revised the single-lined spectroscopic orbit and explored the photometric variability. The orbital light curve displays two minima each orbit that can be partially explained in terms of the tidal distortion of the primary that occurs around the time of periastron. The emission lines in the system are often double-peaked and stationary in their kinematics, indicative of a circumbinary disc. We find that the stellar wind or circumbinary disc is partially resolved in the K\prime-band with the longest baselines of the CHARA Array. From this analysis, we provide a simple, qualitative model in an attempt to explain the observations. From the assumption of Roche Lobe overflow and tidal synchronisation at periastron, we estimate the component masses to be M1 $\approx 5$ M$_\odot$ and M2$\approx 15$ M$_\odot$, which indicates a mass of the LBV that is extremely low. In addition to the orbital modulation, we discovered two pulsational modes with the MOST satellite. These modes are easily supported by a low-mass hydrogen-poor star, but cannot be easily supported by a star with the parameters of an LBV. The combination of these results provides evidence that the primary star was likely never a normal LBV, but rather is the product of binary interactions. As such, this system presents opportunities for studying mass-transfer and binary evolution with many observational techniques.

We have observed and spatially resolved a set of seven A-type stars in the nearby Ursa Major moving group with the Classic, CLIMB, and PAVO beam combiners on the CHARA Array. At least four of these stars have large rotational velocities ($v \sin i$ $\gtrsim$ 170 $\mathrm{km~s^{-1}}$) and are expected to be oblate. These interferometric measurements, the stars' observed photometric energy distributions, and $v \sin i$ values are used to computationally construct model oblate stars from which stellar properties (inclination, rotational velocity, and the radius and effective temperature as a function of latitude, etc.) are determined. The results are compared with MESA stellar evolution models (Paxton et al. 2011, 2013) to determine masses and ages. The value of this new technique is that it enables the estimation of the fundamental properties of rapidly rotating stars without the need to fully image the star. It can thus be applied to stars with sizes comparable to the interferometric resolution limit as opposed to those that are several times larger than the limit. Under the assumption of coevality, the spread in ages can be used as a test of both the prescription presented here and the MESA evolutionary code for rapidly rotating stars. With our validated technique, we combine these age estimates and determine the age of the moving group to be 414 $\pm$ 23 Myr, which is consistent with, but much more precise than previous estimates.

We calculate the abundances of electrons and ions in the hot (> 500 K), dusty parts of protoplanetary disks, treating for the first time the effects of thermionic and ion emission from the dust grains. High-temperature ionization modeling has involved simply assuming that alkali elements such as potassium occur as gas-phase atoms and are collisionally ionized following the Saha equation. We show that the Saha equation often does not hold, because free charges are produced by thermionic and ion emission and destroyed when they stick to grain surfaces. This means the ionization state depends not on the first ionization potential of the alkali atoms, but rather on the grains' work functions. The charged species' abundances typically rise abruptly above about 800 K, with little qualitative dependence on the work function, gas density, or dust-to-gas mass ratio. Applying our results, we find that protoplanetary disks' dead zone, where high diffusivities stifle magnetorotational turbulence, has its inner edge located where the temperature exceeds a threshold value ~1000 K. The threshold is set by ambipolar diffusion except at the highest densities, where it is set by Ohmic resistivity. We find that the disk gas can be diffusively loaded onto the stellar magnetosphere at temperatures below a similar threshold. We investigate whether the "short-circuit" instability of current sheets can operate in disks and find that it cannot, or works only in a narrow range of conditions; it appears not to be the chondrule formation mechanism. We also suggest that thermionic emission is important for determining the rate of Ohmic heating in hot Jupiters.

Presented are the first interferometric images of cool starspots on the chromospherically active giant $\lambda$ Andromedae. These images represent the first model-independent images of cool starspots on a star other than the Sun to date. The interferometric observations, taken with the Michigan Infra-Red Combiner coupled to the Center for High Angular Resolution Astronomy Array, span 26 days from Aug 17$^{th}$, 2008 to Sep 24$^{th}$, 2011. The photometric time series acquired at Fairborn Observatory spanning Sep 20$^{th}$, 2008 to Jan 20$^{th}$, 2011 is also presented. The angular diameter and power law limb-darkening coefficient of this star are 2.759 $\pm$ 0.050 mas and 0.229 $\pm$ 0.111, respectively. Starspot properties are obtained from both modeled and SQUEEZE reconstructed images. The images from 2010 through 2011 show anywhere from one to four starspots. The measured properties of identical starspots identified in both the model and reconstructed images are within two $\sigma$ error bars in 51$\%$ of cases. The cadence in the data for the 2010 and 2011 data sets are sufficient to measure a stellar rotation period based on apparent starspot motion. This leads to estimates of the rotation period (P$_{2010}$ = 60 $\pm$ 13 days, P$_{2011}$ = 54.0 $\pm$ 7.6 days) that are consistent with the photometrically determined period of 54.8 days. In addition, the inclination and position angle of the rotation axis is computed for both the 2010 and 2011 data sets; values ($\bar{\Psi}$ = 21.5$\degree$, $\bar{\emph{i}}$ = 78.0$\degree$) for each are nearly identical between the two years. \end{abstract}

We report on a total of 106 nights of optical interferometric observations of the $\epsilon$ Aurigae system taken during the last 14 years by four beam combiners at three different interferometric facilities. This long sequence of data provides an ideal assessment of the system prior to, during, and after the recent 2009-2011 eclipse. We have reconstructed model-independent images from the 10 in-eclipse epochs which show that a disk-like object is indeed responsible for the eclipse. Using new 3D, time-dependent modeling software, we derive the properties of the F-star (diameter, limb darkening), determine previously unknown orbital elements ($\Omega$, $i$), and access the global structures of the optically thick portion of the eclipsing disk using both geometric models and approximations of astrophysically relevant density distributions. These models may be useful in future hydrodynamical modeling of the system. Lastly, we address several outstanding research questions including mid-eclipse brightening, possible shrinking of the F-type primary, and any warps or sub-features within the disk.

To measure the stellar and orbital properties of the metal-poor RS CVn binary o Draconis (o Dra), we directly detect the companion using interferometric observations obtained with the Michigan InfraRed Combiner at Georgia State University's Center for High Angular Resolution Astronomy (CHARA) Array. The H-band flux ratio between the primary and secondary stars is the highest confirmed flux ratio (370 +/- 40) observed with long-baseline optical interferometry. These detections are combined with radial velocity data of both the primary and secondary stars, including new data obtained with the Tillinghast Reflector Echelle Spectrograph on the Tillinghast Reflector at the Fred Lawrence Whipple Observatory and the 2-m Tennessee State University Automated Spectroscopic Telescope at Fairborn Observatory. We determine an orbit from which we find model-independent masses and ages of the components (M_A = 1.35 +\- 0.05 M_Sun, M_B = 0.99 +\- 0.02 M_Sun, system age = 3.0 -\+ 0.5 Gyr). An average of a 23-year light curve of o Dra from the Tennessee State University Automated Photometric Telescope folded over the orbital period newly reveals eclipses and the quasi-sinusoidal signature of ellipsoidal variations. The modeled light curve for our system's stellar and orbital parameters confirm these ellipsoidal variations due to the primary star partially filling its Roche lobe potential, suggesting most of the photometric variations are not due to stellar activity (starspots). Measuring gravity darkening from the average light curve gives a best-fit of beta = 0.07 +\- 0.03, a value consistent with conventional theory for convective envelope stars. The primary star also exhibits an anomalously short rotation period, which, when taken with other system parameters, suggests the star likely engulfed a low-mass companion that had recently spun-up the star.

Long-baseline interferometry is an important technique to spatially resolve binary or multiple systems in close orbits. By combining several telescopes together and spectrally dispersing the light, it is possible to detect faint components around bright stars. Aims. We provide a rigorous and detailed method to search for high-contrast companions around stars, determine the detection level, and estimate the dynamic range from interferometric observations. We developed the code CANDID (Companion Analysis and Non-Detection in Interferometric Data), a set of Python tools that allows us to search systematically for point-source, high-contrast companions and estimate the detection limit. The search pro- cedure is made on a N x N grid of fit, whose minimum needed resolution is estimated a posteriori. It includes a tool to estimate the detection level of the companion in the number of sigmas. The code CANDID also incorporates a robust method to set a 3{\sigma} detection limit on the flux ratio, which is based on an analytical injection of a fake companion at each point in the grid. We used CANDID to search for the companions around the binary Cepheids V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul. First, we showed that our previous discoveries of the components orbiting V1334 Cyg and AX Cir were detected at > 13 sigmas. The companion around AW Per is detected at more than 15 sigmas with a flux ratio of f = 1.22 +/- 0.30 %. We made a possible detection of the companion orbiting RT Aur with f = 0.22 +/- 0.11 %. It was detected at 3.8{\sigma} using the closure phases only, and so more observations are needed to confirm the detection. We also set the detection limit for possible undetected companions. We found that there is no companion with a spectral type earlier than B7V, A5V, F0V, B9V, A0V, and B9V orbiting V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul, respectively.

A classical nova occurs when material accreting onto the surface of a white dwarf in a close binary system ignites in a thermonuclear runaway. Complex structures observed in the ejecta at late stages could result from interactions with the companion during the common envelope phase. Alternatively, the explosion could be intrinsically bipolar, resulting from a localized ignition on the surface of the white dwarf or as a consequence of rotational distortion. Studying the structure of novae during the earliest phases is challenging because of the high spatial resolution needed to measure their small sizes. Here we report near-infrared interferometric measurements of the angular size of Nova Delphini 2013, starting from one day after the explosion and continuing with extensive time coverage during the first 43 days. Changes in the apparent expansion rate can be explained by an explosion model consisting of an optically thick core surrounded by a diffuse envelope. The optical depth of the ejected material changes as it expands. We detect an ellipticity in the light distribution, suggesting a prolate or bipolar structure that develops as early as the second day. Combining the angular expansion rate with radial velocity measurements, we derive a geometric distance to the nova of 4.54 +/- 0.59 kpc from the Sun.

To measure the properties of both components of the RS CVn binary sigma Geminorum (sigma Gem), we directly detect the faint companion, measure the orbit, obtain model-independent masses and evolutionary histories, detect ellipsoidal variations of the primary caused by the gravity of the companion, and measure gravity darkening. We detect the companion with interferometric observations obtained with the Michigan InfraRed Combiner (MIRC) at Georgia State University's Center for High Angular Resolution Astronomy (CHARA) Array with a primary-to-secondary H-band flux ratio of 270+/-70. A radial velocity curve of the companion was obtained with spectra from the Tillinghast Reflector Echelle Spectrograph (TRES) on the 1.5-m Tillinghast Reflector at Fred Lawrence Whipple Observatory (FLWO). We additionally use new observations from the Tennessee State University Automated Spectroscopic and Photometric Telescopes (AST and APT, respectively). From our orbit, we determine model-independent masses of the components (M_1 = 1.28 +/- 0.07 M_Sun, M_2 = 0.73 +/- 0.03 M_Sun), and estimate a system age of 5 -/+ 1 Gyr. An average of the 27-year APT light curve of sigma Gem folded over the orbital period (P = 19.6027 +/- 0.0005 days) reveals a quasi-sinusoidal signature, which has previously been attributed to active longitudes 180 deg apart on the surface of sigma Gem. With the component masses, diameters, and orbit, we find that the predicted light curve for ellipsoidal variations due to the primary star partially filling its Roche lobe potential matches well with the observed average light curve, offering a compelling alternative explanation to the active longitudes hypothesis. Measuring gravity darkening from the light curve gives beta < 0.1, a value slightly lower than that expected from recent theory.

Recent observations of gaps and non-axisymmetric features in the dust distributions of transition disks have been interpreted as evidence of embedded massive protoplanets. However, comparing the predictions of planet-disk interaction models to the observed features has shown far from perfect agreement. This may be due to the strong approximations used for the predictions. For example, spiral arm fitting typically uses results that are based on low-mass planets in an isothermal gas. In this work, we describe two-dimensional, global, hydrodynamical simulations of disks with embedded protoplanets, with and without the assumption of local isothermality, for a range of planet-to-star mass ratios 1-10 M_jup for a 1 M_sun star. We use the Pencil Code in polar coordinates for our models. We find that the inner and outer spiral wakes of massive protoplanets (M>5 M_jup) produce significant shock heating that can trigger buoyant instabilities. These drive sustained turbulence throughout the disk when they occur. The strength of this effect depends strongly on the mass of the planet and the thermal relaxation timescale; for a 10 M_jup planet embedded in a thin, purely adiabatic disk, the spirals, gaps, and vortices typically associated with planet-disk interactions are disrupted. We find that the effect is only weakly dependent on the initial radial temperature profile. The spirals that form in disks heated by the effects we have described may fit the spiral structures observed in transition disks better than the spirals predicted by linear isothermal theory.

AA Tau, a classical T Tauri star in the Taurus cloud, has been the subject of intensive photometric monitoring for more than two decades due to its quasi-cyclic variation in optical brightness. Beginning in 2011, AA Tau showed another peculiar variation -- its median optical though near-IR flux dimmed significantly, a drop consistent with a 4-mag increase in visual extinction. It has stayed in the faint state since.Here we present 4.7um CO rovibrational spectra of AA Tau over eight epochs, covering an eleven-year time span, that reveal enhanced 12CO and 13CO absorption features in the $J_{\rm low}\leqslant$13 transitions after the dimming. These newly appeared absorptions require molecular gas along the line of sight with T~500 K and a column density of log (N12CO)~18.5 cm^{-2}, with line centers that show a constant 6 km s$^{-1}$ redshift. The properties of the molecular gas confirm an origin in the circumstellar material. We suggest that the dimming and absorption are caused by gas and dust lifted to large heights by a magnetic buoyancy instability. This material is now propagating inward, and on reaching the star within a few years will be observed as an accretion outburst.

The rapidly rotating Be star phi Persei was spun up by mass and angular momentum transfer from a now stripped-down, hot subdwarf companion. Here we present the first high angular resolution images of phi Persei made possible by new capabilities in longbaseline interferometry at near-IR and visible wavelengths. We observed phi Persei with the MIRC and VEGA instruments of the CHARA Array. Additional MIRC-only observations were performed to track the orbital motion of the companion, and these were fit together with new and existing radial velocity measurements of both stars to derive the complete orbital elements and distance. The hot subdwarf companion is clearly detected in the near-IR data at each epoch of observation with a flux contribution of 1.5% in the H band, and restricted fits indicate that its flux contribution rises to 3.3% in the visible. A new binary orbital solution is determined by combining the astrometric and radial velocity measurements. The derived stellar masses are 9.6+-0.3Msol and 1.2+-0.2Msol for the Be primary and subdwarf secondary, respectively. The inferred distance (186 +- 3 pc), kinematical properties, and evolutionary state are consistent with membership of phi Persei in the alpha Per cluster. From the cluster age we deduce significant constraints on the initial masses and evolutionary mass transfer processes that transformed the phi Persei binary system. The interferometric data place strong constraints on the Be disk elongation, orientation, and kinematics, and the disk angular momentum vector is coaligned with and has the same sense of rotation as the orbital angular momentum vector. The VEGA visible continuum data indicate an elongated shape for the Be star itself, due to the combined effects of rapid rotation, partial obscuration of the photosphere by the circumstellar disk, and flux from the bright inner disk.

The classical T Tauri star AA Tau presented photometric variability attributed to an inner disk warp, caused by the interaction between the inner disk and an inclined magnetosphere. Previous studies of NGC 2264 have shown that similar photometric behavior is common among CTTS. The goal of this work is to investigate the main causes of the observed photometric variability of CTTS in NGC 2264 that present AA Tau-like light curves, and verify if an inner disk warp could be responsible for their variability. We investigate veiling variability in their spectra and u-r color variations and estimate parameters of the inner disk warp using an occultation model proposed for AA Tau. We compare infrared and optical light curves to analyze the dust responsible for the occultations. AA Tau-like variability is transient on a timescale of a few years. We ascribe it to stable accretion regimes and aperiodic variability to unstable accretion regimes and show that a transition, and even coexistence, between the two is common. We find evidence of hot spots associated with occultations, indicating that the occulting structures could be located at the base of accretion columns. We find average values of warp maximum height of 0.23 times its radial location, consistent with AA Tau, with variations of on average 11% between rotation cycles. We show that extinction laws in the inner disk indicate the presence of grains larger than interstellar grains. The inner disk warp scenario is consistent with observations for all but one periodic star in our sample. AA Tau-like systems comprise 14% of CTTS observed in NGC 2264, though this increases to 35% among systems of mass 0.7M_sun<M<2.0M_sun. Assuming random inclinations, we estimate that nearly all systems in this mass range likely possess an inner disk warp, possibly because of a change in magnetic field configurations among stars of lower mass.

We identify nine young stellar objects (YSOs) in the NGC 2264 star-forming region with optical {\em CoRoT} light curves exhibiting short-duration, shallow, periodic flux dips. All of these stars have infrared (IR) excesses that are consistent with their having inner disk walls near the Keplerian co-rotation radius. The repeating photometric dips have FWHM generally less than one day, depths almost always less than 15%, and periods (3<P<11 days) consistent with dust near the Keplerian co-rotation period. The flux dips vary considerably in their depth from epoch to epoch, but usually persist for several weeks and, in two cases, were present in data collected on successive years. For several of these stars, we also measure the photospheric rotation period and find that the rotation and dip periods are the same, as predicted by standard "disk-locking" models. We attribute these flux dips to clumps of material in or near the inner disk wall, passing through our line of sight to the stellar photosphere. In some cases, these dips are also present in simultaneous {\em Spitzer} IRAC light curves at 3.6 and 4.5 microns. We characterize the properties of these dips, and compare the stars with light curves exhibiting this behavior to other classes of YSO in NGC 2264. A number of physical mechanisms could locally increase the dust scale height near the inner disk wall, and we discuss several of those mechanisms; the most plausible mechanisms are either a disk warp due to interaction with the stellar magnetic field or dust entrained in funnel-flow accretion columns arising near the inner disk wall.

Protoplanetary disks are believed to accrete onto their central T Tauri star because of magnetic stresses. Recently published shearing box simulations indicate that Ohmic resistivity, ambipolar diffusion and the Hall effect all play important roles in disk evolution. In the presence of a vertical magnetic field, the disk remains laminar between 1-5au, and a magnetocentrifugal disk wind forms that provides an important mechanism for removing angular momentum. Questions remain, however, about the establishment of a true physical wind solution in the shearing box simulations because of the symmetries inherent in the local approximation. We present global MHD simulations of protoplanetary disks that include Ohmic resistivity and ambipolar diffusion, where the time-dependent gas-phase electron and ion fractions are computed under FUV and X-ray ionization with a simplified recombination chemistry. Our results show that the disk remains laminar, and that a physical wind solution arises naturally in global disk models. The wind is sufficiently efficient to explain the observed accretion rates. Furthermore, the ionization fraction at intermediate disk heights is large enough for magneto-rotational channel modes to grow and subsequently develop into belts of horizontal field. Depending on the ionization fraction, these can remain quasi-global, or break-up into discrete islands of coherent field polarity. The disk models we present here show a dramatic departure from our earlier models including Ohmic resistivity only. It will be important to examine how the Hall effect modifies the evolution, and to explore the influence this has on the observational appearance of such systems, and on planet formation and migration.

Since their discovery twenty year ago, transition region bright points (TRBPs) have never been observed spectroscopically. Bright point properties have not been compared with similar transition region and coronal structures. In this work we have investigated three transient quiet Sun brightenings including a TRBP, a coronal BP (CBP) and a blinker. We use time-series observations of the extreme ultraviolet emission lines of a wide range of temperature T (log T = 5.3 - 6.4) from the EUV imaging spectrometer (EIS) onboard the Hinode satellite. We present the EIS temperature maps and Doppler maps, which are compared with magnetograms from the Michelson Doppler Imager (MDI) onboard the SOHO satellite. Doppler velocities of the TR BP and blinker are <,25 km s$^{-1}$, which is typical of transient TR phenomena. The Dopper velocities of the CBP were found to be < 20 km s^{-1} with exception of those measured at log T = 6.2 where a distinct bi-directional jet is observed. From an EM loci analysis we find evidence of single and double isothermal components in the TRBP and CBP, respectively. TRBP and CBP loci curves are characterized by broad distributions suggesting the existence of unresolved structure. By comparing and contrasting the physical characteristics of the events we find the BP phenomena are an indication of multi-scaled self similarity, given similarities in both their underlying magnetic field configuration and evolution in relation to EUV flux changes. In contrast, the blinker phenomena and the TRBP are sufficiently dissimilar in their observed properties as to constitute different event classes. Our work indicates that the measurement of similar characteristics across multiple event types holds class-predictive power, and is a significant step towards automated solar atmospheric multi-class classification of unresolved transient EUV sources.

We carry out radiation hydrodynamical simulations of the formation of massive stars in the super-Eddington regime including both their radiative feedback and protostellar outflows. The calculations start from a prestellar core of dusty gas and continue until the star stops growing. The accretion ends when the remnants of the core are ejected, mostly by the force of the direct stellar radiation in the polar direction and elsewhere by the reradiated thermal infrared radiation. How long the accretion persists depends on whether the protostellar outflows are present. We set the mass outflow rate to 1% of the stellar sink particle's accretion rate. The outflows open a bipolar cavity extending to the core's outer edge, through which the thermal radiation readily escapes. The radiative flux is funneled into the polar directions while the core's collapse proceeds near the equator. The outflow thus extends the "flashlight effect", or anisotropic radiation field, found in previous studies from the few hundred AU scale of the circumstellar disk up to the 0.1 parsec scale of the core. The core's flashlight effect allows core gas to accrete on the disk for longer, in the same way that the disk's flashlight effect allows disk gas to accrete on the star for longer. Thus although the protostellar outflows remove material near the core's poles, causing slower stellar growth over the first few free-fall times, they also enable accretion to go on longer in our calculations. The outflows ultimately lead to stars of somewhat higher mass.

We determined the fundamental properties of HD 140283 by obtaining new interferometric and spectroscopic measurements and combining them with photometry from the literature. The interferometric measurements were obtained using the visible interferometer VEGA on the CHARA array and we determined a 1D limb-darkened angular diameter of 0.353 +/- 0.013 milliarcseconds. Using photometry from the literature we derived the bolometric flux with two solutions: a zero-reddening one of Fbol = 3.890 +/- 0.066 1E-8 erg/s/cm2 and a solution with a maximum of Av = 0.1 mag, Fbol= 4.220 +/- 0.067 1E-8 erg/s/cm2. The interferometric Teff is thus 5534 +/- 103 K or 5647 +/- 105 K and its radius is R = 2.21 +/- 0.08 Rsol. Spectroscopic measurements of HD140283 were obtained using HARPS, NARVAL, and UVES and a 1D LTE analysis of H-alpha line wings yields Teff(Halpha) = 5626 +/- 75 K. Using fine-tuned stellar models including diffusion of elements we then determined the mass M and age t of HD140283. Once the metallicity has been fixed, the age of the star depends on M, initial helium abundance Yi and mixing-length parameter alpha, only two of which are independent. We need to adjust alpha to much lower values than the solar one (~2) in order to fit the observations, and if Av = 0.0 mag then 0.5 < alpha < 1. We give an equation to estimate t from M, Yi (alpha) and Av. Establishing a reference alpha = 1.00 and adopting Yi = 0.245 we derive a mass and age of HD140283: M = 0.780 +/- 0.010 Msol and t = 13.7 +/- 0.7 Gyr (Av = 0.0) or M = 0.805 +/- 0.010 Msol and t = 12.2 +/- 0.6 Gyr (Av=0.1 mag). Our stellar models yield an initial metallicity of [Z/X]i = -1.70 and logg = 3.65 +/- 0.03. Asteroseismic observations are critical for overcoming limitations in our results.

We present direct radii measurements of the well-known transiting exoplanet host stars HD 189733 and HD 209458 using the CHARA Array interferometer. We find the limb-darkened angular diameters to be theta_LD = 0.3848 +/- 0.0055 and 0.2254 +/- 0.0072 milliarcsec for HD 189733 and HD 209458, respectively. HD 189733 and HD 209458 are currently the only two transiting exoplanet systems where detection of the respective planetary companion's orbital motion from high resolution spectroscopy has revealed absolute masses for both star and planet. We use our new measurements together with the orbital information from radial velocity and photometric time series data, Hipparcos distances, and newly measured bolometric fluxes to determine the stellar effective temperatures (T_eff = 4875 +/- 43, 6093 +/- 103 K), stellar linear radii (R_* = 0.805 +/- 0.016, 1.203 +/- 0.061 R_sun), mean stellar densities (rho_* = 1.62 +/- 0.11, 0.58 +/- 0.14 rho_sun), planetary radii (R_p = 1.216 +/- 0.024, 1.451 +/- 0.074 R_Jup), and mean planetary densities (rho_p = 0.605 +/- 0.029, 0.196 +/- 0.033 rho_Jup) for HD 189733 b and HD 209458 b, respectively. The stellar parameters for HD 209458, a F9 dwarf, are consistent with indirect estimates derived from spectroscopic and evolutionary modeling. However, we find that models are unable to reproduce the observational results for the K2 dwarf, HD 189733. We show that, for stellar evolutionary models to match the observed stellar properties of HD 189733, adjustments lowering the solar-calibrated mixing length parameter from 1.83 to 1.34 need to be employed.

Rossby wave instability (RWI) at dead zone boundaries may play an important role in planet formation. Viscous hydrodynamics results suggest RWI is excited only when the viscosity changes over a radial distance less than two density scale heights. However in the disks around Solar-mass T Tauri stars, it is not viscosity but magnetic forces that provide the accretion stress beyond about 10 AU, where surface densities are low enough so stellar X-rays and interstellar cosmic rays can penetrate. Here we aim to explore the conditions for RWI in the smooth transition with increasing distance, from resistive and magnetically-dead to conducting and magnetically-active. We perform 3D unstratified MHD simulations with the Pencil Code, using static resistivity profiles. As a result, we find that in MHD, contrary to viscous models, the RWI is triggered even with a gradual change in resistivity extending from 10 to 40 AU (i.e., spanning 15 scale heights for aspect ratio 0.1). This is because magneto-rotational turbulence sets in abruptly when the resistivity reaches a threshold level. At higher resistivities the longest unstable wavelength is quenched, resulting in a sharp decline of the Maxwell stress towards the star. The sharp gradient in the magnetic forces leads to a localized density bump, that is in turn Rossby wave unstable. We conclude that even weak gradients in the resistivity can lead to sharp transitions in the Maxwell stress. The upshot is that the RWI is more easily activated in the outer disk than previously thought. Rossby vortices at the outer dead zone boundary thus could underlie the dust asymmetries seen in the outer reaches of transition disks.

The aim of this work is to improve the SBC relation for early-type stars in the $-1 \leq V-K \leq 0$ color domain, using optical interferometry. Observations of eight B- and A-type stars were secured with the VEGA/CHARA instrument in the visible. The derived uniform disk angular diameters were converted into limb darkened angular diameters and included in a larger sample of 24 stars, already observed by interferometry, in order to derive a revised empirical relation for O, B, A spectral type stars with a V-K color index ranging from -1 to 0. We also took the opportunity to check the consistency of the SBC relation up to $V-K \simeq 4$ using 100 additional measurements. We determined the uniform disk angular diameter for the eight following stars: $\gamma$ Ori, $\zeta$ Per, $8$ Cyg, $\iota$ Her, $\lambda$ Aql, $\zeta$ Peg, $\gamma$ Lyr, and $\delta$ Cyg with V-K color ranging from -0.70 to 0.02 and typical precision of about $1.5\%$. Using our total sample of 132 stars with $V-K$ colors index ranging from about $-1$ to $4$, we provide a revised SBC relation. For late-type stars ($0 \leq V-K \leq 4$), the results are consistent with previous studies. For early-type stars ($-1 \leq V-K \leq 0$), our new VEGA/CHARA measurements combined with a careful selection of the stars (rejecting stars with environment or stars with a strong variability), allows us to reach an unprecedented precision of about 0.16 magnitude or $\simeq 7\%$ in terms of angular diameter.

The YSOVAR (Young Stellar Object VARiability) Spitzer Space Telescope observing program obtained the first extensive mid-infrared (3.6 & 4.5 um) time-series photometry of the Orion Nebula Cluster plus smaller footprints in eleven other star-forming cores (AFGL490, NGC1333, MonR2, GGD 12-15, NGC2264, L1688, Serpens Main, Serpens South, IRAS 20050+2720, IC1396A, and Ceph C). There are ~29,000 unique objects with light curves in either or both IRAC channels in the YSOVAR data set. We present the data collection and reduction for the Spitzer and ancillary data, and define the "standard sample" on which we calculate statistics, consisting of fast cadence data, with epochs about twice per day for ~40d. We also define a "standard sample of members", consisting of all the IR-selected members and X-ray selected members. We characterize the standard sample in terms of other properties, such as spectral energy distribution shape. We use three mechanisms to identify variables in the fast cadence data--the Stetson index, a chi^2 fit to a flat light curve, and significant periodicity. We also identified variables on the longest timescales possible of ~6 years, by comparing measurements taken early in the Spitzer mission with the mean from our YSOVAR campaign. The fraction of members in each cluster that are variable on these longest timescales is a function of the ratio of Class I/total members in each cluster, such that clusters with a higher fraction of Class I objects also have a higher fraction of long-term variables. For objects with a YSOVAR-determined period and a [3.6]-[8] color, we find that a star with a longer period is more likely than those with shorter periods to have an IR excess. We do not find any evidence for variability that causes [3.6]-[4.5] excesses to appear or vanish within our data; out of members and field objects combined, at most 0.02% may have transient IR excesses.

The ratio of mass and magnetic flux determines the relative importance of magnetic and gravitational forces in the evolution of molecular clouds and their cores. Its measurement is thus central in discriminating between different theories of core formation and evolution. Here we discuss the effect of chemical depletion on measurements of the mass-to-flux ratio using the same molecule (OH) both for Zeeman measurements of the magnetic field and the determination of the mass of the region. The uncertainties entering through the OH abundance in determining separately the magnetic field and the mass of a region have been recognized in the literature. It has been proposed however that, when comparing two regions of the same cloud, the abundance will in both cases be the same. We show that this assumption is invalid. We demonstrate that when comparing regions with different densities, the effect of OH depletion in measuring changes of the mass-to-flux ratio between different parts of the same cloud can even reverse the direction of the underlying trends (for example, the mass-to-flux ratio may appear to decrease as we move to higher density regions). The systematic errors enter primarily through the inadequate estimation of the mass of the region.

We aim at determining the masses of Cepheids in binary systems, as well as their geometric distances and the flux contribution of the companions. The combination of interferometry with spectroscopy will offer a unique and independent estimate of the Cepheid masses. Using long-baseline interferometry at visible and infrared wavelengths, it is possible to spatially resolve binary systems containing a Cepheid down to milliarcsecond separations. Based on the resulting visual orbit and radial velocities, we can then derive the fundamental parameters of these systems, particularly the masses of the components and the geometric distance. We therefore performed interferometric observations of the first-overtone mode Cepheid V1334 Cyg with the CHARA/MIRC combiner. We report the first detection of a Cepheid companion using long-baseline interferometry. We detect the signature of a companion orbiting V1334 Cyg at two epochs. We measure a flux ratio between the companion and the Cepheid f = 3.10+/-0.08%, giving an apparent magnitude mH = 8.47+/-0.15mag. The combination of interferometric and spectroscopic data have enabled the unique determination of the orbital elements: P = 1938.6+/-1.2 days, Tp = 2 443 616.1+/-7.3, a = 8.54+/-0.51mas, i = 124.7+/-1.8{\deg}, e = 0.190+/-0.013, {\omega} = 228.7+/-1.6{\deg}, and {\Omega} = 206.3+/-9.4{\deg}. We derive a minimal distance d ~ 691 pc, a minimum mass for both stars of 3.6 Msol, with a spectral type earlier than B5.5V for the companion star. Our measured flux ratio suggests that radial velocity detection of the companion using spectroscopy is within reach, and would provide an orbital parallax and model-free masses.

When observed with optical long-baseline interferometers (OLBI), components of a binary star which are sufficiently separated produce their own interferometric fringe packets; these are referred to as Separated Fringe Packet (SFP) binaries. These SFP binaries can overlap in angular separation with the regime of systems resolvable by speckle interferometry at single, large-aperture telescopes and can provide additional measurements for preliminary orbits lacking good phase coverage, help constrain elements of already established orbits, and locate new binaries in the undersampled regime between the bounds of spectroscopic surveys and speckle interferometry. In this process, a visibility calibration star is not needed, and the separated fringe packets can provide an accurate vector separation. In this paper, we apply the SFP approach to {\omega} Andromeda, HD 178911, and {\xi} Cephei with the CLIMB three-beam combiner at the CHARA Array. For these systems we determine component masses and parallax of 0.963${\pm}$0.049 $M_{\odot}$ and 0.860${\pm}$0.051 $M_{\odot}$ and 39.54${\pm}$1.85 milliarcseconds (mas) for {\omega} Andromeda, for HD 178911 of 0.802${\pm}$0.055 $M_{\odot}$ and 0.622${\pm}$0.053 $M_{\odot}$ with 28.26${\pm}$1.70 mas, and masses of 1.045${\pm}$0.031 $M_{\odot}$ and 0.408${\pm}$0.066 $M_{\odot}$ and 38.10${\pm}$2.81 mas for {\xi} Cephei.

Two red supergiants of the Per OB1 association, RS Per and T Per, have been observed in H band using the MIRC instrument at the CHARA array. The data show clear evidence of departure from circular symmetry. We present here new techniques specially developed to analyze such cases, based on state-of-the-art statistical frameworks. The stellar surfaces are first modeled as limb-darkened discs based on SATLAS models that fit both MIRC interferometric data and publicly available spectrophotometric data. Bayesian model selection is then used to determine the most probable number of spots. The effective surface temperatures are also determined and give further support to the recently derived hotter temperature scales of red su- pergiants. The stellar surfaces are reconstructed by our model-independent imaging code SQUEEZE, making use of its novel regularizer based on Compressed Sensing theory. We find excellent agreement between the model-selection results and the reconstructions. Our results provide evidence for the presence of near-infrared spots representing about 3-5% of the stellar flux.

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" having 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$\alpha$ 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.

Planets appear to form in environments shaped by the gas flowing through protostellar disks to the central young stars. The flows in turn are governed by orbital angular momentum transfer. In this chapter we summarize current understanding of the transfer processes best able to account for the flows, including magneto-rotational turbulence, magnetically-launched winds, self-gravitational instability and vortices driven by hydrodynamical instabilities. For each in turn we outline the major achievements of the past few years and the outstanding questions. We underscore the requirements for operation, especially ionization for the magnetic processes and heating and cooling for the others. We describe the distribution and strength of the resulting flows and compare with the long-used phenomenological $\alpha$-picture, highlighting issues where the fuller physical picture yields substantially different answers. We also discuss the links between magnetized turbulence and magnetically-launched outflows, and between magnetized turbulence and hydrodynamical vortices. We end with a summary of the status of efforts to detect specific signatures of the flows.

Based on more than four weeks of continuous high cadence photometric monitoring of several hundred members of the young cluster NGC 2264 with two space telescopes, NASA's Spitzer and the CNES CoRoT (Convection, Rotation, and planetary Transits), we provide high quality, multi-wavelength light curves for young stellar objects (YSOs) whose optical variability is dominated by short duration flux bursts, which we infer are due to enhanced mass accretion rates. These light curves show many brief -- several hour to one day -- brightenings at optical and near-infrared (IR) wavelengths with amplitudes generally in the range 5-50% of the quiescent value. Typically, a dozen or more of these bursts occur in a thirty day period. We demonstrate that stars exhibiting this type of variability have large ultraviolet (UV) excesses and dominate the portion of the u-g vs. g-r color-color diagram with the largest UV excesses. These stars also have large Halpha equivalent widths, and either centrally peaked, lumpy Halpha emission profiles or profiles with blue-shifted absorption dips associated with disk or stellar winds. Light curves of this type have been predicted for stars whose accretion is dominated by Rayleigh-Taylor instabilities at the boundary between their magnetosphere and inner circumstellar disk, or where magneto-rotational instabilities modulate the accretion rate from the inner disk. Amongst the stars with the largest UV excesses or largest Halpha equivalent widths, light curves with this type of variability greatly outnumber light curves with relatively smooth sinusoidal variations associated with long-lived hot spots. We provide quantitative statistics for the average duration and strength of the accretion bursts and for the fraction of the accretion luminosity associated with these bursts.

We show that the first order (non co-orbital) corotation torques are significantly modified by entropy gradients in a non-barotropic protoplanetary disk. Such non-barotropic torques can dramatically alter the balance that, for barotropic cases, results in the net eccentricity damping for giant gap-clearing planets embedded in the disk. We demonstrate that stellar illumination can heat the gap enough for the planet's orbital eccentricity to instead be excited. We also discuss the "Eccentricity Valley" noted in the known exoplanet population, where low-metallicity stars have a deficit of eccentric planets between $\sim 0.1$ and $\sim 1$ AU compared to metal-rich systems (Dawson & Murray-Clay 2013). We show that this feature in the planet distribution may be due to the self-shadowing of the disk by a rim located at the dust sublimation radius $\sim 0.1$ AU, which is known to exist for several T Tauri systems. In the shadowed region between $\sim 0.1$ and $\sim 1$ AU lack of gap insolation allows disk interactions to damp eccentricity. Outside such shadowed regions stellar illumination can heat the planetary gaps and drive eccentricity growth for giant planets. We suggest that the self-shadowing does not arise at higher metallicity due to the increased optical depth of the gas interior to the dust sublimation radius.

We examine the conditions under which the disks of gas and dust orbiting young gas giant planets are sufficiently conducting to experience turbulence driven by the magneto-rotational instability. By modeling the ionization and conductivity in the disk around proto-Jupiter, we find that turbulence is possible if the X-rays emitted near the Sun reach the planet's vicinity and either (1) the gas surface densities are in the range of the minimum-mass models constructed by augmenting Jupiter's satellites to Solar composition, while dust is depleted from the disk atmosphere, or (2) the surface densities are much less, and in the range of gas-starved models fed with material from the Solar nebula, but not so low that ambipolar diffusion decouples the neutral gas from the plasma. The results lend support to both minimum-mass and gas-starved models of the protojovian disk: (1) The dusty minimum-mass models have internal conductivities low enough to prevent angular momentum transfer by magnetic forces, as required for the material to remain in place while the satellites form. (2) The gas-starved models have magnetically-active surface layers and a decoupled interior "dead zone". Similar active layers in the Solar nebula yield accretion stresses in the range assumed in constructing the circumjovian gas-starved models. Our results also point to aspects of both classes of models that can be further developed. Non-turbulent minimum-mass models will lose dust from their atmospheres by settling, enabling gas to accrete through a thin surface layer. For the gas-starved models it is crucial to learn whether enough stellar X-ray and ultraviolet photons reach the circumjovian disk. Additionally the stress-to-pressure ratio ought to increase with distance from the planet, likely leading to episodic accretion outbursts.

We use near-infrared interferometric data coupled with trigonometric parallax values and spectral energy distribution fitting to directly determine stellar radii, effective temperatures, and luminosities for the exoplanet host stars 61 Vir, $\rho$ CrB, GJ 176, GJ 614, GJ 649, GJ 876, HD 1461, HD 7924, HD 33564, HD 107383, and HD 210702. Three of these targets are M dwarfs. Statistical uncertainties in the stellar radii and effective temperatures range from 0.5% -- 5% and from 0.2% -- 2%, respectively. For eight of these targets, this work presents the first directly determined values of radius and temperature; for the other three, we provide updates to their properties. The stellar fundamental parameters are used to estimate stellar mass and calculate the location and extent of each system's circumstellar habitable zone. Two of these systems have planets that spend at least parts of their respective orbits in the system habitable zone: two of GJ 876's four planets and the planet that orbits HD 33564. We find that our value for GJ 876's stellar radius is more than 20% larger than previous estimates and frequently used values in the astronomical literature.

Young stars with masses 2-8 Suns, called the Herbig Ae and Be stars, often show a near-infrared excess too large to explain with a hydrostatically-supported circumstellar disk of gas and dust. At the same time the accretion flow carrying the circumstellar gas to the star is thought to be driven by magneto-rotational turbulence, which according to numerical MHD modeling yields an extended low-density atmosphere supported by the magnetic fields. We demonstrate that the base of the atmosphere can be optically-thick to the starlight and that the parts lying near 1 AU are tall enough to double the fraction of the stellar luminosity reprocessed into the near-infrared. We generate synthetic spectral energy distributions (SEDs) using Monte Carlo radiative transfer calculations with opacities for sub-micron silicate and carbonaceous grains. The synthetic SEDs closely follow the median Herbig SED constructed recently by Mulders and Dominik, and in particular match the large near-infrared flux, provided the grains have a mass fraction close to interstellar near the disk's inner rim.

We present global hydrodynamic and magnetohydrodynamic (MHD) simulations with mesh refinement of accreting planets embedded in protoplanetary disks (PPDs). The magnetized disk includes Ohmic resistivity that depends on the overlying mass column, leading to turbulent surface layers and a dead zone near the midplane. The main results are: (i) The accretion flow in the Hill sphere is intrinsically 3D for hydrodynamic and MHD models. Net inflow toward the planet is dominated by high latitude flows. A circumplanetary disk (CPD) forms. Its midplane flows outward in a pattern whose details differ between models. (ii) Gap opening magnetically couples and ignites the dead zone near the planet, leading to stochastic accretion, a quasi-turbulent flow in the Hill sphere and a CPD whose structure displays high levels of variability. (iii) Advection of magnetized gas onto the rotating CPD generates helical fields that launch magnetocentrifugally driven outflows. During one specific epoch a highly collimated, one-sided jet is observed. (iv) The CPD's surface density $\sim30{\rm\,g\,cm^{-2}}$, small enough for significant ionization and turbulence to develop. (v) The accretion rate onto the planet in the MHD simulation reaches a steady value $8 \times 10^{-3} {\rm M_\oplus yr^{-1}}$, and is similar in the viscous hydrodynamic runs. Our results suggest that gas accretion onto a forming giant planet within a magnetized PPD with dead zone allows rapid growth from Saturnian to Jovian masses. As well as being relevant for giant planet formation, these results have important implications for the formation of regular satellites around gas giant planets.

Due to the strong magnetic field and related abnormal surface layers existing in rapidly oscillating Ap stars, systematic errors are likely to be present when determining their effective temperatures, which potentially compromises asteroseismic studies of these pulsators. Using long-baseline interferometry, our goal is to determine accurate angular diameters of a number of roAp targets to provide a temperature calibration for these stars. We obtained interferometric observations of 10 Aql with the visible spectrograph VEGA at the CHARA array. We determined a limb-darkened angular diameter of 0.275+/-0.009 mas and deduced a linear radius of 2.32+/-0.09 R_sun. We estimated the star's bolometric flux and used it, in combination with its parallax and angular diameter, to determine the star's luminosity and effective temperature. For two data sets of bolometric flux we derived an effective temperature of 7800+/-170 K and a luminosity of 18+/-1 L_sun or of 8000+/-210 K and 19+/-2 L_sun. We used these fundamental parameters together with the large frequency separation to constrain the mass and the age of 10 Aql, using the CESAM stellar evolution code. Assuming a solar chemical composition and ignoring all kinds of diffusion and settling of elements, we obtained a mass of 1.92 M_sun and an age of 780 Gy or a mass of 1.95 M_sun and an age of 740 Gy, depending on the considered bolometric flux. For the first time, we managed to determine an accurate angular diameter for a star smaller than 0.3 mas and to derive its fundamental parameters. In particular, by only combining our interferometric data and the bolometric flux, we derived an effective temperature that can be compared to those derived from atmosphere models. Such fundamental parameters can help for testing the mechanism responsible for the excitation of the oscillations observed in the magnetic pulsating stars.

(Abridged) Dust is expected to be ubiquitous in extrasolar planetary systems owing to the dynamical activity of minor bodies. Inner dust populations are, however, still poorly known because of the high contrast and small angular separation with respect to their host star. We aim to determine the level of near-infrared exozodiacal dust emission around a sample of 42 nearby main sequence stars with spectral types ranging from A to K and to investigate its correlation with various stellar parameters and with the presence of cold dust belts. We use high-precision K-band visibilities obtained with the FLUOR interferometer on the shortest baseline of the CHARA array. The calibrated visibilities are compared with the expected visibility of the stellar photosphere to assess whether there is an additional, fully resolved circumstellar emission. Near-infrared circumstellar emission amounting to about 1% of the stellar flux is detected around 13 of our 42 target stars. Follow-up observations showed that one of them (eps Cep) is associated with a stellar companion, while another one was detected around what turned out to be a giant star (kap CrB). The remaining 11 excesses found around single main sequence stars are most probably associated with hot circumstellar dust, yielding an overall occurrence rate of 28+8-6% for our (biased) sample. We show that the occurrence rate of bright exozodiacal discs correlates with spectral type, K-band excesses being more frequent around A-type stars. It also correlates with the presence of detectable far-infrared excess emission in the case of solar-type stars. This study provides new insight into the phenomenon of bright exozodiacal discs, showing that hot dust populations are probably linked to outer dust reservoirs in the case of solar-type stars. For A-type stars, no clear conclusion can be made regarding the origin of the detected near-infrared excesses.

Based on CHARA Array measurements, we present the angular diameters of 23 nearby, main- sequence stars, ranging from spectral type A7 to K0, five of which are exoplanet host stars. We derive linear radii, effective temperatures, and absolute luminosities of the stars using HIPPARCOS parallaxes and measured bolometric fluxes. The new data are combined with previously published values to create an Angular Diameter Anthology of measured angular diameters to main-sequence stars (luminosity class V and IV). This compilation consists of 125 stars with diameter uncertainties of less than 5%, ranging in spectral types from A to M. The large quantity of empirical data are used to derive color-temperature relations to an assortment of color indices in the Johnson (BVRIJHK), Cousins (RI), Kron (RI), Sloan (griz), and WISE (W3W4) photometric systems. These relations have an average standard deviation of ~3% and are valid for stars with spectral types A0 to M4. To derive even more accurate relations for Sun-like stars, we also determined these temperature relations omitting early-type stars (Teff > 6750 K) that may have biased luminosity estimates because of rapid rotation; for this subset the dispersion is only ~2.5%. We find effective temperatures in agreement within a couple percent for the interferometrically characterized sample of main sequence stars compared to those derived via the infrared-flux method and spectroscopic analysis.

We present the fundamental parameters of HR 2582, a high-mass red giant star whose evolutionary state is a mystery. We used the CHARA Array interferometer to directly measure the star's limb-darkened angular diameter (1.006+/-0.020 mas) and combined our measurement with parallax and photometry from the literature to calculate its physical radius (35.76+/-5.31 R_Sun), luminosity (517.8+/-17.5 L_Sun), bolometric flux (14.8+/-0.5 e-8 erg s-1 cm-2) and effective temperature (4577+/-60 K). We then determined the star's mass (5.6+/-1.7 M_Sun) using our new values with stellar oscillation results from Baudin et al. Finally, using the Yonsei-Yale evolutionary models, we estimated HR 2582's age to be 165 +20/-15 Myr. While our measurements do not provide the precision required to definitively state where the star is in its evolution, it remains an excellent test case for evaluating stellar interior models.

Young close binaries clear central cavities in their surrounding circumbinary disk from which the stars can still accrete material. This process takes place within the very first astronomical units, and is still not well constrained as the observational evidence has been gathered, until now, only by means of spectroscopy. The young object HD200775 (MWC361) is a massive spectroscopic binary (separation of ~15.9mas, ~5.0~AU), with uncertain classification (early/late Be), that shows a strong and variable Halpha emission. We aim to study the mechanisms that produce the Halpha line at the AU-scale. Combining the radial velocity measurements and astrometric data available in the literature, we determined new orbital parameters. With the VEGA instrument on the CHARA array, we spatially and spectrally resolved the Halpha emission of HD200775, at low and medium spectral resolutions (R~1600 and 5000) over a full orbital period (~3.6 years). We observe that the Halpha equivalent width varies with the orbital phase, and increases close to periastron, as expected from theoretical models that predict an increase of the mass transfer from the circumbinary disk to the primary disk. In addition, using spectral visibilities and differential phases, we find marginal variations of the typical extent of the Halpha emission (at 1 to 2-sigma level) and location (at 1 to 5-sigma level). The spatial extent of the Halpha emission, as probed by a Gaussian FWHM, is minimum at the ascending node (0.67+/-0.20 mas, i.e., 0.22+/-0.06 AU), and more than doubles at periastron. In addition, the Gaussian photocenter is slightly displaced in the direction opposite to the secondary, ruling out the scenario in which all or most of the Halpha emission is due to accretion onto the secondary. These findings, together with the wide Halpha line profile, may be due to a non-spherical wind enhanced at periastron.

We examine the observational appearance of partial gaps being opened by planets in protoplanetary disks, considering the effects of the inclination relative to the line of sight. The gap's trough is darkened by both shadowing and cooling, relative to the uninterrupted disk. The gap's outer wall is brightened by direct illumination and also by heating, which puffs it up so that it intercepts more starlight. In this paper, we examine the effects of inclination on resolved images of disks with and without gaps at a wide range of wavelengths. The scattering surface's offset from the disk midplane creates a brightness asymmetry along the axis of inclination, making the disk's near side appear brighter than the far side in scattered light. Finite disk thickness also causes the projected distances of equidistant points on the disk surface to be smaller on the near side of the disk as compared to the far side. Consequently, the gap shoulder on the near side of the disk should appear brighter and closer to the star than on the far side. However, if the angular resolution of the observation is coarser than the width of the brightened gap shoulder, then the gap shoulder on the far side may appear brighter because of its larger apparent size. We present a formula to recover the scale height and inclination angle of an imaged disk using simple geometric arguments and measuring disk asymmetries. Resolved images of circumstellar disks have revealed clearings and gaps, such as the transitional disk in LkCa 15. Models created using our synthetic imaging attempting to match the morphology of observed scattered light images of LkCa 15 indicate that the H-band flux deficit in the inner $\sim0.5\arcsec$ of the disk can be explained with a planet of mass greater than 0.5 Jupiter mass.

We present the results of long-baseline optical interferometry observations using the Precision Astronomical Visual Observations (PAVO) beam combiner at the Center for High Angular Resolution Astronomy (CHARA) Array to measure the angular sizes of three bright Kepler stars: {\theta} Cygni, and both components of the binary system 16 Cygni. Supporting infrared observations were made with the Michigan Infrared Combiner (MIRC) and Classic beam combiner, also at the CHARA Array. We find limb-darkened angular diameters of 0.753+/-0.009 mas for {\theta} Cyg, 0.539+/-0.007 mas for 16 Cyg A and 0.490+/-0.006 mas for 16 Cyg B. The Kepler Mission has observed these stars with outstanding photometric precision, revealing the presence of solar-like oscillations. Due to the brightness of these stars the oscillations have exceptional signal-to-noise, allowing for detailed study through asteroseismology, and are well constrained by other observations. We have combined our interferometric diameters with Hipparcos parallaxes, spectrophotometric bolometric fluxes and the asteroseismic large frequency separation to measure linear radii ({\theta} Cyg: 1.48+/-0.02 Rsun, 16 Cyg A: 1.22+/-0.02 Rsun, 16 Cyg B: 1.12+/-0.02 Rsun), effective temperatures ({\theta} Cyg: 6749+/-44 K, 16 Cyg A: 5839+/-42 K, 16 Cyg B: 5809+/-39 K), and masses ({\theta} Cyg: 1.37+/-0.04 Msun, 16 Cyg A: 1.07+/-0.05 Msun, 16 Cyg B: 1.05+/-0.04 Msun) for each star with very little model dependence. The measurements presented here will provide strong constraints for future stellar modelling efforts.

Magnetorotational instability (MRI) is the most promising mechanism behind accretion in low-mass protostellar disks. Here we present the first analysis of the global structure and evolution of non-ideal MRI-driven T-Tauri disks on million-year timescales. We accomplish this in a 1+1D simulation by calculating magnetic diffusivities and utilizing turbulence activity criteria to determine thermal structure and accretion rate without resorting to a 3-D magnetohydrodynamical (MHD) simulation. Our major findings are as follows. First, even for modest surface densities of just a few times the minimum-mass solar nebula, the dead zone encompasses the giant planet-forming region, preserving any compositional gradients. Second, the surface density of the active layer is nearly constant in time at roughly 10 g/cm2, which we use to derive a simple prescription for viscous heating in MRI-active disks for those who wish to avoid detailed MHD computations. Furthermore, unlike a standard disk with constant-alpha viscosity, the disk midplane does not cool off over time, though the surface cools as the star evolves along the Hayashi track. The ice line is firmly in the terrestrial planet-forming region throughout disk evolution and can move either inward or outward with time, depending on whether pileups form near the star. Finally, steady-state mass transport is a poor description of flow through an MRI-active disk. We caution that MRI activity is sensitive to many parameters, including stellar X-ray flux, grain size, gas/small grain mass ratio and magnetic field strength, and we have not performed an exhaustive parameter study here.

We present the first high angular resolution observations in the nearinfrared H-band (1.6 microns) of the Luminous Blue Variable star P Cygni. We obtained six-telescope interferometric observations with the CHARA Array and the MIRC beam combiner. These show that the spatial flux distribution is larger than expected for the stellar photosphere. A two component model for the star (uniform disk) plus a halo (two-dimensional Gaussian) yields an excellent fit of the observations, and we suggest that the halo corresponds to flux emitted from the base of the stellar wind. This wind component contributes about 45% of the H-band flux and has an angular FWHM = 0.96 mas, compared to the predicted stellar diameter of 0.41 mas. We show several images reconstructed from the interferometric visibilities and closure phases, and they indicate a generally spherical geometry for the wind. We also obtained near-infrared spectrophotometry of P Cygni from which we derive the flux excess compared to a purely photospheric spectral energy distribution. The H-band flux excess matches that from the wind flux fraction derived from the two component fits to the interferometry. We find evidence of significant near-infrared flux variability over the period from 2006 to 2010 that appears similar to the variations in the H-alpha emission flux from the wind. Future interferometric observations may be capable of recording the spatial variations associated with temporal changes in the wind structure.

We report on a high angular resolution survey of circumstellar disks around 24 northern sky Be stars. The K-band continuum survey was made using the CHARA Array long baseline interferometer (baselines of 30 to 331 m). The interferometric visibilities were corrected for the flux contribution of stellar companions in those cases where the Be star is a member of a known binary or multiple system. For those targets with good uv coverage, we used a four parameter Gaussian elliptical disk model to fit the visibilities and to determine the axial ratio, position angle, K-band photospheric flux contribution, and angular diameter of the disk major axis. For the other targets with relatively limited uv coverage, we constrained the axial ratio, inclination angle, and or disk position angle where necessary in order to resolve the degeneracy between possible model solutions. We also made fits of the ultraviolet and infrared spectral energy distributions to estimate the stellar angular diameter and infrared flux excess of each target. The mean ratio of the disk diameter (measured in K-band emission) to stellar diameter (from SED modeling) is 4.4 among the 14 cases where we reliably resolved the disk emission, a value which is generally lower than the disk size ratio measured in the higher opacity Halpha emission line. We estimated the equatorial rotational velocity from the projected rotational velocity and disk inclination for 12 stars, and most of these stars rotate close to or at the critical rotational velocity.

The edges of magnetically-dead zones in protostellar disks have been proposed as locations where density bumps may arise, trapping planetesimals and helping form planets. Magneto-rotational turbulence in magnetically-active zones provides both accretion of gas on the star and transport of mass to the dead zone. We investigate the location of the magnetically-active regions in a protostellar disk around a solar-type star, varying the disk temperature, surface density profile, and dust-to-gas ratio. We also consider stellar masses between 0.4 and 2 $M_\odot$, with corresponding adjustments in the disk mass and temperature. The dead zone's size and shape are found using the Elsasser number criterion with conductivities including the contributions from ions, electrons, and charged fractal dust aggregates. The charged species' abundances are found using the approach proposed by S. Okuzumi. The dead zone is in most cases defined by the ambipolar diffusion. In our maps, the dead zone takes a variety of shapes, including a fish-tail pointing away from the star and islands located on and off the midplane. The corresponding accretion rates vary with radius, indicating locations where the surface density will increase over time, and others where it will decrease. We show that density bumps do not readily grow near the dead zone's outer edge, independently of the disk parameters and the dust properties. Instead, the accretion rate peaks at the radius where the gas-phase metals freeze out. This could lead to clearing a valley in the surface density, and to a trap for pebbles located just outside the metal freeze-out line.

Context. V2492 Cyg is a young eruptive star that went into outburst in 2010. The near-infrared color changes observed since the outburst peak suggest that the source belongs to a newly defined sub-class of young eruptive stars, where time-dependent accretion and variable line-of-sight extinction play a combined role in the flux changes. Aims. In order to learn about the origin of the light variations and to explore the circumstellar and interstellar environment of V2492 Cyg, we monitored the source at ten different wavelengths, between 0.55 \mu m and 2.2 \mu m from the ground and between 3.6 \mu m and 160 \mu m from space. Methods. We analyze the light curves and study the color-color diagrams via comparison with the standard reddening path. We examine the structure of the molecular cloud hosting V2492 Cyg by computing temperature and optical depth maps from the far-infrared data. Results. We find that the shapes of the light curves at different wavelengths are strictly self-similar and that the observed variability is related to a single physical process, most likely variable extinction. We suggest that the central source is episodically occulted by a dense dust cloud in the inner disk, and, based on the invariability of the far-infrared fluxes, we propose that it is a long-lived rather than a transient structure. In some respects, V2492 Cyg can be regarded as a young, embedded analog of UX Orionis-type stars. Conclusions. The example of V2492 Cyg demonstrates that the light variations of young eruptive stars are not exclusively related to changing accretion. The variability provided information on an azimuthally asymmetric structural element in the inner disk. Such an asymmetric density distribution in the terrestrial zone may also have consequences for the initial conditions of planet formation.

We investigate the viability of the magnetorotational instability (MRI) in accretion disks around both solar-type stars and very low mass stars. In particular, we determine the disk regions where the MRI can be shut off either by Ohmic resistivity (the so-called Dead and Undead Zones) or by ampipolar diffusion (a region we term the Zombie Zone). We consider 2 stellar masses: Mstar = 0.7 and 0.1 Msun. In each case, we assume that: the disk surface density profile is that of a scaled Minimum Mass Solar Nebula, with Mdisk/Mstar ~ 0.01 as currently estimated; disk ionisation is driven primarily by stellar X-rays, complemented by cosmic rays and radionuclides; and the stellar X-ray luminosity scales with bolometric luminosity as Lx/Lstar ~ 10^-3.5, as observed. Ionization rates are calculated with the MOCCASIN code, and ionisation balance determined using a simplified chemical network, including well-mixed 0.1 um grains at various levels of depletion. We find that (1) ambipolar diffusion is the primary factor controlling MRI activity in disks around both solar-type and very low mass stars. Assuming that the MRI yields the maximum possible field strength at each radius, we further find that: (2) the MRI-active layer constitutes only ~ 5-10% of the total disk mass; (3) the accretion rate (Mdot) varies radially in both magnitude and sign (inward or outward), implying time-variable accretion as well as the creation of disk gaps and overdensities, with consequences for planet formation and migration; (4) achieving the empirical accretion rates in solar-type and very low mass stars requires a depletion of well-mixed small grains by a factor of 10-1000 relative to the standard dust-to-gas mass ratio of 10^-2; and (5) the current non-detection of polarized emission from field-aligned grains in the outer disk regions is consistent with active MRI at those radii.

Optical and infrared interferometers definitively established that the photometric standard Vega (alpha Lyrae) is a rapidly rotating star viewed nearly pole-on. Recent independent spectroscopic analyses could not reconcile the inferred inclination angle with the observed line profiles, preferring a larger inclination. In order to resolve this controversy, we observed Vega using the six-beam Michigan Infrared Combiner on the Center for High Angular Resolution Astronomy Array. With our greater angular resolution and dense (u,v)-coverage, we find Vega is rotating less rapidly and with a smaller gravity darkening coefficient than previous interferometric results. Our models are compatible with low photospheric macroturbulence and also consistent with the possible rotational period of ~0.71 days recently reported based on magnetic field observations. Our updated evolutionary analysis explicitly incorporates rapid rotation, finding Vega to have a mass of 2.15+0.10_-0.15 Msun and an age 700-75+150 Myrs, substantially older than previous estimates with errors dominated by lingering metallicity uncertainties (Z=0.006+0.003-0.002).

HR 8799 is an hF0 mA5 gamma Doradus, lambda Bootis, Vega-type star best known for hosting four directly imaged candidate planetary companions. Using the CHARA Array interferometer, we measure HR 8799's limb-darkened angular diameter to be 0.342 +/- 0.008 mas; this is the smallest interferometrically measured stellar diameter to date, with an error of only 2%. By combining our measurement with the star's parallax and photometry from the literature, we greatly improve upon previous estimates of its fundamental parameters, including stellar radius (1.44 +/- 0.06 R_Sun), effective temperature (7193 +/- 87 K, consistent with F0), luminosity (5.05 +/- 0.29 L_Sun), and the extent of the habitable zone (1.62 AU to 3.32 AU). These improved stellar properties permit much more precise comparisons with stellar evolutionary models, from which a mass and age can be determined, once the metallicity of the star is known. Considering the observational properties of other lambda Bootis stars and the indirect evidence for youth of HR 8799, we argue that the internal abundance, and what we refer to as the effective abundance, is most likely near-solar. Finally, using the Yonsei-Yale evolutionary models with uniformly scaled solar-like abundances, we estimate HR 8799's mass and age considering two possibilities: 1.516 +0.038/-0.024 M_Sun and 33 +7/-13 Myr if the star is contracting toward the zero age main-sequence or 1.513 +0.023/-0.024 M_Sun and 90 +381/-50 Myr if it is expanding from it. This improved estimate of HR 8799's age with realistic uncertainties provides the best constraints to date on the masses of its orbiting companions, and strongly suggests they are indeed planets. They nevertheless all appear to orbit well outside the habitable zone of this young star.

We present results of a long-baseline interferometry campaign using the PAVO beam combiner at the CHARA Array to measure the angular sizes of five main-sequence stars, one subgiant and four red giant stars for which solar-like oscillations have been detected by either Kepler or CoRoT. By combining interferometric angular diameters, Hipparcos parallaxes, asteroseismic densities, bolometric fluxes and high-resolution spectroscopy we derive a full set of near model-independent fundamental properties for the sample. We first use these properties to test asteroseismic scaling relations for the frequency of maximum power (nu_max) and the large frequency separation (Delta_nu). We find excellent agreement within the observational uncertainties, and empirically show that simple estimates of asteroseismic radii for main-sequence stars are accurate to <~4%. We furthermore find good agreement of our measured effective temperatures with spectroscopic and photometric estimates with mean deviations for stars between T_eff = 4600-6200 K of -22+/-32 K (with a scatter of 97K) and -58+/-31 K (with a scatter of 93 K), respectively. Finally we present a first comparison with evolutionary models, and find differences between observed and theoretical properties for the metal-rich main-sequence star HD173701. We conclude that the constraints presented in this study will have strong potential for testing stellar model physics, in particular when combined with detailed modelling of individual oscillation frequencies.

The growing interest in solar twins is motivated by the possibility of comparing them directly to the Sun. To carry on this kind of analysis, we need to know their physical characteristics with precision. Our first objective is to use asteroseismology and interferometry on the brightest of them: 18 Sco. We observed the star during 12 nights with HARPS for seismology and used the PAVO beam-combiner at CHARA for interferometry. An average large frequency separation $134.4\pm0.3$ $\mu$Hz and angular and linear radiuses of $0.6759 \pm 0.0062$ mas and $1.010\pm0.009$ R$_{\odot}$ were estimated. We used these values to derive the mass of the star, $1.02\pm0.03$ M$_{\odot}$.

We present interferometric diameter measurements of 21 K- and M- dwarfs made with the CHARA Array. This sample is enhanced by literature radii measurements to form a data set of 33 K-M dwarfs with diameters measured to better than 5%. For all 33 stars, we compute absolute luminosities, linear radii, and effective temperatures (Teff). We develop empirical relations for \simK0 to M4 main- sequence stars between the stellar Teff, radius, and luminosity to broad-band color indices and metallicity. These relations are valid for metallicities between [Fe/H] = -0.5 and +0.1 dex, and are accurate to ~2%, ~5%, and ~4% for Teff, radius, and luminosity, respectively. Our results show that it is necessary to use metallicity dependent transformations to convert colors into stellar Teffs, radii, and luminosities. We find no sensitivity to metallicity on relations between global stellar properties, e.g., Teff-radius and Teff-luminosity. Robust examinations of single star Teffs and radii compared to evolutionary model predictions on the luminosity-Teff and luminosity-radius planes reveals that models overestimate the Teffs of stars with Teff < 5000 K by ~3%, and underestimate the radii of stars with radii < 0.7 R\odot by ~5%. These conclusions additionally suggest that the models overestimate the effects that the stellar metallicity may have on the astrophysical properties of an object. By comparing the interferometrically measured radii for single stars to those of eclipsing binaries, we find that single and binary star radii are consistent. However, the literature Teffs for binary stars are systematically lower compared to Teffs of single stars by ~ 200 to 300 K. Lastly, we present a empirically determined HR diagram for a total of 74 nearby, main-sequence, A- to M-type stars, and define regions of habitability for the potential existence of sub-stellar mass companions in each system. [abridged]

Studying exoplanet host stars is of the utmost importance to establish the link between the presence of exoplanets around various types of stars and to understand the respective evolution of stars and exoplanets. Using the limb-darkened diameter (LDD) obtained from interferometric data, we determine the fundamental parameters of four exoplanet host stars. We are particularly interested in the F4 main-sequence star, {\theta} Cyg, for which Kepler has recently revealed solar-like oscillations that are unexpected for this type of star. Furthermore, recent photometric and spectroscopic measurements with SOPHIE and ELODIE (OHP) show evidence of a quasi-periodic radial velocity of \sim150 days. Models of this periodic change in radial velocity predict either a complex planetary system orbiting the star, or a new and unidentified stellar pulsation mode. We performed interferometric observations of {\theta} Cyg, 14 Andromedae, {\upsilon} Andromedae and 42 Draconis for two years with VEGA/CHARA (Mount Wilson, California) in several three-telescope configurations. We measured accurate limb darkened diameters and derived their radius, mass and temperature using empirical laws. We obtain new accurate fundamental parameters for stars 14 And, {\upsilon} And and 42 Dra. We also obtained limb darkened diameters with a minimum precision of \sim 1.3%, leading to minimum planet masses of Msini=5.33\pm 0.57, 0.62 \pm 0.09 and 3.79\pm0.29 MJup for 14 And b, {\upsilon} And b and 42 Dra b, respectively. The interferometric measurements of {\theta} Cyg show a significant diameter variability that remains unexplained up to now. We propose that the presence of these discrepancies in the interferometric data is caused by either an intrinsic variation of the star or an unknown close companion orbiting around it.

We have determined the angular diameters of two metal-poor stars, HD 122563 and Gmb 1830, using CHARA and Palomar Testbed Interferometer observations. For the giant star HD 122563, we derive an angular diameter theta_3D = 0.940 +- 0.011 milliarcseconds (mas) using limb-darkening from 3D convection simulations and for the dwarf star Gmb 1830 (HD 103095) we obtain a 1D limb-darkened angular diameter theta_1D = 0.679 +- 0.007 mas. Coupling the angular diameters with photometry yields effective temperatures with precisions better than 55 K (Teff = 4598 +- 41 K and 4818 +- 54 K --- for the giant and the dwarf star, respectively). Including their distances results in very well-determined luminosities and radii (L = 230 +- 6 L_sun, R = 23.9 +- 1.9 R_sun and L = 0.213 +- 0.002 L_sun, R = 0.664 +- 0.015 R_sun, respectively). We used the CESAM2k stellar structure and evolution code in order to produce models that fit the observational data. We found values of the mixing-length parameter alpha (which describes 1D convection) that depend on the mass of the star. The masses were determined from the models with precisions of <3% and with the well-measured radii excellent constraints on the surface gravity are obtained (log g = 1.60 +- 0.04, 4.59 +- 0.02, respectively). The very small errors on both log g and Teff provide stringent constraints for spectroscopic analyses given the sensitivity of abundances to both of these values. The precise determination of Teff for the two stars brings into question the photometric scales for metal-poor stars.

We study the effect that non-equilibrium chemistry in dynamical models of collapsing molecular cloud cores has on measurements of the magnetic field in these cores, the degree of ionization, and the mean molecular weight of ions. We find that OH and CN, usually used in Zeeman observations of the line-of-sight magnetic field, have an abundance that decreases toward the center of the core much faster than the density increases. As a result, Zeeman observations tend to sample the outer layers of the core and consistently underestimate the core magnetic field. The degree of ionization follows a complicated dependence on the number density at central densities up to 10^5 cm^{-3} for magnetic models and 10^6 cm^{-3} in non-magnetic models. At higher central densities the scaling approaches a power-law with a slope of -0.6 and a normalization which depends on the cosmic-ray ionization rate {\zeta} and the temperature T as ({\zeta}T)^1/2. The mean molecular weight of ions is systematically lower than the usually assumed value of 20 - 30, and, at high densities, approaches a value of 3 due to the asymptotic dominance of the H3+ ion. This significantly lower value implies that ambipolar diffusion operates faster.

The late-type dwarf GJ 436 is known to host a transiting Neptune-mass planet in a 2.6-day orbit. We present results of our interferometric measurements to directly determine the stellar diameter ($R_{\star} = 0.455 \pm 0.018 R_{\odot}$) and effective temperature ($T_{\rm EFF} = 3416 \pm 54$ K). We combine our stellar parameters with literature time-series data, which allows us to calculate physical and orbital system parameters, including GJ 436's stellar mass ($M_{\star} = 0.507^{+ 0.071}_{- 0.062} M_{\odot}$) and density ($\rho_* = 5.37^{+ 0.30}_{- 0.27} \rho_\odot$), planetary radius ($R_{p} = 0.369^{+ 0.015}_{- 0.015} R_{Jupiter}$), planetary mass ($M_{p} = 0.078^{+ 0.007}_{- 0.008} M_{Jupiter}$), implying a mean planetary density of $\rho_{p} = 1.55^{+ 0.12}_{- 0.10} \rho_{Jupiter}$. These values are generally in good agreement with previous literature estimates based on assumed stellar mass and photometric light curve fitting. Finally, we examine the expected phase curves of the hot Neptune GJ 436b, based on various assumptions concerning the efficiency of energy redistribution in the planetary atmosphere, and find that it could be constrained with {\it Spitzer} monitoring observations.

Algol (Beta Per) is an extensively studied hierarchical triple system whose inner pair is a prototype semi-detached binary with mass transfer occurring from the sub-giant secondary to the main-sequence primary. We present here the results of our Algol observations made between 2006 and 2010 at the CHARA interferometer with the Michigan Infrared Combiner in the H band. The use of four telescopes with long baselines allows us to achieve better than 0.5 mas resolution and to unambiguously resolve the three stars. The inner and outer orbital elements, as well as the angular sizes and mass ratios for the three components are determined independently from previous studies. We report a significantly improved orbit for the inner stellar pair with the consequence of a 15% change in the primary mass compared to previous studies. We also determine the mutual inclination of the orbits to be much closer to perpendicularity than previously established. State-of-the-art image reconstruction algorithms are used to image the full triple system. In particular an image sequence of 55 distinct phases of the inner pair orbit is reconstructed, clearly showing the Roche-lobe-filling secondary revolving around the primary, with several epochs corresponding to the primary and secondary eclipses.

We combine dynamical and non-equilibrium chemical modeling of evolving prestellar molecular cloud cores, and explore the evolution of molecular abundances in the contracting core. We model both magnetic cores, with varying degrees of initial magnetic support, and non-magnetic cores, with varying collapse delay times. We explore, through a parameter study, the competing effects of various model parameters in the evolving molecular abundances, including the elemental C/O ratio, the temperature, and the cosmic-ray ionization rate. We find that different models show their largest quantitative differences at the center of the core, whereas the outer layers, which evolve slower, have abundances which are severely degenerate among different dynamical models. There is a large range of possible abundance values for different models at a fixed evolutionary stage (central density), which demonstrates the large potential of chemical differentiation in prestellar cores. However, degeneracies among different models, compounded with uncertainties induced by other model parameters, make it difficult to discriminate among dynamical models. To address these difficulties, we identify abundance ratios between particular molecules, the measurement of which would have maximal potential for discrimination among the different models examined here. In particular, we find that the ratios between NH3 and CO; NH2 and CO; NH3 and HCO+ are sensitive to the evolutionary timescale, and that the ratio between HCN and OH is sensitive to the C/O ratio. Finally, we demonstrate that measurements of the central deviation (central depletion or enhancement) of abundances of certain molecules are good indicators of the dynamics of the core.

Optical and near-infrared variability is a well-known property of young stellar objects. However, a growing number of recent studies claim that a considerable fraction of them also exhibit mid-infrared flux changes. With the aim of studying and interpreting variability on a decadal timescale, here we present a mid-infrared spectral atlas containing observations of 68 low- and intermediate mass young stellar objects. The atlas consists of 2.5-11.6 um low-resolution spectra obtained with the ISOPHOT-S instrument on-board the Infrared Space Observatory (ISO) between 1996 and 1998, as well as 5.2-14.5 um low-resolution spectra obtained with the IRS instrument on-board the Spitzer Space Telescope between 2004 and 2007. The observations were retrieved from the ISO and Spitzer archives and were post-processed interactively by our own routines. For those 47 objects where multi-epoch spectra were available, we analyze mid-infrared spectral variability on annual and/or decadal timescales. We identify 37 variable candidate sources. Many stars show wavelength-independent flux changes, possibly due to variable accretion rate. In several systems, all exhibiting 10 um silicate emission, the variability of the 6-8 um continuum and the silicate feature exhibit different amplitudes. A possible explanation is variable shadowing of the silicate emitting region by an inner disk structure of changing height or extra silicate emission from dust clouds in the disk atmosphere. Our results suggest that mid-infrared variability, in particular the wavelength-dependent changes, are more ubiquitous than was known before. Interpreting this variability is a new possibility to explore the structure of the disk and its dynamical processes.

We report the first angular diameter measurements of two classical Cepheids, FF Aql and T Vul, that we have obtained with the FLUOR instrument installed at the CHARA interferometric array. We obtain average limb-darkened angular diameters of \theta_LD = 0.878 +/- 0.013 mas and \theta_LD = 0.629 +/- 0.013 mas, respectively for FF Aql and T Vul. Combining these angular diameters with the HST-FGS trigonometric parallaxes leads to linear radii R = 33.6 +/- 2.2 Rsol and R = 35.6 +/- 4.4 Rsol, respectively. The comparison with empirical and theoretical Period-Radius relations leads to the conclusion that these Cepheids are pulsating in their fundamental mode. The knowledge of the pulsation mode is of prime importance to calibrate the Period-Luminosity relation with a uniform sample of fundamental mode Cepheids.

We present long-baseline interferometry of the Kepler exoplanet host star HD179070 (Kepler-21) using the PAVO beam combiner at the CHARA Array. The visibility data are consistent with a single star and exclude stellar companions at separations ~1-1000 mas (~ 0.1-113 AU) and contrasts < 3.5 magnitudes. This result supports the validation of the 1.6 R_{earth} exoplanet Kepler-21b by Howell et al. (2012) and complements the constraints set by adaptive optics imaging, speckle interferometry, and radial velocity observations to rule out false-positives due to stellar companions. We conclude that long-baseline interferometry has strong potential to validate transiting extrasolar planets, particularly for future projects aimed at brighter stars and for host stars where radial velocity follow-up is not available.

We examine the observational consequences of partial gaps being opened by planets in protoplanetary disks. We model the disk using a static alpha-disk model with detailed radiative transfer, parametrizing the shape and size of the partially cleared gaps based on the results of hydrodynamic simulations. Shadowing and illumination by stellar irradiation at the surface of the gap leads to increased contrast as the gap trough is deepened by shadowing and cooling and the far gap wall is puffed up by illumination and heating. In calculating observables, we find that multiple scattering is important and derive an approximation to include these effects. A gap produced by a 200 M_Earth (70 M_Earth) planet at 10 AU can lower/raise the midplane temperature of the disk by up to ~-25/+29% (~-11/+19%) by shadowing in the gap trough and illumination on the far shoulder of the gap. At the distance of Taurus, this gap would be resolvable with ~0.01" angular resolution. The gap contrast is most significant in scattered light and at thermal continuum wavelengths characteristic of the surface temperature, reducing or raising the surface brightness by up to order of magnitude. Since gaps sizes are correlated to planet mass, this is a promising way of finding and determining the masses of planets embedded in protoplanetary disks.

Hierarchical triple systems comprise a close binary and a more distant component. They are important for testing theories of star formation and of stellar evolution in the presence of nearby companions. We obtained 218 days of Kepler photometry of HD 181068 (magnitude of 7.1), supplemented by groundbased spectroscopy and interferometry, which show it to be a hierarchical triple with two types of mutual eclipses. The primary is a red giant that is in a 45-day orbit with a pair of red dwarfs in a close 0.9-day orbit. The red giant shows evidence for tidally-induced oscillations that are driven by the orbital motion of the close pair. HD 181068 is an ideal target for studies of dynamical evolution and testing tidal friction theories in hierarchical triple systems.

(Abridged) Planetesimals embedded in a protoplanetary disc are stirred by gravitational torques exerted by density fluctuations in the surrounding turbulence. In particular, planetesimals in a disc supporting fully developed magneto-rotational turbulence are readily excited to velocity dispersions above the threshold for catastrophic disruption, halting planet formation. We aim to examine the stirring of planetesimals lying instead in a magnetically-decoupled midplane dead zone, stirred only by spiral density waves propagating out of the disc's magnetically-coupled turbulent surface layers. We extend previous studies to include a wider range of disc models, and explore the effects of varying the disc column density and external magnetic field strength. [...] The strength of the stirring is found to be independent of the gas surface density, which is contrary to the increase with disc mass expected from a simple linear wave picture. The discrepancy arises from the shearing out of density waves as they propagate into the dead zone, resulting in density structures near the midplane that exert weaker stochastic torques on average. We provide a simple analytic fit to our numerically obtained torque amplitudes that accounts for this effect. The stirring on the other hand depends sensitively on the net vertical magnetic flux, up to a saturation level above which magnetic forces dominate in the turbulent layers. For the majority of our models, the equilibrium planetesimal velocity dispersions lie between the thresholds for disrupting strong and weak aggregates, suggesting that collision outcomes will depend on material properties. However, discs with relatively weak magnetic fields yield reduced stirring, and their dead zones provide safe-havens even for the weakest planetesimals against collisional destruction.

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

The Sun was an order of magnitude more luminous during the first few hundred thousand years of its existence, due in part to the gravitational energy released by material accreting from the Solar nebula. If Jupiter was already near its present mass, the planet's tides opened an optically-thin gap in the nebula. We show using Monte Carlo radiative transfer calculations that sunlight absorbed by the nebula and re-radiated into the gap raised temperatures well above the sublimation threshold for water ice, with potentially drastic consequences for the icy bodies in Jupiter's feeding zone. Bodies up to a meter in size were vaporized within a single orbit if the planet was near its present location during this early epoch. Dust particles lost their ice mantles, and planetesimals were partially to fully devolatilized, depending on their size. Scenarios in which Jupiter formed promptly, such as those involving a gravitational instability of the massive early nebula, must cope with the high temperatures. Enriching Jupiter in the noble gases through delivery trapped in clathrate hydrates will be more difficult, but might be achieved by either forming the planet much further from the star, or capturing planetesimals at later epochs. The hot gap resulting from an early origin for Jupiter also would affect the surface compositions of any primordial Trojan asteroids.

We have executed a survey of nearby, main sequence A, F, and G-type stars with the CHARA Array, successfully measuring the angular diameters of fortyfour stars with an average precision of ~ 1.5%. We present new measures of the bolometric flux, which in turn leads to an empirical determination of the effective temperature for the stars observed. In addition, these CHARA-determined temperatures, radii, and luminosities are fit to Yonsei-Yale model isochrones to constrain the masses and ages of the stars. These results are compared to indirect estimates of these quantities obtained by collecting photometry of the stars and applying them to model atmospheres and evolutionary isochrones. We find that for most cases, the models overestimate the effective temperature by ~ 1.5-4%, when compared to our directly measured values. The overestimated temperatures and underestimated radii in these works appear to cause an additional offset in the star's surface gravity measurements, which consequently yield higher masses and younger ages, in particular for stars with masses greater than ~1.3 Msun. Additionally, we compare our measurements to a large sample of eclipsing binary stars, and excellent agreement is seen within both data sets. Finally, we present temperature relations with respect to (B-V) and (V-K) color as well as spectral type showing that calibration of effective temperatures with errors ~ 1% is now possible from interferometric angular diameters of stars.

Wolf-Rayet stars represent one of the final stages of massive stellar evolution. Relatively little is known about this short-lived phase and we currently lack reliable mass, distance, and binarity determinations for a representative sample. Here we report the first visual orbit for WR 140(=HD193793), a WC7+O5 binary system known for its periodic dust production episodes triggered by intense colliding winds near periastron passage. The IOTA and CHARA interferometers resolved the pair of stars in each year from 2003--2009, covering most of the highly-eccentric, 7.9 year orbit. Combining our results with the recent improved double-line spectroscopic orbit of Fahed et al. (2011), we find the WR 140 system is located at a distance of 1.67 +/- 0.03 kpc, composed of a WR star with M_WR = 14.9 +/- 0.5 Msun and an O star with M_O = 35.9 +/- 1.3 Msun. Our precision orbit yields key parameters with uncertainties times 6 smaller than previous work and paves the way for detailed modeling of the system. Our newly measured flux ratios at the near-infrared H and Ks bands allow an SED decomposition and analysis of the component evolutionary states.

We investigate the significance of large scale azimuthal, magnetic and velocity modes for the MRI turbulence in accretion disks. We perform 3D global ideal MHD simulations of global stratified proto-planetary disk models. Our domains span azimuthal angles of \pi/4, \pi/2, \pi and 2\pi. We observe up to 100% stronger magnetic fields and stronger turbulence for the restricted azimuthal domain models \pi/2 and \pi/4 compared to the full 2\pi model. We show that for those models, the Maxwell Stress is larger due to strong axisymmetric magnetic fields, generated by the \alpha \Omega dynamo. Large radial extended axisymmetric toroidal fields trigger temporal magnification of accretion stress. All models display a positive dynamo-\alpha in the northern hemisphere (upper disk). The parity is distinct in each model and changes on timescales of 40 local orbits. In model 2\pi, the toroidal field is mostly antisymmetric in respect to the midplane. The eddies of the MRI turbulence are highly anisotropic. The major wavelengths of the turbulent velocity and magnetic fields are between one and two disk scale heights. At the midplane, we find magnetic tilt angles around 8-9 degree increasing up to 12-13 degree in the corona. We conclude that an azimuthal extent of \pi is sufficient to reproduce most turbulent properties in 3D global stratified simulations of magnetised accretion disks.

Context. The interpretation of stellar pulsations in terms of internal structure depends on the knowledge of the fundamental stellar parameters. Long-base interferometers permit us to determine very accurate stellar radii, which are independent constraints for stellar models that help us to locate the star in the HR diagram. Aims: Using a direct interferometric determination of the angular diameter and advanced three-dimensional (3D) modeling, we derive the radius of the CoRoT target HD 49933 and reduce the global stellar parameter space compatible with seismic data. Methods: The VEGA/CHARA spectro-interferometer is used to measure the angular diameter of the star. A 3D radiative hydrodynamical simulation of the surface is performed to compute the limb darkening and derive a reliable diameter from visibility curves. The other fundamental stellar parameters (mass, age, and Teff) are found by fitting the large and small p-mode frequency separations using a stellar evolution model that includes microscopic diffusion. Results: We obtain a limb-darkened angular diameter of {\theta}LD = 0.445 \pm 0.012 mas. With the Hipparcos parallax, we obtain a radius of R = 1.42 \pm 0.04 Rsun. The corresponding stellar evolution model that fits both large and small frequency separations has a mass of 1.20 \pm 0.08 Msun and an age of 2.7 Gy. The atmospheric parameters are Teff = 6640 \pm 100 K, log g = 4.21 \pm 0.14, and [Fe/H] = -0.38.

We aim to constrain the geometrical extent of the chromosphere of non-binary K giant stars and detect any spatial structures in the chromosphere. We performed observations with the CHARA interferometer and the VEGA beam combiner at optical wavelengths. We observed seven non-binary K giant stars. We measured the ratio of the radii of the photosphere to the chromosphere using the interferometric measurements in the Halpha and the Ca II infrared triplet line cores. For beta Ceti, spectro-interferometric observations are compared to an non-local thermal equilibrium (NLTE) semi-empirical model atmosphere including a chromosphere. The NLTE computations provide line intensities and contribution functions that indicate the relative locations where the line cores are formed and can constrain the size of the limb-darkened disk of the stars with chromospheres. We measured the angular diameter of seven K giant stars and deduced their fundamental parameters: effective temperatures, radii, luminosities, and masses. We determined the geometrical extent of the chromosphere for four giant stars. The chromosphere extents obtained range between 16% to 47% of the stellar radius. The NLTE computations confirm that the Ca II/849 nm line core is deeper in the chromosphere of ? Cet than either of the Ca II/854 nm and Ca II/866 nm line cores. We present a modified version of a semi-empirical model atmosphere derived by fitting the Ca II triplet line cores of this star. In four of our targets, we also detect the signature of a differential signal showing the presence of asymmetries in the chromospheres. Conclusions. It is the first time that geometrical extents and structure in the chromospheres of non-binary K giant stars are determined by interferometry. These observations provide strong constrains on stellar atmosphere models.

Direct detection of thermal emission from nearby hot Jupiters has greatly advanced our knowledge of extrasolar planets in recent years. Since hot Jupiter systems can be regarded as analogs of high contrast binaries, ground-based infrared long baseline interferometers have the potential to resolve them and detect their thermal emission with precision closure phase - a method that is immune to the systematic errors induced by the Earth's atmosphere. In this work, we present closure phase studies toward direct detection of nearby hot Jupiters using the CHARA interferometer array outfitted with the MIRC instrument. We carry out closure phase simulations and conduct a large number of observations for the best candidate {\upsion} And. Our experiments suggest the method is feasible with highly stable and precise closure phases. However, we also find much larger systematic errors than expected in the observations, most likely caused by dispersion across different wavelengths. We find that using higher spectral resolution modes (e.g., R=150) can significantly reduce the systematics. By combining all calibrators in an observing run together, we are able to roughly recalibrate the lower spectral resolution data, allowing us to obtain upper limits of the star-planet contrast ratios of {\upsion} And b across the H band. The data also allow us to get a refined stellar radius of 1.625\pm0.011 R\odot. Our best upper limit corresponds to a contrast ratio of 2.1\times10^3:1 with 90% confidence level at 1.52{\mu}m, suggesting that we are starting to have the capability of constraining atmospheric models of hot Jupiters with interferometry. With recent and upcoming improvements of CHARA/MIRC, the prospect of detecting emission from hot Jupiters with closure phases is promising.

We measured the angular diameter of the exoplanet host star iota Dra with Georgia State University's Center for High Angular Resolution Astronomy (CHARA) Array interferometer, and, using the star's parallax and photometry from the literature, calculated its physical radius and effective temperature. We then combined our results with stellar oscillation frequencies from Zechmeister et al. (2008) and orbital elements from Kane et al. (2010) to determine the masses for the star and exoplanet. Our value for the central star's mass is 1.82 +/- 0.23 M_Sun, which means the exoplanet's minimum mass is 12.6 +/- 1.1 M_Jupiter. Using our new effective temperature, we recalculated the habitable zone for the system, though it is well outside the star-planet separation.

A study of the host stars to exoplanets is important to understanding their environment. To that end, we report new speckle observations of a sample of exoplanet host primaries. The bright exoplanet host HD 8673 (= HIP 6702) is revealed to have a companion, although at this time we cannot definitively establish the companion as physical or optical. The observing lists for planet searches and for these observations have for the most part been pre-screened for known duplicity, so the detected binary fraction is lower than what would otherwise be expected. Therefore, a large number of double stars were observed contemporaneously for verification and quality control purposes, to ensure the lack of detection of companions for exoplanet hosts was valid. In these additional observations, ten pairs are resolved for the first time and sixty pairs are confirmed. These observations were obtained with the USNO speckle camera on the NOAO 4m telescopes at both KPNO and CTIO from 2001 to 2010.

The results of an adaptive optics survey of exoplanet host stars for stellar companions is presented. We used the AEOS Telescope and its adaptive optics system to collect deep images of the stars in $I$-band. Sixty-two exoplanet host stars were observed and fifteen multiple star systems were resolved. Of these eight are known multiples, while seven are new candidate binaries. For all binaries, we measured the relative astrometry of the pair and the differential magnitude in $I$-band. We improved the orbits of HD 19994 and $\tau$ Boo. These observations will provide improved statistics on the duplicity of exoplanet hosts stars and provide an increased understanding of the dynamics of known binary star exoplanet hosts.

Using CHARA and VLTI near-infrared spectro-interferometry with hectometric baseline lengths (up to 330m) and with high spectral resolution (up to 12000), we studied the gas distribution and kinematics around two classical Be stars. The combination of high spatial and spectral resolution achieved allows us to constrain the gas velocity field on scales of a few stellar radii and to obtain, for the first time in optical interferometry, a dynamical mass estimate using the position-velocity analysis technique known from radio astronomy. For our first target star, Beta Canis Minoris, we model the H+K-band continuum and Br-gamma line geometry with a near-critical rotating stellar photosphere and a geometrically thin equatorial disk. Testing different disk rotation laws, we find that the disk is in Keplerian rotation (v(r)=r^(-0.5+/-0.1)) and derive the disk position angle (140+/-1.7 deg), inclination (38.5+/-1 deg), and the mass of the central star (3.5+/-0.2 M_sun). As a second target star, we observed the prototypical Be star Zeta Tauri and spatially resolved the Br{\gamma} emission as well as nine transitions from the hydrogen Pfund series (Pf14-Pf22). Comparing the spatial origin of the different line transitions, we find that the Brackett (Br-gamma), Pfund (Pf14-17), and Balmer (H-alpha) lines originate from different stellocentric radii (R_cont < R_Pf < R_Br-gamma = R_H-alpha), which we can reproduce with an LTE line radiative transfer computation. Discussing different disk-formation scenarios, we conclude that our constraints are inconsistent with wind compression models predicting a strong outflowing velocity component, but support viscous decretion disk models, where the Keplerian-rotating disk is replenished with material from the near-critical rotating star.

The bright star 55 Cancri is known to host five planets, including a transiting super-Earth. The study presented here yields directly determined values for 55 Cnc's stellar astrophysical parameters based on improved interferometry: $R=0.943 \pm 0.010 R_{\odot}$, $T_{\rm EFF} = 5196 \pm 24$ K. We use isochrone fitting to determine 55 Cnc's age to be 10.2 $\pm$ 2.5 Gyr, implying a stellar mass of $0.905 \pm 0.015 M_{\odot}$. Our analysis of the location and extent of the system's habitable zone (0.67--1.32 AU) shows that planet f, with period $\sim$ 260 days and $M \sin i = 0.155 M_{Jupiter}$, spends the majority of the duration of its elliptical orbit in the circumstellar habitable zone. Though planet f is too massive to harbor liquid water on any planetary surface, we elaborate on the potential of alternative low-mass objects in planet f's vicinity: a large moon, and a low-mass planet on a dynamically stable orbit within the habitable zone. Finally, our direct value for 55 Cancri's stellar radius allows for a model-independent calculation of the physical diameter of the transiting super-Earth 55 Cnc e ($\sim 2.05 \pm 0.15 R_{\earth}$), which, depending on the planetary mass assumed, implies a bulk density of 0.76 $\rho_{\earth}$ or 1.07 $\rho_{\earth}$.

The CHARA Array is a six 1-m telescope optical and near infrared interferometer located at the Mount Wilson Observatory in southern California and operated by the Center for High Angular Resolution Astronomy of Georgia State University. The CHARA Array has been in regular scientific operation since 2005 and now has over 55 publications in the refereed literature, including two in Science and one in Nature. The Array now supports seven beam combiners ranging from 0.5 microns up to 2.3 microns and combing from 2 to 4 beams at a time. An upgrade to a full 6 beam combiner is now underway and fringes with all six telescopes were achieved soon after the meeting.We present some of the more recent results from the CHARA-Array.

The bright star 55 Cancri is known to host five planets, including a transiting super-Earth. We use the CHARA Array to directly determine the following of 55 Cnc's stellar astrophysical parameters: $R=0.943 \pm 0.010 R_{\odot}$, $T_{\rm EFF} = 5196 \pm 24$ K. Planet 55 Cnc f ($M \sin i = 0.155 M_{Jupiter}$) spends the majority of the duration of its elliptical orbit in the circumstellar habitable zone (0.67--1.32 AU) where, with moderate greenhouse heating, it could harbor liquid water. Our determination of 55 Cancri's stellar radius allows for a model-independent calculation of the physical diameter of the transiting super-Earth 55 Cnc e ($\simeq 2.1 R_{\earth}$), which, depending on the assumed literature value of planetary mass, implies a bulk density of 0.76 $\rho_{\earth}$ or 1.07 $\rho_{\earth}$.

(abridged) Accretion in protoplanetary discs is thought to be driven by [...] turbulence via the magnetorotational instability (MRI). Recent work has shown that a planetesimal swarm embedded in a fully turbulent disc is subject to strong excitation of the velocity dispersion, leading to collisional destruction of bodies with radii R_p < 100 km. Significant diffusion of planetesimal semimajor axes also arises, leading to large-scale spreading of the planetesimal population throughout the inner regions of the protoplanetary disc, in apparent contradiction of constraints provided by the distribution of asteroids within the asteroid belt. In this paper, we examine the dynamics of planetesimals embedded in vertically stratified turbulent discs, with and without dead zones. Our main aims are to examine the turbulent excitation of the velocity dispersion, and the radial diffusion, of planetesimals in these discs. We employ three dimensional MHD simulations [...], along with an equilibrium chemistry model [...] We find that planetesimals in fully turbulent discs develop large random velocities that will lead to collisional destruction/erosion for bodies with sizes below 100 km, and undergo radial diffusion on a scale \sim 2.5 au over a 5 Myr disc life time. But planetesimals in a dead zone experience a much reduced excitation of their random velocities, and equilibrium velocity dispersions lie between the disruption thresholds for weak and strong aggregates for sizes R_p < 100 km. We also find that radial diffusion occurs over a much reduced length scale \sim 0.25 au over the disc life time, this being consistent with solar system constraints. We conclude that planetesimal growth via mutual collisions between smaller bodies cannot occur in a fully turbulent disc. By contrast, a dead zone may provide a safe haven in which km-sized planetesimals can avoid mutual destruction through collisions.

Classical Be stars are hot non-supergiant stars surrounded by a gaseous circumstellar disk that is responsible for the observed IR-excess and emission lines. The influence of binarity on these phenomena remains controversial. delta Sco is a binary system whose primary suddently began to exhibit the Be phenomenon at the last periastron in 2000. We want to constrain the geometry and kinematics of its circumstellar environment. We observed the star between 2007 and 2010 using spectrally-resolved interferometry with the VLTI/AMBER and CHARA/VEGA instruments. We found orbital elements that are compatible with previous estimates. The next periastron should take place around July 5, 2011 (+- 4,days). We resolved the circumstellar disk in the HAlpha (FWHM = 4.8+-1.5mas), BrGamma (FWHM = 2.9 0.,mas), and the 2.06$ \mu$m HeI (FWHM = 2.4+-0.3mas) lines as well as in the K band continuum (FWHM ~2.4mas). The disk kinematics are dominated by the rotation, with a disk expansion velocity on the order of 0.2km/s. The rotation law within the disk is compatible with Keplerian rotation. As the star probably rotates at about 70% of its critical velocity the ejection of matter doesn't seems to be dominated by rotation. However, the disk geometry and kinematics are similar to that of the previously studied quasi-critically rotating Be stars, namely Alpha Ara, Psi Per and 48 Per.

Near-infrared interferometers have recently imaged a number of rapidly rotating A-type stars, finding levels of gravity darkening inconsistent with theoretical expectations. Here, we present new imaging of both a cooler star {\beta} Cas (F2IV) and a hotter one {\alpha} Leo (B7V) using the CHARA array and the MIRC instrument at the H band. Adopting a solid-body rotation model with a simple gravity darkening prescription, we modeled the stellar geometric properties and surface temperature distributions, confirming both stars are rapidly rotating and show gravity darkening anomalies. We estimate the masses and ages of these rapid rotators on L-Rpol and HR diagrams constructed for non-rotating stars by tracking their non-rotating equivalents. The unexpected fast rotation of the evolved subgiant {\beta} Cas offers a unique test of the stellar core-envelope coupling, revealing quite efficient coupling over the past ~ 0.5 Gyr. Lastly we summarize all our interferometric determinations of the gravity darkening coefficient for rapid rotators, finding none match the expectations from the widely used von Zeipel gravity darkening laws. Since the conditions of the von Zeipel law are known to be violated for rapidly rotating stars, we recommend using the empirically-derived {\beta} = 0.19 for such stars with radiation-dominated envelopes. Furthermore, we note that no paradigm exists for self-consistently modeling heavily gravity-darkened stars that show hot radiative poles with cool convective equators.

The star HD 193322 is a remarkable multiple system of massive stars that lies at the heart of the cluster Collinder 419. Here we report on new spectroscopic observations and radial velocities of the narrow-lined component Ab1 that we use to determine its orbital motion around a close companion Ab2 ($P = 312$ d) and around a distant third star Aa ($P = 35$ y).We have also obtained long baseline interferometry of the target in the $K^\prime$-band with the CHARA Array that we use in two ways. First, we combine published speckle interferometric measurements with CHARA separated fringe packet measurements to improve the visual orbit for the wide Aa,Ab binary. Second, we use measurements of the fringe packet from Aa to calibrate the visibility of the fringes of the Ab1,Ab2 binary, and we analyze these fringe visibilities to determine the visual orbit of the close system. The two most massive stars, Aa and Ab1, have masses of approximately 21 and $23 M_\odot$, respectively, and their spectral line broadening indicates that they represent extremes of fast and slow projected rotational velocity, respectively.

We present full 2 Pi global 3-D stratified MHD simulations of accretion disks. We interpret our results in the context of proto-planetary disks. We investigate the turbulence driven by the magneto-rotational instability (MRI) using the PLUTO Godunov code in spherical coordinates with the accurate and robust HLLD Riemann solver. We follow the turbulence for more than 1500 orbits at the innermost radius of the domain to measure the overall strength of turbulent motions and the detailed accretion flow pattern. We find that regions within two scale heights of the midplane have a turbulent Mach number of about 0.1 and a magnetic pressure two to three orders of magnitude less than the gas pressure, while outside three scale heights the magnetic pressure equals or exceeds the gas pressure and the turbulence is transonic, leading to large density fluctuations. The strongest large-scale density disturbances are spiral density waves, and the strongest of these waves has m=5. No clear meridional circulation appears in the calculations because fluctuating radial pressure gradients lead to changes in the orbital frequency, comparable in importance to the stress gradients that drive the meridional flows in viscous models. The net mass flow rate is well-reproduced by a viscous model using the mean stress distribution taken from the MHD calculation. The strength of the mean turbulent magnetic field is inversely proportional to the radius, so the fields are approximately force-free on the largest scales. Consequently the accretion stress falls off as the inverse square of the radius.

We examine heating and cooling in protostellar disks using 3-D radiation-MHD calculations of a patch of the Solar nebula at 1 AU, employing the shearing-box and flux-limited radiation diffusion approximations. The disk atmosphere is ionized by stellar X-rays, well-coupled to magnetic fields, and sustains a turbulent accretion flow driven by magneto-rotational instability, while the interior is resistive and magnetically dead. The turbulent layers heat by absorbing the light from the central star and by dissipating the magnetic fields. They are optically-thin to their own radiation and cool inefficiently. The optically-thick interior in contrast is heated only weakly, by re-emission from the atmosphere. The interior is colder than a classical viscous model, and isothermal. The magnetic fields support an extended atmosphere that absorbs the starlight 1.5 times higher than the hydrostatic viscous model. The disk thickness thus measures not the internal temperature, but the magnetic field strength. Fluctuations in the fields move the starlight-absorbing surface up and down. The height ranges between 13% and 24% of the radius over timescales of several orbits, with implications for infrared variability. The fields are buoyant, so the accretion heating occurs higher in the atmosphere than the stresses. The heating is localized around current sheets, caused by magneto-rotational instability at lower elevations and by Parker instability at higher elevations. Gas in the sheets is heated above the stellar irradiation temperature, even though accretion is much less than irradiation power when volume-averaged. The hot optically-thin current sheets might be detectable through their line emission.

We measured the angular diameter of the lithium-rich K giant star HD 148293 using Georgia State University's Center for High Angular Resolution Astronomy (CHARA) Array interferometer. We used our measurement to calculate the star's effective temperature, which allowed us to place it on an H-R diagram to compare it with other Li-rich giants. Its placement supports the evidence presented by Charbonnel & Balachandran that it is undergoing a brief stage in its evolution where Li is being created.

GJ 581 is an M dwarf host of a multiplanet system. We use long-baseline interferometric measurements from the CHARA Array, coupled with trigonometric parallax information, to directly determine its physical radius to be $0.299 \pm 0.010 R_{\odot}$. Literature photometry data are used to perform spectral energy distribution fitting in order to determine GJ 581's effective surface temperature $T_{\rm EFF}=3498 \pm 56$ K and its luminosity $L=0.01205 \pm 0.00024 L_{\odot}$. From these measurements, we recompute the location and extent of the system's habitable zone and conclude that two of the planets orbiting GJ 581, planets d and g, spend all or part of their orbit within or just on the edge of the habitable zone.

Variability is a defining characteristic of young stellar systems, and optical variability has been heavily studied to select and characterize the photospheric properties of young stars. In recent years, multi-epoch observations sampling a wider range of wavelengths and time-scales have revealed a wealth of time-variable phenomena at work during the star formation process. This splinter session was convened to summarize recent progress in providing improved coverage and understanding of time-variable processes in young stars and circumstellar disks. We begin by summarizing results from several multi-epoch Spitzer campaigns, which have demonstrated that many young stellar objects evidence significant mid-IR variability. While some of these variations can be attributed to processes in the stellar photosphere, others appear to trace short time-scale changes in the circumstellar disk which can be successfully modeled with axisymmetric or non-axisymmetric structures. We also review recent studies probing variability at shorter wavelengths that provide evidence for high frequency pulsations associated with accretion outbursts, correlated optical/X-ray variability in Classical T Tauri stars, and magnetic reversals in young solar analogs.

We present measurements of fundamental astrophysical properties of nearby, low-mass, K- and M-dwarfs from our DISCOS survey (DIameterS of COol Stars). The principal goal of our study is the determination of linear radii and effective temperatures for these stars. We calculate their radii from angular diameter measurements using the CHARA Array and Hipparcos distances. Combined with bolometric flux measurements based on literature photometry, we use our angular diameter results to calculate their effective surface temperatures. We present preliminary results established on an assortment of empirical relations to the stellar effective temperature and radius that are based upon these measurements. We elaborate on the discrepancy seen between theoretical and observed stellar radii, previously claimed to be related to stellar activity and/or metallicity. Our preliminary conclusion, however, is that convection plays a larger role in the determination of radii of these late-type stars. Understanding the source of the radius disagreement is likely to impact other areas of study for low-mass stars, such as the detection and characterization of extrasolar planets in the habitable zones.

Physical processes working in the stellar interiors as well as the evolution of stars depend on some fundamental stellar properties, such as mass, radius, luminosity, and chemical abundances. A classical way to test stellar interior models is to compare the predicted and observed location of a star on theoretical evolutionary tracks in a H-R diagram. This requires the best possible determinations of stellar mass, radius, luminosity and abundances. To derive its fundamental parameters, we observed the well-known rapidly oscillating Ap star, $\gamma$ Equ, using the visible spectro-interferometer VEGA installed on the optical CHARA array. We computed the calibrated squared visibility and derived the limb-darkened diameter. We used the whole energy flux distribution, the parallax and this angular diameter to determine the luminosity and the effective temperature of the star. We obtained a limb-darkened angular diameter of 0.564~$\pm$~0.017~mas and deduced a radius of $R$~=~2.20~$\pm$~0.12~${\rm R_{\odot}}$. Without considering the multiple nature of the system, we derived a bolometric flux of $(3.12\pm 0.21)\times 10^{-7}$ erg~cm$^{-2}$~s$^{-1}$ and an effective temperature of 7364~$\pm$~235~K, which is below the effective temperature that has been previously determined. Under the same conditions we found a luminosity of $L$~=~12.8~$\pm$~1.4~${\rm L_{\odot}}$. When the contribution of the closest companion to the bolometric flux is considered, we found that the effective temperature and luminosity of the primary star can be, respectively, up to $\sim$~100~K and up to $\sim$~0.8~L$_\odot$ smaller than the values mentioned above.These new values of the radius and effective temperature should bring further constraints on the asteroseismic modelling of the star.

High-precision interferometric measurements of pulsating stars help to characterize their close environment. In 1974, a close companion was discovered around the pulsating star beta Cep using the speckle interferometry technique and features at the limit of resolution (20 milli-arcsecond or mas) of the instrument were mentioned that may be due to circumstellar material. Beta Cep has a magnetic field that might be responsible for a spherical shell or ring-like structure around the star as described by the MHD models. Using the visible recombiner VEGA installed on the CHARA long-baseline interferometer at Mt. Wilson, we aim to determine the angular diameter of beta Cep and resolve its close environment with a spatial resolution up to 1 mas level. Medium spectral resolution (R=6000) observations of beta Cep were secured with the VEGA instrument over the years 2008 and 2009. These observations were performed with the S1S2 (30m) and W1W2 (100m) baselines of the array. We investigated several models to reproduce our observations. A large-scale structure of a few mas is clearly detected around the star with a typical flux relative contribution of 0.23 +- 0.02. Our best model is a co-rotational geometrical thin ring around the star as predicted by magnetically-confined wind shock models. The ring inner diameter is 8.2 +- 0.8 mas and the width is 0.6 +- 0.7 mas. The orientation of the rotation axis on the plane of the sky is PA = 60 +- 1 deg, while the best fit of the mean angular diameter of beta Cep gives UD[V] = 0.22 +- 0.05 mas. Our data are compatible with the predicted position of the close companion of beta Cep. These results bring additional constraints on the fundamental parameters and on the future MHD and asteroseismological models of the star.

We present interferometric observations of the Be star Zeta Tau obtained using the MIRC beam combiner at the CHARA Array. We resolved the disk during four epochs in 2007-2009. We fit the data with a geometric model to characterize the circumstellar disk as a skewed elliptical Gaussian and the central Be star as a uniform disk. The visibilities reveal a nearly edge-on disk with a FWHM major axis of ~ 1.8 mas in the H-band. The non-zero closure phases indicate an asymmetry within the disk. Interestingly, when combining our results with previously published interferometric observations of Zeta Tau, we find a correlation between the position angle of the disk and the spectroscopic V/R ratio, suggesting that the tilt of the disk is precessing. This work is part of a multi-year monitoring campaign to investigate the development and outward motion of asymmetric structures in the disks of Be stars.

BA-type supergiants are amongst the most optically-bright stars. They are observable in extragalactic environments, hence potential accurate distance indicators. Emission activity in the Halpha line of the BA supergiants Rigel (B8Ia) and Deneb (A2Ia) is indicative of presence of localized time-dependent mass ejections. Here, we employ optical interferometry to study the Halpha line-formation region in these stellar environments. High spatial- (0.001 arcsec) and spectral- (R=30 000) resolution observations of Halpha were obtained with the visible recombiner VEGA installed on the CHARA interferometer, using the S1S2 array-baseline (34m). Six independent observations were done on Deneb over the years 2008 and 2009, and two on Rigel in 2009. We analyze this dataset with the 1D non-LTE radiative-transfer code CMFGEN, and assess the impact of the wind on the visible and near-IR interferometric signatures, using both Balmer-line and continuum photons. We observe a visibility decrease in Halpha for both Rigel and Deneb, suggesting that the line-formation region is extended (1.5-1.75 R*). We observe a significant visibility decrease for Deneb in the SiII6371 line. We witness time variations in the differential phase for Deneb, implying an inhomogeneous and unsteady circumstellar environment, while no such variability is seen in differential visibilities. Radiative-transfer modeling of Deneb, with allowance for stellar-wind mass loss, accounts fairly well for the observed decrease in the Halpha visibility. Based on the observed differential visibilities, we estimate that the mass-loss rate of Deneb has changed by less than 5%.

The young open cluster Collinder 419 surrounds the massive O star, HD 193322, that is itself a remarkable multiple star system containing at least four components. Here we present a discussion of the cluster distance based upon new spectral classifications of the brighter members, UBV photometry, and an analysis of astrometric and photometric data from the UCAC3 and 2MASS catalogs. We determine an average cluster reddening of E(B-V)=0.37 +- 0.05 mag and a cluster distance of 741 +- 36 pc. The cluster probably contains some very young stars that may include a reddened M3 III star, IRAS~20161+4035.

A crucial issue in star formation is to understand the physical mechanism by which mass is accreted onto and ejected by a young star. The visible spectrometer VEGA on the CHARA array can be an efficient means of probing the structure and the kinematics of the hot circumstellar gas at sub-AU. For the first time, we observed the Herbig Ae star AB Aur in the H$\alpha$ emission line, using the VEGA low spectral resolution on two baselines of the array. We computed and calibrated the spectral visibilities between 610 nm and 700 nm. To simultaneously reproduce the line profile and the visibility, we used a 1-D radiative transfer code that calculates level populations for hydrogen atoms in a spherical geometry and synthetic spectro-interferometric observables. We clearly resolved AB Aur in the H$\alpha$ line and in a part of the continuum, even at the smallest baseline of 34 m. The small P-Cygni absorption feature is indicative of an outflow but could not be explained by a spherical stellar wind model. Instead, it favors a magneto-centrifugal X-disk or disk-wind geometry. The fit of the spectral visibilities could not be accounted for by a wind alone, so we considered a brightness asymmetry possibly caused by large-scale nebulosity or by the known spiral structures, inducing a visibility modulation around H$\alpha$. Thanks to the unique capabilities of VEGA, we managed to simultaneously record for the first time a spectrum at a resolution of 1700 and spectral visibilities in the visible range on a target as faint as $m_{V}$ = 7.1. It was possible to rule out a spherical geometry for the wind of AB Aur and provide realistic solutions to account for the H$\alpha$ emission compatible with magneto-centrifugal acceleration. The study illustrates the advantages of optical interferometry and motivates observations of other bright young stars to shed light on the accretion/ejection processes.

Of the over 450 exoplanets known to date, more than 420 of them have been discovered using radial velocity studies, a method that tells nothing about the inclination of the planet's orbit. Because it is more likely that the companion is a planetary-mass object in a moderate- to high-inclination orbit than a low-mass stellar object in a nearly face-on orbit, the secondary bodies are presumed to be planets. Interferometric observations allow us to inspect the angular diameter fit residuals to calibrated visibilities in order to rule out the possibility of a low-mass stellar companion in a very low-inclination orbit. We used the Center for High Angular Resolution Astronomy (CHARA) Array interferometer to observe 20 exoplanet host stars and considered five potential secondary spectral types: G5 V, K0 V, K5 V, M0 V, and M5 V. If a secondary star is present and is sufficiently bright, the effects of the added light will appear in interferometric observations where the planet will not. All secondary types could be eliminated from consideration for 7 host stars and no secondary stars of any spectral type could be ruled out for 7 more. The remaining 6 host stars showed a range of possible secondary types.

Eclipses of the single-line spectroscopic binary star, epsilon Aurigae, provide an opportunity to study the poorly-defined companion. We used the MIRC beam combiner on the CHARA array to create interferometric images during eclipse ingress. Our results demonstrate that the eclipsing body is a dark disk that is opaque and tilted, and therefore exclude alternative models for the system. These data constrain the geometry and masses of the components, providing evidence that the F-star is not a massive supergiant star.

The poorly-ionized interior of the protoplanetary disk is the location where dust coagulation processes may be most efficient. However even here, planetesimal formation may be limited by the loss of solid material through radial drift, and by collisional fragmentation of the particles. Our aim is to investigate the possibility that solid particles are trapped at local pressure maxima in the dynamically evolving disk. We perform the first 3-D global non-ideal MHD calculations of the disk treating the turbulence driven by the magneto-rotational instability. The domain contains an inner MRI-active region near the young star and an outer midplane dead zone, with the transition between the two modeled by a sharp increase in the magnetic diffusivity. The azimuthal magnetic fields generated in the active zone oscillate over time, changing sign about every 150 years. We thus observe the radial structure of the `butterfly pattern' seen previously in local shearing-box simulations. The mean magnetic field diffuses from the active zone into the dead zone, where the Reynolds stress nevertheless dominates. The greater total accretion stress in the active zone leads to a net reduction in the surface density, so that after 800 years an approximate steady state is reached in which a local radial maximum in the midplane pressure lies near the transition radius. We also observe the formation of density ridges within the active zone. The dead zone in our models possesses a mean magnetic field, significant Reynolds stresses and a steady local pressure maximum at the inner edge, where the outward migration of planetary embryos and the efficient trapping of solid material are possible.

Using Georgia State University's CHARA Array interferometer, we measured angular diameters for 25 giant stars, six of which host exoplanets. The combination of these measurements and Hipparcos parallaxes produce physical linear radii for the sample. Except for two outliers, our values match angular diameters and physical radii estimated using photometric methods to within the associated errors with the advantage that our uncertainties are significantly lower. We also calculated the effective temperatures for the stars using the newly-measured diameters. Our values do not match those derived from spectroscopic observations as well, perhaps due to the inherent properties of the methods used or because of a missing source of extinction in the stellar models that would affect the spectroscopic temperatures.

The prospects for using asteroseismology of rapidly oscillating Ap (roAp) stars are hampered by the large uncertainty in fundamental stellar parameters. Results in the literature for the effective temperature (Teff) often span a range of 1000 K. Our goal is to reduce systematic errors and improve the Teff calibration of Ap stars based on new interferometric measurements. We obtained long-baseline interferometric observations of beta CrB using the CHARA/FLUOR instrument. To disentangle the flux contributions of the two components of this binary star, we obtained VLT/NACO adaptive optics images. We determined limb darkened angular diameters of 0.699+-0.017 mas for beta CrB A (from interferometry) and 0.415+-0.017 mas for beta CrB B (from surface brightness- color relations), corresponding to radii of 2.63+-0.09 Rsun (3.4 percent uncertainty) and 1.56+-0.07 Rsun (4.5 percent). The combined bolometric flux of the A and B components was determined from satellite UV data, spectrophotometry in the visible and broadband data in the infrared. The flux from the B component constitutes 16+-4 percent of the total flux and was determined by fitting an ATLAS9 model atmosphere to the broad-band NACO J and K magnitudes. Combining the flux of the A component with its measured angular diameter, we determine the effective temperature Teff(A) = 7980+-180 K (2.3 percent). Our new interferometric and imaging data enable a nearly model-independent determination of the effective temperature of beta CrB A. Including our recent study of alpha Cir, we now have direct Teff measurements of two of the brightest roAp stars, providing a strong benchmark for an improved calibration of the Teff scale for Ap stars. This will support the use of potentially strong constraints imposed by asteroseismic studies of roAp stars.

We apply ionization balance and MHD calculations to investigate whether magnetic activity moderated by recombination on dust can account for the mass accretion rates and the mid-infrared spectra and variability of protostellar disks. The MHD calculations use the stratified shearing-box approach and include grain settling and the feedback from the changing dust abundance on the resistivity of the gas. The two-decade spread in accretion rates among T Tauri stars is too large to result solely from variety in the grain size and stellar X-ray luminosity, but can be produced by varying these together with the disk magnetic flux. The diversity in the silicate bands can come from the coupling of grain settling to the distribution of the magneto-rotational turbulence, through three effects: (1) Recombination on grains yields a magnetically inactive dead zone extending above two scale heights, while turbulence in the magnetically active disk atmosphere overshoots the dead zone boundary by only about one scale height. (2) Grains deep in the dead zone oscillate vertically in waves driven by the turbulent layer above, but on average settle at the laminar rates, so the interior of the dead zone is a particle sink and the disk atmosphere becomes dust-depleted. (3) With sufficient depletion, the dead zone is thinner and mixing dredges grains off the midplane. The MHD results also show that the magnetic activity intermittently lifts clouds of dust into the atmosphere. The photosphere height changes by up to one-third over a few orbits, while the extinction along lines of sight grazing the disk surface varies by factors of two over times down to 0.1 orbit. We suggest that the changing shadows cast by the dust clouds on the outer disk are a cause of the daily to monthly mid-infrared variability in some young stars. (Abridged.)

We determined the spatial orientation of the Algol AB close pair orbital plane using optical interferometry with the CHARA Array, and radio interferometry with the European VLBI Network (EVN). We found the longitude of the line of nodes for the close pair being $\Omega_1=48\deg\pm2\deg$ and the mutual inclination of the orbital planes of the close and the wide pairs being $95\deg\pm3\deg$. This latter value differs by $5\deg$ from the formerly known $100\deg$ which would imply a very fast inclination variation of the system, not supported by the photometric observations. We also investigated the dynamics of the system with numerical integration of the equations of motions using our result as an initial condition. We found large variations in the inclination of the close pair (its amplitude $\sim 170\deg$) with a period of about 20 millenia. This result is in good agreement with the photometrically observed change of amplitude in Algol's primary minimum.

Algol is a triple stellar system consisting of a close semidetached binary orbited by a third object. Due to the disputed spatial orientation of the close pair, the third body perturbation of this pair is a subject of much research. In this study, we determine the spatial orientation of the close pair orbital plane using the CHARA Array, a six-element optical/IR interferometer located on Mount Wilson, and state-of-the-art e-EVN interferometric techniques. We find that the longitude of the line of nodes for the close pair is $\Omega_1=48\degr\pm2\degr$ and the mutual inclination of the orbital planes of the close and the wide pairs is $95\degr\pm3\degr$. This latter value differs by $5\degr$ from the formerly known $100\degr$ which would imply a very fast inclination variation of the system, not borne out by the photometric observations. We also investigated the dynamics of the system with numerical integration of the equations of motions using our result as an initial condition. We found large variations in the inclination of the close pair (its amplitude $\sim 170\degr$) with a period of about 20 millennia. This result is in good agreement with the photometrically observed change of amplitude in Algol's primary minimum.

We calculate the location of the magnetically-inactive dead zone in the minimum-mass protosolar disk, under ionization scenarios including stellar X-rays, long- or short-lived radionuclide decay, and energetic protons arriving from the general interstellar medium, from a nearby supernova explosion, from the disk corona, or from the corona of the young star. The disk contains a dead zone in all scenarios except those with small dust grains removed and a fraction of the short-lived radionuclides remaining in the gas. All the cases without exception have an "undead zone" where intermediate resistivities prevent magneto-rotational turbulence while allowing shear-generated large-scale magnetic fields. The mass column in the undead zone is typically greater than the column in the turbulent surface layers. The results support the idea that the dead and undead zones are robust consequences of cold, dusty gas with mass columns exceeding 1000 g/cm^2.

The GIII red giant star epsilon Oph has been found to exhibit several modes of oscillation by the MOST mission. We interpret the observed frequencies of oscillation in terms of theoretical radial p-mode frequencies of stellar models. Evolutionary models of this star, in both shell H-burning and core He-burning phases of evolution, are constructed using as constraints a combination of measurements from classical ground-based observations (for luminosity, temperature, and chemical composition) and seismic observations from MOST. Radial frequencies of models in either evolutionary phase can reproduce the observed frequency spectrum of epsilon Oph almost equally well. The best-fit models indicate a mass in the range of 1.85 +/- 0.05 Msun with radius of 10.55 +/- 0.15 Rsun. We also obtain an independent estimate of the radius of epsilon Oph using high accuracy interferometric observations in the infrared K' band, using the CHARA/FLUOR instrument. The measured limb darkened disk angular diameter of epsilon Oph is 2.961 +/- 0.007 mas. Together with the Hipparcos parallax, this translates into a photospheric radius of 10.39 +/- 0.07 Rsun. The radius obtained from the asteroseismic analysis matches the interferometric value quite closely even though the radius was not constrained during the modelling.

We directly measured the angular diameters for 11 exoplanet host stars using Georgia State University's CHARA Array interferometer and calculated their linear radii and effective temperatures. The sample tends towards evolving or evolved stars and includes one dwarf, four subgiants, and six giants. We then estimated masses and ages for the stars using our effective temperatures combined with metallicity measurements from the literature.

We present sub-milliarcseond resolution imaging and modeling of two nearby rapid rotators Alpha Cephei and Alpha Ophiuchi, obtained with the CHARA array - the largest optical/IR interferometer in the world. Incorporating a gravity darkening model, we are able to determine the inclination, the polar and equatorial radius and temperature, as well as the fractional rotation speed of the two stars with unprecedented precision. The polar and equatorial regions of the two stars have ~2000K temperature gradient, causing their apparent temperatures and luminosities to be dependent on their viewing angles. Our modeling allow us to determine the true effective temperatures and luminosities of Alpha Cep and Alpha Oph, permitting us to investigate their true locations on the H-R diagram. These properties in turn give us estimates of the masses and ages of the two stars within a few percent of error using stellar evolution models. Also, based on our gravity darkening modeling, we propose a new method to estimate the masses of single stars in a more direct way through Vsin(i) measurements and precise geometrical constraint. Lastly, we investigate the degeneracy between the inclination and the gravity darkening coefficient, which especially affects the modeling of Alpha Oph. Although incorporating Vsin(i) has lifted the degeneracy to some extent, higher resolution observations are still needed to further constrain the parameters independently.

Almost all stars in the 1-8 Msun range evolve through the Asymptotic Giant Branch (AGB), preplanetary nebula (PPN) and planetary nebula (PN) evolutionary phases. Most stars that leave the main sequence in a Hubble time will end their lives in this way. The heavy mass loss which occurs during the AGB phase is important across astrophysics, and the particulate matter crucial for the birth of new solar systems is made and ejected by AGB stars. Yet stellar evolution from the beginning of the AGB phase to the PN phase remains poorly understood. We do not understand how the mass-loss (rate, geometry, temporal history) depends on fundamental stellar parameters or the presence of a binary companion. While the study of evolved non-massive stars has maintained a relatively modest profile in recent decades, we are nonetheless in the midst of a quiet but exciting revolution in this area, driven by new observational results, such as the discovery of jets and disks in stellar environments where these were never expected, and by the recognition of new symmetries such as multipolarity and point-symmetry occuring frequently in the nebulae resulting from the outflows. In this paper we summarise the major unsolved problems in this field, and specify the areas where allocation of effort and resources is most likely to help make significant progress.

To date, there is no core accretion simulation that can successfully account for the formation of Uranus or Neptune within the observed 2-3 Myr lifetimes of protoplanetary disks. Since solid accretion rate is directly proportional to the available planetesimal surface density, one way to speed up planet formation is to take a full accounting of all the planetesimal-forming solids present in the solar nebula. By combining a viscously evolving protostellar disk with a kinetic model of ice formation, we calculate the solid surface density in the solar nebula as a function of heliocentric distance and time. We find three effects that strongly favor giant planet formation: (1) a decretion flow that brings mass from the inner solar nebula to the giant planet-forming region, (2) recent lab results (Collings et al. 2004) showing that the ammonia and water ice lines should coincide, and (3) the presence of a substantial amount of methane ice in the trans-Saturnian region. Our results show higher solid surface densities than assumed in the core accretion models of Pollack et al. (1996) by a factor of 3 to 4 throughout the trans-Saturnian region. We also discuss the location of ice lines and their movement through the solar nebula, and provide new constraints on the possible initial disk configurations from gravitational stability arguments.

We present infrared interferometric observations of the inner regions of two A-star debris disks, beta Leo and zeta Lep, using the FLUOR instrument at the CHARA interferometer on both short (30 m) and long (>200 m) baselines. For the target stars, the short baseline visibilities are lower than expected for the stellar photosphere alone, while those of a check star, delta Leo, are not. We interpret this visibility offset of a few percent as a near-infrared excess arising from dust grains which, due to the instrumental field of view, must be located within several AU of the central star. For beta Leo, the near-infrared excess producing grains are spatially distinct from the dust which produces the previously known mid-infrared excess. For zeta Lep, the near-infrared excess may be spatially associated with the mid-infrared excess producing material. We present simple geometric models which are consistent with the near and mid-infrared excess and show that for both objects, the near-infrared producing material is most consistent with a thin ring of dust near the sublimation radius with typical grain sizes smaller than the nominal radiation pressure blowout radius. Finally, we discuss possible origins of the near-infrared emitting dust in the context of debris disk evolution models.

We present angular diameters of the Hyades giants, gamma, delta^1, epsilon, and theta^1 Tau from interferometric measurements with the CHARA Array. Our errors in the limb-darkened angular diameters for these stars are all less than 2%, and in combination with additional observable quantities, we determine the effective temperatures, linear radii and absolute luminosities for each of these stars. Additionally, stellar masses are inferred from model isochrones to determine the surface gravities. These data show that a new calibration of effective temperatures with errors well under 100K is now possible from interferometric angular diameters of stars.

We present two new in situ core accretion simulations of Saturn with planet formation timescales of 3.37 Myr (model S0) and 3.48 Myr (model S1), consistent with observed protostellar disk lifetimes. In model S0, we assume rapid grain settling reduces opacity due to grains from full interstellar values (Podolak 2003). In model S1, we do not invoke grain settling, instead assigning full interstellar opacities to grains in the envelope. Surprisingly, the two models produce nearly identical formation timescales and core/atmosphere mass ratios. We therefore observe a new manifestation of core accretion theory: at large heliocentric distances, the solid core growth rate (limited by Keplerian orbital velocity) controls the planet formation timescale. We argue that this paradigm should apply to Uranus and Neptune as well.

We present an updated spectroscopic orbit and a new visual orbit for the double-lined spectroscopic binary \sigma^2 Coronae Borealis based on radial velocity measurements at the Oak Ridge Observatory in Harvard, Massachusetts and interferometric visibility measurements at the CHARA Array on Mount Wilson. \sigma^2 CrB is composed of two Sun-like stars of roughly equal mass in a circularized orbit with a period of 1.14 days. The long baselines of the CHARA Array have allowed us to resolve the visual orbit for this pair, the shortest period binary yet resolved interferometrically, enabling us to determine component masses of 1.137 \pm 0.037 M_sun and 1.090 \pm 0.036 M_sun. We have also estimated absolute V-band magnitudes of MV (primary) = 4.35 \pm 0.02 and MV(secondary) = 4.74 \pm 0.02. A comparison with stellar evolution models indicates a relatively young age of 1-3 Gyr, consistent with the high Li abundance measured previously. This pair is the central component of a quintuple system, along with another similar-mass star, \sigma^1 CrB, in a ~ 730-year visual orbit, and a distant M-dwarf binary, \sigma CrB C, at a projected separation of ~ 10 arcmin. We also present differential proper motion evidence to show that components C & D (ADS 9979C & D) listed for this system in the Washington Double Star Catalog are optical alignments that are not gravitationally bound to the \sigma CrB system.

We present comprehensive models for the Herbig Ae stars MWC275 and AB Aur that aim to explain their spectral energy distribution (from UV to millimeter) and long baseline interferometry (from near-infrared to millimeter) simultaneously. Data from the literature, combined with new mid-infrared (MIR) interferometry from the Keck Segment Tilting Experiment, are modeled using an axisymmetric Monte Carlo radiative transfer code. Models in which most of the near-infrared (NIR) emission arises from a dust rim fail to fit the NIR spectral energy distribution (SED) and sub-milli-arcsecond NIR CHARA interferometry. Following recent work, we include an additional gas emission component with similar size scale to the dust rim, inside the sublimation radius, to fit the NIR SED and long-baseline NIR interferometry on MWC275 and AB Aur. In the absence of shielding of star light by gas, we show that the gas-dust transition region in these YSOs will have to contain highly refractory dust, sublimating at ~1850K. Despite having nearly identical structure in the thermal NIR, the outer disks of MWC275 and AB Aur differ substantially. In contrast to the AB Aur disk, MWC275 lacks small grains in the disk atmosphere capable of producing significant 10-20micron emission beyond ~7AU, forcing the outer regions into the "shadow" of the inner disk

We present the first resolved images of the eclipsing binary Beta Lyrae, obtained with the CHARA Array interferometer and the MIRC combiner in the H band. The images clearly show the mass donor and the thick disk surrounding the mass gainer at all six epochs of observation. The donor is brighter and generally appears elongated in the images, the first direct detection of photospheric tidal distortion due to Roche-lobe filling. We also confirm expectations that the disk component is more elongated than the donor and is relatively fainter at this wavelength. Image analysis and model fitting for each epoch were used for calculating the first astrometric orbital solution for Beta Lyrae, yielding precise values for the orbital inclination and position angle. The derived semi-ma jor axis also allows us to estimate the distance of Beta Lyrae; however, systematic differences between the models and the images limit the accuracy of our distance estimate to about 15%. To address these issues, we will need a more physical, self-consistent model to account for all epochs as well as the multi-wavelength information from the eclipsing light curves.

High-precision interferometric observations of six early-type main sequence stars known to harbour cold debris discs have been obtained in the near-infrared K band with the FLUOR instrument at the CHARA Array. The measured squared visibilities are compared to the expected visibility of the stellar photospheres based on theoretical photospheric models taking into account rotational distortion, searching for potential visibility reduction at short baselines due to circumstellar emission. Our observations bring to light the presence of resolved circumstellar emission around one of the six target stars (zeta Aql) at the 5 sigma level. The morphology of the emission source cannot be directly constrained because of the sparse spatial frequency sampling of our interferometric data. Using complementary adaptive optics observations and radial velocity measurements, we find that the presence of a low-mass companion is a likely origin for the excess emission. The potential companion has a K-band contrast of four magnitudes, a most probable mass of about 0.6 Msun, and is expected to orbit between about 5.5 AU and 8 AU from its host star assuming a purely circular orbit. Nevertheless, by adjusting a physical debris disc model to the observed Spectral Energy Distribution of the zeta Aql system, we also show that the presence of hot dust within 10 AU from zeta Aql, producing a total thermal emission equal to 1.69 +- 0.31% of the photospheric flux in the K band, is another viable explanation for the observed near-infrared excess. Our re-interpretation of archival near- to far-infrared photometric measurements shows however that cold dust is not present around zeta Aql at the sensitivity limit of the IRS and MIPS instruments onboard Spitzer, and urges us to remove zeta Aql from the category of bona fide debris disc stars.

Context: The main sequence binary star 61 Cyg (K5V+K7V) is our nearest stellar neighbour in the northern hemisphere. This proximity makes it a particularly well suited system for very high accuracy interferometric radius measurements. Aims: Our goal is to constrain the poorly known evolutionary status and age of this bright binary star. Methods: We obtained high accuracy interferometric observations in the infrared K' band, using the CHARA/FLUOR instrument. We then computed evolutionary models of 61 Cyg A & B with the CESAM2k code. As model constraints, we used a combination of observational parameters from classical observation methods (photometry, spectroscopy) as well as our new interferometric radii. Results: The measured limb darkened disk angular diameters are theta_LD(A) = 1.775 +/- 0.013 mas and theta_LD(B) = 1.581 +/- 0.022 mas, respectively for 61 Cyg A and B. Considering the high accuracy parallaxes available, these values translate into photospheric radii of R(A) = 0.665 +/- 0.005 Rsun and R(B) = 0.595 +/- 0.008 Rsun. The new radii constrain efficiently the physical parameters adopted for the modeling of both stars, allowing us to predict asteroseismic frequencies based on our best-fit models. Conclusions: The CESAM2k evolutionary models indicate an age around 6 Gyrs and are compatible with small values of the mixing length parameter. The measurement of asteroseismic oscillation frequencies in 61 Cyg A & B would be of great value to improve the modeling of this important fiducial stellar system, in particular to better constrain the masses.

We examine a physical process that leads to the efficient formation of gas giant planets around intermediate mass stars. In the gaseous protoplanetary disks surrounding rapidly-accreting intermediate-mass stars we show that the midplane temperature (heated primarily by turbulent dissipation) can reach > 1000 K out to 1 AU. Thermal ionization of this hot gas couples the disk to the magnetic field, allowing the magneto-rotational instability (MRI) to generate turbulence and transport angular momentum. Further from the central star the ionization fraction decreases, decoupling the disk from the magnetic field and reducing the efficiency of angular momentum transport. As the disk evolves towards a quasi-steady state, a local maximum in the surface density and in the midplane pressure both develop at the inner edge of the MRI-dead zone, trapping inwardly migrating solid bodies. Small particles accumulate and coagulate into planetesimals which grow rapidly until they reach isolation mass. In contrast to the situation around solar type stars, we show that the isolation mass for cores at this critical radius around the more massive stars is large enough to promote the accretion of significant amounts of gas prior to disk depletion. Through this process, we anticipate a prolific production of gas giants at ~1 AU around intermediate-mass stars.

We present the results of an adaptive optics survey for faint companions among Galactic O-type star systems (V < 8) using the Advanced Electro-Optical System (AEOS) 3.6-meter telescope on Haleakala. We surveyed these O star systems in I-band, typically being able to detect a companion with a magnitude difference of about 6 in the projected separation range of 0.5 to 1.0 arcseconds, and of about 9.5 in the range of 1.0 to 5.0 arcseconds. In the course of the survey, we discovered 40 new companions among 31 of the 116 objects examined and made astrometric and differential magnitude measurements of 24 additional known pairs, several of them being confirmation detections. We present new astrometric orbits for two binaries, BU 1032AB (WDS 05387-0236; sig Ori AB) and SEE 322 (WDS 17158-3344; HD 155889AB). We lack magnitude differences for other filter bands, so it is difficult to determine physical from line-of-sight companions, but we present empirical arguments for the limiting magnitude difference where field contamination is significant.

Using the longest baselines of the CHARA Array, we have measured the angular diameter of the G5 V subdwarf $\mu$ Cas A, the first such determination for a halo population star. We compare this result to new diameters for the higher metallicity K0 V stars, $\sigma$ Dra and HR 511, and find that the metal-poor star, $\mu$ Cas A, has an effective temperature ($T_{\rm eff}=5297\pm32$ K), radius ($R=0.791\pm0.008 R_{\rm \odot}$), and absolute luminosity ($L=0.442\pm0.014 L_{\rm \odot}$) comparable to the other two stars with later spectral types. We show that stellar models show a discrepancy in the predicted temperature and radius for $\mu$ Cas A, and we discuss these results and how they provide a key to understanding the fundamental relationships for stars with low metallicity.

Planets form inside protostellar disks in a dead zone where the electrical resistivity of the gas is too high for magnetic forces to drive turbulence. We show that much of the dead zone nevertheless is active and flows toward the star while smooth, large-scale magnetic fields transfer the orbital angular momentum radially outward. Stellar X-ray and radionuclide ionization sustain a weak coupling of the dead zone gas to the magnetic fields, despite the rapid recombination of free charges on dust grains. Net radial magnetic fields are generated in the magneto-rotational turbulence in the electrically conducting top and bottom surface layers of the disk, and reach the midplane by Ohmic diffusion. A toroidal component to the fields is produced near the midplane by the orbital shear. The process is similar to the magnetization of the Solar tachocline. The result is a laminar, magnetically-driven accretion flow in the region where the planets form.

Most exoplanets have been discovered via radial velocity studies, which are inherently insensitive to orbital inclination. Interferometric observations will show evidence of a stellar companion if it sufficiently bright, regardless of the inclination. Using the CHARA Array, we observed 22 exoplanet host stars to search for stellar companions in low-inclination orbits that may be masquerading as planetary systems. While no definitive stellar companions were discovered, it was possible to rule out certain secondary spectral types for each exoplanet system observed by studying the errors in the diameter fit to calibrated visibilities and by searching for separated fringe packets.

Using the longest optical-interferometeric baselines currently available, we have detected strong near-infrared (NIR) emission from inside the dust-destruction radius of Herbig Ae stars MWC275 and AB Aur. Our sub-milli-arcsecond resolution observations unambiguously place the emission between the dust-destruction radius and the magnetospheric co-rotation radius. We argue that this new component corresponds to hot gas inside the dust-sublimation radius, confirming recent claims based on spectrally-resolved interferometry and dust evaporation front modeling.

We have measured the angular diameters for a sample of 24 exoplanet host stars using Georgia State University's CHARA Array interferometer. We use these improved angular diameters together with Hipparcos parallax measurements to derive linear radii and to estimate the stars' evolutionary states.

We perform numerical simulations of solid particle motion in a shearing box model of a protoplanetary disc. The accretion flow is turbulent due to the action of the magnetorotational instability. Aerodynamic drag on the particles is modelled using the Epstein law with the gas velocity interpolated to the particle position. The effect of the magnetohydrodynamic turbulence on particle velocity dispersions is quantified for solids of different stopping times t_s, or equivalently, different sizes. The anisotropy of the turbulence is reflected upon the dispersions of the particle velocity components, with the radial component larger than both the azimuthal and vertical components for particles larger than ~ 10 cm (assuming minimum-mass solar nebula conditions at 5 AU). The dispersion of the particle velocity magnitude, as well as that of the radial and azimuthal components, as functions of stopping time, agree with previous analytical results for isotropic turbulence. The relative speed between pairs of particles with the same value of t_s decays faster with decreasing separation than in the case of solids with different stopping time. Correlations in the particle number density introduce a non-uniform spatial distribution of solids in the 10 to 100 cm size range. Any clump of particles is disrupted by the turbulence in less than one tenth on an orbital period, and the maximally concentrated clumps are stable against self-gravitational collapse.

We show that space-based telescopes such as the proposed Terrestrial Planet Finder Coronagraph will be able to detect the light scattered by the interstellar grains along lines of sight passing near stars in our Galaxy. The relative flux of the scattered light within one arcsecond of a star at 100 pc in a uniform interstellar medium of 0.1 H atoms cm^-3 is about 10^-7. The halo increases in strength with the distance to the star and is unlikely to limit the coronagraphic detection of planets around the nearest stars. Grains passing within 100 AU of Sun-like stars are deflected by radiation, gravity and magnetic forces, leading to features in the scattered light that can potentially reveal the strength of the stellar wind, the orientation of the stellar magnetic field and the relative motion between the star and the surrounding interstellar medium.

We probed the first 3AU around tau Ceti and epsilon Eridani with the CHARA array (Mt Wilson, USA) in order to gauge the 2micron excess flux emanating from possible hot dust grains in the debris disks and to also resolve the stellar photospheres. High precision visibility amplitude measurements were performed with the FLUOR single mode fiber instrument and telescope pairs on baselines ranging from 22 to 241m of projected length. The short baseline observations allow us to disentangle the contribution of an extended structure from the photospheric emission, while the long baselines constrain the stellar diameter. We have detected a resolved emission around tau Cet, corresponding to a spatially integrated, fractional excess flux of 0.98 +/- 0.21 x 10^{-2} with respect to the photospheric flux in the K'-band. Around eps Eri, our measurements can exclude a fractional excess of greater than 0.6x10^{-2} (3sigma). We interpret the photometric excess around tau Cet as a possible signature of hot grains in the inner debris disk and demonstrate that a faint, physical or background, companion can be safely excluded. In addition, we measured both stellar angular diameters with an unprecedented accuracy: Theta_LD(tau Cet)= 2.015 +/- 0.011 mas and Theta_LD(eps Eri)=2.126 +/- 0.014 mas.

Spatially resolving the surfaces of nearby stars promises to advance our knowledge of stellar physics. Using optical long-baseline interferometry, we present here a near-infrared image of the rapidly rotating hot star Altair with <1 milliarcsecond resolution. The image clearly reveals the strong effect of gravity darkening on the highly-distorted stellar photosphere. Standard models for a uniformly rotating star can not explain our results, requiring differential rotation, alternative gravity darkening laws, or both.

We have measured the angular diameter of the transiting extrasolar planet host star HD 189733 using the CHARA O/IR interferometric array. Combining our new angular diameter of 0.377+/-0.024 mas with the Hipparcos parallax leads to a linear radius for the host star of 0.779+/-0.052 Rsol and a radius for the planet of 1.19+/-0.08 RJup. Adopting the mass of the planet as derived by its discoverers, we derive a mean density of the planet of 0.91+/-0.18 g cm-3. This is the first determination of the diameter of an extrasolar planet through purely direct means.

We show that the optically-thick dusty envelopes surrounding young high-mass stars are subject to the photon bubble instability. The infrared radiation passing through the envelope amplifies magnetosonic disturbances, with growth rates in our local numerical radiation MHD calculations that are consistent with a linear analysis. Modes with wavelengths comparable to the gas pressure scale height grow by more than two orders of magnitude in a thousand years, reaching non-linear amplitudes within the envelope lifetime. If the magnetic pressure in the envelope exceeds the gas pressure, the instability develops into trains of propagating shocks. Radiation escapes readily through the low-density material between the shocks, enabling accretion to continue despite the Eddington limit imposed by the dust opacity. The supersonic motions arising from the photon bubble instability can help explain the large velocity dispersions of hot molecular cores, while conditions in the shocked gas are suitable for maser emission. We conclude that the photon bubble instability may play a key role in the formation of massive stars.

Using the CHARA Array and the Palomar Testbed Interferometer, the chemically peculiar star $\lambda$ Bo\"{o}tis has been spatially resolved. We have measured the limb darkened angular diameter to be $\theta_{LD} = 0.533\pm0.029$ mas, corresponding to a linear radius of $R_{\star} = 1.70 \pm 0.10 R_\odot$. The measured angular diameter yields an effective temperature for $\lambda$ Boo of $T_{eff} = 8887 \pm 242$ K. Based upon literature surface gravity estimates spanning $\log{(g)} = 4.0-4.2$ $[\rm{cm s}^{-\rm{2}}]$, we have derived a stellar mass range of $M_{\star} = 1.1 - 1.7$ $M_\odot$. For a given surface gravity, the linear radius uncertainty contributes approximately $\sigma(M_\star) = 0.1-0.2 M_\odot$ to the total mass uncertainty. The uncertainty in the mass (i.e., the range of derived masses) is primarily a result of the uncertainty in the surface gravity. The upper bound of our derived mass range ($\log(g)=4.2, M_\star = 1.7\pm0.2 M_\odot$) is consistent with 100-300 MYr solar-metallicity evolutionary models. The mid-range of our derived masses ($\log(g)=4.1, M_\star = 1.3\pm0.2 M_\odot$) is consistent with 2-3 GYr metal-poor evolutionary models. A more definitive surface gravity determination is required to determine a more precise mass for $\lambda$ Boo.

We investigate the conditions for the presence of a magnetically inactive dead zone in protostellar disks, using 3-D shearing-box MHD calculations including vertical stratification, Ohmic resistivity and time-dependent ionization chemistry. Activity driven by the magnetorotational instability fills the whole thickness of the disk at 5 AU, provided cosmic ray ionization is present, small grains are absent and the gas-phase metal abundance is sufficiently high. At 1 AU the larger column density of 1700 g/cm^2 means the midplane is shielded from ionizing particles and remains magnetorotationally stable even under the most favorable conditions considered. Nevertheless the dead zone is effectively eliminated. Turbulence mixes free charges into the interior as they recombine, leading to a slight coupling of the midplane gas to the magnetic fields. Weak, large-scale radial fields diffuse to the midplane where they are sheared out to produce stronger azimuthal fields. The resulting midplane accretion stresses are just a few times less than in the surface layers on average.

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.

Only a handful of debris disks have been imaged up to now. Due to the need for high dynamic range and high angular resolution, very little is known about the inner planetary region, where small amounts of warm dust are expected to be found. We investigate the close neighbourhood of Vega with the help of infrared stellar interferometry and estimate the integrated K-band flux originating from the central 8 AU of the debris disk. We performed precise visibility measurements at both short (~30 m) and long (~150 m) baselines with the FLUOR beam-combiner installed at the CHARA Array (Mt Wilson, California) in order to separately resolve the emissions from the extended debris disk (short baselines) and from the stellar photosphere (long baselines). After revising Vega's K-band angular diameter (3.202+/-0.005 mas), we show that a significant deficit in squared visibility (1.88+/-0.34%) is detected at short baselines with respect to the best-fit uniform disk stellar model. This deficit can be either attributed to the presence of a low-mass stellar companion around Vega, or as the signature of the thermal and scattered emissions from the debris disk. We show that the presence of a close companion is highly unlikely, as well as other possible perturbations (stellar morphology, calibration), and deduce that we have most probably detected the presence of dust in the close neighbourhood of Vega. The resulting flux ratio between the stellar photosphere and the debris disk amounts to 1.29+/-0.19% within the FLUOR field-of-view (~7.8 AU). Finally, we complement our K-band study with archival photometric and interferometric data in order to evaluate the main physical properties of the inner dust disk. The inferred properties suggest that the Vega system could be currently undergoing major dynamical perturbations.

With the MIPS camera on Spitzer, we have searched for far-infrared emission from dust in the planetary system orbiting pulsar PSR 1257+12. With accuracies of 0.05 mJy at 24 um and 1.5 mJy at 70 um, photometric measurements find no evidence for emission at these wavelengths. These observations place new upper limits on the luminosity of dust with temperatures between 20 and 1000 K. They are particularly sensitive to dust temperatures of 100-200 K, for which they limit the dust luminosity to below $3 \times 10^{-5}$ of the pulsar's spin-down luminosity, three orders of magnitude better than previous limits. Despite these improved constraints on dust emission, an asteroid belt similar to the Solar System's cannot be ruled out.

We have obtained high-precision interferometric measurements of Vega with the CHARA Array and FLUOR beam combiner in the K' band at projected baselines between 103m and 273m. The measured visibility amplitudes beyond the first lobe are significantly weaker than expected for a slowly rotating star characterized by a single effective temperature and surface gravity. Our measurements, when compared to synthetic visibilities and synthetic spectrophotometry from a Roche-von Zeipel gravity-darkened model atmosphere, provide strong evidence for the model of Vega as a rapidly rotating star viewed very nearly pole-on. Our best fitting model indicates that Vega is rotating at ~91% of its angular break-up rate with an equatorial velocity of 275 km/s. Together with the measured vsin(i), this velocity yields an inclination for the rotation axis of 5 degrees. For this model the pole-to-equator effective temperature difference is 2250 K, a value much larger than previously derived from spectral line analyses. The derived equatorial T_eff of 7900 K indicates Vega's equatorial atmosphere may be convective and provides a possible explanation for the discrepancy. The model has a luminosity of ~37 Lsun, a value 35% lower than Vega's apparent luminosity based on its bolometric flux and parallax, assuming a slowly rotating star. The model luminosity is consistent with the mean absolute magnitude of A0V stars. Our model predicts the spectral energy distribution of Vega as viewed from its equatorial plane; a model which may be employed in radiative models for the surrounding debris disk.

Quasi-periodic oscillations (QPOs) in the X-ray lightcurves of accreting neutron star and black hole binaries have been widely interpreted as being due to standing wave modes in accretion disks. These disks are thought to be highly turbulent due to the magnetorotational instability (MRI). We study wave excitation by MRI turbulence in the shearing box geometry. We demonstrate that axisymmetric sound waves and radial epicyclic motions driven by MRI turbulence give rise to narrow, distinct peaks in the temporal power spectrum. Inertial waves, on the other hand, do not give rise to distinct peaks which rise significantly above the continuum noise spectrum set by MRI turbulence, even when the fluid motions are projected onto the eigenfunctions of the modes. This is a serious problem for QPO models based on inertial waves.

We have measured the angular diameters of six M dwarfs with the CHARA Array, a long-baseline optical interferometer located at Mount Wilson Observatory. Spectral types range from M1.0 V to M3.0 V and linear radii from 0.38 to 0.69 Rsun. These results are consistent with the seven other M-dwarf radii measurements from optical interferometry and with those for sixteen stars in eclipsing binary systems. We compare all directly measured M dwarf radii to model predictions and find that current models underestimate the true stellar radii by up to 15-20%. The differences are small among the metal-poor stars but become significantly larger with increasing metallicity. This suggests that theoretical models for low mass stars may be missing some opacity source that alters the computed stellar radii.

The effects of turbulence on the mixing of gases and dust in the outer Solar nebula are examined using 3-D MHD calculations in the shearing-box approximation with vertical stratification. The turbulence is driven by the magneto-rotational instability. The magnetic and hydrodynamic stresses in the turbulence correspond to an accretion time at the midplane about equal to the lifetimes of T Tauri disks, while accretion in the surface layers is thirty times faster. The mixing resulting from the turbulence is also fastest in the surface layers. The mixing rate is similar to the rate of radial exchange of orbital angular momentum, so that the Schmidt number is near unity. The vertical spreading of a trace species is well-matched by solutions of a damped wave equation when the flow is horizontally-averaged. The damped wave description can be used to inexpensively treat mixing in 1-D chemical models. However, even in calculations reaching a statistical steady state, the concentration at any given time varies substantially over horizontal planes, due to fluctuations in the rate and direction of the transport. In addition to mixing species that are formed under widely varying conditions, the turbulence intermittently forces the nebula away from local chemical equilibrium. The different transport rates in the surface layers and interior may affect estimates of the grain evolution and molecular abundances during the formation of the Solar system.

We present observations of the A7IV-V star Alderamin ($\alpha$ Cep, HR 8162, HD 203280) from the Georgia State University CHARA Array. These infrared interferometric angular size measurements indicate a non-circular projected disk brightness distribution for this known rapid rotator. The interferometric observations are modeled as arising from an elongated rigid atmosphere, with apparent polar and equatorial radii of $r_p=0.6753^{+0.0119}_{-0.0135}$ milliarcseconds (mas) and $r_e=0.8767^{+0.0293}_{-0.0183}$ mas, respectively, for a difference of $201\pm 32$ microarcseconds ($\mu$as), and with an axial ratio of $r_e/r_p =1.298\pm 0.051$. Using the Hipparcos distance of $14.96\pm 0.11$ pc, these angular measures translate to $2.18\pm0.05$ and $2.82\pm0.10$ $R_\odot$. The inclination of Alderamin to the line of sight indicated by this modeling is effectively edge-on ($i=88.2^{+1.8}_{-13.3}$). The star has a true rotational velocity of $283\pm 10$ km/s ($\sim83$% of breakup velocity), and a polar temperature of roughly 8400 K. Significantly, a necessary aspect of this modeling is a determination of the gravity darkening coefficient, which at a value of $\beta=0.084^{+0.026}_{-0.049}$ is consistent with a convective photosphere, as expected for an A7IV-V star. Our detailed characterization of this object allows us to investigate various scenarios for the angular momentum history of Alderamin and the appropriateness of certain stellar evolution models.

Cepheids play a key role in astronomy as standard candles for measuring intergalactic distances. Their distance is usually inferred from the Period-Luminosity relationship, calibrated using the semi-empirical Baade-Wesselink method. Using this method, the distance is known to a multiplicative factor, called the projection factor. Presently, this factor is computed using numerical models - it has hitherto never been measured directly. Based on our new interferometric measurements obtained with the CHARA Array and the already published parallax, we present a geometrical measurement of the projection factor of a Cepheid, delta Cep. The value we determined, p = 1.27$\pm$0.06, confirms the generally adopted value of p = 1.36 within 1.5 sigmas. Our value is in line with recent theoretical predictions of Nardetto et al. (2004).

The CHARA Array is a six 1-m telescope optical/IR interferometric array located on Mount Wilson California, designed and built by the Center for High Angular Resolution Astronomy of Georgia State University. In this paper we describe the main elements of the Array hardware and software control systems as well as the data reduction methods currently being used. Our plans for upgrades in the near future are also described.

We report on K-band interferometric observations of the bright, rapidly rotating star Regulus (type B7 V) made with the CHARA Array on Mount Wilson, California. Through a combination of interferometric and spectroscopic measurements, we have determined for Regulus the equatorial and polar diameters and temperatures, the rotational velocity and period, the inclination and position angle of the spin axis, and the gravity darkening coefficient. These first results from the CHARA Array provide the first interferometric measurement of gravity darkening in a rapidly rotating star and represent the first detection of gravity darkening in a star that is not a member of an eclipsing binary system.

We examine the effects of photon bubble instability in radiation-dominated accretion disks such as those found around black holes in active galactic nuclei and X-ray binary star systems. Two- and 3-D numerical radiation MHD calculations of small patches of disk are used. Modes with wavelengths shorter than the gas pressure scale height grow faster than the orbital frequency in the surface layers. The fastest growth rate observed is five times the orbital frequency and occurs on nearly-vertical magnetic fields. The spectrum of linear modes agrees with a WKB analysis indicating still faster growth at unresolved scales, with a maximum proportional to the gravity and inversely proportional to the gas sound speed. Disturbances reaching non-linear amplitudes steepen into trains of shocks similar to a 1-D periodic non-linear analytic solution. Variations in propagation speed result in merging of adjacent fronts, and over time the shock spacing and amplitude increase. Growth is limited by the strength of the field, and the structure is disrupted when the ram pressure exceeds the magnetic pressure. The largest horizontal density variations are similar to the ratio of magnetic to gas pressure, and in our calculations are more than 100. Under the conditions considered, radiation diffuses through the inhomogeneneous flow 5x faster than through the initial hydrostatic equilibrium, and the net cooling rate is several times greater than in a similar calculation without magnetic fields showing the effects of convection. These results indicate photon bubbles may be important in cooling radiation-dominated disks.

One of the principal scientific objectives of the upcoming Constellation-X mission is to attempt to map the inner regions of accretion disks around black holes in Seyfert galaxies by reverberation mapping of the Fe K fluorescence line. This area of the disk is likely radiation pressure dominated and subject to various dynamical instabilities. Here, we show that density inhomogeneities in the disk atmosphere resulting from the photon bubble instability (PBI) can cause rapid changes in the X-ray reflection features, even when the illuminating flux is constant. Using a simulation of the development of the PBI, we find that, for the disk parameters chosen, the Fe K and O VIII Ly\alpha lines vary on timescales as short as a few hundredths of an orbital time. In response to the changes in accretion disk structure, the Fe K equivalent width (EW) shows variations as large as ~100 eV. The magnitude and direction (positive or negative) of the changes depends on the ionization state of the atmosphere. The largest changes are found when the disk is moderately ionized. The O VIII EW varies by tens of eV, as well as exhibiting plenty of rapid, low-amplitude changes. This effect provides a natural explanation for some observed instances of short timescale Fe K variability which was uncorrelated with the continuum (e.g., Mrk 841). New predictions for Fe K reverberation mapping should be made which include the effects of this accretion disk driven line variability and a variable ionization state. Reflection spectra averaged over the evolution of the instability are well fit by constant density models in the 2-10 keV region.

We present calculations of non-LTE, relativistic accretion disk models applicable to the high/soft state of black hole X-ray binaries. We include the effects of thermal Comptonization and bound-free and free-free opacities of all abundant ion species. We present spectra calculated for a variety of accretion rates, black hole spin parameters, disk inclinations, and stress prescriptions. We also consider nonzero inner torques on the disk, and explore different vertical dissipation profiles, including some which are motivated by recent radiation MHD simulations of magnetorotational turbulence. Bound-free metal opacity generally produces significantly less spectral hardening than previous models which only considered Compton scattering and free-free opacity. It also tends to keep the effective photosphere near the surface, resulting in spectra which are remarkably independent of the stress prescription and vertical dissipation profile, provided little dissipation occurs above the effective photosphere. We provide detailed comparisons between our models and the widely used multicolor disk model. Frequency dependent discrepancies exist that may affect the parameters of other spectral components when this simpler disk model is used to fit modern X-ray data. For a given source, our models predict that the luminosity in the high/soft state should approximately scale with the fourth power of the empirically inferred maximum temperature, but with a slight hardening at high luminosities. This is in good agreement with observations. (abridged)

Spectroscopic studies indicate that gas in the photospheres of young O stars moves at speeds up to the sound speed. We show, using two-dimensional radiation MHD calculations and results from a local linear analysis, that the motions may be due to photon bubble instability if young O stars have magnetic fields.

The vertical structure of black hole accretion disks in which radiation dominates the total pressure is investigated using a three-dimensional radiation-MHD calculation. The domain is a small patch of disk centered 100 Schwarzschild radii from a black hole of 10^8 Solar masses, and the stratified shearing-box approximation is used. Magneto-rotational instability converts gravitational energy to turbulent magnetic and kinetic energy. The gas is heated by magnetic dissipation and by radiation damping of the turbulence, and cooled by diffusion and advection of radiation through the vertical boundaries. The resulting structure differs in several fundamental ways from the standard Shakura-Sunyaev picture. The disk consists of three layers. At the midplane, the density is large, and the magnetic pressure and total accretion stress are less than the gas pressure. In lower-density surface layers that are optically thick, the magnetic pressure and stress are greater than the gas pressure but less than the radiation pressure. Horizontal density variations in the surface layers exceed an order of magnitude. Magnetic fields in the regions of greatest stress are buoyant, and dissipate as they rise, so the heating rate declines more slowly with height than the stress. Much of the dissipation occurs at low column depth, and the interior is cooler and less radiation-dominated than in the Shakura-Sunyaev model with the same surface mass density and flux. The mean structure is convectively stable.

The saturation level of the magnetorotational instability (MRI) is investigated using three-dimensional MHD simulations. The shearing box approximation is adopted and the vertical component of gravity is ignored, so that the evolution of the MRI is followed in a small local part of the disk. We focus on the dependence of the saturation level of the stress on the gas pressure, which is a key assumption in the standard alpha disk model. From our numerical experiments it is found that there is a weak power-law relation between the saturation level of the Maxwell stress and the gas pressure in the nonlinear regime; the higher the gas pressure, the larger the stress. Although the power-law index depends slightly on the initial field geometry, the relationship between stress and gas pressure is independent of the initial field strength, and is unaffected by Ohmic dissipation if the magnetic Reynolds number is at least 10. The relationship is the same in adiabatic calculations, where pressure increases over time, and nearly-isothermal calculations, where pressure varies little with time. Our numerical results are qualitatively consistent with an idea that the saturation level of the MRI is determined by a balance between the growth of the MRI and the dissipation of the field through reconnection. The quantitative interpretation of the pressure-stress relation, however, may require advances in the theoretical understanding of non-steady magnetic reconnection.

Numerical calculations and linear theory of radiation magnetohydrodynamic flows indicate that the photon bubble and magnetorotational instability (MRI) may produce large density inhomogeneities in radiation pressure supported media. We study the effects of the photon bubble instability on accretion disk spectra using 2-D Monte Carlo (MC) and 1-D Feautrier radiative transfer calculations on a snapshot of a 2-D numerical simulation domain. We find an enhancement in the thermalization of the MC spectra over that of the Feautrier calculation. In the inner-most regions of these disks, the turbulent magnetic pressure may greatly exceed that of the gas. It is then possible for bulk turbulent Alfvenic motions driven by the MRI to exceed the thermal velocity making turbulent Comptonization the dominant radiative process. We estimate the spectral distortion due to turbulent Comptonization utilizing a 1-D MC calculation.

Numerical simulations of the interiors of radiation dominated accretion disks show that significant density inhomogeneities can be generated in the gas. Here, we present the first results of our study on X-ray reflection spectra from such heterogeneous density structures. We consider two cases: first, we produce a number of toy models where a sharp increase or decrease in density of variable width is placed at different depths in a uniform slab. Comparing the resulting reflection spectra to those from an unaltered slab shows that the inhomogeneity can affect the emission features, in particular the Fe K and O VIII Lyalpha lines. The magnitude of any differences depends on both the parameters of the density change and the ionizing power of the illuminating radiation, but the inhomogeneity is required to be within ~2 Thomson depths of the surface to cause an effect. However, only relatively small variations in density (on the order of a few) are necessary for significant changes in the reflection features to be possible. Our second test was to compute reflection spectra from the density structure predicted by a simulation of the non-linear outcome of the photon bubble instability. The resulting spectra also exhibited differences from the constant density models, caused primarily by a strong 6.7 keV iron line. Nevertheless, constant density models can provide a good fit to simulated spectra, albeit with a low reflection fraction, between 2 and 10 keV. Below 2 keV, differences in the predicted soft X-ray line emission result in very poor fits with a constant density ionized disk model. The results indicate that density inhomogeneities may further complicate the relationship between the Fe K equivalent width and the X-ray continuum.

We examine the small-scale dynamics of black hole accretion disks in which radiation pressure exceeds gas pressure. Local patches of disk are modeled by numerically integrating the equations of radiation MHD in the flux-limited diffusion approximation. The shearing-box approximation is used, and the vertical component of gravity is neglected. Magneto-rotational instability (MRI) leads to turbulence in which accretion stresses are due primarily to magnetic torques. When radiation is locked to gas over the length and time scales of fluctuations in the turbulence, the accretion stress, density contrast, and dissipation differ little from those in the corresponding calculations with radiation replaced by extra gas pressure. However, when radiation diffuses each orbit a distance that is comparable to the RMS vertical wavelength of the MRI, radiation pressure is less effective in resisting squeezing. Large density fluctuations occur, and radiation damping of compressive motions converts PdV work into photon energy. The accretion stress in calculations having a net vertical magnetic field is found to be independent of opacity over the range explored, and approximately proportional to the square of the net field. In calculations with zero net magnetic flux, the accretion stress depends on the portion of the total pressure that is effective in resisting compression. The stress is lower when radiation diffuses rapidly with respect to the gas. We show that radiation-supported Shakura-Sunyaev disks accreting via internal magnetic stresses are likely in their interiors to have radiation marginally coupled to turbulent gas motions.

We perform numerical simulations of magneto-rotational instability in a local patch of accretion disk in which radiation pressure exceeds gas pressure. Such conditions may occur in the central regions of disks surrounding compact objects in active galactic nuclei and Galactic X-ray sources. We assume axisymmetry, and neglect vertical stratification. The growth rates of the instability on initially uniform magnetic fields are consistent with the linear analysis of Blaes & Socrates 2001. As is the case when radiation effects are neglected, the non-linear development of the instability leads to transitory turbulence when the initial magnetic field has no net vertical flux. During the turbulent phase, angular momentum is transported outwards. The Maxwell stress is a few times the Reynolds stress, and their sum is about four times the mean pressure in the vertical component of the magnetic field. For magnetic pressure exceeding gas pressure, turbulent fluctuations in the field produce density contrasts about equal to the ratio of magnetic to gas pressure. These are many times larger than in the corresponding gas pressure dominated situation, and may have profound implications for the steady-state vertical structure of radiation-dominated disks. Diffusion of radiation from compressed regions damps turbulent motions, converting kinetic energy into photon energy.

The standard equilibrium for radiation-dominated accretion disks has long been known to be viscously, thermally, and convectively unstable, but the nonlinear development of these instabilities---hence the actual state of such disks---has not yet been identified. By performing local two-dimensional hydrodynamic simulations of disks, we demonstrate that convective motions can release heat sufficiently rapidly as to substantially alter the vertical structure of the disk. If the dissipation rate within a vertical column is proportional to its mass, the disk settles into a new configuration thinner by a factor of two than the standard radiation-supported equilibrium. If, on the other hand, the vertically-integrated dissipation rate is proportional to the vertically-integrated total pressure, the disk is subject to the well-known thermal instability. Convection, however, biases the development of this instability toward collapse. The end result of such a collapse is a gas pressure-dominated equilibrium at the original column density.

A module for the ZEUS-2D code is described which may be used to solve the equations of radiation hydrodynamics to order unity in v/c, in the flux-limited diffusion (FLD) approximation. In this approximation, the tensor Eddington factor f which closes the radiation moment equations is chosen to be an empirical function of radiation energy density. This is easier to implement and faster than full-transport techniques, in which f is computed by solving the transfer equation. However, FLD is less accurate when the flux has a component perpendicular to the gradient in radiation energy density, and in optically thin regions when the radiation field depends strongly on angle. The material component of the fluid is here assumed to be in local thermodynamic equilibrium. The energy equations are operator-split, with transport terms, radiation diffusion term, and other source terms evolved separately. Transport terms are applied using the same consistent transport algorithm as in ZEUS-2D. The radiation diffusion term is updated using an alternating-direction implicit method with convergence checking. Remaining source terms are advanced together implicitly using numerical root-finding. However when absorption opacity is zero, accuracy is improved by treating compression and expansion source terms using time-centered differencing. Results are discussed for test problems including radiation-damped linear waves, radiation fronts propagating in optically-thin media, subcritical and supercritical radiating shocks, and an optically-thick shock in which radiation dominates downstream pressure.