High resolution spectra of typical wind diagnostics ([OI] 6300 \r{A} and other forbidden emission lines) can often be decomposed into multiple components: high-velocity components with blueshifts up to several 100 km/s are usually attributed to fast jets, while narrow (NLVC) and broad (BLVC) low-velocity components are believed to trace slower disc winds. Under the assumption that the line-broadening is dominated by Keplerian rotation, several studies have found that the BLVCs should trace gas launched between 0.05 and 0.5 au and correlations between the properties of BLVCs and NLVCs have been interpreted as evidence for the emission tracing an extended MHD wind and not a photoevaporative wind. We calculated synthetic line profiles obtained from detailed photoionisation calculations of an X-ray photoevaporation model and a simple MHD wind model and analyzed the emission regions of different diagnostic lines and the resulting spectral profiles. The photoevaporation model reproduces most of the observed NLVCs but not the BLVCs or HVCs. The MHD model is able to reproduce all components but produces Keplerian double peaks at average inclinations that are rarely observed. The combination of MHD and photoevaporative winds could solve this problem. Our results suggest that the Gaussian decomposition does not allow for a clear distinction of flux from different wind regions and that the line broadening is often dominated by the velocity gradient in the outflow rather than by Keplerian rotation. We show that observed correlations between BLVC and NLVC do not necessarily imply a common origin in an extended MHD wind.

FU Orionis objects (FUors) are rapidly-accreting, pre-main sequence objects that are known to exhibit large outbursts at optical and near-infrared wavelengths, with post-eruption, small-scale photometric variability superimposed on longer-term trends. In contrast, little is known about the variability of FUors at longer wavelengths. To explore this further, we observed six FUor objects using the NOrthern Extended Millimeter Array (NOEMA) and for a subset of three objects we obtained coordinated observations with NOEMA and the Lowell Discovery Telescope (LDT). In combination with previously published NOEMA observations from 2014, our 2017 observations of V1735 Cyg provide the first detection of variability in an FUor object at 2.7 mm. In the absence of significant optical variability, we discount the possibility that the mm flux density changed as a result of irradiation from the central disk. In addition, a change in the dust mass due to infall is highly unlikely. A plausible explanation for the change in 2.7 mm flux density is variability in free-free emission due to changes in the object's jet/wind. Thus, it may be that free-free emission in some FUor objects is significant at $\sim$3 mm and must be considered when deriving disk masses in order to help constrain the mechanism responsible for triggering FUor outbursts.

In young multiple stellar systems the gravitational influence of the stars shapes the circumstellar disk, controlling accretion and the material available for planet formation. Our observations of the triple system GW Orionis provide direct evidence for disk tearing, where the gravitational torques due to misaligned disk/orbit planes warp the disk and break the inner disk into precessing rings. We image an eccentric ring that is misaligned with respect to the orbital planes and outer disk, and casts shadows on the strongly warped intermediate disk. Our ring/warp geometry constraints and the fully characterized perturber orbits make the system a potential Rosetta Stone for studying disk hydrodynamics. The ring might offer suitable conditions for planet formation, providing a mechanism for forming wide-separation planets on highly oblique orbits.

Building on our previous hydrodynamic study of the angular momenta of cloud cores formed during gravitational collapse of star-forming molecular gas in our previous work, we now examine core properties assuming ideal magnetohydrodynamics (MHD). Using the same sink-patch implementation for the \emph{Athena} MHD code, we characterize the statistical properties of cores, including the mass accretion rates, specific angular momenta, and alignments between the magnetic field and the spin axis of the core on the $0.1 \ \mathrm{pc}$ scale. Our simulations, which reproduce the observed relation between magnetic field strength and gas density, show that magnetic fields can help collimate low density flows and help seed the locations of filamentary structures. Consistent with our previous purely hydrodynamic simulations, stars (sinks) form within the heterogeneous environments of filaments, such that accretion onto cores is highly episodic leading to short-term variability but no long-term monotonic growth of the specific angular momenta. With statistical characterization of protostellar cores properties and behaviors, we aim to provide a starting point for building more realistic and self-consistent disk formation models, helping to address whether magnetic fields can prevent the development of (large) circumstellar disks in the ideal MHD limit.

We are undertaking a large survey of over thirty disks using the Gemini Planet Imager (GPI) to see whether the observed dust structures match spectral energy distribution (SED) predictions and have any correlation with stellar properties. GPI can observe near-infrared light scattered from dust in circumstellar environments using high-resolution Polarimetric Differential Imaging (PDI) with coronagraphy and adaptive optics. The data have been taken in J and H bands over two years, with inner working angles of 0.08'' and 0.11'' respectively. Ahead of the release of the complete survey results, here we present five objects with extended and irregular dust structures within 2'' of the central star. These objects are: FU Ori; MWC 789; HD 45677; Hen 3-365; and HD 139614. The observed structures are consistent with each object being a pre-main-sequence star with protoplanetary dust. The five objects' circumstellar environments could result from extreme youth and complex initial conditions, from asymmetric scattering patterns due to shadows cast by misaligned disks, or in some cases from interactions with companions. We see complex U_phi structures in most objects that could indicate multiple scattering or result from the illumination of companions. Specific key findings include the first high-contrast observation of MWC 789 revealing a newly-discovered companion candidate and arc, and two faint companion candidates around Hen 3-365. These two objects should be observed further to confirm whether the companion candidates are co-moving. Further observations and modeling are required to determine the causes of the structures.

This textbook gives a perspective of heliophysics in a way that emphasizes universal processes from a perspective that draws attention to what provides Earth (and similar (exo-)planets) with a relatively stable setting in which life as we know it can thrive. The book is intended for students in physical sciences in later years of their university training and for beginning graduate students in fields of solar, stellar, (exo-)planetary, and planetary-system sciences.

We investigate the deuteration of methanol towards the high-mass star forming region NGC 7538-IRS1. We have carried out a multi-transition study of CH$_3$OH, $^{13}$CH$_3$OH and of the deuterated fllavors, CH$_2$DOH and CH$_3$OD, between 1.0--1.4 mm with the IRAM-30~m antenna. In total, 34 $^{13}$CH$_3$OH, 13 CH$_2$DOH lines and 20 CH$_3$OD lines spanning a wide range of upper-state energies (E$_{up}$) were detected. From the detected transitions, we estimate that the measured D/H does not exceed 1$\%$, with a measured CH$_2$DOH/CH$_3$OH and CH$_3$OD/CH$_3$OH of about (32$\pm$8)$\times$10$^{-4}$ and (10$\pm$4)$\times$10$^{-4}$, respectively. This finding is consistent with the hypothesis of a short-time scale formation during the pre-stellar phase. We find a relative abundance ratio CH$_2$DOH/CH$_3$OD of 3.2 $\pm$ 1.5. This result is consistent with a statistical deuteration. We cannot exclude H/D exchanges between water and methanol if water deuteration is of the order 0.1$\%$, as suggested by recent Herschel observations.

We present radiation-magnetohydrodynamic simulations aimed at studying evolutionary properties of H\,{\normalsize II} regions in turbulent, magnetised, and collapsing molecular clouds formed by converging flows in the warm neutral medium. We focus on the structure, dynamics and expansion laws of these regions. Once a massive star forms in our highly structured clouds, its ionising radiation eventually stops the accretion (through filaments) toward the massive star-forming regions. The new over-pressured H\,{\normalsize II} regions push away the dense gas, thus disrupting the more massive collapse centres. Also, because of the complex density structure in the cloud, the H\,{\normalsize II} regions expand in a hybrid manner: they virtually do not expand toward the densest regions (cores), while they expand according to the classical analytical result towards the rest of the cloud, and in an accelerated way, as a blister region, towards the diffuse medium. Thus, the ionised regions grow anisotropically, and the ionising stars generally appear off-centre of the regions. Finally, we find that the hypotheses assumed in standard H\,{\normalsize II}-region expansion models (fully embedded region, blister-type, or expansion in a density gradient) apply simultaneously in different parts of our simulated H\,{\normalsize II} regions, producing a net expansion law ($R \propto t^\alpha$, with $\alpha$ in the range of 0.93-1.47 and a mean value of $1.2 \pm 0.17$) that differs from any of those of the standard models.

We present the results of a suite of numerical simulations designed to explore the origin of the angular momenta of protostellar cores. Using the hydrodynamic grid code \emph{Athena} with a sink implementation, we follow the formation of protostellar cores and protostars (sinks) from the subvirial collapse of molecular clouds on larger scales to investigate the range and relative distribution of core properties. We find that the core angular momenta are relatively unaffected by large-scale rotation of the parent cloud; instead, we infer that angular momenta are mainly imparted by torques between neighboring mass concentrations and exhibit a log-normal distribution. Our current simulation results are limited to size scales $\sim 0.05$~pc ($\sim 10^4 \rm AU$), but serve as first steps toward the ultimate goal of providing initial conditions for higher-resolution studies of core collapse to form protoplanetary disks.

We present multi-epoch continuum observations of the Class 0 protostellar system IRAS 16293-2422 taken with the Very Large Array (VLA) at multiple wavelengths between 7 mm and 15 cm (41 GHz down to 2 GHz), as well as single-epoch Atacama Large Millimeter/submillimeter Array (ALMA) continuum observations covering the range from 0.4 to 1.3 mm (700 GHz down to 230 GHz). The new VLA observations confirm that source A2 is a protostar driving episodic mass ejections, and reveal the complex relative motion between A2 and A1. The spectrum of component B can be described by a single power law ($S_\nu \propto \nu^{2.28}$) over the entire range from 3 to 700 GHz (10 cm down to 0.4 mm), suggesting that the emission is entirely dominated by dust even at $\lambda$ = 10 cm. Finally, the size of source B appears to increase with frequency up to 41 GHz, remaining roughly constant (at $0''.39$ $\equiv$ 55 AU) at higher frequencies. We interpret this as evidence that source B is a dusty structure of finite size that becomes increasingly optically thick at higher frequencies until, in the millimeter regime, the source becomes entirely optically thick. The lack of excess free-free emission at long wavelengths, combined with the absence of high-velocity molecular emission indicates that source B does not drive a powerful outflow, and might indicate that source B is at a particularly early stage of its evolution.

We present the first images of the transition disk around the close binary system HD 34700A in polarized scattered light using the Gemini Planet Imager instrument on Gemini South. The J and H band images reveal multiple spiral-arm structures outside a large (R = 0.49" = 175 au) cavity along with a bluish spiral structure inside the cavity. The cavity wall shows a strong discontinuity and we clearly see significant non-azimuthal polarization Uphi consistent with multiple scattering within a disk at an inferred inclination ~42deg. Radiative transfer modeling along with a new Gaia distance suggest HD 37400A is a young (~5 Myr) system consisting of two intermediate-mass (~2Msun) stars surrounded by a transitional disk and not a solar-mass binary with a debris disk as previously classified. Conventional assumptions of the dust-to-gas ratio would rule out a gravitational instability origin to the spirals while hydrodynamical models using the known external companion or a hypothetical massive protoplanet in the cavity both have trouble reproducing the relatively large spiral arm pitch angles (~30deg) without fine tuning of gas temperature. We explore the possibility that material surrounding a massive protoplanet could explain the rim discontinuity after also considering effects of shadowing by an inner disk. Analysis of archival Hubble Space Telescope data suggests the disk is rotating counterclockwise as expected from the spiral arm structure and revealed a new low-mass companion at 6.45" separation. We include an appendix which sets out clear definitions of Q, U, Qphi, Uphi, correcting some confusion and errors in the literature.

We present Gaia-DR2 astrometry of a sample of YSO candidates in Ophiuchus, Serpens Main and Serpens South/W40 in the Aquila Rift, which had been mainly identified by their infrared excess with Spitzer. We compare the Gaia-DR2 parallaxes against published and new parallaxes obtained from our Very Long Baseline Array (VLBA) program GOBELINS. We obtain consistent results between Gaia and the VLBA for the mean parallaxes in each of the regions analyzed here. We see small offsets, when comparing mean values, of a few tens of micro-arcseconds in the parallaxes, which are either introduced by the Gaia zero-point error or due to a selection effect by Gaia toward the brightest, less obscured stars. Gaia-DR2 data alone conclusively places Serpens Main and Serpens South at the same distance, as we first inferred from VLBA data alone in a previous publication. Thus, Serpens Main, Serpens South and W40 are all part of the same complex of molecular clouds, located at a mean distance of 436+/-9 pc. In Ophiuchus, both Gaia and VLBA suggest a small parallax gradient across the cloud, and the distance changes from 144.2+/-1.3 pc to 138.4+/-2.6 pc when going from L1689 to L1688.

We present disc wind models aimed at reproducing the main features of the strong Na I resonance line P-Cygni profiles in the rapidly-accreting pre-main sequence FU Ori objects. We conducted Monte Carlo radiative transfer simulations for a standard magnetocentrifugally driven wind (MHD) model and our own "Genwind" models, which allows for a more flexible wind parameterisation. We find that the fiducial MHD wind and similar Genwind models, which have flows emerging outward from the inner disc edge, and thus have polar cavities with no absorbing gas, cannot reproduce the deep, wide Na I absorption lines in FU Ori objects viewed at low inclination. We find that it is necessary to include an "inner wind" to fill this polar cavity to reproduce observations. In addition, our models assuming pure scattering source functions in the Sobolev approximation at intermediate viewing angles ($30^{\circ} \lesssim i \lesssim 60^{\circ}$) do not yield sufficiently deep line profiles. Assuming complete absorption yields better agreement with observations, but simple estimates strongly suggest that pure scattering should be a much better approximation. The discrepancy may indicate that the Sobolev approximation is not applicable, possibly due to turbulence or non-monotonic velocity fields; there is some observational evidence for the latter. Our results provide guidance for future attempts to constrain FU Ori wind properties using full MHD wind simulations, by pointing to the importance of the boundary conditions necessary to give rise to an inner wind, and by suggesting that the winds must be turbulent to produce sufficiently deep line profiles.

Episodic accretion may be a common occurrence in the evolution of young pre-main sequence stars and has important implications for our understanding of star and planet formation. Many fundamental aspects of what drives the accretion physics, however, are still unknown. The ngVLA will be a key tool in understanding the nature of these events. The high spatial resolution, broad spectral coverage, and unprecedented sensitivity will allow for the detailed analysis of outburst systems. The proposed frequency range of the ngVLA allows for observations of the gas, dust, and non-thermal emission from the star and disk.

Numerical simulations of star formation have found that a power-law mass function can develop at high masses. In a previous paper, we employed isothermal simulations which created large numbers of sinks over a large range in masses to show that the power law exponent of the mass function, $dN/d\log M \propto M^{\Gamma}$, asymptotically and accurately approaches $\Gamma = -1.$ Simple analytic models show that such a power law can develop if the mass accretion rate $\dot{M} \propto M^2$, as in Bondi-Hoyle accretion; however, the sink mass accretion rates in the simulations show significant departures from this relation. In this paper we show that the expected accretion rate dependence is more closely realized provided the gravitating mass is taken to be the sum of the sink mass and the mass in the near environment. This reconciles the observed mass functions with the accretion rate dependencies, and demonstrates that power-law upper mass functions are essentially the result of gravitational focusing, a mechanism present in, for example, the competitive accretion model.

We use numerical simulations of molecular cloud formation in the colliding flow scenario to investigate the reliability of observational estimates of the angular momenta of early-state, low-mass protostellar cores. We show that, with suitable corrections for projection factors, molecular line observations of velocity gradients in NH$_3$ can be used to provide reasonable estimates of core angular momenta within a factor of two to three, with a few large underestimates due to unfavorable viewing angles. Our results differ from previous investigations which suggested that observations might overestimate true angular momenta by as much as an order of magnitude; the difference is probably due to the much smoother velocity field on small scales in our simulations, which result from allowing turbulence to decay and gravitational infall to dominate. The results emphasize the importance of understanding the nature of 'turbulent' velocities, with implications for the formation of protostellar disks during core collapse.

We derive the distance and structure of the Perseus molecular cloud by combining trigonometric parallaxes from Very Long Baseline Array (VLBA) observations, taken as part of the GOBELINS survey, and Gaia Data Release 2. Based on our VLBA astrometry, we obtain a distance of 321+/-10 pc for IC 348. This is fully consistent with the mean distance of 320+/-26 measured by Gaia. The VLBA observations toward NGC 1333 are insufficient to claim a successful distance measurement to this cluster. Gaia parallaxes, on the other hand, yield a mean distance of 293+/-22 pc. Hence, the distance along the line of sight between the eastern and western edges of the cloud is ~30 pc, which is significantly smaller than previously inferred. We use Gaia proper motions and published radial velocities to derive the spatial velocities of a selected sample of stars. The average velocity vectors with respect to the LSR are (u,v,w) = (-6.1+/-1.6, 6.8+/-1.1, -0.9+/-1.2) and (-6.4+/-1.0, 2.1+/-1.4, -2.4+/-1.0) km/s for IC 348 and NGC 1333, respectively. Finally, our analysis of the kinematics of the stars has shown that there is no clear evidence of expansion, contraction, or rotational motions within the clusters.

We present new trigonometric parallaxes and proper motions of young stellar objects in the Taurus molecular cloud complex from observations collected with the Very Long Baseline Array as part of the Gould's Belt Distances Survey (GOBELINS). We detected 26 young stellar objects and derived trigonometric parallaxes for 18 stars with an accuracy of 0.3$\%$ to a few percent. We modeled the orbits of six binaries and determined the dynamical masses of the individual components in four of these systems (V1023 Tau, T Tau S, V807 Tau and V1000 Tau). Our results are consistent with the first trigonometric parallaxes delivered by the Gaia satellite and reveal the existence of significant depth effects. We find that the central portion of the dark cloud Lynds 1495 is located at $d=129.5\pm 0.3$ pc while the B 216 clump in the filamentary structure connected to it is at $d=158.1\pm 1.2$ pc. The closest and remotest stars in our sample are located at $d=126.6\pm 1.7$ pc and $d=162.7\pm 0.8$ pc yielding a distance difference of about 36 pc. We also provide a new distance estimate for HL Tau that was recently imaged. Finally, we compute the spatial velocity of the stars with published radial velocity and investigate the kinematic properties of the various clouds and gas structures in this region.

We present results of our large scale, optical, multi-epoch photometric survey of ~180 square degrees across the Orion OB1 association, complemented with extensive follow up spectroscopy. We map and characterize in an uniform way the off-cloud, low-mass, pre-main sequence populations. We report 2064, mostly K and M-type, confirmed T Tauri members. Most (59%) are located in the OB1a subassociation, 27% in OB1b, and 14% within the confines of the A and B molecular clouds. There is significant structure in the spatial distribution of the young stars. We characterize two new clusterings of T Tauri stars, HD 35762 and HR 1833 in the OB1a subassociation, and two stellar overdensities in OB1b. There is indication of two populations of young stars in the OB1b region, located at two different distances, which may be due to the OB1a subassociation overlapping on front of OB1b. The various groups and regions can be ordered in an age sequence that agrees with the long standing picture of star formation starting in Orion OB1a some 10-15 Myr ago. We define a new type of T Tauri star, the C/W class, objects we propose may be nearing the end of their accretion phase. We detect the observational signature of Li depletion in young K and M stars with a timescale of 8.5 Myr. The decline of the accretion fraction from ~2 - 10 Myr, implies an accretion e-folding timescale of 2.1 Myr. Finally, the median amplitude of the V-band variability shows the decline of stellar activity, from accreting Classical T Tauri stars to the least active field dwarfs.

We present observations of VLBA 20, a radio source found towards the edge of the Orion Nebula Cluster (ONC). Nonthermal emission dominates the spectral energy distribution of this object from radio to mid-infrared regime, suggesting that VLBA 20 is extragalactic. This source is heavily scattered in the radio regime. Very Long Baseline Array observations resolve it to ~34x19 mas at 5 GHz, and the wavelength dependence of the scattering disk is consistent with $\nu^{-2}$ at other frequencies. The origin of the scattering is most likely the ionized X-ray emitting gas from the winds of the most massive stars of the ONC. The scattering is highly anisotropic, with the axis ratio of 2:1, higher than what is typically observed towards other sources.

As spiral waves driven by a planet in a gaseous disk steepen into a shock, they deposit angular momentum, opening a gap in the disk. This has been well studied using both linear theory and numerical simulations, but so far, only for the primary spiral arm -- the one directly attached to the planet. Using two-dimensional hydrodynamic simulations, we show that the secondary and tertiary arms driven by a planet can also open gaps as they steepen into shocks. The depths of the secondary/tertiary gaps in surface density grow with time in a low viscosity disk ($\alpha = 5\times10^{-5}$), so even low-mass planets (e.g., super-Earth or mini-Neptune-mass) embedded in the disk can open multiple observable gaps, provided that sufficient time has passed. Applying our results to the HL Tau disk, we show that a single 30 Earth-mass planet embedded in the ring at 68.8~au (B5) can reasonably well reproduce the positions of the two major gaps at 13.2 and 32.3 au (D1 and D2), and roughly reproduce two other major gaps at 64.2 and 74.7 au (D5 and D6) seen in the mm continuum. The positions of secondary/tertiary gaps are found to be sensitive to the planetary mass and the disk temperature profile, so with accurate observational measurements of the temperature structure the positions of multiple gaps can be used to constrain the mass of the planet. We also comment on the gaps seen in the TW Hya and HD~163296 disk.

We conjecture that observed protoplanetary disc accretion rates may be explained with low viscosities which could be the result of hydrodynamic turbulence. We show that viscosities parameterized in the usual way with $\alpha \gtrsim 10^{-4}$, comparable to values suggested for hydrodynamic turbulence, can explain the observed accretion rates and lifetimes with plausible inner disc surface densities. Our models are also in better agreement with surface density estimates of the minimum mass solar nebula than models with rapid transport for a given mass accretion rate, such as recent models of accretion driven by magnetic winds. The required surface densities are a natural result of the protostellar infall phase, as long as non-gravitational transport is limited. We argue that, in addition to possible non-ideal magnetic transport due to disc winds possibly modified by the Hall effect, the effects of low-viscosity hydrodynamic accretion deserve more consideration.

The origin of the observed morphological and kinematic substructure of young star forming regions is a matter of debate. We offer a new analysis of data from simulations of globally gravitationally collapsing clouds of progenitor gas to answer questions about sub-structured star formation in the context of cold collapse. As a specific example, we compare our models to recent radial velocity survey data from the IN-SYNC survey of Orion and new observations of dense gas kinematics, and offer possible interpretations of kinematic and morphological signatures in the region. In the context of our model, we find the frequently-observed hub-filament morphology of the gas naturally arises during gravitational evolution, as well as the dynamically-distinct kinematic substructure of stars. We emphasize that the global and not just the local gravitational potential plays an important role in determining the dynamics of both clusters and filaments.

We analyse the morphology and kinematics of dense filamentary structures produced in a numerical simulation of a star--forming cloud of $1.4 \times 10^4 \, \Msun$ evolving under their own self--gravity in a magnetized medium. This study is motivated by recent observations of velocity--coherent substructures ("fibres") in star-forming filaments. We find such "fibres" ubiquitously in our simulated filament. We found that a fibre in one projection is not necessarily a fibre in another projection, and thus, caution should be taken into account when considering them as real objects. We found that only the densest parts of the filament ($\sim$30\% of the densest volume, which contains $\sim$70\% of the mass) belong to fibres in 2 projections. Moreover, it is quite common that they are formed by separated density enhancements superposed along the line of sight. Observations of fibres can yield insight into the level of turbulent substructure driven by gravity, but care should be taken in interpreting the results given the problem of line of sight superposition. We also studied the morphology and kinematics of the 3D skeleton (spine), finding that subfilaments accrete structured material mainly along the magnetic field lines, which are preferentially perpendicular to the skeleton. The magnetic field is at the same time dragged by the velocity field due to the gravitational collapse.

We review the general theoretical concepts and observational constraints on the distribution and evolution of water vapor and ice in protoplanetary disks, with a focus on the Solar System. Water is expected to freeze out at distances greater than 1-3 AU from solar-type central stars; more precise estimates are difficult to obtain due to uncertainties in the complex processes involved in disk evolution, including dust growth, settling, and radial drift, and the level of turbulence and viscous dissipation within disks. Interferometric observations are now providing constraints on the positions of CO snow lines, but extrapolation to the unresolved regions where water ice sublimates will require much better theoretical understanding of mass and angular momentum transport in disks as well as more refined comparison of observations with sophisticated disk models.

We present results from high resolution optical spectra toward 66 young stars in the Orion B molecular cloud to study their kinematics and other properties. Observations of the H$\alpha$ and Li I 6707 {\AA} lines are used to check membership and accretion properties. While the stellar radial velocities of in NGC 2068 and L1622 show good agreement with that of the molecular gas, many of the stars in NGC 2024 show a considerable offset. This could be a signature of either expansion of the cluster, high degree of the ejection of the stars from the cluster through the dynamical interaction, or the acceleration of the gas due to stellar feedback.

NGC 1980 is a young cluster that is located about 0.5 degrees south of the Orion Nebula Cluster (ONC). Recent studies by Bouy et al. and Pillitteri et al. have suggested that NGC 1980 contains an older population of stars compared to a much younger ONC, and that it belongs to a foreground population that may be located in front of the Orion A molecular gas by as much as 40 pc. In this work we present low-resolution spectra towards 148 young stars found towards the NGC 1980 region. We determine the spectral types of these stars, examine accretion signatures and measure the extinction towards them. We determine that based on these observations, the age of the population of NGC 1980 is indistinguishable from L1641, estimated to be ~3 Myr, comparable with the study by Fang et al.

The Herschel Orion Protostar Survey obtained well-sampled 1.2 - 870 micron spectral energy distributions (SEDs) of over 300 protostars in the Orion molecular clouds, home to most of the young stellar objects (YSOs) in the nearest 500 pc. We plot the bolometric luminosities and temperatures for 330 Orion YSOs, 315 of which have bolometric temperatures characteristic of protostars. The histogram of bolometric temperature is roughly flat; 29% of the protostars are in Class 0. The median luminosity decreases by a factor of four with increasing bolometric temperature; consequently, the Class 0 protostars are systematically brighter than the Class I protostars, with a median luminosity of 2.3 L_sun as opposed to 0.87 L_sun. At a given bolometric temperature, the scatter in luminosities is three orders of magnitude. Using fits to the SEDs, we analyze how the luminosities corrected for inclination and foreground reddening relate to the mass in the inner 2500 AU of the best-fit model envelopes. The histogram of envelope mass is roughly flat, while the median corrected luminosity peaks at 15 L_sun for young envelopes and falls to 1.7 L_sun for late-stage protostars with remnant envelopes. The spread in luminosity at each envelope mass is three orders of magnitude. Envelope masses that decline exponentially with time explain the flat mass histogram and the decrease in luminosity, while the formation of a range of stellar masses explains the dispersion in luminosity.

In order to look for signs of on-going planet formation in young disks, we carried out the first J-band polarized emission imaging of the Herbig Ae/Be stars HD 150193, HD 163296, and HD 169142 using the Gemini Planet Imager (GPI), along with new H band observations of HD 144432. We confirm the complex "double ring" structure for the nearly face-on system HD 169142 first seen in H-band, finding the outer ring to be substantially redder than the inner one in polarized intensity. Using radiative transfer modeling, we developed a physical model that explains the full spectral energy distribution (SED) and J- and H-band surface brightness profiles, suggesting that the differential color of the two rings could come from reddened starlight traversing the inner wall and may not require differences in grain properties. In addition, we clearly detect an elongated, off-center ring in HD 163296 (MWC 275), locating the scattering surface to be 18 AU above the midplane at a radial distance of 77 AU, co-spatial with a ring seen at 1.3mm by ALMA linked to the CO snow line. Lastly, we report a weak tentative detection of scattered light for HD 150193 (MWC 863) and a non-detection for HD 144432; the stellar companion known for each of these targets has likely disrupted the material in the outer disk of the primary star. For HD 163296 and HD 169142, the prominent outer rings we detect could be evidence for giant planet formation in the outer disk or a manifestation of large-scale dust growth processes possibly related to snow-line chemistry.

We discuss the possibility that gravitational focusing, is responsible for the power-law mass function of star clusters $N(\log M) \propto M^{-1}$. This power law can be produced asymptotically when the mass accretion rate of an object depends upon the mass of the accreting body as $\dot{M} \propto M^2$. While Bondi-Hoyle-Littleton accretion formally produces this dependence on mass in a uniform medium, realistic environments are much more complicated. However, numerical simulations in SPH allowing for sink formation yield such an asymptotic power-law mass function. We perform pure N-body simulations to isolate the effects of gravity from those of gas physics and to show that clusters naturally result with the power-law mass distribution. We also consider the physical conditions necessary to produce clusters on appropriate timescales. Our results help support the idea that gravitationally-dominated accretion is the most likely mechanism for producing the cluster mass function.

We present the first results of the Gould's Belt Distances Survey (GOBELINS), a project aimed at measuring the proper motion and trigonometric parallax of a large sample of young stars in nearby regions using multi-epoch Very Long Baseline Array (VLBA) radio observations. Enough VLBA detections have now been obtained for 16 stellar systems in Ophiuchus to derive their parallax and proper motion. This leads to distance determinations for individual stars with an accuracy of 0.3 to a few percent. In addition, the orbits of 6 multiple systems were modelled by combining absolute positions with VLBA (and in some cases, near infrared) angular separations. Twelve stellar systems are located in the dark cloud Lynds 1688, the individual distances for this sample are highly consistent with one another, and yield a mean parallax for Lynds 1688 of $\varpi=7.28\pm0.06$ mas, corresponding to a distance $d=137.3\pm1.2$ pc. This represents an accuracy better than 1%. Three systems for which astrometric elements could be measured are located in the eastern streamer (Lynds 1689) and yield an estimate of $\varpi=6.79\pm0.16$ mas, corresponding to a distance $d=147.3\pm3.4$ pc. This suggests that the eastern streamer is located about 10 pc farther than the core, but this conclusion needs to be confirmed by observations (currently being collected) of additional sources in the eastern streamer. From the measured proper motions, we estimate the one-dimensional velocity dispersion in Lynds 1688 to be 2.8$\pm$1.8 and 3.0$\pm$2.0 ${\rm km~s}^{-1}$, in R.A. and DEC., respectively, these are larger than, but still consistent within $1\sigma$, with those found in other studies.

Using multi-epoch VLA observations, covering a time baseline of 29.1 years, we have measured the proper motions of 88 young stars with compact radio emission in the core of the Orion Nebula Cluster (ONC) and the neighboring BN/KL region. Our work increases the number of young stars with measured proper motion at radio frequencies by a factor of 2.5 and enables us to perform a better statistical analysis of the kinematics of the region than was previously possible. Most stars (79 out of 88) have proper motions consistent with a Gaussian distribution centered on $\overline{\mu_{\alpha}\cos{\delta}}=1.07\pm0.09\quad{\rm mas\,yr^{-1}}$, and $\overline{\mu_{\delta}}=-0.84\pm0.16\quad{\rm mas\,yr^{-1}}$, with velocity dispersions of $\sigma_{\alpha}=1.08\pm0.07\quad{\rm mas\,\,yr^{-1}},$ $\sigma_{\delta}=1.27\pm0.15\quad{\rm mas\,\,yr^{-1}}$. We looked for organized movements of these stars but found no clear indication of radial expansion/contraction or rotation. The remaining nine stars in our sample show peculiar proper motions that differ from the mean proper motions of the ONC by more than 3-$\sigma$. One of these stars, V 1326 Ori, could have been expelled from the Orion Trapezium 7,000 years ago. Two could be related to the multi-stellar disintegration in the BN/KL region, in addition to the previously known sources BN, I and n. The others either have high uncertainties (so their anomalous proper motions are not firmly established) or could be foreground objects.

We have recently shown that spiral density waves propagating in accretion disks can undergo a parametric instability by resonantly coupling with and transferring energy into pairs of inertial waves (or inertial-gravity waves when buoyancy is important). In this paper, we perform inviscid three-dimensional global hydrodynamic simulations to examine the growth and consequence of this instability operating on the spiral waves driven by a Jupiter-mass planet in a protoplanetary disk. We find that the spiral waves are destabilized via the spiral wave instability (SWI), generating hydrodynamic turbulence and sustained radially-alternating vertical flows that appear to be associated with long wavelength inertial modes. In the interval $0.3~R_p \leq R \leq 0.7~R_p$, where $R_p$ denotes the semi-major axis of the planetary orbit (assumed to be 5~au), the estimated vertical diffusion rate associated with the turbulence is characterized by $\alpha_{\rm diff} \sim (0.2-1.2) \times 10^{-2}$. For the disk model considered here, the diffusion rate is such that particles with sizes up to several centimeters are vertically mixed within the first pressure scale height. This suggests that the instability of spiral waves launched by a giant planet can significantly disperse solid particles and trace chemical species from the midplane. In planet formation models where the continuous local production of chondrules/pebbles occurs over Myr time scales to provide a feedstock for pebble accretion onto these bodies, this stirring of solid particles may add a time constraint: planetary embryos and large asteroids have to form before a gas giant forms in the outer disk, otherwise the SWI will significantly decrease the chondrule/pebble accretion efficiency.

We report on new distances and proper motions to seven stars across the Serpens/Aquila complex. The observations were obtained as part of the Gould's Belt Distances Survey (GOBELINS) project between September 2013 and April 2016 with the Very Long Baseline Array (VLBA). One of our targets is the proto-Herbig AeBe object EC 95, which is a binary system embedded in the Serpens Core. For this system, we combined the GOBELINS observations with previous VLBA data to cover a total period of ~8 years, and derive the orbital elements and an updated source distance. The individual distances to sources in the complex are fully consistent with each other, and the mean value corresponds to a distance of $436.0\pm9.2$~pc for the Serpens/W40 complex. Given this new evidence, we argue that Serpens Main, W40 and Serpens South are physically associated and form a single cloud structure.

We present the results of the Gould's Belt Distances Survey (GOBELINS) of young star forming regions towards the Orion Molecular Cloud Complex. We detected 36 YSOs with the Very Large Baseline Array (VLBA), 27 of which have been observed in at least 3 epochs over the course of 2 years. At least half of these YSOs belong to multiple systems. We obtained parallax and proper motions towards these stars to study the structure and kinematics of the Complex. We measured a distance of 388$\pm$5 pc towards the Orion Nebula Cluster, 428$\pm$10 pc towards the southern portion L1641, 388$\pm$10 pc towards NGC 2068, and roughly $\sim$420 pc towards NGC 2024. Finally, we observed a strong degree of plasma radio scattering towards $\lambda$ Ori.

We present results from a suite of three-dimensional global hydrodynamic simulations which show that spiral density waves propagating in circumstellar disks are unstable to the growth of a parametric instability that leads to break-down of the flow into turbulence. This spiral wave instability (SWI) arises from a resonant interaction between pairs of inertial waves, or inertial-gravity waves, and the background spiral wave. The development of the instability in the linear regime involves the growth of a broad spectrum of inertial modes, with growth rates on the order of the orbital time, and results in a nonlinear saturated state in which turbulent velocity perturbations are of a similar magnitude to those induced by the spiral wave. The turbulence induces angular momentum transport, and vertical mixing, at a rate that depends locally on the amplitude of the spiral wave (we obtain a stress parameter $\alpha \sim 5 \times 10^{-4}$ in our reference model). The instability is found to operate in a wide-range of disk models, including those with isothermal or adiabatic equations of state, and in viscous disks where the dimensionless kinematic viscosity $\nu\le 10^{-5}$. This robustness suggests that the instability will have applications to a broad range of astrophysical disk-related phenomena, including those in close binary systems, planets embedded in protoplanetary disks (including Jupiter in our own Solar System) and FU Orionis outburst models. Further work is required to determine the nature of the instability, and to evaluate its observational consequences, in physically more complete disk models than we have considered in this paper.

We present key results from the Herschel Orion Protostar Survey (HOPS): spectral energy distributions (SEDs) and model fits of 330 young stellar objects, predominantly protostars, in the Orion molecular clouds. This is the largest sample of protostars studied in a single, nearby star-formation complex. With near-infrared photometry from 2MASS, mid- and far-infrared data from Spitzer and Herschel, and sub-millimeter photometry from APEX, our SEDs cover 1.2-870 $\mu$m and sample the peak of the protostellar envelope emission at ~100 $\mu$m. Using mid-IR spectral indices and bolometric temperatures, we classify our sample into 92 Class 0 protostars, 125 Class I protostars, 102 flat-spectrum sources, and 11 Class II pre-main-sequence stars. We implement a simple protostellar model (including a disk in an infalling envelope with outflow cavities) to generate a grid of 30400 model SEDs and use it to determine the best-fit model parameters for each protostar. We argue that far-IR data are essential for accurate constraints on protostellar envelope properties. We find that most protostars, and in particular the flat-spectrum sources, are well-fit. The median envelope density and median inclination angle decrease from Class 0 to Class I to flat-spectrum protostars, despite the broad range in best-fit parameters in each of the three categories. We also discuss degeneracies in our model parameters. Our results confirm that the different protostellar classes generally correspond to an evolutionary sequence with a decreasing envelope infall rate, but the inclination angle also plays a role in the appearance, and thus interpretation, of the SEDs.

We present a near-IR survey for the visual multiples in the Orion molecular clouds region at separations between 100 and 1000 AU. These data were acquired at 1.6~$\mu$m with the NICMOS and WFC3 cameras on the Hubble Space Telescope. Additional photometry was obtained for some of the sources at 2.05~$\mu$m with NICMOS and in the $L'$-band with NSFCAM2 on the IRTF. Towards 129 protostars and 197 pre-main sequence stars with disks observed with WFC3, we detect 21 and 28 candidate companions between the projected separations of 100---1000 AU, of which less than 5 and 8, respectively, are chance line of sight coincidences. The resulting companion fraction ($CF$) after the correction for the line of sight contamination is 14.4$^{+1.1}_{-1.3}$% for protostars and 12.5$^{+1.2}_{-0.8}$% for the pre-main sequence stars. These values are similar to those found for main sequence stars, suggesting that there is little variation in the $CF$ with evolution, although several observational biases may mask a decrease in the $CF$ from protostars to the main sequence stars. After segregating the sample into two populations based on the surrounding surface density of YSOs, we find that the $CF$ in the high stellar density regions ($\Sigma_{YSO} > 45$~pc$^{-2}$) is approximately 50% higher than that found in the low stellar density regions ($\Sigma_{YSO} < 45$~pc$^{-2}$). We interpret this as evidence for the elevated formation of companions at 100 to 1000 AU in the denser environments of Orion. We discuss possible reasons for this elevated formation.

The disk surrounding SAO 206462, an 8 Myr-old Herbig Ae star, has recently been reported to exhibit spiral arms, an asymmetric dust continuum, and a dust-depleted inner cavity. By carrying out two-dimensional, two-fluid hydrodynamic calculations, we find that a planetary-mass companion located at the outer disk could be responsible for these observed structures. In this model, the planet excites primary and secondary arms interior to its orbit. It also carves a gap and generates a local pressure bump at the inner gap edge where a vortex forms through Rossby wave instability. The vortex traps radially drifting dust particles, forming a dust-depleted cavity in the inner disk. We propose that the vortex is responsible for the brightest southwestern peak seen in infrared scattered light and sub-millimeter dust continuum emission. In particular, it is possible that the scattered light is boosted as one of the spiral arms passes through the high density vortex region, although the vortex alone may be able to explain the peak. We suggest that a planetary companion with a mass of 10-15 M_J is orbiting SAO 206462 at 100-120 au. Monitoring of the brightest peak over the next few years will help reveal its origin because the spiral arms and vortex will show distinguishable displacement.

We examine the spectroscopic binary population for two massive nearby regions of clustered star formation, the Orion Nebula Cluster and NGC 2264, supplementing the data presented by Tobin et al. (2009, 2015) with more recent observations and more extensive analysis. The inferred multiplicity fraction up to 10 AU based on these observations is $5.3\pm 1.2$% for NGC 2264 and $5.8\pm 1.1$% for the ONC; they are consistent with the distribution of binaries in the field in the relevant parameter range. Eight of the multiple systems in the sample have enough epochs to make an initial fit for the orbital parameters. Two of these sources are double-lined spectroscopic binaries; for them we determine the mass ratio. Our reanalysis of the distribution of stellar radial velocities towards these clusters presents a significantly better agreement between stellar and gas kinematics than was previously thought.

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.

Sub-virial gravitational collapse is one mechanism by which star clusters may form. Here we investigate whether this mechanism can be inferred from observations of young clusters. To address this question, we have computed SPH simulations of the initial formation and evolution of a dynamically young star cluster through cold (sub-virial) collapse, starting with an ellipsoidal, turbulently seeded distribution of gas, and forming sink particles representing (proto)stars. While the initial density distributions of the clouds do not have large initial mass concentrations, gravitational focusing due to the global morphology leads to cluster formation. We use the resulting structures to extract observable morphological and kinematic signatures for the case of sub-virial collapse. We find that the signatures of the initial conditions can be erased rapidly as the gas and stars collapse, suggesting that kinematic observations need to be made either early in cluster formation and/or at larger scales, away from the growing cluster core. Our results emphasize that a dynamically young system is inherently evolving on short timescales, so that it can be highly misleading to use current-epoch conditions to study aspects such as star formation rates as a function of local density. Our simulations serve as a starting point for further studies of collapse including other factors such as magnetic fields and stellar feedback.

We analyze the spatial distribution of dusty young stellar objects (YSOs) identified in the Spitzer Survey of the Orion Molecular clouds, augmenting these data with Chandra X-ray observations to correct for incompleteness in dense clustered regions. We also devise a scheme to correct for spatially varying incompleteness when X-ray data are not available. The local surface densities of the YSOs range from 1 pc$^{-2}$ to over 10,000 pc$^{-2}$, with protostars tending to be in higher density regions. This range of densities is similar to other surveyed molecular clouds with clusters, but broader than clouds without clusters. By identifying clusters and groups as continuous regions with surface densities $\ge10$ pc$^{-2}$, we find that 59% of the YSOs are in the largest cluster, the Orion Nebular Cluster (ONC), while 13% of the YSOs are found in a distributed population. A lower fraction of protostars in the distributed population is evidence that it is somewhat older than the groups and clusters. An examination of the structural properties of the clusters and groups show that the peak surface densities of the clusters increase approximately linearly with the number of members. Furthermore, all clusters with more than 70 members exhibit asymmetric and/or highly elongated structures. The ONC becomes azimuthally symmetric in the inner 0.1 pc, suggesting that the cluster is only $\sim 2$ Myr in age. We find the star formation efficiency (SFE) of the Orion B cloud is unusually low, and that the SFEs of individual groups and clusters are an order of magnitude higher than those of the clouds. Finally, we discuss the relationship between the young low mass stars in the Orion clouds and the Orion OB 1 association, and we determine upper limits to the fraction of disks that may be affected by UV radiation from OB stars or by dynamical interactions in dense, clustered regions.

We report on a series of numerical simulations of gas clouds with self-gravity forming sink particles, adopting an isothermal equation of state to isolate the effects of gravity from thermal physics on the resulting sink mass distributions. Simulations starting with supersonic velocity fluctuations develop sink mass functions with a high-mass power-law tail $dN/d\log M \propto M^{\Gamma}$, $\Gamma = -1 \pm 0.1$, independent of the initial Mach number of the velocity field. Similar results but with weaker statistical significance hold for a simulation starting with initial density fluctuations. This mass function power-law dependence agrees with the asymptotic limit found by Zinnecker assuming Bondi-Hoyle-Littleton (BHL) accretion, even though the mass accretion rates of individual sinks show significant departures from the predicted $\mdot \propto M^2$ behavior. While BHL accretion is not strictly applicable due to the complexity of the environment, we argue that the final mass functions are the result of a {\em relative} $M^2$ dependence resulting from gravitationally-focused accretion. Our simulations may show the power-law mass function particularly clearly compared with others because our adoption of an isothermal equation of state limits the effects of thermal physics in producing a broad initial fragmentation spectrum; $\Gamma \rightarrow -1$ is an asymptotic limit found only when sink masses grow well beyond their initial values. While we have purposely eliminated many additional physical processes (radiative transfer, feedback) which can affect the stellar mass function, our results emphasize the importance of gravitational focusing for massive star formation.

We present a study of 15 new brown dwarfs belonging to the $\sim7$ Myr old 25 Orionis group and Orion OB1a sub-association with spectral types between M6 and M9 and estimated masses between $\sim0.07$M$_\odot$ and $\sim0.01$ M$_\odot$. By comparing them through a Bayesian method with low mass stars ($0.8\lesssim$ M/M$_\odot\lesssim0.1$) from previous works in the 25 Orionis group, we found statistically significant differences in the number fraction of classical T Tauri stars, weak T Tauri stars, class II, evolved discs and purely photospheric emitters at both sides of the sub-stellar mass limit. Particularly we found a fraction of $3.9^{+2.4}_{-1.6}~\%$ low mass stars classified as CTTS and class II or evolved discs, against a fraction of $33.3^{+10.8}_{-9.8}~\%$ in the sub-stellar mass domain. Our results support the suggested scenario in which the dissipation of discs is less efficient for decreasing mass of the central object.

We present deep ($\sim 17~\mu$Jy) radio continuum observations of the Serpens molecular cloud, the Serpens south cluster, and the W40 region obtained using the Very Large Array in its A configuration. We detect a total of 146 sources, 29 of which are young stellar objects (YSOs), 2 are BV stars and 5 more are associated with phenomena related to YSOs. Based on their radio variability and spectral index, we propose that about 16 of the remaining 110 unclassified sources are also YSOs. For approximately 65% of the known YSOs detected here as radio sources, the emission is most likely non-thermal, and related to stellar coronal activity. As also recently observed in Ophiuchus, our sample of YSOs with X-ray counterparts lies below the fiducial G\"udel & Benz relation. Finally, we analyze the proper motions of 9 sources in the W40 region. This allows us to better constrain the membership of the radio sources in the region.

We carry out two-fluid, two-dimensional global hydrodynamic simulations to test whether protostellar infall can trigger Rossby wave instability (RWI) in protoplanetry disks. Our results show that infall can trigger the RWI and generate vortices near the outer edge of the mass landing on the disk (i.e. centrifugal radius). We find that the RWI is triggered under a variety of conditions, although the details depend on the disk parameters and the infall pattern. The common key feature of triggering the RWI is the steep radial gradient of the azimuthal velocity induced by the local increase in density at the outer edge of the infall region. Vortices form when the instability enters the nonlinear regime. In our standard model where self-gravity is neglected, vortices merge together to a single vortex within $\sim 20$ local orbital times, and the merged vortex survives for the remaining duration of the calculation ($> 170$ local orbital times). The vortex takes part in outward angular momentum transport, with a Reynolds stress of $\lesssim10^{-2}$. Our two-fluid calculations show that vortices efficiently trap dust particles with stopping times of the order of the orbital time, locally enhancing the dust to gas ratio for particles of the appropriate size by a factor of $\sim 40$ in our standard model. When self-gravity is considered, however, vortices tend to be impeded from merging and may eventually dissipate. We conclude it may well have that protoplanetary disks have favorable conditions for vortex formation during the protostellar infall phase, which might enhance early planetary core formation.

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.

We present an expanded kinematic study of the young cluster NGC 2264 based upon optical radial velocities measured using multi-fiber echelle spectroscopy at the 6.5 meter MMT and Magellan telescopes. We report radial velocities for 695 stars, of which approximately 407 stars are confirmed or very likely members. Our results more than double the number of members with radial velocities from F{\H u}r{\'e}sz et al., resulting in a much better defined kinematic relationship between the stellar population and the associated molecular gas. In particular, we find that there is a significant subset of stars that are systematically blueshifted with respect to the molecular ($^{13}$CO) gas. The detection of Lithium absorption and/or infrared excesses in this blue-shifted population suggests that at least some of these stars are cluster members; we suggest some speculative scenarios to explain their kinematics. Our results also more clearly define the redshifted population of stars in the northern end of the cluster; we suggest that the stellar and gas kinematics of this region are the result of a bubble driven by the wind from O7 star S Mon. Our results emphasize the complexity of the spatial and kinematic structure of NGC 2264, important for eventually building up a comprehensive picture of cluster formation.

We present a multi-epoch radio study of the Taurus-Auriga star-forming complex made with the Karl G. Jansky Very Large Array at frequencies of 4.5 GHz and 7.5 GHz. We detect a total of 610 sources, 59 of which are related to young stellar objects and 18 to field stars. The properties of 56\% of the young stars are compatible with non-thermal radio emission. We also show that the radio emission of more evolved young stellar objects tends to be more non-thermal in origin and, in general, that their radio properties are compatible with those found in other star forming regions. By comparing our results with previously reported X-ray observations, we notice that young stellar objects in Taurus-Auriga follow a G\"{u}del-Benz relation with $\kappa$=0.03, as we previously suggested for other regions of star formation. In general, young stellar objects in Taurus-Auriga and in all the previous studied regions seem to follow this relation with a dispersion of $\sim1$ dex. Finally, we propose that most of the remaining sources are related with extragalactic objects but provide a list of 46 unidentified radio sources whose radio properties are compatible with a YSO nature.

We present CARMA 2.9 mm dust continuum emission observations of a sample of 14 Herschel-detected Class 0 protostars in the Orion A and B molecular clouds, drawn from the PACS Bright Red Sources (PBRS) sample (Stutz et al.). These objects are characterized by very red 24 \micron\ to 70 \micron\ colors and prominent submillimeter emission, suggesting that they are very young Class 0 protostars embedded in dense envelopes. We detect all of the PBRS in 2.9 mm continuum emission and emission from 4 protostars and 1 starless core in the fields toward the PBRS; we also report 1 new PBRS source. The ratio of 2.9 mm luminosity to bolometric luminosity is higher by a factor of $\sim$5 on average, compared to other well-studied protostars in the Perseus and Ophiuchus clouds. The 2.9 mm visibility amplitudes for 6 of the 14 PBRS are very flat as a function of uv-distance, with more than 50\% of the source emission arising from radii $<$ 1500 AU. These flat visibility amplitudes are most consistent with spherically symmetric envelope density profiles with $\rho$~$\propto$~R$^{-2.5}$. Alternatively, there could be a massive unresolved structure like a disk or a high-density inner envelope departing from a smooth power-law. The large amount of mass on scales $<$ 1500 AU (implying high average central densities) leads us to suggest that that the PBRS with flat visibility amplitude profiles are the youngest PBRS and may be undergoing a brief phase of high mass infall/accretion and are possibly among the youngest Class 0 protostars. The PBRS with more rapidly declining visibility amplitudes still have large envelope masses, but could be slightly more evolved.

We improve on our previous treatments of long-term evolution of protostellar disks by explicitly solving disk self-gravity in two dimensions. The current model is an extension of the one-dimensional layered accretion disk model of Bae et al. We find that gravitational instability (GI)-induced spiral density waves heat disks via compressional heating (i.e. $P\rm{d}V$ work), and can trigger accretion outbursts by activating the magnetorotational instability (MRI) in the magnetically inert disk dead-zone. The GI-induced spiral waves propagate well inside of gravitationally unstable region before they trigger outbursts at $R \lesssim 1$ AU where GI cannot be sustained. This long-range propagation of waves cannot be reproduced with the previously used local $\alpha$ treatments for GI. In our standard model where zero dead-zone residual viscosity ($\alpha_{\rm rd}$) is assumed, the GI-induced stress measured at the onset of outbursts is locally as large as $0.01$ in terms of the generic $\alpha$ parameter. However, as suggested in our previous one-dimensional calculations, we confirm that the presence of a small but finite $\alpha_{\rm rd}$ triggers thermally-driven bursts of accretion instead of the GI + MRI-driven outbursts that are observed when $\alpha_{\rm rd}=0$. The inclusion of non-zero residual viscosity in the dead-zone decreases the importance of GI soon after mass feeding from the envelope cloud ceases. During the infall phase while the central protostar is still embedded, our models stay in a quiescent accretion phase with $\dot{M}_{acc}\sim10^{-8}-10^{-7}~M_{\odot}~{\rm yr^{-1}}$ over $60~\%$ of the time and spend less than $15~\%$ of the infall phase in accretion outbursts.

We present a spectroscopic survey of the stellar population of the Sigma Orionis cluster. We have obtained spectral types for 340 stars. Spectroscopic data for spectral typing come from several spectrographs with similar spectroscopic coverage and resolution. More than a half of stars of our sample are members confirmed by the presence of lithium in absorption, strong H$\alpha$ in emission or weak gravity-sensitive features. In addition, we have obtained high resolution (R~34000) spectra in the H$\alpha$ region for 169 stars in the region. Radial velocities were calculated from this data set. The radial velocity distribution for members of the cluster is in agreement with previous work. Analysis of the profile of the H$\alpha$ line and infrared observations reveals two binary systems or fast rotators that mimic the H$\alpha$ width expected in stars with accretion disks. On the other hand there are stars with optically thick disks and narrow H$\alpha$ profile not expected in stars with accretion disks. This contribution constitutes the largest homogeneous spectroscopic data set of the Sigma Orionis cluster to date.

We present the results of a survey of the low mass star and brown dwarf population of the 25 Orionis group. Using optical photometry from the CIDA Deep Survey of Orion, near IR photometry from the Visible and Infrared Survey Telescope for Astronomy and low resolution spectroscopy obtained with Hectospec at the MMT, we selected 1246 photometric candidates to low mass stars and brown dwarfs with estimated masses within $0.02 \lesssim M/M_\odot \lesssim 0.8$ and spectroscopically confirmed a sample of 77 low mass stars as new members of the cluster with a mean age of $\sim$7 Myr. We have obtained a system initial mass function of the group that can be well described by either a Kroupa power-law function with indices $\alpha_3=-1.73\pm0.31$ and $\alpha_2=0.68\pm0.41$ in the mass ranges $0.03\leq M/M_\odot\leq 0.08$ and $0.08\leq M/M_\odot\leq0.5$ respectively, or a Scalo log-normal function with coefficients $m_c=0.21^{+0.02}_{-0.02}$ and $\sigma=0.36\pm0.03$ in the mass range $0.03\leq M/M_\odot\leq0.8$. From the analysis of the spatial distribution of this numerous candidate sample, we have confirmed the East-West elongation of the 25 Orionis group observed in previous works, and rule out a possible southern extension of the group. We find that the spatial distributions of low mass stars and brown dwarfs in 25 Orionis are statistically indistinguishable. Finally, we found that the fraction of brown dwarfs showing IR excesses is higher than for low mass stars, supporting the scenario in which the evolution of circumstellar discs around the least massive objects could be more prolonged.

A comparison of accretion and (turbulent) magnetic diffusion timescales for sheets and filaments demonstrates that dense star-forming clouds generally will -- under realistic conditions -- become supercritical due to mass accretion on timescales at least an order of magnitude shorter than ambipolar and/or turbulent diffusion timescales. Thus, ambipolar or turbulent diffusion -- while present -- is unlikely to control the formation of cores and stars.

We present results from a high-sensitivity (60 $\mu$Jy), large-scale (2.26 square degree) survey obtained with the Karl G. Jansky Very Large Array as part of the Gould's Belt Survey program. We detected 374 and 354 sources at 4.5 and 7.5 GHz, respectively. Of these, 148 are associated with previously known Young Stellar Objects (YSOs). Another 86 sources previously unclassified at either optical or infrared wavelengths exhibit radio properties that are consistent with those of young stars. The overall properties of our sources at radio wavelengths such as their variability and radio to X-ray luminosity relation are consistent with previous results from the Gould's Belt Survey. Our detections provide target lists for followup VLBA radio observations to determine their distances as YSOs are located in regions of high nebulosity and extinction, making it difficult to measure optical parallaxes.

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.

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.

Edge-on protostars are valuable for understanding the disk and envelope properties of embedded young stellar objects, since the disk, envelope, and envelope cavities are all distinctly visible in resolved images and well constrained in modeling. Comparing 2MASS, WISE, Spitzer, Herschel, APEX, and IRAM photometry and limits from 1.25 to 1200 microns, Spitzer spectroscopy from 5 to 40 microns, and high-resolution Hubble imaging at 1.60 and 2.05 microns to radiative transfer modeling, we determine envelope and disk properties for the Class I protostar HOPS 136, an edge-on source in Orion's Lynds 1641 region. The source has a bolometric luminosity of 0.8 L_sun, a bolometric temperature of 170 K, and a ratio of submillimeter to bolometric luminosity of 0.8%. Via modeling, we find a total luminosity of 4.7 L_sun (larger than the observed luminosity due to extinction by the disk), an envelope mass of 0.06 M_sun, and a disk radius and mass of 450 AU and 0.002 M_sun. The stellar mass is highly uncertain but is estimated to fall between 0.4 and 0.5 M_sun. To reproduce the flux and wavelength of the near-infrared scattered-light peak in the spectral energy distribution, we require 5.4 x 10^-5 M_sun of gas and dust in each cavity. The disk has a large radius and a mass typical of more evolved T Tauri stars in spite of the significant remaining envelope. HOPS 136 appears to be a key link between the protostellar and optically revealed stages of star formation.

We present observations of three Class 0/I protostars (L1157-mm, CB230 IRS1, and L1165-SMM1) using the Karl G. Jansky Very Large Array (VLA) and observations of two (L1165-SMM1 and CB230 IRS1) with the Combined Array for Research in Millimeter-wave Astronomy (CARMA). The VLA observations were taken at wavelengths of $\lambda = 7.3$ mm, 1.4 cm, 3.3 cm, 4.0 cm, and 6.5 cm with a best resolution of $\sim$0\farcs06 (18 AU) at 7.3 mm. The L1165-SMM1 CARMA observations were taken at $\lambda = 1.3$ mm with a best resolution of $\sim0\farcs3$ (100 AU), and the CB230 IRS1 observations were taken at $\lambda = 3.4$ mm with a best resolution of $\sim$3\arcsec\ (900 AU). We find that L1165-SMM1 and CB230 IRS1 have probable binary companions at separations of $\sim$0\farcs3 (100 AU) from detections of secondary peaks at multiple wavelengths. The position angles of these companions are nearly orthogonal to the direction of the observed bipolar outflows, consistent with the expected protostellar disk orientations. We suggest that these companions may have formed from disk fragmentation; turbulent fragmentation would not preferentially arrange the binary companions to be orthogonal to the outflow direction. For L1165-SMM1, both the 7.3 mm and 1.3 mm emission show evidence of a large (R $>$ 100 AU) disk. For the L1165-SMM1 primary protostar and the CB230 IRS1 secondary protostar, the 7.3 mm emission is resolved into structures consistent with $\sim20$ AU radius disks. For the other protostars, including L1157-mm, the emission is unresolved, suggesting disks with radii $< 20$ AU.

The dust sublimation walls of disks around T Tauri stars represent a directly observable cross-section through the disk atmosphere and midplane. Their emission properties can probe the grain size distribution and composition of the innermost regions of the disk, where terrestrial planets form. Here we calculate the inner dust sublimation wall properties for four classical T Tauri stars with a narrow range of spectral types and inclination angles and a wide range of mass accretion rates to determine the extent to which the walls are radially curved. Best-fits to the near- and mid-IR excesses are found for curved, 2-layer walls in which the lower layer contains larger, hotter, amorphous pyroxene grains with Mg/(Mg+Fe)=0.6 and the upper layer contains submicron, cooler, mixed amorphous olivine and forsterite grains. As the mass accretion rates decrease from 10^(-8) to 10^(-10) Msol/yr, the maximum grain size in the lower layer decreases from 3 to 0.5 microns. We attribute this to a decrease in fragmentation and turbulent support for micron-sized grains with decreasing viscous heating. The atmosphere of these disks is depleted of dust with dust-gas mass ratios 1x10^(-4) of the ISM value, while the midplane is enhanced to 8 times the ISM value. For all accretion rates, the wall contributes at least half of the flux in the optically thin 10 micron silicate feature. Finally, we find evidence for an iron gradient in the disk, suggestive of that found in our solar system.

We present large-scale ($\sim$ 2000 square arcminutes), deep ($\sim$ 20 $\mu$Jy), high-resolution ($\sim$ 1$''$) radio observations of the Ophiuchus star-forming complex obtained with the Karl G.\ Jansky Very Large Array at $\lambda$ = 4 and 6 cm. In total, 189 sources were detected, 56 of them associated with known young stellar sources, and 4 with known extragalactic objects; the other 129 remain unclassified, but most of them are most probably background quasars. The vast majority of the young stars detected at radio wavelengths have spectral types K or M, although we also detect 4 objects of A/F/B types and 2 brown dwarf candidates. At least half of these young stars are non-thermal (gyrosynchrotron) sources, with active coronas characterized by high levels of variability, negative spectral indices, and (in some cases) significant circular polarization. As expected, there is a clear tendency for the fraction of non-thermal sources to increase from the younger (Class 0/I or flat spectrum) to the more evolved (Class III or weak line T Tauri) stars. The young stars detected both in X-rays and at radio wavelengths broadly follow a G\"udel-Benz relation, but with a different normalization than the most radio-active types of stars. Finally, we detect a $\sim$ 70 mJy compact extragalactic source near the center of the Ophiuchus core, which should be used as gain calibrator for any future radio observations of this region.

We perform calculations of our one-dimensional, two-zone disk model to study the long-term evolution of the circumstellar disk. In particular, we adopt published photoevaporation prescriptions and examine whether the photoevaporative loss alone, coupled with a range of initial angular momenta of the protostellar cloud, can explain the observed decline of the frequency of optically-thick dusty disks with increasing age. In the parameter space we explore, disks have accreting and/or non-accreting transitional phases lasting of $\lesssim20 %$ of their lifetime, which is in reasonable agreement with observed statistics. Assuming that photoevaporation controls disk clearing, we find that initial angular momentum distribution of clouds needs to be weighted in favor of slowly rotating protostellar cloud cores. Again, assuming inner disk dispersal by photoevaporation, we conjecture that this skewed angular momentum distribution is a result of fragmentation into binary or multiple stellar systems in rapidly-rotating cores. Accreting and non-accreting transitional disks show different evolutionary paths on the $\dot{M}-R_{\rm wall}$ plane, which possibly explains the different observed properties between the two populations. However, we further find that scaling the photoevaporation rates downward by a factor of 10 makes it difficult to clear the disks on the observed timescales, showing that the precise value of the photoevaporative loss is crucial to setting the clearing times. While our results apply only to pure photoevaporative loss (plus disk accretion), there may be implications for models in which planets clear disks preferentially at radii of order 10 AU.

Studies by Lada (2010) and Heiderman (2010) have suggested that star formation mostly occurs above a threshold in gas surface density Sigma of Sigma_c = 120 Msun pc^{-2} (A_K = 0.8). Heiderman infer a threshold by combining low-mass star-forming regions, which show a steep increase in the star formation rate per unit area Sigma_SFR with increasing Sigma, and massive cores forming luminous stars which show a linear relation. We argue that these observations do not require a particular density threshold. The steep dependence of Sigma_SFR, approaching unity at protostellar core densities, is a natural result of the increasing importance of self-gravity at high densities along with the corresponding decrease in evolutionary timescales. The linear behavior of Sigma_SFR vs. Sigma in massive cores is consistent with probing dense gas in gravitational collapse, forming stars at a characteristic free-fall timescale given by the use of a particular molecular tracer. The low-mass and high-mass regions show different correlations between gas surface density and the area A spanned at that density, with A=Sigma^{-3} for low-mass regions and A=Sigma^{-1} for the massive cores; this difference, along with the use of differing techniques to measure gas surface density and star formation, suggests that connecting the low-mass regions with massive cores is problematic. We show that the approximately linear relationship between dense gas mass and stellar mass used by Lada similarly does not demand a particular threshold for star formation, and requires continuing formation of dense gas. Our results are consistent with molecular clouds forming by galactic hydrodynamic flows with subsequent gravitational collapse.

We present high-resolution sub/millimeter interferometric imaging of the Class 0 protostar L1527 IRS (IRAS 04368+2557) at 870 micron and 3.4 mm from the Submillimeter Array (SMA) and Combined Array for Research in Millimeter Astronomy (CARMA). We detect the signature of an edge-on disk surrounding the protostar with an observed diameter of 180 AU in the sub/millimeter images. The mass of the disk is estimated to be 0.007 M_sun, assuming optically thin, isothermal dust emission. The millimeter spectral index is observed to be quite shallow at all the spatial scales probed; alpha ~ 2, implying a dust opacity spectral index beta ~ 0. We model the emission from the disk and surrounding envelope using Monte Carlo radiative transfer codes, simultaneously fitting the sub/millimeter visibility amplitudes, sub/millimeter images, resolved L\arcmin\ image, spectral energy distribution, and mid-infrared spectrum. The best fitting model has a disk radius of R = 125 AU, is highly flared (H ~ R^1.3), has a radial density profile rho ~ R^-2.5, and has a mass of 0.0075 M_sun. The scale height at 100 AU is 48 AU, about a factor of two greater than vertical hydrostatic equilibrium. The resolved millimeter observations indicate that disks may grow rapidly throughout the Class 0 phase. The mass and radius of the young disk around L1527 is comparable to disks around pre-main sequence stars; however, the disk is considerably more vertically extended, possibly due to a combination of lower protostellar mass, infall onto the disk upper layers, and little settling of ~1 micron-sized dust grains.

Using NASA IRTF SpeX data from 0.8 to 4.5 $\mu$m, we determine self-consistently the stellar properties and excess emission above the photosphere for a sample of classical T Tauri stars (CTTS) in the Taurus molecular cloud with varying degrees of accretion. This process uses a combination of techniques from the recent literature as well as observations of weak-line T Tauri stars (WTTS) to account for the differences in surface gravity and chromospheric activity between the TTS and dwarfs, which are typically used as photospheric templates for CTTS. Our improved veiling and extinction estimates for our targets allow us to extract flux-calibrated spectra of the excess in the near-infrared. We find that we are able to produce an acceptable parametric fit to the near-infrared excesses using a combination of up to three blackbodies. In half of our sample, two blackbodies at temperatures of 8000 K and 1600 K suffice. These temperatures and the corresponding solid angles are consistent with emission from the accretion shock on the stellar surface and the inner dust sublimation rim of the disk, respectively. In contrast, the other half requires three blackbodies at 8000, 1800, and 800 K, to describe the excess. We interpret the combined two cooler blackbodies as the dust sublimation wall with either a contribution from the disk surface beyond the wall or curvature of the wall itself, neither of which should have single-temperature blackbody emission. In these fits, we find no evidence of a contribution from optically thick gas inside the inner dust rim.

We perform a census of the reddest, and potentially youngest, protostars in the Orion molecular clouds using data obtained with the PACS instrument onboard the Herschel Space Observatory and the LABOCA and SABOCA instruments on APEX as part of the Herschel Orion Protostar Survey (HOPS). A total of 55 new protostar candidates are detected at 70 um and 160 um that are either too faint (m24 > 7 mag) to be reliably classified as protostars or undetected in the Spitzer/MIPS 24 um band. We find that the 11 reddest protostar candidates with log (lambda F_lambda 70) / (lambda F_lambda 24) > 1.65 are free of contamination and can thus be reliably explained as protostars. The remaining 44 sources have less extreme 70/24 colors, fainter 70 um fluxes, and higher levels of contamination. Taking the previously known sample of Spitzer protostars and the new sample together, we find 18 sources that have log (lambda F_lambda 70) / (lambda F_lambda 24) > 1.65; we name these sources "PACS Bright Red sources", or PBRs. Our analysis reveals that the PBRs sample is composed of Class 0 like sources characterized by very red SEDs (T_bol < 45 K) and large values of sub-millimeter fluxes (L_smm/L_bol > 0.6%). Modified black-body fits to the SEDs provide lower limits to the envelope masses of 0.2 M_sun to 2 M_sun and luminosities of 0.7 L_sun to 10 L_sun. Based on these properties, and a comparison of the SEDs with radiative transfer models of protostars, we conclude that the PBRs are most likely extreme Class 0 objects distinguished by higher than typical envelope densities and hence, high mass infall rates.

We present a study on the spatial distribution of N2D+ and N2H+ in thirteen protostellar systems. Eight of thirteen objects observed with the IRAM 30m telescope show relative offsets between the peak N2D+ (J=2-1) and N2H+ (J=1-0) emission. We highlight the case of L1157 using interferometric observations from the Submillimeter Array and Plateau de Bure Interferometer of the N2D+ (J=3-2) and N2H+ (J=1-0) transitions respectively. Depletion of N2D+ in L1157 is clearly observed inside a radius of ~2000 AU (7") and the N2H+ emission is resolved into two peaks at radii of ~1000 AU (3.5"), inside the depletion region of N2D+. Chemical models predict a depletion zone in N2H+ and N2D+ due to destruction of H2D+ at T ~ 20 K and the evaporation of CO off dust grains at the same temperature. However, the abundance offsets of 1000 AU between the two species are not reproduced by chemical models, including a model that follows the infall of the protostellar envelope. The average abundance ratios of N2D+ to N2H+ have been shown to decrease as protostars evolve by Emprechtinger et al., but this is the first time depletion zones of N2D+ have been spatially resolved. We suggest that the difference in depletion zone radii for N2H+ and N2D+ is caused by either the CO evaporation temperature being above 20 K or an H2 ortho-to-para ratio gradient in the inner envelope.

We extend the one-dimensional, two-zone models of long-term protostellar disk evolution with infall of Zhu et al. to consider the potential effects of a finite viscosity in regions where the ionization is too low for the magnetorotational instability (MRI) to operate (the "dead zone"). We find that the presence of a small but finite dead zone viscosity, as suggested by simulations of stratified disks with MRI-active outer layers, can trigger inside-out bursts of accretion, starting at or near the inner edge of the disk, instead of the previously-found outside-in bursts with zero dead zone viscosity, which originate at a few AU in radius. These inside-out bursts of accretion bear a qualitative resemblance to the outburst behavior of one FU Ori object, V1515 Cyg, in contrast to the outside-in burst models which more closely resemble the accretion events in FU Ori and V1057 Cyg. Our results suggest that the type and frequency of outbursts are potentially a probe of transport efficiency in the dead zone. Simulations must treat the inner disk regions, $R \lesssim 0.5$ AU, to show the detailed time evolution of accretion outbursts in general and to observe the inside-out bursts in particular.

We extend our previous study of the stellar population of L1641, the lower-density star-forming region of the Orion A cloud south of the dense Orion Nebula Cluster (ONC), with the goal of testing whether there is a statistically significant deficiency of high-mass stars in low-density regions. Previously, we compared the observed ratio of low-mass stars to high-mass stars with theoretical models of the stellar initial mass function (IMF) to infer a deficiency of the highest-mass stars in L1641. We expand our population study to identify the intermediate mass (late B to G) L1641 members in an attempt to make a more direct comparison with the mass function of the nearby ONC. The spectral type distribution and the K-band luminosity function of L1641 are similar to those of the ONC (Hillenbrand 1997; Muench et al. 2002), but problems of incompleteness and contamination prevent us from making a detailed test for differences. We limit our analysis to statistical tests of the ratio of high-mass to low-mass stars, which indicate a probability of only 3% that the ONC and the southern region of L1641 were drawn from the same population, supporting the hypothesis that the upper mass end of the IMF is dependent on environmental density.

In their earliest stages, protostars accrete mass from their surrounding envelopes through circumstellar disks. Until now, the smallest observed protostar/envelope mass ratio was ~2.1. The protostar L1527 IRS is thought to be in the earliest stages of star formation. Its envelope contains ~1 solar mass of material within a ~0.05 pc radius, and earlier observations suggested the presence of an edge-on disk. Here we report observations of dust continuum emission and 13CO (J=2-1) line emission from the disk around L1527, from which we determine a protostellar mass of M = 0.19 +/- 0.04 solar masses and a protostar/envelope mass ratio of ~0.2. We conclude that most of the luminosity is generated through the accretion process, with an accretion rate of ~6.6 x 10^-7 solar masses per year. If it has been accreting at that rate through much of its life, its age is ~300,000 yr, though theory suggests larger accretion rates earlier, so it may be younger. The presence of a rotationally--supported disk is confirmed and significantly more mass may be added to its planet-forming region as well as the protostar itself.

We present ALMA and VLA observations of the molecular and ionized gas at 0.1-0.3 arcsec resolution in the Class 0 protostellar system IRAS 16293-2422. These data clarify the origins of the protostellar outflows from the deeply embedded sources in this complex region. Source A2 is confirmed to be at the origin of the well known large scale north-east--south-west flow. The most recent VLA observations reveal a new ejection from that protostar, demonstrating that it drives an episodic jet. The central compact part of the other known large scale flow in the system, oriented roughly east-west, is well delineated by the CO(6-5) emission imaged with ALMA and is confirmed to be driven from within component A. Finally, a one-sided blueshifted bubble-like outflow structure is detected here for the first time from source B to the north-west of the system. Its very short dynamical timescale (~ 200 yr), low velocity, and moderate collimation support the idea that source B is the youngest object in the system, and possibly one of the youngest protostars known.

High-energy radiation from T Tauri stars (TTS) influences the amount and longevity of gas in disks, thereby playing a crucial role in the creation of gas giant planets. Here we probe the high-energy ionizing radiation from TTS using high-resolution mid-infrared (MIR) Spitzer IRS Neon forbidden line detections in a sample of disks from IC 348, NGC 2068, and Chamaeleon. We report three new detections of [Ne III] from CS Cha, SZ Cha, and T 54, doubling the known number of [Ne III] detections from TTS. Using [Ne III]-to-[Ne II] ratios in conjunction with X-ray emission measurements, we probe high-energy radiation from TTS. The majority of previously inferred [Ne III]/[Ne II] ratios based on [Ne III] line upper limits are significantly less than 1, pointing to the dominance of either X-ray radiation or soft Extreme-Ultraviolet (EUV) radiation in producing these lines. Here we report the first observational evidence for hard EUV dominated Ne forbidden line production in a T Tauri disk: [Ne III]/[Ne II]~1 in SZ Cha. Our results provide a unique insight into the EUV emission from TTS, by suggesting that EUV radiation may dominate the creation of Ne forbidden lines, albeit in a minority of cases.

We present a survey of the Orion A and B molecular clouds undertaken with the IRAC and MIPS instruments onboard Spitzer. In total, five distinct fields were mapped covering 9 sq. degrees in five mid-IR bands spanning 3-24 microns. The survey includes the Orion Nebula Cluster, the Lynds 1641, 1630 and 1622 dark clouds, and the NGC 2023, 2024, 2068 and 2071 nebulae. These data are merged with the 2MASS point source catalog to generate a catalog of eight band photometry. We identify 3479 dusty young stellar objects (YSOs) in the Orion molecular clouds by searching for point sources with mid-IR colors indicative of reprocessed light from dusty disks or infalling envelopes. The YSOs are subsequently classified on the basis of their mid-IR colors and their spatial distributions are presented. We classify 2991 of the YSOs as pre-main sequence stars with disks and 488 as likely protostars. Most of the sources were observed with IRAC in 2-3 epochs over 6 months; we search for variability between the epochs by looking for correlated variability in the 3.6 and 4.5 micron bands. We find that 50% of the dusty YSOs show variability. The variations are typically small (0.2 mag.) with the protostars showing a higher incidence of variability and larger variations. The observed correlations between the 3.6, 4.5, 5.8 and 8 micron variability suggests that we are observing variations in the heating of the inner disk due to changes in the accretion luminosity or rotating accretion hot spots.

By carrying out two-dimensional two-fluid global simulations, we have studied the response of dust to gap formation by a single planet in the gaseous component of a protoplanetary disk - the so-called "dust filtration" mechanism. We have found that a gap opened by a giant planet at 20 AU in a \alpha=0.01, \dot{M}=10^{-8} Msun/yr disk can effectively stop dust particles larger than 0.1 mm drifting inwards, leaving a sub-millimeter dust cavity/hole. However, smaller particles are difficult to filter by a planet-induced gap due to 1) dust diffusion, and 2) a high gas accretion velocity at the gap edge. An analytic model is also derived to understand what size particles can be filtered by the gap edge. Finally, with our updated understanding of dust filtration, we have computed Monte-Carlo radiative transfer models with variable dust size distributions to generate the spectral energy distributions (SEDs) of disks with gaps. By comparing with transitional disk observations (e.g. GM Aur), we have found that dust filtration alone has difficulties to deplete small particles sufficiently to explain the near-IR deficit of transitional disks, except under some extreme circumstances. The scenario of gap opening by multiple planets studied previously suffers the same difficulty. One possible solution is by invoking both dust filtration and dust growth in the inner disk. In this scenario, a planet induced gap filters large dust particles in the disk, and the remaining small dust particles passing to the inner disk can grow efficiently without replenishment from fragmentation of large grains. Predictions for ALMA have also been made based on all these scenarios. We conclude that dust filtration with planet(s) in the disk is a promising mechanism to explain submm observations of transitional disks but it may need to be combined with other processes (e.g. dust growth) to explain the near-IR deficit.

V921 Scorpii is a close binary system (separation 0.025") showing the B[e]-phenomenon. The system is surrounded by an enigmatic bipolar nebula, which might have been shaped by episodic mass-loss events, possibly triggered by dynamical interactions between the companion and the circumprimary disk (Kraus et al. 2012a). In this paper, we investigate the spatial structure and kinematics of the circumprimary disk, with the aim to obtain new insights into the still strongly debated evolutionary stage. For this purpose, we combine, for the first time, infrared spectro-interferometry (VLTI/AMBER, R=12,000) and spectro-astrometry (VLT/CRIRES, R=100,000), which allows us to study the AU-scale distribution of circumstellar gas and dust with an unprecedented velocity resolution of 3 km*s^-1. Using a model-independent photocenter analysis technique, we find that the Br-gamma-line emission rotates in the same plane as the dust disk. We can reproduce the wavelength-differential visibilities and phases and the double-peaked line profile using a Keplerian-rotating disk model. The derived mass of the central star is 5.4+/-0.4 M_sun*(d/1150 pc), which is considerably lower than expected from the spectral classification, suggesting that V921 Sco might be more distant (d approx 2kpc) than commonly assumed. Using the geometric information provided by our Br-gamma spectro-interferometric data and Paschen, Brackett, and Pfund line decrement measurements in 61 hydrogen recombination line transitions, we derive the density of the line-emitting gas (N_e=2...6*10^19 m^-3). Given that our measurements can be reproduced with a Keplerian velocity field without outflowing velocity component and the non-detection of age-indicating spectroscopic diagnostics, our study provides new evidence for the pre-main-sequence nature of V921 Sco.

We present results from an optical photometric and spectroscopic survey of the young stellar population in L1641, the low-density star-forming region of the Orion A cloud south of the Orion Nebula Cluster (ONC). Our goal is to determine whether L1641 has a large enough low-mass population to make the known lack of high-mass stars a statistically-significant demonstration of environmental dependence of the upper mass stellar initial mass function (IMF). Our spectroscopic sample consists of IR-excess objects selected from the Spitzer/IRAC survey and non-excess objects selected from optical photometry. We have spectral confirmation of 864 members, with another 98 probable members; of the confirmed members, 406 have infrared excesses and 458 do not. Assuming the same ratio of stars with and without IR excesses in the highly-extincted regions, L1641 may contain as many as ~1600 stars down to ~0.1 solar mass, comparable within a factor of two to the the ONC. Compared to the standard models of the IMF, L1641 is deficient in O and early B stars to a 3-4 sigma significance level, assuming that we know of all the massive stars in L1641. With a forthcoming survey of the intermediate-mass stars, we will be in a better position to make a direct comparison with the neighboring, dense ONC, which should yield a stronger test of the dependence of the high-mass end of the stellar initial mass function upon environment.

Transitional circumstellar disks around young stellar objects have a distinctive infrared deficit around 10 microns in their Spectral Energy Distributions (SED), recently measured by the Spitzer Infrared Spectrograph (IRS), suggesting dust depletion in the inner regions. These disks have been confirmed to have giant central cavities by imaging of the submillimeter (sub-mm) continuum emission using the Submillimeter Array (SMA). However, the polarized near-infrared scattered light images for most objects in a systematic IRS/SMA cross sample, obtained by HiCIAO on the Subaru telescope, show no evidence for the cavity, in clear contrast with SMA and Spitzer observations. Radiative transfer modeling indicates that many of these scattered light images are consistent with a smooth spatial distribution for micron-sized grains, with little discontinuity in the surface density of the micron-sized grains at the cavity edge. Here we present a generic disk model that can simultaneously account for the general features in IRS, SMA, and Subaru observations. Particularly, the scattered light images for this model are computed, which agree with the general trend seen in Subaru data. Decoupling between the spatial distributions of the micron-sized dust and mm-sized dust inside the cavity is suggested by the model, which, if confirmed, necessitates a mechanism, such as dust filtration, for differentiating the small and big dust in the cavity clearing process. Our model also suggests an inwardly increasing gas-to-dust-ratio in the inner disk, and different spatial distributions for the small dust inside and outside the cavity, echoing the predictions in grain coagulation and growth models.

Belonging to the group of B[e] stars, V921 Scorpii is associated with a strong infrared excess and permitted and forbidden line emission, indicating the presence of low- and high-density circumstellar gas and dust. Many aspects of V921 Sco and other B[e] stars still remain mysterious, including their evolutionary state and the physical conditions resulting in the class-defining characteristics. In this paper, we employ VLTI/AMBER spectro-interferometry in order to reconstruct high-resolution (lambda/2B=0.0013") model-independent interferometric images for three wavelength bands around 1.65, 2.0, and 2.3 micrometer. In our images, we discover a close (25.0+/-0.8 milliarcsecond, corresponding to 29+/-0.9 AU at 1.15 kpc) companion around V921 Sco. Between two epochs in 2008 and 2009, we measure orbital motion of 7 degrees, implying an orbital period of about 35 years (for a circular orbit). Around the primary star, we detect a disk-like structure with indications for a radial temperature gradient. The polar axis of this AU-scale disk is aligned with the arcminute-scale bipolar nebula in which V921 Sco is embedded. Using Magellan/IMACS imaging, we detect multi-layered arc-shaped sub-structure in the nebula, suggesting episodic outflow activity from the system with a period of about 25 years, roughly matching the estimated orbital period of the companion. Our study supports the hypothesis that the B[e] phenomenon is related to dynamical interaction in a close binary system.

We present an interferometric kinematic study of morphologically complex protostellar envelopes based on observations of the dense gas tracers N2H+ and NH3. The strong asymmetric nature of most envelopes in our sample leads us to question the common interpretation of velocity gradients as rotation, given the possibility of projection effects in the observed velocities. Several "idealized" sources with well-ordered velocity fields and envelope structures are now analyzed in more detail. We compare the interferometric data to position-velocity diagrams of kinematic models for spherical rotating collapse and filamentary rotating collapse. For this purpose, we developed a filamentary parametrization of the rotating collapse model to explore the effects of geometric projection on the observed velocity structures. We find that most envelopes in our sample have PV structures that can be reproduced by an infalling filamentary envelope projected at different angles within the plane of the sky. The infalling filament produces velocity shifts across the envelope that can mimic rotation, especially when viewed at single-dish resolutions and the axisymmetric rotating collapse model does not uniquely describe any dataset. Furthermore, if the velocities are assumed to reflect rotation, then the inferred centrifugal radii are quite large in most cases, indicating significant fragmentation potential or more likely another component to the line-center velocity. We conclude that ordered velocity gradients cannot be interpreted as rotation alone when envelopes are non-axisymmetric and that projected infall velocities likely dominate the velocity field on scales larger than 1000 AU.

We present two-dimensional hydrodynamic simulations of self-gravitating protostellar disks subject to axisymmetric infall from envelopes and irradiation from the central star, to explore disk fragmentation due to gravitational instability (GI), and the fragmented clump evolution. We assume that the disk is built gradually and smoothly by the infall, resulting in good numerical convergence. We confirm that for disks around solar-mass stars, infall at high rates at radii beyond ~50 AU leads to disk fragmentation. At lower infall rates <1e-5 Msun/yr, however, irradiation suppresses fragmentation. We find that, once formed, the fragments or clumps migrate inward on typical type-I time scales of ~2e3 yr initially, but later migration deviates from the type-I time scale when the clump becomes more massive than the local disk mass, and/or when they starts to open gaps. As they migrate, the clumps accrete from the disk at a rate 1e-3 to 1e-1 MJupiter/yr, consistent with analytic estimates that assume a 1-2 Hill radii cross section. Most clumps can grow to their isolation masses >0.1 Msun quickly. The eventual fates of these clumps, however, diverges depending on the migration speed: 3 out of 13 clumps become massive enough (brown dwarf mass) to open gaps in the disk and essentially stop migrating; 4 are tidally destroyed during inward migration; 6 migrate across the inner simulated disk boundary. A simple analytic model for clump evolution is derived to explain these different fates. Overall, our results indicate that fast migration, accretion, and tidal destruction of the clumps pose challenges to the scenario of giant planet formation by GI in situ, but may provide a formation mechanism for close binary systems.

Most young stars in nearby molecular clouds have estimated ages of 1-2 Myr, suggesting that star formation is rapid. However, small numbers of stars in these regions with inferred ages of >= 5-10 Myr have been cited to argue that star formation is instead a slow, quasi-static process. When considering these alternative pictures it is important to recognize that the age spread in a given star-forming cloud is necessarily an upper limit to the timescales of local collapse, as not all spatially-distinct regions will start contracting at precisely the same instant. Moreover, star-forming clouds may dynamically evolve on timescales of a few Myr; in particular, global gravitational contraction will tend to yield increasing star formation rates with time due to generally increasing local gas densities. We show that two different numerical simulations of dynamic, flow-driven molecular cloud formation and evolution 1) predict age spreads for the main stellar population roughly consistent with observations, and 2) raise the possibility of forming small numbers of stars early in cloud evolution, before global contraction concentrates the gas and the bulk of the stellar population is produced. In general, the existence of a small number of older stars among a generally much-younger population is consistent with the picture of dynamic star formation, and may even provide clues to the time evolution of star-forming clouds.

We present a study of dense molecular gas kinematics in seventeen nearby protostellar systems using single-dish and interferometric molecular line observations. The non-axisymmetric envelopes around a sample of Class 0/I protostars were mapped in the N2H+ (J=1-0) tracer with the IRAM 30m, CARMA and PdBI as well as NH3 (1,1) with the VLA. The molecular line emission is used to construct line-center velocity and linewidth maps for all sources to examine the kinematic structure in the envelopes on spatial scales from 0.1 pc to ~1000 AU. The direction of the large-scale velocity gradients from single-dish mapping is within 45 degrees of normal to the outflow axis in more than half the sample. Furthermore, the velocity gradients are often quite substantial, the average being ~2.3 km\s\pc. The interferometric data often reveal small-scale velocity structure, departing from the more gradual large-scale velocity gradients. In some cases, this likely indicates accelerating infall and/or rotational spin-up in the inner envelope; the median velocity gradient from the interferometric data is ~10.7 km/s/pc. In two systems, we detect high-velocity HCO+ (J=1-0) emission inside the highest-velocity \nthp\ emission. This enables us to study the infall and rotation close to the disk and estimate the central object masses. The velocity fields observed on large and small-scales are more complex than would be expected from rotation alone, suggesting that complex envelope structure enables other dynamical processes (i.e. infall) to affect the velocity field.

Low-mass protostars may accrete most of their material through short-lived episodes of rapid disk accretion; yet until recently evolutionary tracks for these protostars assumed only constant or slowly-varying accretion. Important initial steps toward examining the potential effects of rapid accretion were recently made by Baraffe, Chabrier, & Gallardo, who showed that in the limit of low-temperature ("cold") accretion, protostars may have much smaller radii than found in previous treatments. Such small radii at the end of protostellar accretion would have the effect of making some young stars appear much older - perhaps as much as 10 Myr - than they really are. However, we argue that very rapid disk accretion is unlikely to be cold, because observations of the best-studied pre-main sequence disks with rapid disk accretion outbursts - the FU Ori objects - have spectral energy distributions which imply large, not small, protostellar radii. In addition, theory indicates models that at high accretion rates, protostellar disks become internally hot and geometrically thick, making it much more likely that hot material is added to the star. In addition, the very large luminosity of the accretion disk is likely to irradiate the central star strongly, heating up the outer layers and potentially expanding them. Nevertheless, the Baraffe et al. calculations emphasize the importance of initial protostellar radii for subsequent evolution.

It has been recently shown that molecular clouds do not exhibit a unique shape for the column density probability distribution function (Npdf). Instead, clouds without star formation seem to possess a lognormal distribution, while clouds with active star formation develope a power-law tail at high column densities. The lognormal behavior of the Npdf has been interpreted in terms of turbulent motions dominating the dynamics of the clouds, while the power-law behavior occurs when the cloud is dominated by gravity. In the present contribution we use thermally bi-stable numerical simulations of cloud formation and evolution to show that, indeed, these two regimes can be understood in terms of the formation and evolution of molecular clouds: a very narrow lognormal regime appears when the cloud is being assembled. However, as the global gravitational contraction occurs, the initial density fluctuations are enhanced, resulting, first, in a wider lognormal Npdf, and later, in a power-law Npdf. We thus suggest that the observed Npdf of molecular clouds are a manifestation of their global gravitationally contracting state. We also show that, contrary to recent suggestions, the exact value of the power-law slope is not unique, as it depends on the projection in which the cloud is being observed.

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

We present new X-ray and Far Ultraviolet (FUV) observations of T Tauri stars covering the age range 1 to 10 Myr. Our goals are to observationally constrain the intensity of radiation fields responsible for evaporating gas from the circumstellar disk and to assess the feasibility of current photoevaporation models, focusing on X-ray and UV radiation. We greatly increase the number of 7-10 Myr old T Tauri stars observed in the X-rays by including observations of the well populated 25 Ori aggregate in the Orion OB1a subassociation. With these new 7-10 Myr objects, we confirm that X-ray emission remains constant from 1-10 Myr. We also show, for the first time, observational evidence for the evolution of FUV radiation fields with a sample of 56 accreting and non-accreting young stars spanning 1 Myr to 1 Gyr. We find that the FUV emission decreases on timescales consistent with the decline of accretion in classical T Tauri stars until reaching the chromospheric level in weak T Tauri stars and debris disks. Overall, we find that the observed strength of high energy radiation is consistent with that required by photoevaporation models to dissipate the disks in timescales of approximately 10 Myr. Finally, we find that the high energy fields that affect gas evolution are not similarly affecting dust evolution; in particular, we find that disks with inner clearings, the transitional disks, have similar levels of FUV emission as full disks.

We use 2D hydrodynamic simulations of viscous disks to examine whether dynamically-interacting multiple giant planets can explain the large gaps (spanning over one order of magnitude in radius) inferred for the transitional and pre-transitional disks around T Tauri stars. In the absence of inner disk dust depletion, we find that it requires three to four giant planets to open up large enough gaps to be consistent with inferences from spectral energy distributions, because the gap width is limited by the tendency of the planets to be driven together into 2:1 resonances. With very strong tidal torques and/or rapid planetary accretion, fewer planets can also generate a large cavity interior to the locally formed gap(s) by preventing outer disk material from moving in. In these cases, however, the reduction of surface density produces a corresponding reduction in the inner disk accretion rate onto the star; this makes it difficult to explain the observed accretion rates of the pre/transitional disks. We find that even with four planets in disks, additional substantial dust depletion is required to explain observed disk gaps/holes. Substantial dust settling and growth, with consequent significant reductions in optical depths, is inferred for typical T Tauri disks in any case, and an earlier history of dust growth is consistent with the hypothesis that pre/transitional disks are explained by the presence of giant planets. We conclude that the depths and widths of gaps, and disk accretion rates in pre/transitional disks cannot be reproduced by a planet-induced gap opening scenario alone. Significant dust depletion is also required within the gaps/holes. Order of magnitude estimates suggest the mass of small dust particles (<1 micron) relative to the gas must be depleted to 1e-5 -- 1e-2 of the interstellar medium value, implying a very efficient mechanism of small dust removal or dust growth.

Using the Spitzer Infrared Spectrograph, we have performed mid-infrared spectroscopy on the young binary brown dwarf 2MASS J04414489+2301513 (15 AU) in the Taurus star-forming region. The spectrum exhibits excess continuum emission that likely arises from a circumstellar disk around the primary. Silicate emission is not detected in these data, indicating the presence of significant grain growth. This is one of the few brown dwarf disks at such a young age (~1 Myr) that has been found to lack silicate emission. To quantitatively constrain the properties of the disk, we have compared the spectral energy distribution of 2MASS J04414489+2301513 to the predictions of our vertical structure codes for irradiated accretion disks. Our models suggest that the remaining atmospheric grains of moderately depleted layers may have grown to a size of $\gtrsim5$ micron. In addition, our model fits indicate an outer radius of 0.2-0.3 AU for the disk. The small size of this circumprimary disk could be due to truncation by the secondary. The absence of an outer disk containing a reservoir of small, primordial grains, combined with a weak turbulent mechanism, may be responsible for the advanced grain growth in this disk.

We discuss the nature of the velocity dispersion vs. size relation for molecular clouds. In particular, we add to previous observational results showing that the velocity dispersions in molecular clouds and cores are not purely functions of spatial scale but involve surface gas densities as well. We emphasize that hydrodynamic turbulence is required to produce the first condensations in the progenitor medium. However, as the cloud is forming, it also becomes bound, and gravitational accelerations dominate the motions. Energy conservation in this case implies $|E_g| \sim E_k$, in agreement with observational data, and providing an interpretation for two recent observational results: the scatter in the $\delta v-R$ plane, and the dependence of the velocity dispersion on the surface density ${\delta v^2/ R} \propto \Sigma$. We argue that the observational data are consistent with molecular clouds in a state of hierarchical gravitational collapse, i.e., developing local centers of collapse throughout the whole cloud while the cloud itself is collapsing, and making equilibrium unnecessary at all stages prior to the formation of actual stars. Finally, we discuss how this mechanism need not be in conflict with the observed star formation rate.

We present IRAC/MIPS Spitzer Space Telescope observations of the solar type and the low mass stellar population of the young (~5 Myr) Lambda Orionis cluster. Combining optical and 2MASS photometry, we identify 436 stars as probable members of the cluster. Given the distance (450 pc) and the age of the cluster, our sample ranges in mass from 2 solar mass to objects below the substellar limit. With the addition of the Spitzer mid-infrared data, we have identified 49 stars bearing disks in the stellar cluster. Using spectral energy distribution (SED) slopes, we place objects in several classes: non-excess stars (diskless), stars with optically thick disks, stars with ``evolved disks''( with smaller excesses than optically thick disk systems), and ``transitional disks'' candidates (in which the inner disk is partially or fully cleared). The disk fraction depends on the stellar mass, ranging from ~6% for K type stars (Rc-J<2) to ~27% for stars with spectral type M5 or later (Rc-J>4). We confirm the dependence of disk fraction on stellar mass in this age range found in other studies. Regarding clustering levels, the overall fraction of disks in the Lambda Orionis cluster is similar to those reported in other stellar groups with ages normally quoted as ~5 Myr.

We present high-resolution L'-band imaging of the inner scattered light structure of Class 0 protostar L1527 IRS (IRAS 04368+2557) taken with the Gemini North telescope. The central point-source like feature seen in Spitzer Space Telescope IRAC images is resolved in the Gemini image into a compact bipolar structure with a narrow dark lane in the center. Two scattered light lobes are extended ~1.8" (200 AU) perpendicular to the direction of the outflow and ~2.5" (350 AU) along the outflow axis; the narrow dark lane between the scattered light lobes is ~0.45" (60 AU) thick. The observations are consistent with our initial modeling of a bright inner cavity separated by a dark lane due to extinction along the line of sight of the central protostar by the disk (Tobin et al. 2008). The bright, compact scattered light might be due to complex inner structure generated by the outflow, as suggested in our first paper, or it may more likely be the upper layers of the disk forming from infalling matter.

Extinction maps at 8 micron from the Spitzer Space Telescope show that many Class 0 protostars exhibit complex, irregular, and on-axisymmetric structure within the densest regions of their dusty envelopes. Many of the systems have highly irregular and on-axisymmetric morphologies on scales $\sim$1000 AU, with a quarter of the sample exhibiting filamentary or flattened dense structures. Complex envelope structure is observed in regions spatially distinct from outflow cavities, and the densest structures often show no systematic alignment perpendicular to the cavities. We suggest that the observed envelope complexity is the result of collapse from protostellar cores with initially non-equilibrium structures. The striking non-axisymmetry in many envelopes could provide favorable conditions for the formation of binary systems. We then show that the kinematics around L1165 as probed with N2H+ are indicative of asymmetric infall; the velocity gradient is not perpendicular to the outflow.

We report a set of numerical experiments aimed at addressing the applicability of competitive accretion to explain the high-mass end of the stellar initial mass function in a sheet geometry with shallow gravitational potential, in contrast to most previous simulations which have assumed formation in a cluster gravitational potential. Our flat cloud geometry is motivated by models of molecular cloud formation due to large-scale flows in the interstellar medium. The experiments consisted of SPH simulations of gas accretion onto sink particles formed rapidly from Jeans-unstable dense clumps placed randomly in the finite sheet. These simplifications allow us to study accretion with a minimum of free parameters, and to develop better statistics on the resulting mass spectra. We considered both clumps of equal mass and gaussian distributions of masses, and either uniform or spatially-varying gas densities. In all cases, the sink mass function develops a power law tail at high masses, with $dN/dlog M \propto M^{-\Gamma}$. The accretion rates of individual sinks follow $\dot{M} \propto M^2$ at high masses; this results in a continual flattening of the slope of the mass function towards an asymptotic form $\Gamma \sim 1$ (where the Salpeter slope is $\Gamma = 1.35$). The asymptotic limit is most rapidly reached when starting from a relatively broad distribution of initial sink masses. In general the resulting upper mass slope is correlated with the maximum sink mass; higher sink masses are found in simulations with flatter upper mass slopes. Although these simulations are of a highly idealized situation, the results suggest that competitive accretion may be relevant in a wider variety of environments than previously considered, and in particular that the upper mass distribution may generally evolve towards a limiting value of $\Gamma \sim 1$.

Spectral observations of classical T Tauri stars show a wide range of line profiles, many of which reveal signs of matter inflow and outflow. Halpha is the most commonly observed line profile due to its intensity, and it is highly dependent on the characteristics of the surrounding environment of these stars. Our aim is to analyze how the Halpha line profile is affected by the various parameters of our model which contains both the magnetospheric and disk wind contributions to the Halpha flux. We used a dipolar axisymmetric stellar magnetic field to model the stellar magnetosphere and a modified Blandford & Payne model was used in our disk wind region. A three-level atom with continuum was used to calculate the required Hydrogen level populations. We use the Sobolev approximation and a ray-by-ray method to calculate the integrated line profile. Through an extensive study of the model parameter space, we have investigated the contribution of many of the model parameters on the calculated line profiles. Our results show that the Halpha line is strongly dependent on the densities and temperatures inside the magnetosphere and the disk wind region. The bulk of the flux comes, most of the time, from the magnetospheric component for standard classical T Tauri stars parameters, but the disk wind contribution becomes more important as the mass accretion rate, the temperatures and densities inside the disk wind increase. We have also found that most of the disk wind contribution to the Halpha line is emitted at the innermost region of the disk wind. Models that take into consideration both inflow and outflow of matter are a necessity to fully understand and describe classical T Tauri stars.

We have analyzed nearly all images of the Taurus star-forming region at 3.6-24um that were obtained during the cryogenic mission of the Spitzer Space Telescope (46 deg^2) and have measured photometry for all known members of the region that are within these data, corresponding to 348 sources. We have classified the members of Taurus according to whether they show evidence of disks and envelopes (classes I, II, and III). The disk fraction in Taurus is 75% for solar-mass stars and declines to 45% for low-mass stars and brown dwarfs (0.01-0.3 M_sun). This dependence on stellar mass is similar to that measured for Cha I, although the disk fraction in Taurus is slightly higher overall, probably because of its younger age (1 vs. 2-3 Myr). In comparison, the disk fraction for solar-mass stars is much lower (20%) in IC 348 and Sigma Ori, which are denser than Taurus and Cha I and are roughly coeval with the latter. These data indicate that disk lifetimes for solar-mass stars are longer in regions that have lower stellar densities. Through an analysis of multiple epochs of photometry that are available for ~200 Taurus members, we find that stars with disks exhibit significantly greater mid-IR variability than diskless stars. Finally, we have used our data in Taurus to refine the criteria for primordial, evolved, and transitional disks. The number ratio of evolved and transitional disks to primordial disks in Taurus is 15/98 for K5-M5, indicating a timescale of 0.15 x tau(primordial)=0.45 Myr for the clearing of the inner regions of optically thick disks. After applying the same criteria to older clusters (2-10 Myr), we find that the proportions of evolved and transitional disks in those populations are consistent with the measurements in Taurus when their star formation histories are properly taken into account. ERRATUM: In Table 7, we inadvertently omitted the spectral type bins in which class II sources were placed in Table 8 based on their bolometric luminosities (applies only to stars that lack spectroscopic classifications). The bins were K6-M3.5 for FT Tau, DK Tau B, and IRAS 04370+2559, M3.5-M6 for IRAS 04200+2759, IT Tau B, and ITG 1, and M6-M8 for IRAS 04325+2402 C. In addition, the values of K_s-[3.6] in Table 13 and Figure 26 for spectral types of M4-M9 are incorrect. We present corrected versions of Table 13 and Figure 26.

Upcoming surveys with the Herschel Space Observatory will yield far-IR photometry of large samples of young stellar objects, which will require careful interpretation. We investigate the color and luminosity diagnostics based on Herschel broad-band filters to identify and discern the properties of low-mass protostars. We compute a grid of 2,016 protostars in various physical congurations, present the expected flux densities and flux density ratios for this grid of protostars, and compare Herschel observations of three protostars to the model results. These provide useful constraints on the range of colors and fluxes of protostar in the Herschel filters. We find that Herschel data alone is likely a useful diagnostic of the envelope properties of young stars

The NGC 1999 reflection nebula features a dark patch with a size of ~10,000 AU, which has been interpreted as a small, dense foreground globule and possible site of imminent star formation. We present Herschel PACS far-infrared 70 and 160mum maps, which reveal a flux deficit at the location of the globule. We estimate the globule mass needed to produce such an absorption feature to be a few tenths to a few Msun. Inspired by this Herschel observation, we obtained APEX LABOCA and SABOCA submillimeter continuum maps, and Magellan PANIC near-infrared images of the region. We do not detect a submillimer source at the location of the Herschel flux decrement; furthermore our observations place an upper limit on the mass of the globule of ~2.4x10^-2 Msun. Indeed, the submillimeter maps appear to show a flux depression as well. Furthermore, the near-infrared images detect faint background stars that are less affected by extinction inside the dark patch than in its surroundings. We suggest that the dark patch is in fact a hole or cavity in the material producing the NGC 1999 reflection nebula, excavated by protostellar jets from the V 380 Ori multiple system.

We present 70 and 160 micron Herschel science demonstration images of a field in the Orion A molecular cloud that contains the prototypical Herbig-Haro objects HH 1 and 2, obtained with the Photodetector Array Camera and Spectrometer (PACS). These observations demonstrate Herschel's unprecedented ability to study the rich population of protostars in the Orion molecular clouds at the wavelengths where they emit most of their luminosity. The four protostars previously identified by Spitzer 3.6-40 micron imaging and spectroscopy are detected in the 70 micron band, and three are clearly detected at 160 microns. We measure photometry of the protostars in the PACS bands and assemble their spectral energy distributions (SEDs) from 1 to 870 microns with these data, Spitzer spectra and photometry, 2MASS data, and APEX sub-mm data. The SEDs are fit to models generated with radiative transfer codes. From these fits we can constrain the fundamental properties of the protostars. We find luminosities in the range 12-84 L_sun and envelope densities spanning over two orders of magnitude. This implies that the four protostars have a wide range of envelope infall rates and evolutionary states: two have dense, infalling envelopes, while the other two have only residual envelopes. We also show the highly irregular and filamentary structure of the cold dust and gas surrounding the protostars as traced at 160 microns.

As an initial investigation into the long-term evolution of protostellar disks, we explore the conditions required to explain the large outbursts of disk accretion seen in some young stellar objects. We use one-dimensional time-dependent disk models with a phenomenological treatment of the magnetorotational instability (MRI) and gravitational torques to follow disk evolution over long timescales. Comparison with our previous two-dimensional disk model calculations (Zhu et al. 2009b, Z2009b) indicates that the neglect of radial effects and two-dimensional disk structure in the one-dimensional case makes only modest differences in the results; this allows us to use the simpler models to explore parameter space efficiently. We find that the mass infall rates typically estimated for low-mass protostars generally result in AU-scale disk accretion outbursts, as predicted by our previous analysis (Zhu et al. 2009a,Z2009a). We also confirm quasi-steady accretion behavior for high mass infall rates if the values of $\alpha$-parameter for the magnetorotational instability is small, while at this high accretion rate convection from the thermal instability may lead to some variations. We further constrain the combinations of the $\alpha$-parameter and the MRI critical temperature, which can reproduce observed outburst behavior. Our results suggest that dust sublimation may be connected with full activation of the MRI. This is consistent with the idea that small dust captures ions and electrons to suppress the MRI. In a later paper we will explore both long-term outburst and disk evolution with this model, allowing for infall from protostellar envelopes with differing angular momenta.

We use one-dimensional two-zone time-dependent accretion disk models to study the long-term evolution of protostellar disks subject to mass addition from the collapse of a rotating cloud core. Our model consists of a constant surface density magnetically coupled active layer, with transport and dissipation in inactive regions only via gravitational instability. We start our simulations after a central protostar has formed, containing ~ 10% of the mass of the protostellar cloud. Subsequent evolution depends on the angular momentum of the accreting envelope. We find that disk accretion matches the infall rate early in the disk evolution because much of the inner disk is hot enough to couple to the magnetic field. Later infall reaches the disk beyond ~10 AU, and the disk undergoes outbursts of accretion in FU Ori-like events as described in Zhu et al. 2009c. If the initial cloud core is moderately rotating most of the central star's mass is built up by these outburst events. Our results suggest that the protostellar "luminosity problem" is eased by accretion during these FU Ori-like outbursts. After infall stops the disk enters the T Tauri phase. An outer, viscously evolving disk has structure that is in reasonable agreement with recent submillimeter studies and its surface density evolves from $\Sigma \propto R^{-1}$ to $R^{-1.5}$. An inner, massive belt of material-- the "dead zone" -- would not have been observed yet but should be seen in future high angular resolution observations by EVLA and ALMA. This high surface density belt is a generic consequence of low angular momentum transport efficiency at radii where the disk is magnetically decoupled, and would strongly affect planet formation and migration.

We use archived IRAC images from the Spitzer Space Telescope to show that many Class 0 protostars exhibit complex, irregular, and non-axisymmetric structure within their dusty envelopes. Our 8 $\mu$m extinction maps probe some of the densest regions in these protostellar envelopes. Many of the systems are observed to have highly irregular and non-axisymmetric morphologies on scales >= 1000 AU, with a quarter of the sample exhibiting filamentary or flattened dense structures. Complex envelope structure is observed in regions spatially distinct from outflow cavities, and the densest structures often show no systematic alignment perpendicular to the cavities. These results indicate that mass ejection is not responsible for much of the irregular morphologies we detect; rather, we suggest that the observed envelope complexity is mostly the result of collapse from protostellar cores with initially non-equilibrium structures. The striking non-axisymmetry in many envelopes could provide favorable conditions for the formation of binary systems. We also note that protostars in the sample appear to be formed preferentially near the edges of clouds or bends in filaments, suggesting formation by gravitational focusing.

We present accretion rates for a large number of solar-type stars in the Cep OB2 region, based on U band observations. Our study comprises 95 members of the ~4 Myr-old cluster Tr 37 (including 20 "transition" objects; TO), as well as the only CTTS in the ~12 Myr-old cluster NGC 7160. The stars show different disk morphologies, with the majority of them having evolved and flattened disks. The typical accretion rates are about one order of magnitude lower than in regions aged 1-2 Myr, and we find no strong correlation between disk morphology and accretion rates. Although half of the TO are not accreting, the median accretion rates of normal CTTS and accreting "transition" disks are similar (~3 10^{-9} and 2 10^{-9} Msun/yr, respectively). Comparison with other regions suggests that the TO observed at different ages do not necessarily represent the same type of objects, which is consistent with the fact that the different processes that can lead to reduced IR excess/inner disk clearing (e.g., binarity, dust coagulation/settling, photoevaporation, giant planet formation) do not operate on the same timescales. Accreting TO in Tr 37 are probably suffering strong dust coagulation/settling. Regarding the equally large number of non-accreting TO in the region, other processes, like photoevaporation, the presence of stellar/substellar companions, and/or giant planet formation may account for their "transitional" SEDs and negligible accretion rates.

We present observations of the Class 0 protostar L1157-mm using the Combined Array for Research in Millimeter-wave Astronomy (CARMA) in 3 mm dust continuum and N2H+ line emission. In the N2H+ line, we detect a large-scale envelope extended over a linear size of ~20,000AU flattened in the direction perpendicular to the outflow. This N2H+ feature coincides with the outer envelope seen in the 8 micron extinction by Looney et al. Meanwhile, the dust continuum traces the compact, nearly spherical structure of the inner envelope, where N2H+ becomes depleted. This highly flattened N2H+ envelope also shows dynamical signatures consistent with gravitational infall in the inner region, but a slow, solid-body rotation at large scales. This flattened structure is not a rotationally supported circumstellar disk; instead, it resembles a prestellar core both morphologically and kinematically, representing the early phase of a Class 0 system. In this paper, we construct a simple model to interpret both the dust continuum and N2H+ emission and suggest a possible dynamical scenario for the overall properties of the envelope.

We present IRAC/MIPS Spitzer observations of intermediate-mass stars in the 5 Myr old Lambda Orionis cluster. In a representative sample of stars earlier than F5 (29 stars), we find a population of 9 stars with a varying degree of moderate 24um excess comparable to those produced by debris disks in older stellar groups. As expected in debris disks systems, those stars do not exhibit emission lines in their optical spectra. We also include in our study the star HD 245185, a known Herbig Ae object which displays excesses in all Spitzer bands and shows emission lines in its spectrum. We compare the disk population in the Lambda Orionis cluster with the disk census in other stellar groups studied using similar methods to detect and characterize their disks and spanning a range of ages from 3 Myr to 10 Myr. We find that for stellar groups of 5 Myr or older the observed disk frequency in intermediate mass stars (with spectral types from late B to early F) is higher than in low mass stars (with spectral types K and M). This is in contradiction with the observed trend for primordial disks evolution, in which stars with higher stellar masses dissipate their primordial disks faster. At 3 Myr the observed disk frequency in intermediate mass stars is still lower than for low mass stars indicating that second generation dusty disks start to dominate the disk population at 5 Myr for intermediate mass stars. This result agrees with recent models of evolution of solids in the region of the disk where icy objects form (>30 AU), which suggest that at 5-10 Myr collisions start to produce large amount of dust during the transition from runaway to oligarchic growth (reaching sizes of ~500 km) and then dust production peaks at 10-30 Myr, when objects reach their maximum sizes (>1000 km)

Recent models of molecular cloud formation and evolution suggest that such clouds are dynamic and generally exhibit gravitational collapse. We present a simple analytic model of global collapse onto a filament and compare this with our numerical simulations of the flow-driven formation of an isolated molecular cloud to illustrate the supersonic motions and infall ram pressures expected in models of gravity-driven cloud evolution. We apply our results to observations of the Pipe Nebula, an especially suitable object for our purposes as its low star formation activity implies insignifcant perturbations from stellar feedback. We show that our collapsing cloud model can explain the magnitude of the velocity dispersions seen in the $^{13}$CO filamentary structure by Onishi et al. and the ram pressures required by Lada et al. to confine the lower-mass cores in the Pipe nebula. We further conjecture that higher-resolution simulations will show small velocity dispersions in the densest core gas, as observed, but which are infall motions and not supporting turbulence. Our results point out the inevitability of ram pressures as boundary conditions for molecular cloud filaments, and the possibility that especially lower-mass cores still can be accreting mass at significant rates, as suggested by observations.

We analyze samples of Spitzer Infrared Spectrograph (IRS) spectra of T Tauri stars in the Ophiuchus, Taurus, and Chamaeleon I star-forming regions, whose median ages lie in the <1 to 2 Myr range. The median mid-infrared spectra of objects in these three regions are similar in shape, suggesting, on average, similar disk structures. When normalized to the same stellar luminosity, the medians follow each other closely, implying comparable mid-infrared excess emission from the circumstellar disks. We use the spectral index between 13 and 31 micron and the equivalent width of the 10 micron silicate emission feature to identify objects whose disk configuration departs from that of a continuous, optically thick accretion disk. Transitional disks, whose steep 13-31 micron spectral slope and near-IR flux deficit reveal inner disk clearing, occur with about the same frequency of a few percent in all three regions. Objects with unusually large 10 micron equivalent widths are more common (20-30%); they could reveal the presence of disk gaps filled with optically thin dust. Based on their medians and fraction of evolved disks, T Tauri stars in Taurus and Chamaeleon I are very alike. Disk evolution sets in early, since already the youngest region, the Ophiuchus core (L1688), has more settled disks with larger grains. Our results indicate that protoplanetary disks show clear signs of dust evolution at an age of a few Myr, even as early as ~1 Myr, but age is not the only factor determining the degree of evolution during the first few million years of a disk's lifetime.

We analyze the far-ultraviolet (FUV) spectra of 33 classical T Tauri stars (CTTS), including 20 new spectra obtained with the Advanced Camera for Surveys Solar Blind Channel (ACS/SBC) on the Hubble Space Telescope. Of the sources, 28 are in the ~1 Myr old Taurus-Auriga complex or Orion Molecular Cloud, 4 in the 8-10 Myr old Orion OB1a complex and one, TW Hya, in the 10 Myr old TW Hydrae Association. We also obtained FUV ACS/SBC spectra of 10 non-accreting sources surrounded by debris disks with ages between 10 and 125 Myr. We use a feature in the FUV spectra due mostly to electron impact excitation of \h2 to study the evolution of the gas in the inner disk. We find that the \h2 feature is absent in non-accreting sources, but is detected in the spectra of CTTS and correlates with accretion luminosity. Since all young stars have active chromospheres which produce strong X-ray and UV emission capable of exciting \h2 in the disk, the fact that the non-accreting sources show no \h2 emission implies that the \h2 gas in the inner disk has dissipated in the non-accreting sources, although dust (and possibly gas) remains at larger radii. Using the flux at 1600 {\AA}, we estimate that the column density of \h2 left in the inner regions of the debris disks in our sample is less than ~ 3x10^-6 g cm^-2, nine orders of magnitude below the surface density of the minimum mass solar nebula at 1 AU.

We have used Spitzer/IRAC to conduct a photometric monitoring program of the IC1396A dark globule in order to study the mid-IR (3.6 - 8 micron) variability of the heavily embedded Young Stellar Objects (YSOs) present in that area. We obtained light curves covering a 14 day timespan with a twice daily cadence for 69 YSOs, and continuous light curves with approximately 12 second cadence over 7 hours for 38 YSOs. Typical accuracies for our relative photometry were 1-2% for the long timespan data and a few mmag, corresponding to less than 0.5%, for the 7 hour continuous "staring-mode" data. More than half of the YSOs showed detectable variability, with amplitudes from ~0.05 mag to ~0.2 mag. About thirty percent of the YSOs showed quasi-sinusoidal light curve shapes with apparent periods from 5-12 days and light curve amplitudes approximately independent of wavelength over the IRAC bandpasses. We have constructed models which simulate the time dependent spectral energy distributions of Class I and I I YSOs in order to attempt to explain these light curves. Based on these models, the apparently periodic light curves are best explained by YSO models where one or two high latitude photospheric spots heat the inner wall of the circumstellar disk, and where we view the disk at fairly large inclination angle. Disk inhomogeneities, such as increasing the height where the accretion funnel flows to the stellar hotspot, enhances the light curve modulations. The other YSOs in our sample show a range of light curve shapes, some of which are probably due to varying accretion rate or disk shadowing events. One star, IC1396A-47, shows a 3.5 hour periodic light curve; this object may be a PMS Delta Scuti star.

We have developed time-dependent models of FU Ori accretion outbursts to explore the physical properties of protostellar disks. Our two-dimensional, axisymmetric models incorporate full vertical structure with a new treatment of the radiative boundary condition for the disk photosphere. We find that FU Ori-type outbursts can be explained by a slow accumulation of matter due to gravitational instability. Eventually this triggers the magnetorotational instability, which leads to rapid accretion. The thermal instability is triggered in the inner disk but this instability is not necessary for the outburst. An accurate disk vertical structure, including convection, is important for understanding the outburst behavior. Large convective eddies develop during the high state in the inner disk. The models are in agreement with Spitzer IRS spectra and also with peak accretion rates and decay timescales of observed outbursts, though some objects show faster rise timescale. We also propose that convection may account for the observed mild-supersonic turbulence and the short-timescale variations of FU Orionis objects.

Motivated by kinematic observations of young embedded clusters, this paper explores possible kinematic signatures produced by asphericity and departures from initial virial equilibrium in these systems. Specifically, the kinematic quantity that is measured and calculated in this study is the distribution of the line-of-sight velocities as a function of position along the cluster. Although clusters are found within a wide range of sizes, we focus on the regime with stellar membership $N \sim 10^3$. The gravitational potential of these young clusters is dominated by the gas, and the geometry of the gas distribution is generalized to include axisymmetric (and triaxial) forms. With this loss of symmetry, the kinematic results thus depend on viewing angle. This work also considers the stars to begin their trajectories with subvirial speeds, as indicated by observations of core motions in such clusters. Our results determine the conditions necessary for the kinematic signature to display interesting structure, i.e., a non-spherical potential, a viewing angle that is not along one of the principal axes, and subvirial starting conditions. We characterize the effects on this signature due to projection angle, initial stellar velocities, cluster elongation, and star formation efficiency. Finally, we compare our theoretical results to recent kinematic observations of the Orion Nebula Cluster; we find that the observations can be explained provided that the cluster is non-spherical, starts with subvirial initial velocities, and is not viewed along a principal axis.

We present a kinematic study of the Orion Nebula Cluster based upon radial velocities measured by multi-fiber echelle spectroscopy at the 6.5 meter MMT and Magellan telescopes. Velocities are reported for 1613 stars, with multi-epoch data for 727 objects as part of our continuing effort to detect and analyze spectroscopic binaries. We confirm and extend the results of Furesz et al. showing that the ONC is not relaxed, consistent with its youth, and that the stars generally follow the position-velocity structure of the moderate density gas in the region, traced by $^{13}$CO. The additional radial velocities we have measured enable us to probe some discrepancies between stellar and gaseous structure which can be attributed to binary motion and the inclusion of non-members in our kinematic sample. Our multi-epoch data allow us to identify 89 spectroscopic binaries; more will be found as we continue monitoring. Our results reinforce the idea that the ONC is a cluster in formation, and thus provides a valuable testing ground for theory. In particular, our observations are not consistent with the quasi-equilibrium or slow contraction models of cluster formation, but are consistent with cold collapse models.

The emission of FU Orionis objects in outburst has been identified as arising in rapidly accreting protoplanetary disks, based on a number of observational properties. A fundamental test of the accretion disk scenario is that the differentially rotating disk spectrum should produce a variation of rotational velocity with the wavelength of observation, as spectra taken at longer wavelengths probe outer, more slowly rotating disk regions. Previous observations of FU Ori have shown smaller rotation at near-infrared (~ 2.2 micron) wavelengths than observed at optical (~ 0.6 micron) wavelengths consistent with the assumption of Keplerian rotation. Here we report a spectrum from the Phoenix instrument on Gemini South which shows that differential (slower) rotation continues to be observed out to ~ 5 micron. The observed spectrum is well matched by the prediction of our accretion disk model previously constructed to match the observed spectral energy distribution and the differential rotation at wavelengths < 2.2 micron. This kinematic result allows us to confirm our previous inference of a large outer radius (~ 1 AU) for the rapidly accreting region of the FU Ori disk, which presents difficulties for outburst models relying purely on thermal instability. While some optical spectra have been interpreted to pose problems for the disk interpretation of FU Ori, we show that the adjustment of the maximum effective temperature of the disk model, proposed in a previous paper, greatly reduces these difficulties.

We present a set of numerical simulations addressing the effects of magnetic field strength and orientation on the flow-driven formation of molecular clouds. Fields perpendicular to the flows sweeping up the cloud can efficiently prevent the formation of massive clouds but permit the build-up of cold, diffuse filaments. Fields aligned with the flows lead to substantial clouds, whose degree of fragmentation and turbulence strongly depends on the background field strength. Adding a random field component leads to a "selection effect" for molecular cloud formation: high column densities are only reached at locations where the field component perpendicular to the flows is vanishing. Searching for signatures of colliding flows should focus on the diffuse, warm gas, since the cold gas phase making up the cloud will have lost the information about the original flow direction because the magnetic fields redistribute the kinetic energy of the inflows.

Observations indicate that mass accretion rates onto low-mass protostars are generally lower than the rates of infall to their disks; this suggests that much of the protostellar mass must be accreted during rare, short outbursts of rapid accretion. We explore when protostellar disk accretion is likely to be highly variable. While constant $\alpha$ disks can in principle adjust their accretion rates to match infall rates, protostellar disks are unlikely to have constant $\alpha$. In particular we show that neither models with angular momentum ransport due solely to the magnetorotational instability (MRI) nor ravitational instability (GI) are likely to transport disk mass at rotostellar infall rates over the large range of radii needed to move infalling envelope material down to the central protostar. We show that the MRI and GI are likely to combine to produce outbursts of rapid accretion starting at a few AU. Our analysis is consistent with the time-dependent models of Armitage, Livio, & Pringle (2001) and agrees with our observational study of the outbursting object FU Ori.

Here we present the Spitzer IRS spectrum of CVSO 224, the sole transitional disk located within the ~10 Myr old 25 Orionis group in Orion OB1a. A model fit to the spectral energy distribution of this object indicates a ~7 AU inner disk hole that contains a small amount of optically thin dust. In previous studies, CVSO 224 had been classified as a weak-line T Tauri star based on its Halpha equivalent width, but here we find an accretion rate of 7x10^-11 Msun/yr based on high-resolution Hectochelle data. CVSO 224's low mass accretion rate is in line with photoevaporative clearing theories. However, the Spitzer IRS spectrum of CVSO 224 has a substantial mid-infrared excess beyond 20microns which indicates that it is surrounded by a massive outer disk. Millimeter measurements are necessary to constrain the mass of the outer disk around CVSO 224 in order to confirm that photoevaporation is not the mechanism behind creating its inner disk hole.

We present the Spitzer Infrared Spectrograph (IRS) spectrum of SR20, a 5--10 AU binary T Tauri system in the $\rho$ Ophiuchi star forming region. The spectrum has features consistent with the presence of a disk; however, the continuum slope is steeper than the $\lambda^{-4/3}$ slope of an infinite geometrically thin, optically thick disk, indicating that the disk is outwardly truncated. Comparison with photometry from the literature shows a large increase in the mid-infrared flux from 1993 to 1996. We model the spectral energy distribution and IRS spectrum with a wall + optically thick irradiated disk, yielding an outer radius of 0.39$_{+0.03}^{-0.01}$ AU, much smaller than predicted by models of binary orbits. Using a two temperature $\chi^2$ minimization model to fit the dust composition of the IRS spectrum, we find the disk has experienced significant grain growth: its spectrum is well-fit using opacities of grains larger than 1 $\mu$m. We conclude that the system experienced a significant gravitational perturbation in the 1990s.

We present spectra of a sample of Herbig Ae and Be (HAeBe) stars obtained with the Infrared Spectrograph on the Spitzer Space Telescope. All but one of the Herbig stars show emission from polycyclic aromatic hydrocarbons (PAHs) and seven of the spectra show PAH emission, but no silicate emission at 10 microns. The central wavelengths of the 6.2, 7.7--8.2, and 11.3 micron emission features decrease with stellar temperature, indicating that the PAHs are less photo-processed in cooler radiation fields. The apparent low level of photo processing in HAeBe stars, relative to other PAH emission sources, implies that the PAHs are newly exposed to the UV-optical radiation fields from their host stars. HAeBe stars show a variety of PAH emission intensities and ionization fractions, but a narrow range of PAH spectral classifications based on positions of major PAH feature centers. This may indicate that, regardless of their locations relative to the stars, the PAH molecules are altered by the same physical processes in the proto-planetary disks of intermediate-mass stars. Analysis of the mid-IR spectral energy distributions indicates that our sample likely includes both radially flared and more flattened/settled disk systems, but we do not see the expected correlation of overall PAH emission with disk geometry. We suggest that the strength of PAH emission from HAeBe stars may depend not only on the degree of radial flaring, but also on the abundance of PAHs in illuminated regions of the disks and possibly on the vertical structure of the inner disk as well.

We characterize the crystalline silicate content and spatial distribution of small dust grains in a large sample of protoplanetary disks in the Taurus-Auriga young cluster, using Spitzer Space Telescope mid-infrared spectra. In turn we use the results to analyze the evolution of structure and composition of these 1-2 Myr-old disks around Solar- and later-type young stars, and test the standard models of dust processing which result in the conversion of originally amorphous dust into minerals. We find strong evidence of evolution of the dust crystalline mass fraction in parallel with that of the structure of the disks, in the sense that increasing crystalline mass fraction is strongly linked to dust settling to the disk midplane. We also confirm that the crystalline silicates are confined to small radii, r < 10 AU. However, we see no significant correlation of crystalline mass fraction with stellar mass or luminosity, stellar accretion rate, disk mass, or disk/star mass ratio, as would be expected in the standard models of dust processing based upon photo-evaporation and condensation close to the central star, accretion-heating-driven annealing at r < 1 AU, or spiral-shock heating at r < 10 AU, with or without effective radial mixing mechanisms. Either another grain-crystallizing mechanism dominates over these, or another process must be at work within the disks to erase the correlations they produce. We propose one of each sort that seem to be worth further investigation, namely X-ray heating and annealing of dust grains, and modulation of disk structure by giant-planetary formation and migration.

We demonstrate that the observationally inferred rapid onset of star formation after parental molecular clouds have assembled can be achieved by flow-driven cloud formation of atomic gas, using our previous three-dimensional numerical simulations. We post-process these simulations to approximate CO formation, which allows us to investigate the times at which CO becomes abundant relative to the onset of cloud collapse. We find that global gravity in a finite cloud has two crucial effects on cloud evolution. (a) Lateral collapse (perpendicular to the flows sweeping up the cloud) leads to rapidly increasing column densities above the accumulation from the one-dimensional flow. This in turn allows fast formation of CO, allowing the molecular cloud to ``appear'' rapidly. (b) Global gravity is required to drive the dense gas to the high pressures necessary to form solar-mass cores, in support of recent analytical models of cloud fragmentation. While the clouds still appear ``supersonically turbulent'', this turbulence is relegated to playing a secondary role, in that it is to some extent a consequence of gravitational forces.

We have performed a census of disks around brown dwarfs in the Sigma Ori cluster using all available images from IRAC onboard the Spitzer Space Telescope. To search for new low-mass cluster members with disks, we have measured photometry for all sources in the Spitzer images and have identified the ones that have red colors that are indicative of disks. We present 5 promising candidates, which may consist of 2 brown dwarfs, 2 stars with edge-on disks, and a low-mass protostar if they are bona fide members. Spectroscopy is needed to verify the nature of these sources. We have also used the Spitzer data to determine which of the previously known probable members of Sigma Ori are likely to have disks. By doing so, we measure disk fractions of ~40% and ~60% for low-mass stars and brown dwarfs, respectively. These results are similar to previous estimates of disk fractions in IC 348 and Cha I, which have roughly the same median ages as Sigma Ori (3 Myr). Finally, we note that our photometric measurements and the sources that we identify as having disks differ significantly from those of other recent studies that analyzed the same Spitzer images. For instance, previous work has suggested that the T dwarf S Ori 70 is redder than typical field dwarfs, which has been cited as possible evidence of youth and cluster membership. However, we find that this object is only slightly redder than the reddest field dwarfs in [3.6]-[4.5] (1.56+/-0.07 vs. 0.93-1.46). We measure a larger excess in [3.6]-[5.8] (1.75+/-0.21 vs. 0.87-1.19), but the flux at 5.8um may be overestimated because of the low signal-to-noise ratio of the detection. Thus, the Spitzer data do not offer strong evidence of youth and membership for this object, which is the faintest and coolest candidate member of Sigma Ori that has been identified to date.

The mid- to far-infrared emission of the outbursting FU Orionis objects has been attributed either to a flared outer disk or to an infalling envelope. We revisit this issue using detailed radiative transfer calculations to model the recent, high signal-to-noise data from the IRS instrument on the {Spitzer Space Telescope}. In the case of FU Ori, we find that a physically-plausible flared disk irradiated by the central accretion disk matches the observations. Building on our previous work, our accretion disk model with outer disk irradiation by the inner disk reproduces the spectral energy distribution between ~4000 angstroms to ~40 microns. Our model is consistent with near-infrared interferometry but there are some inconsistencies with mid-infared interferometric results. Including the outer disk allows us to refine our estimate of the outer radius of the outbursting, high mass accretion rate disk in FU Ori as ~ 0.5 AU, which is a crucial parameter in assessing theories of the FU Orionis phenomenon. We are able to place an upper limit on the mass infall rate of any remnant envelope infall rate to ~ 7e-7 Msun/yr assuming a centrifugal radius of 200 AU. The FUor BBW 76 is also well modelled by a 0.6 AU inner disk and a flared outer disk. However, V1515 Cyg requires an envelope with an outflow cavity to adequately reproduce the IRS spectrum. In contrast with the suggestion by Green et al., we do not require a flattened envelope to match the observations; the inferred cavity shape is qualitatively consistent with typical protostellar envelopes. This variety of dusty structures suggests that the FU Orionis phase can be present at either early or late stages of protostellar evolution.

We present new Spitzer Space Telescope observations of stars in the young ~5 Myr gamma Velorum stellar cluster. Combining optical and 2MASS photometry, we have selected 579 stars as candidate members of the cluster. With the addition of the Spitzer mid-infrared data, we have identified 5 debris disks around A-type stars, and 5-6 debris disks around solar-type stars, indicating that the strong radiation field in the cluster does not completely suppress the production of planetesimals in the disks of cluster members. However, we find some evidence that the frequency of circumstellar primordial disks is lower, and the IR flux excesses are smaller than for disks around stellar populations with similar ages. This could be evidence for a relatively fast dissipation of circumstellar dust by the strong radiation field from the highest mass star(s) in the cluster. Another possibility is that gamma Velorum stellar cluster is slightly older than reported ages and the the low frequency of primordial disks reflects the fast disk dissipation observed at ~5 Myr.

The observed rapid onset of star formation in molecular clouds requires rapid formation of dense fragments which can collapse individually before being overtaken by global gravitationally-driven flows. Many previous investigations have suggested that supersonic turbulence produces the necessary fragmentation, without addressing however the source of this turbulence. Motivated by our previous (numerical) work on the flow-driven formation of molecular clouds, we investigate the expected timescales of the dynamical and thermal instabilities leading to the rapid fragmentation of gas swept up by large-scale flows, and compare them with global gravitational collapse timescales. We identify parameter regimes in gas density, temperature and spatial scale within which a given instability will dominate. We find that dynamical instabilities disrupt large-scale coherent flows via generation of turbulence, while strong thermal fragmentation amplifies the resulting low-amplitude density perturbations, thus leading to small-scale, high-density fragments as seeds for {\em local} gravity to act upon. Global gravity dominates only on the largest scales; large-scale gravitationally-driven flows promote the formation of groups and clusters of stars formed by turbulence, thermal fragmentation, and rapid cooling.

We report the initial results of a large-scale optical-near infrared survey to extend the known young population of the entire Orion star-forming region down to the substellar domain. Using deep optical I-band photometry and data from the 2MASS survey, we selected candidates across ~14.8 deg^2 in the ~8 Myr old Ori OB1a subassociation and over ~6.7 deg^2 in the Ori OB1b subassociation (age ~3), with completeness down to 0.05Mo and 0.072Mo respectively. We obtained low resolution optical spectra for a subsample of 4 candidates in Ori OB1a and 26 in Ori OB1b; as a result we confirmed 3 new members in Ori OB1a, one of which is substellar, and 19 new members in Ori OB1b, out of which 7 are at the substellar limit and 5 are substellar. We looked into the presence of accretion signatures by measuring the strength of the Ha line in emission. Accordingly, we classified the new members as having Classical T-Tauri star (CTTS) or Weak Lined T Tauri star-like (WTTS) nature. We found that all the new members confirmed in Ori OB1a are WTTSs, while 39 +25/-22 % of the new members in Ori OB1b exhibit CTTS-like behavior, suggestive of ongoing accretion from a circum(sub)stellar disk. Additionally we found that none of the members confirmed in OB1a show near-IR color excess while 38 +26/-21 % of OB1b members show H-K color excess. These results are consistent with recent findings for low mass young stars in Orion OB1. The similarity in CTTS-like properties and near-IR excess across the substellar boundary gives support to the idea of a common formation mechanism for low mass stars and at least the most massive brown dwarfs. Finally, we remark the discovery of two new members classified as CTTSs, both exhibiting W(Ha) < -140 A, suggesting significant ongoing accretion.

We present a census of circumstellar disks in the Chamaeleon I star-forming region. Using the Infrared Array Camera and the Multiband Imaging Photometer onboard the Spitzer Space Telescope, we have obtained images of Chamaeleon I at 3.6, 4.5, 5.8, 8.0, and 24 um. To search for new disk-bearing members of the cluster, we have performed spectroscopy on objects that have red colors in these data. Through this work, we have discovered four new members of Chamaeleon I with spectral types of M4, M6, M7.5, and L0. The first three objects are highly embedded (A_J~5) and reside near known protostars, indicating that they may be among the youngest low-mass sources in the cluster (<1 Myr). The L0 source is the coolest known member of Chamaeleon I. Its luminosity implies a mass of 0.004-0.01 M_sun, making it the least massive brown dwarf for which a circumstellar disk has been reliably detected. To characterize the disk population in Chamaeleon I, we have classified the infrared spectral energy distributions of the 203 known members that are encompassed by the Spitzer images. Through these classifications, we find that the disk fraction in Chamaeleon I is roughly constant at ~50% from 0.01 to 0.3 M_sun. These data are similar to the disk fraction of IC 348, which is a denser cluster at the same age as Chamaeleon I. However, the disk fraction at M>1 M_sun is significantly higher in Chamaeleon I than in IC 348 (65% vs. 20%), indicating longer disk lifetimes in Chamaeleon I for this mass range. Thus, low-density star-forming regions like Chamaeleon I may offer more time for planet formation around solar-type stars than denser clusters.

We model Spitzer Space Telescope observations of the Taurus Class 0 protostar L1527 IRS (IRAS 04368+2557) to provide constraints on its protostellar envelope structure. The nearly edge-on inclination of L1527 IRS, coupled with the highly spatially-resolved near to mid-infrared images of this object and the detailed IRS spectrum, enable us to constrain the outflow cavity geometry quite well, reducing uncertainties in the other derived parameters. The mid-infrared scattered light image shows a bright central source within a dark lane; the aspect ratio of this dark lane is such that it appears highly unlikely to be a disk shadow. In modeling this dark lane, we conclude that L1527 IRS is probably not described by a standard TSC envelope with simple bipolar cavities. We find it necessary to model the dark lane and central source as a modified inner envelope structure. This structure may be due either to a complex wind-envelope interaction or induced by the central binary. To fit the overall SED, we require the central source to have a large near to mid-infrared excess, suggesting substantial disk accretion. Our model reproduces the overall morphology and surface brightness distribution of L1527 IRS fairly well, given the limitations of using axisymmetric models to fit the non-axisymmetric real object, and the derived envelope infall rates are in reasonable agreement with some other investigations. IRAC observations of L1527 IRS taken 12 months apart show variability in total flux and variability in the opposing bipolar cavities, suggesting asymmetric variations in accretion. We also provide model images at high resolution for comparison to future observations with current ground-based instrumentation and future space-based telescopes.

Using SpeX at the NASA Infrared Telescope Facility and the Spitzer Infrared Spectrograph, we have obtained infrared spectra from 0.7 to 30um for three young brown dwarfs in the TW Hydra Association (~10 Myr), 2M J1207-3932, 2M J1139-3159, and SS J1102-3431. The spectral energy distribution for 2M J1139-3159 is consistent with a stellar photosphere for the entire wavelength range of our data while the other two objects exhibit significant excess emission at >5um. We are able to reproduce the excess emission from each brown dwarf using our models of irradiated accretion disks. According to our model fits, both disks have experienced a high degree of dust settling. We also find that silicate emission at 10 and 20um is absent from the spectra of these disks, indicating that grains in the upper disk layers have grown to sizes larger than ~5um. Both of these characteristics are consistent with previous observations of decreasing silicate emission with lower stellar masses and older ages. These trends suggest that either 1) the growth of dust grains, and perhaps planetesimal formation, occurs faster in disks around brown dwarfs than in disks around stars, or 2) the radii of the mid-IR-emitting regions of disks are smaller for brown dwarfs than for stars, and grains grow faster at smaller disk radii. Finally, we note the possible detection of an unexplained emission feature near 14um in the spectra of both of the disk-bearing brown dwarfs.

We study the formation of the calcium II infrared triplet lines 8498\AA, 8542\AA and 8662\AA, in the accreting magnetospheric flows of Classical T Tauri stars (CTTS), and present a grid of models for a large range of magnetospheric conditions. We apply our models to the interpretation of multi epoch observations of the CTTS DI Cep. We find that these lines form in the magnetospheric infall and that the variability of the CaII triplet lines in DI Cep can be explained in the context of changes in the mass accretion rate/temperature of the accretion column gas.

We present Spitzer Infrared Spectrograph spectra of 28 Class I protostars in the Taurus star-forming region. The 5 to 36 micron spectra reveal excess emission from the inner regions of the envelope and accretion disk surrounding these predecessors of low-mass stars, as well as absorption features due to silicates and ices. Together with shorter- and longer-wavelength data from the literature, we construct spectral energy distributions and fit envelope models to 22 protostars of our sample, most of which are well-constrained due to the availability of the IRS spectra. We infer that the envelopes of the Class I objects in our sample cover a wide range in parameter space, particularly in density and centrifugal radius, implying different initial conditions for the collapse of protostellar cores.

We present results from 1078 high resolution spectra of 990 stars in the young open cluster NGC 2264, obtained with the Hectochelle multiobject echelle spectrograph on the 6.5m MMT. We confirm 471 stars as members, based on their radial velocity and/or H-alpha emission. The radial velocity distribution of cluster members is non-Gaussian with a dispersion of approx 3.5 km/s. We find a substantial north-south velocity gradient and spatially coherent structure inthe radial velocity distribution, similar to that seen in the molecular gas in the region. Our results suggest that there are at least three distinguishable subclusters in NGC 2264, correlated with similar structure seen in 13CO emission, which is likely to be a remnant of initial structure in this very young cluster. We propose that this substructure is the result of gravitational amplification of initial inhomogeneities during overall collapse to a filamentary distribution of gas and stars, as found in simulations by Burkert & Hartman (2004).

We present results from 1351 high resolution spectra of 1215 stars in the Orion Nebula Cluster (ONC) and the surrounding Orion 1c association, obtained with the Hectochelle multiobject echelle spectrograph on the 6.5m MMT. We confirmed 1111 stars as members, based on their radial velocity and/or H-alpha emission. The radial velocity distribution of members shows a dispersion of 3.1 km/s. We found a substantial north-south velocity gradient and spatially coherent structure in the radial velocity distribution, similar to that seen in the molecular gas in the region. We also identified several binary and high velocity stars, a region exhibiting signs of triggered star formation, and a possible foreground population of stars somewhat older than the ONC. Stars without infrared excesses (as detected with the IRAC instrument on the Spitzer Space Telescope) exhibit a wider spread in radial velocity than the infrared excess stars; this spread is mostly due to a blue-shifted population of stars that may constitute a foreground population. We also identify some accreting stars, based on H-alpha, that do not have detectable infrared excesses with IRAC, and thus are potential transitional disk systems (objects with inner disk holes). We propose that the substructure seen both in stellar and gaseous component is the result of non-uniform gravitational collapse to a filamentary distribution of gas. The spatial and kinematic correlation between the stellar and gaseous components suggests the region is very young, probably only ~1 crossing time old or less to avoid shock dissipation and gravitational interactions which would tend to destroy the correlation between stars and gas.

The recently recognized class of "transitional disk" systems consists of young starswith optically-thick outer disks but inner disks which are mostly devoid of small dust. Here we introduce a further class of "pre-transitional disks" with significant near-infrared excesses which indicate the presence of an optically thick inner disk separated from an optically thick outer disk; thus, the spectral energy distributions of pre-transitional disks suggest the incipient development of disk gaps rather than inner holes. In UX Tau A, our analysis of the Spitzer IRS spectrum finds that the near-infrared excess is produced by an inner optically thick disk and a gap of ~56 AU is present. The Spitzer IRS spectrum of LkCa 15 is suggestive of a gap of ~46 AU, confirming previous millimeter imaging. In addition, UX Tau A contains crystalline silicates in its disk at radii >~ 56 AU which poses a challenge to our understanding of the production of this crystalline material. In contrast, LkCa 15's silicates are amorphous and pristine. UX Tau A and LkCa 15 increase our knowledge of the diversity of dust clearing in low-mass star formation.

We present optical, IR and millimeter observations of the solar-type star 13-277, also known as GM Cep, in the 4 Myr-old cluster Tr 37. GM Cep experiences rapid magnitude variations of more than 2 mag at optical wavelengths. We explore the causes of the variability, which seem to be dominated by strong increases in the accretion, being similar to EX-or episodes. The star shows high, variable accretion rates (up to ~10$^{-6}$ Msun/yr), signs of powerful winds, and it is a very fast rotator (Vsini~43 km/s). Its strong mid-IR excesses reveal a very flared disk and/or a remnant envelope, most likely out of hydrostatic equilibrium. The 1.3 millimeter fluxes suggest a relatively massive disk (Mdisk~0.1 Msun). Nevertheless, the millimeter mass is not enough to sustain increased accretion episodes over large timescales, unless the mass is underestimated due to significant grain growth. We finally explore the possibility of GM Cep having a binary companion, which could trigger disk instabilities producing the enhanced accretion episodes.

We study numerically the formation of molecular clouds in large-scale colliding flows including self-gravity. The models emphasize the competition between the effects of gravity on global and local scales in an isolated cloud. Global gravity builds up large-scale filaments, while local gravity -- triggered by a combination of strong thermal and dynamical instabilities -- causes cores to form. The dynamical instabilities give rise to a local focusing of the colliding flows, facilitating the rapid formation of massive protostellar cores of a few 100 M$_\odot$. The forming clouds do not reach an equilibrium state, though the motions within the clouds appear comparable to ``virial''. The self-similar core mass distributions derived from models with and without self-gravity indicate that the core mass distribution is set very early on during the cloud formation process, predominantly by a combination of thermal and dynamical instabilities rather than by self-gravity.

We present new Spitzer Space Telescope observations of two fields in the Orion OB1 association. We report here IRAC/MIPS observations for 115 confirmed members and 41 photometric candidates of the ~10 Myr 25 Orionis aggregate in the OB1a subassociation, and 106 confirmed members and 65 photometric candidates of the 5 Myr region located in the OB1b subassociation. The 25 Orionis aggregate shows a disk frequency of 6% while the field in the OB1b subassociation shows a disk frequency of 13%. Combining IRAC, MIPS and 2MASS photometry we place stars bearing disks in several classes: stars with optically thick disks (class II systems), stars with an inner transitional disks (transitional disk candidates) and stars with "evolved disks"; the last exhibit smaller IRAC/MIPS excesses than class II systems. In all, we identify 1 transitional disk candidate in the 25 Orionis aggregate and 3 in the OB1b field; this represents ~10% of the disk bearing stars, indicating that the transitional disk phase can be relatively fast. We find that the frequency of disks is a function of the stellar mass, suggesting a maximum around stars with spectral type M0. Comparing the infrared excess in the IRAC bands among several stellar groups we find that inner disk emission decays with stellar age, showing a correlation with the respective disk frequencies. The disk emission at the IRAC and MIPS bands in several stellar groups indicates that disk dissipation takes place faster in the inner region of the disks. Comparison with models of irradiated accretion disks, computed with several degrees of settling, suggests that the decrease in the overall accretion rate observed in young stellar groups is not sufficient to explain the weak disk emission observed in the IRAC bands for disk bearing stars with ages 5 Myr or older.

We have constructed a detailed radiative transfer disk model which reproduces the main features of the spectrum of the outbursting young stellar object FU Orionis from ~ 4000 angstrom, to ~ 8 micron. Using an estimated visual extinction Av~1.5, a steady disk model with a central star mass ~0.3 Msun and a mass accretion rate ~ 2e-4 Msun/yr, we can reproduce the spectral energy distribution of FU Ori quite well. With the mid-infrared spectrum obtained by the Infrared Spectrograph (IRS) on board the Spitzer Space Telescope, we estimate that the outer radius of the hot, rapidly accreting inner disk is ~ 1 AU using disk models truncated at this outer radius. Inclusion of radiation from a cooler irradiated outer disk might reduce the outer limit of the hot inner disk to ~ 0.5 AU. In either case, the radius is inconsistent with a pure thermal instability model for the outburst. Our radiative transfer model implies that the central disk temperature Tc > 1000 K out to ~ 0.5 - 1 AU, suggesting that the magnetorotational instability (MRI) can be supported out to that distance. Assuming that the ~ 100 yr decay timescale in brightness of FU Ori represents the viscous timescale of the hot inner disk, we estimate the viscosity parameter (alpha) to be ~ 0.2 - 0.02 in the outburst state, consistent with numerical simulations of MRI in disks. The radial extent of the high mass accretion region is inconsistent with the model of Bell & Lin, but may be consistent with theories incorporating both gravitational instability and MRI.

Here we present the Spitzer IRS spectrum of CS Cha, a member of the ~2 Myr old Chamaeleon star-forming region, which reveals an optically thick circumstellar disk truncated at ~43 AU, the largest hole modeled in a transitional disk to date. Within this inner hole, ~5x10^-5 lunar masses of dust are located in a small optically thin inner region which extends from 0.1 to 1 AU. In addition, the disk of CS Cha has bigger grain sizes and more settling than the previously modeled transitional disks DM Tau, GM Aur, and CoKu Tau/4, suggesting that CS Cha is in a more advanced state of dust evolution. The Spitzer IRS spectrum also shows [Ne II] 12.81 micron fine-structure emission with a luminosity of 1.3x10^29 ergs s^-1, indicating that optically thin gas is present in this ~43 AU hole, in agreement with H_alpha measurements and a UV excess which indicate that CS Cha is still accreting 1.2x10^-8 M_sun yr^-1. We do not find a correlation of the [Ne II] flux with L_X, however, there is a possible correlation with mass accretion rate, which if confirmed would suggest that EUV fluxes due to accretion are the main agent for formation of the [Ne II] line.

We present observations of a circumstellar disk that is inclined close to edge-on around a young brown dwarf in the Taurus star-forming region. Using data obtained with SpeX at the NASA Infrared Telescope Facility, we find that the slope of the 0.8-2.5 um spectrum of the brown dwarf 2MASS J04381486+2611399 cannot be reproduced with a photosphere reddened by normal extinction. Instead, the slope is consistent with scattered light, indicating that circumstellar material is occulting the brown dwarf. By combining the SpeX data with mid-IR photometry and spectroscopy from the Spitzer Space Telescope and previously published millimeter data from Scholz and coworkers, we construct the spectral energy distribution for 2MASS J04381486+2611399 and model it in terms of a young brown dwarf surrounded by an irradiated accretion disk. The presence of both silicate absorption at 10 um and silicate emission at 11 um constrains the inclination of the disk to be ~70 deg, i.e. ~20 deg from edge-on. Additional evidence of the high inclination of this disk is provided by our detection of asymmetric bipolar extended emission surrounding 2MASS J04381486+2611399 in high-resolution optical images obtained with the Hubble Space Telescope. According to our modeling for the SED and images of this system, the disk contains a large inner hole that is indicative of a transition disk (R_in~58 R_star~0.275 AU) and is somewhat larger than expected from embryo ejection models (R_out=20-40 AU vs. R_out<10-20 AU).

We present the mid-infrared spectrum, obtained with the Spitzer Infrared Spectrograph (IRS), of HD 98800, a quadruple star system located in the 10-Myr-old TW Hydrae association. It has a known mid-infrared excess that arises from a circumbinary disk around the B components of the system. The IRS spectrum confirms that the disk around HD 98800 B displays no excess emission below about 5.5 micron, implying an optically thick disk wall at 5.9 AU and an inner, cleared-out region; however, some optically thin dust, consisting mainly of 3-micron-sized silicate dust grains, orbits the binary in a ring between 1.5 and 2 AU. The peculiar structure and apparent lack of gas in the HD 98800 B disk suggests that this system is likely already at the debris disks stage, with a tidally truncated circumbinary disk of larger dust particles and an inner, second-generation dust ring, possibly held up by the resonances of a planet. The unusually large infrared excess can be explained by gravitational perturbations of the Aa+Ab pair puffing up the outer dust ring and causing frequent collisions among the larger particles.

We present multi-wavelength optical and infrared photometry of 170 previously known low mass stars and brown dwarfs of the 5 Myr Collinder 69 cluster (Lambda Orionis). The new photometry supports cluster membership for most of them, with less than 15% of the previous candidates identified as probable non-members. The near infrared photometry allows us to identify stars with IR excesses, and we find that the Class II population is very large, around 25% for stars (in the spectral range M0 - M6.5) and 40% for brown dwarfs, down to 0.04 Msun, despite the fact that the H(alpha) equivalent width is low for a significant fraction of them. In addition, there are a number of substellar objects, classified as Class III, that have optically thin disks. The Class II members are distributed in an inhomogeneous way, lying preferentially in a filament running toward the south-east. The IR excesses for the Collinder 69 members range from pure Class II (flat or nearly flat spectra longward of 1 micron), to transition disks with no near-IR excess but excesses beginning within the IRAC wavelength range, to two stars with excess only detected at 24 micron. Collinder 69 thus appears to be at an age where it provides a natural laboratory for the study of primordial disks and their dissipation.

We make use of new near and mid-IR photometry of the Pleiades cluster in order to help identify proposed cluster members. We also use the new photometry with previously published photometry to define the single-star main sequence locus at the age of the Pleiades in a variety of color-magnitude planes. The new near and mid-IR photometry extend effectively two magnitudes deeper than the 2MASS All-Sky Point Source catalog, and hence allow us to select a new set of candidate very low mass and sub-stellar mass members of the Pleiades in the central square degree of the cluster. We identify 42 new candidate members fainter than Ks =14 (corresponding to 0.1 Mo). These candidate members should eventually allow a better estimate of the cluster mass function to be made down to of order 0.04 solar masses. We also use new IRAC data, in particular the images obtained at 8 um, in order to comment briefly on interstellar dust in and near the Pleiades. We confirm, as expected, that -- with one exception -- a sample of low mass stars recently identified as having 24 um excesses due to debris disks do not have significant excesses at IRAC wavelengths. However, evidence is also presented that several of the Pleiades high mass stars are found to be impacting with local condensations of the molecular cloud that is passing through the Pleiades at the current epoch.

Using the Infrared Array Camera and the Multiband Imaging Photometer aboard the {\it Spitzer Space Telescope}, we have obtained images of the Lupus 3 star-forming cloud at 3.6, 4.5, 5.8, 8.0, and 24 \micron. We present photometry in these bands for the 41 previously known members that are within our images. In addition, we have identified 19 possible new members of the cloud based on red 3.6-8.0 \micron colors that are indicative of circumstellar disks. We have performed optical spectroscopy on 6 of these candidates, all of which are confirmed as young low-mass members of Lupus 3. The spectral types of these new members range from M4.75 to M8, corresponding to masses of 0.2-0.03 $M_\odot$ for ages of $\sim1$ Myr according to theoretical evolutionary models. We also present optical spectroscopy of a candidate disk-bearing object in the vicinity of the Lupus 1 cloud, 2M 1541-3345, which Jayawardhana & Ivanov recently classified as a young brown dwarf ($M\sim0.03$ $M_\odot$) with a spectral type of M8. In contrast to their results, we measure an earlier spectral type of M5.75$\pm$0.25 for this object, indicating that it is probably a low-mass star ($M\sim0.1$ $M_\odot$). In fact, according to its gravity-sensitive absorption lines and its luminosity, 2M 1541-3345 is older than members of the Lupus clouds ($\tau\sim1$ Myr) and instead is probably a more evolved pre-main-sequence star that is not directly related to the current generation of star formation in Lupus.

We report here on the photometric and kinematic properties of a well defined group of nearly 200 low-mass pre-main sequence stars, concentrated within ~ 1 deg of the early-B star 25 Ori, in the Orion OB1a sub-association. We refer to this stellar aggregate as the 25 Orionis group. The group also harbors the Herbig Ae/Be star V346 Ori and a dozen other early type stars with photometry, parallaxes, and some with IR excess emission, consistent with group membership. The number of high and low-mass stars is in agreement with expectations from a standard Initial Mass Function. The velocity distribution for the young stars in 25 Ori shows a narrow peak centered at 19.7 km/s, very close to the velocity of the star 25 Ori. Our results provide new and compelling evidence that the 25 Ori group is a distinct kinematic entity, and that considerable space and velocity structure is present in the Ori OB1a sub-association. The low-mass members follow a well defined band in the color-magnitude diagram, consistent with an isochronal age of ~ 7-10 Myr, depending on the assumed evolutionary model. The highest density of members is located near the star 25 Ori, but the actual extent of the cluster cannot be well constrained with our present data. In a simple-minded kinematic evolution scenario, the 25 Ori group may represent the evolved counterpart of a younger aggregate like the sigma Ori cluster. The 25 Ori stellar aggregate is the most populous ~ 10 Myr sample yet known within 500 pc, setting it as an excellent laboratory to study the evolution of solar-like stars and protoplanetary disks.

We report new Spitzer Space Telescope observations from the IRAC and MIPS instruments of the young (~ 3 Myr) sigma Orionis cluster. We identify 336 stars as members of the cluster using optical and near-infrared color magnitude diagrams. Using the spectral energy distribution (SED) slopes in the IRAC spectral range, we place objects in several classes: non-excess stars, stars with optically thick disks(like classical T Tauri stars), class I (protostellar) candidates, and stars with ``evolved disks''; the last exhibit smaller IRAC excesses than optically thick disk systems. In general, this classification agrees with the location expected in IRAC-MIPS color-color diagrams for these objects. We find that the evolved disk systems are mostly a combination of objects with optically thick but non-flared disks, suggesting grain growth and/or settling, and transition disks, systems in which the inner disk is partially or fully cleared of small dust. In all, we identify 7 transition disk candidates and 3 possible debris disk systems. As in other young stellar populations, the fraction of disks depends on the stellar mass, ranging from ~10% for stars in the Herbig Ae/Be mass range (>2 msun) to ~35% in the T Tauri mass range (1-0.1 msun). We find that the disk fraction does not decrease significantly toward the brown dwarf candidates (<0.1 msun). The IRAC infrared excesses found in stellar clusters and associations with and without central high mass stars are similar, suggesting that external photoevaporation is not very important in many clusters. Finally, we find no correlation between the X-ray luminosity and the disk infrared excess, suggesting that the X-rays are not strongly affected by disk accretion.

We present the Spitzer IRS spectra for 33 young stars in Tr 37 and NGC 7160. The sample includes the high- and intermediate-mass stars with MIPS 24 microns excess, the only known active accretor in the 12 Myr-old cluster NGC 7160, and 19 low-mass stars with disks in the 4 Myr-old cluster Tr 37. We examine the 10 microns silicate feature, present in the whole sample of low-mass star and in 3 of the high- and intermediate-mass targets, and we find that PAH emission is detectable only in the Herbig Be star. We analyze the composition and size of the warm photospheric silicate grains by fitting the 10 microns silicate feature, and study the possible correlations between the silicate characteristics and the stellar and disk properties (age, SED slope, accretion rate, spectral type). We find indications of dust settling with age and of the effect of turbulent enrichment of the disk atmosphere with large grains. Crystalline grains are only small contributors to the total silicate mass in all disks, and do not seem to correlate with any other property, except maybe binarity. We also observe that spectra with very weak silicate emission are at least 3 times more frequent among M stars than among earlier spectral types, which may be an evidence of inner disk evolution. Finally, we find that 5 of the high- and intermediate-mass stars have SEDs and IRS spectra consistent with debris disk models involving planet formation, which could indicate debris disk formation at ages as early as 4 Myr.

We revisit the analysis of the Non-linear Thin Shell Instability (NTSI) numerically, including magnetic fields. The magnetic tension force is expected to work against the main driver of the NTSI -- namely transverse momentum transport. However, depending on the field strength and orientation, the instability may grow. For fields aligned with the inflow, we find that the NTSI is suppressed only when the Alfv\'en speed surpasses the (supersonic) velocities generated along the collision interface. Even for fields perpendicular to the inflow, which are the most effective at preventing the NTSI from developing, internal structures form within the expanding slab interface, probably leading to fragmentation in the presence of self-gravity or thermal instabilities. High Reynolds numbers result in local turbulence within the perturbed slab, which in turn triggers reconnection and dissipation of the excess magnetic flux. We find that when the magnetic field is initially aligned with the flow, there exists a (weak) correlation between field strength and gas density. However, for transverse fields, this correlation essentially vanishes. In light of these results, our general conclusion is that instabilities are unlikely to be erased unless the magnetic energy in clouds is much larger than the turbulent energy. Finally, while our study is motivated by the scenario of molecular cloud formation in colliding flows, our results span a larger range of applicability, from supernovae shells to colliding stellar winds.

We discuss problems with some observational estimates indicating long protostellar core lifetimes and large stellar age spreads in molecular clouds. We also point out some additional observational constraints which suggest that protostellar cores do not have long lifetimes before collapsing. For external galaxies, we argue that the widths of spiral arms does not imply a long star-formation process, since the formation of massive stars will disrupt molecular clouds, move material around, compress it in other regions which produce new star-forming clouds. Thus, it seems unavoidable that this cyclical process will result in an extended period of enhanced star formation, which does not represent the survival time of any individual molecular cloud. We argue that the rapid star formation indicated observationally is also easier to understand theoretically than the traditional scenario of slow quasi-static contraction with ambipolar diffusion.

We suggest that the Orion A cloud is gravitationally collapsing on large scales, and is producing the Orion Nebula Cluster due to the focusing effects of gravity acting within a finite cloud geometry. In support of this suggestion, we show how an elliptical rotating sheet of gas with a modest density gradient along the major axis can collapse to produce a structure qualitatively resembling Orion A, with a fan-shaped structure at one end, ridges or filaments along the fan, and a narrow curved filament at the other end reminiscent of the famous integral-shaped filament. The model produces a local concentration of mass within the narrow filament which in principle could form a dense cluster of stars like that of the Orion Nebula. We suggest that global gravitational contraction might be a more common feature of molecular clouds than previously recognized, and that the formation of star clusters is a dynamic process resulting from the focusing effects of gravity acting upon the geometry of finite clouds.

We present mid-infrared spectra of T Tauri stars in the Taurus star-forming region obtained with the Spitzer Infrared Spectrograph (IRS). For the first time, the 5-36 micron spectra of a large sample of T Tauri stars belonging to the same star-forming region is studied, revealing details of the mid-infrared excess due to dust in circumstellar disks. We analyze common features and differences in the mid-IR spectra based on disk structure, dust grain properties, and the presence of companions. Our analysis encompasses spectral energy distributions from the optical to the far-infrared, a morphological sequence based on the IRS spectra, and spectral indices in IRS wave bands representative of continuum emission. By comparing the observed spectra to a grid of accretion disk models, we infer some basic disk properties for our sample of T Tauri stars, and find additional evidence for dust settling.

We report new Spitzer observations of intermediate mass stars in two regions of the Orion OB1 association located in the subassociations OB1a ($\sim$10 Myr) and OB1b ($\sim$5 Myr). In a representative sample of stars earlier than F5 of both stellar groups, we find a population of stars surrounded of debris disks, without excess in the IRAC bands and without emission lines in their optical spectra, but with a varying degree of 24{\micron} excess. Comparing our samples with 24{\micron} observations of intermediate mass stars in other stellar groups, spanning a range of ages from 2.5 Myr to 150 Myr, we find that debris disks are more frequent and have larger 24{\micron} excess at 10 Myr (OB1a). This trend agrees with predictions of models of evolution of solids in the outer regions of disks ($>$30 AU), where large icy objects ($\sim$1000 Km) begin to form at $\sim$10 Myr; the presence of these objects in the disk initiates a collisional cascade, producing enough dust particles to explain the relatively large 24 {\micron} excess observed in OB1a. The dust luminosity observed in the stellar groups older than 10 Myr declines roughly as predicted by collisional cascade models. Combining Spitzer observations, optical spectra and 2MASS data, we found a new Herbig Ae/Be star (HD290543) and a star (HD36444) with a large 24 {\micron} excess, both in OB1b. This last object could be explained as a intermediate stage between HAeBe and true debris systems or as a massive debris disk produced by a collision between two large objects ($>$1000 Km).

Using the Hectochelle multifiber spectrograph, we have obtained high-resolution (R~34,000) spectra in the Halpha region for a large number of stars in the 4 Myr-old cluster Tr 37, containing 146 previously known members and 26 newly identified ones. We present the Halpha line profiles of all members, compare them to our IR observations of dusty disks (2MASS/JHK + IRAC + MIPS 24 micron), use the radial velocities as a membership criterion, and calculate the rotational velocities. We find a good correlation between the accretion-broadened profiles and the presence of protoplanetary disks, noting that a small fraction of the accreting stars presents broad profiles with Halpha equivalent widths smaller than the canonical limit separating CTTS and WTTS. The number of strong accretors appears to be lower than in younger regions, and a large number of CTTS have very small accretion rates (dM/dt<10^{-9} Msun/yr). Taking into account that the spectral energy distributions are consistent with dust evolution (grain growth/settling) in the innermost disk, this suggests a parallel evolution of the dusty and gaseous components. We also observe that about half of the "transition objects" (stars with no IR excesses at wavelengths shorter than ~6 micron) do not show any signs of active accretion, whereas the other half is accreting with accretion rates <10^{-9} Msun/yr. These zero or very low accretion rates reveal important gas evolution and/or gas depletion in the innermost disk, which could be related to grain growth up to planetesimal or even planet sizes. Finally, we examine the rotational velocities of accreting and non accreting stars, finding no significant differences that could indicate disk locking at these ages.

We present the first Spitzer Infrared Spectrograph (IRS; The IRS was a collaborative venture between Cornell University and Ball Aerospace Corporation funded by NASA through the Jet Propulsion Laboratory and the Ames Research Center.) observations of the disks around classical T Tauri stars: spectra in the 5.2-30 micron range of six stars. The spectra are dominated by emission features from amorphous silicate dust, and a continuous component from 5 to 8 microns that in most cases comprises an excess above the photosphere throughout our spectral range. There is considerable variation in the silicate feature/continuum ratio, which implies variations of inclination, disk flaring, and stellar mass accretion rate. In most of our stars, structure in the silicate feature suggests the presence of a crystalline component. In one, CoKu Tau/4, no excess above the photosphere appears at wavelengths shortward of the silicate features, similar to 10 Myr old TW Hya, Hen 3-600, and HR 4796A. This indicates the optically thick inner disk is largely absent. The silicate emission features with peaks at 9.7 and 18 microns indicate small dust grains are present. The extremely low 10-20 micron color temperature of the dust excess, 135 K, indicates these grains are located more than 10 AU from the star. These features are suggestive of gravitational influence by planets or close stellar companions and grain growth in the region within 10 AU of the star, somewhat surprising for a star this young (1 Myr).

Molecular Cloud Complexes (MCCs) are highly structured and ``turbulent''. Observational evidence suggests that MCCs are dynamically dominated systems, rather than quasi-equilibrium entities. The observed structure is more likely a consequence of the formation process rather than something that is imprinted after the formation of the MCC. Converging flows provide a natural mechanism to generate MCC structure. We present a detailed numerical analysis of this scenario. Our study addresses the evolution of a MCC from its birth in colliding atomic hydrogen flows up until the point when H$_2$ may begin to form. A combination of dynamical and thermal instabilities breaks up coherent flows efficiently, seeding the small-scale non-linear density perturbations necessary for local gravitational collapse and thus allowing (close to) instantaneous star formation. Many observed properties of MCCs come as a natural consequence of this formation scenario. Since converging flows are omnipresent in the ISM, we discuss the general applicability of this mechanism, from local star formation regions to galaxy mergers.

The 8-14 micron emission spectra of 12 T Tauri stars in the Taurus/Auriga dark clouds and in the TW Hydrae association obtained with the Infrared Spectrograph (IRS; The IRS is a collaborative venture between Cornell University and Ball Aerospace Corporation funded by NASA through the Jet Propulsion Laboratory and the Ames Research Center.) on board Spitzer are analyzed. Assuming the 10 micron features originate from silicate grains in the optically thin surface layers of T Tauri disks, the 8-14 micron dust emissivity for each object is derived from its Spitzer spectrum. The emissivities are fit with the opacities of laboratory analogs of cosmic dust. The fits include small nonspherical grains of amorphous silicates (pyroxene and olivine), crystalline silicates (forsterite and pyroxene), and quartz, together with large fluffy amorphous silicate grains. A wide range in the fraction of crystalline silicate grains as well as large silicate grains among these stars are found. The dust in the transitional-disk objects CoKu Tau/4, GM Aur, and DM Tau has the simplest form of silicates, with almost no hint of crystalline components and modest amounts of large grains. This indicates that the dust grains in these objects have been modified little from their origin in the interstellar medium. Other stars show various amounts of crystalline silicates, similar to the wide dispersion of the degree of crystallinity reported for Herbig Ae/Be stars of mass <2.5 solar masses. Late spectral type, low-mass stars can have significant fractions of crystalline silicate grains. Higher quartz mass fractions often accompany low amorphous olivine-to-amorphous pyroxene ratios. It is also found that lower contrast of the 10 micron feature accompanies greater crystallinity.

We present 5-35 $\mu$m spectra, taken with the Infrared Spectrograph (IRS) on the Spitzer Space Telescope, of five FU Orionis objects: FU Ori, V1515 Cyg, V1057 Cyg, BBW 76, and V346 Nor. All but V346 Nor reveal amorphous silicate grains in emission at 10 $\mu$m and 20 $\mu$m, and show water-vapor absorption bands at 5.8 and 6.8 $\mu$m and SiO or possibly methane absorption at 8 $\mu$m. These absorption features closely match these bands in model stellar photospheres -- signs of the gaseous photospheres of the inner regions of these objects' accretion disks. The continuum emission at 5-8 $\mu$m is also consistent with such disks, and, for FU Orionis and BBW 76, longer-wavelength emission may be fit by a model which includes moderate disk flaring. V1057 Cyg and V1515 Cyg have much more emission at longer wavelengths than the others, perhaps evidence of substantial remnant of their natal, infalling envelopes.

Observations of T Tauri stars and young brown dwarfs suggest that the accretion rates of their disks scale strongly with the central stellar mass, approximately $\mdot \propto M_*^2$. No dependence of accretion rate on stellar mass is predicted by the simplest version of the layered disk model of Gammie (1996), in which non-thermal ionization of upper disk layers allows accretion to occur via the magnetorotational instability. We show that a minor modification of Gammie's model to include heating by irradiation from the central star yields a modest dependence of $\mdot$ upon the mass of the central star. A purely viscous disk model could provide a strong dependence of accretion rate on stellar mass if the initial disk radius (before much viscous evolution has occurred) has a strong dependence on stellar mass. However, it is far from clear that at least the most massive pre-main sequence disks can be totally magnetically activated by X-rays or cosmic rays. We suggest that a combination of effects are responsible for the observed dependence, with the lowest-mass stars having the lowest mass disks, which can be thoroughly magnetically active, while the higher-mass stars have higher mass disks which have layered accretion and relatively inactive or ``dead'' central zones at some radii. In such dead zones, we suggest that gravitational instabilities may play a role in allowing accretion to proceed. In this connection, we emphasize the uncertainty in disk masses derived from dust emission, and argue that T Tauri disk masses have been systematically underestimated by conventional analyses. Further study of accretion rates, especially in the lowest-mass stars, would help to clarify the mechanisms of accretion in T Tauri stars.

Early observations of T Tauri stars suggested that stars with evidence of circumstellar accretion disks rotated slower than stars without such evidence, but more recent results are not as clear. Near-IR circumstellar disk indicators, though the most widely available, are subject to uncertainties that can result from inner disk holes and/or the system inclination. Mid-infrared observations are less sensitive to such effects, but until now, these observations have been difficult to obtain. The Spitzer Space Telescope now easily enables mid-infrared measurements of large samples of PMS stars covering a broad mass range in nearby star-forming regions. Megeath and collaborators surveyed the Orion Molecular Clouds (~1 Myr) with the IRAC instrument (3.6, 4.5, 5.8, 8 microns) as part of a joint IRAC and MIPS GTO program. We examine the relationship between rotation and Spitzer mid-IR fluxes for ~900 stars in Orion for stars between 3 and 0.1 Msun. We find in these Spitzer data the clearest indication to date that stars with longer periods are more likely than those with short periods to have IR excesses suggestive of disks.

We present observations of six late-type members of the young cluster IC 348 detected at 24 microns with the Multiband Imaging Photometer for Spitzer(MIPS). At least four of the objects are probably substellar. Combining these data with ground-based optical and near-infrared photometry and complementary observations with the Infrared Array Camera (IRAC), we have modeled the spectral energy distributions using detailed models of irradiated accretion disks. We are able to fit the observations with models using a range of maximum grain sizes from ISM-type dust to grains as large as 1 millimeter. Two objects show a lack of excess emission at wavelengths shortward of 5.8-8 microns but significant excess at longer wavelengths, indicative of large optically thin or evacuated inner holes. Our models indicate a inner hole of radius ~ 0.5-0.9 AU for the brown dwarf L316; this is the first brown dwarf with evidence for an AU-scale inner disk hole. We examine several possible mechanisms for the inner disk clearing in this case, including photoevaporation and planet formation.

We present new K-band long baseline interferometer observations of three young stellar objects of the FU Orionis class, V1057 Cyg, V1515 Cyg and Z CMa-SE, obtained at the Keck Interferometer during its commissioning science period. The interferometer clearly resolves the source of near-infrared emission in all three objects. Using simple geometrical models we derive size scales (0.5-4.5 AU) for this emission. All three objects appear significantly more resolved than expected from simple models of accretion disks tuned to fit the broadband optical and infrared spectro-photometry. We explore variations in the key parameters that are able to lower the predicted visibility amplitudes to the measured levels, and conclude that accretion disks alone do not reproduce the spectral energy distributions and K-band visibilities simultaneously. We conclude that either disk models are inadequate to describe the near-infrared emission, or additional source components are needed. We hypothesize that large scale emission (10s of AU) in the interferometer field of view is responsible for the surprisingly low visibilities. This emission may arise in scattering by large envelopes believed to surround these objects.

Using the Hubble Space Telescope, the 4 m Blanco telescope at the Cerro Tololo Inter-American Observatory, and the Spitzer Space Telescope, we have performed deep imaging from 0.8 to 8 um of the southern subcluster in the Chamaeleon I star-forming region. In these data, we have discovered an object, Cha 110913-773444, whose colors and magnitudes are indicative of a very low-mass brown dwarf with a circumstellar disk. In a near-infrared spectrum of this source obtained with the Gemini Near-Infrared Spectrograph, the presence of strong steam absorption confirms its late-type nature (>=M9.5) while the shapes of the H- and K-band continua and the strengths of the Na I and K I lines demonstrate that it is a young, pre-main-sequence object rather than a field dwarf. A comparison of the bolometric luminosity of Cha 110913-773444 to the luminosities predicted by the evolutionary models of Chabrier and Baraffe and Burrows and coworkers indicates a mass of 8+7/-3 M_Jup, placing it fully within the mass range observed for extrasolar planetary companions (M<=15 M_Jup). The spectral energy distribution of this object exhibits mid-infrared excess emission at >5 um, which we have successfully modeled in terms of an irradiated viscous accretion disk with M'<=10e-12 M_sun/year. Cha 110913-773444 is now the least massive brown dwarf observed to have a circumstellar disk, and indeed is one of the least massive free-floating objects found to date. These results demonstrate that the raw materials for planet formation exist around free-floating planetary-mass bodies.

We present new Spitzer Space Telescope observations of the young cluster NGC2264. Observations at 24 micron with the Multiband Imaging Photometer has enabled us to identify the most highly embedded and youngest objects in NGC2264. This letter reports on one particular region of NGC2264 where bright 24 micron sources are spatially configured in curious linear structures with quasi-uniform separations. The majority of these sources (~60% are found to be protostellar in nature with Class I spectral energy distributions. Comparison of their spatial distribution with sub-millimeter data from Wolf-Chase (2003) and millimeter data from Peretto et al. (2005) shows a close correlation between the dust filaments and the linear spatial configurations of the protostars, indicating that star formation is occurring primarily within dense dusty filaments. Finally, the quasi-uniform separations of the protostars are found to be comparable in magnitude to the expected Jeans length suggesting thermal fragmentation of the dense filamentary material.

We present near and mid-infrared photometry obtained with the Spitzer Space Telescope of 300 known members of the IC348 cluster. We merge this photometry with existing ground-based optical and near-infrared photometry in order to construct optical-infrared spectral energy distributions (SEDs) for all the cluster members and present a complete atlas of these SEDs. We employ these observations to both investigate the frequency and nature of the circumstellar disk population in the cluster. The observations are sufficiently sensitive to enable the first detailed measurement of the disk frequency for very low mass stars at the peak of the stellar IMF. Using measurements of infrared excess between 3.6 and 8 microns we find the total frequency of disk-bearing stars in the cluster to be 50 +/- 6%. However, only 30 +/- 4% of the member stars are surrounded by optically thick, primordial disks, while the remaining disk-bearing stars are surrounded by what appear to be optically thin,anemic disks. The disk fraction appears to be a function of spectral type and stellar mass. The disk longevity and thus conditions for planet formation appear to be most favorable for stars which are of comparable mass to the sun. The optically thick disks around later type (> M4) stars appear to be less flared than the disks around earlier type stars. This may indicate a greater degree of dust settling and a more advanced evolutionary state for the late M disk population. Finally we find that the presence of an optically thick dust disk is correlated with gaseous accretion as measured by the strength of Halpha emission. These results suggest that it is more likely for dust disks to persist in the absence of active gaseous accretion than for active accretion to persist in the absence of dusty disks.

We present self-consistent disk models for T Tauri stars which include a parameterized treatment of dust settling and grain growth, building on techniques developed in a series of papers by D'Alessio etal. The models incorporate depleted distributions of dust in upper disk layers along with larger-sized particles near the disk midplane, as expected theoretically and as we suggested earlier is necessary to account for mm-wave emission, SEDs, scattered light images, and silicate emission features simultaneously. By comparing the models with recent mid- and near-IR observations, we find that the dust to gas mass ratio of small grains at the upper layers should be < 10 % of the standard value. The grains that have disappeared from the upper layers increase the dust to gas mass ratio of the disk interior; if those grains grow to maximum sizes of the order of mm during the settling process, then both the millimeter-wave fluxes and spectral slopes can be consistently explained. Depletion and growth of grains can also enhance the ionization of upper layers, enhancing the possibility of the magnetorotational instability for driving disk accretion.

We present IRAC 3.6, 4.5, 5.8 and 8 micron photometry for the 17 A, K and M type members of the Eta Chameleontis association. These data show infrared excesses toward six of the 15 K and M stars, indicating the presence of circumstellar disks around 40% of the stars with masses of 0.1-1 solar mass. The two A-stars show no infrared excesses. The excess emission around one of the stars is comparable to the median excess for classical T Tauri stars in the Taurus association; the remaining five show comparatively weak excess emission. Taking into account published Halpha spectroscopy that shows that five of the six stars are accreting, we argue that the disks with weak mid-infrared excesses are disks in which the inner disks have been largely depleted of small grains by grain growth, or, in one case, the small grains have settled to the midplane. This suggests that Eta Cha has a much higher fraction of disks caught in the act of transitioning into optically thin disks than that measured in younger clusters and associations.

Using the Infrared Array Camera (IRAC) aboard the Spitzer Space Telescope, we have obtained mid-infrared photometry for 25 and 18 low-mass members of the IC 348 and Chamaeleon I star-forming clusters, respectively (>M6, M<=0.08 Msun). We find that 42+/-13% and 50+/-17% of the two samples exhibit excess emission indicative of circumstellar disks. In comparison, the disk fractions for stellar members of these clusters are 33+/-4% and 45+/-7% (M0-M6, 0.7 Msun>=M>=0.1 Msun). The similarity of the disk fractions of stars and brown dwarfs is consistent with a common formation mechanism and indicates that the raw materials for planet formation are available around brown dwarfs as often as around stars.

We present observations of the T Tauri stars BP Tau, DG Tau, DI Tau, GM Aur, LkCa 15, RW Aur and V830 Tau, using long baseline infrared interferometry at K band (2.2 microns) from the Keck Interferometer. The target sources have a range of mass accretion rates and excess near-infrared emission. The interferometer is most sensitive to extended emission on characteristic size scales of 1 to 5 millarcseconds. All sources show evidence for resolved K band emission on these scales, although a few of the sources are marginally consistent with being unresolved. We calculate the infrared excess based on fitting stellar photosphere models to the optical photometry and estimate the physical size of the emission region using simple geometric models for the sources with a significant infrared excess. Assuming that the K band resolved emission traces the inner edge of the dust disk, we compare the measured characteristic sizes to predicted dust sublimation radii and find that the models require a range of dust sublimation temperatures and possibly optical depths within the inner rim to match the measured radii.

Molecular clouds (MCs) are highly structured and ``turbulent''. Colliding gas streams of atomic hydrogen have been suggested as a possible source of MCs, imprinting the filamentary structure as a consequence of dynamical and thermal instabilities. We present a 2D numerical analysis of MC formation via converging HI flows. Even with modest flow speeds and completely uniform inflows, non-linear density perturbations as possible precursors of MCs arise. Thus, we suggest that MCs are inevitably formed with substantial structure, e.g., strong density and velocity fluctuations, which provide the initial conditions for subsequent gravitational collapse and star formation in a variety of galactic and extragalactic environments.

We present initial results from a survey of the Orion A and B molecular clouds made with the InfraRed Array Camera (IRAC) onboard the Spitzer Space Telescope. This survey encompasses a total of 5.6 square degrees with the sensitivity to detect objects below the hydrogen burning limit at an age of 1 Myr. These observations cover a number of known star forming regions, from the massive star forming clusters in the Orion Nebula and NGC 2024, to small groups of low mass stars in the L1641. We combine the IRAC photometry with photometry from the 2MASS point source catalog and use the resulting seven band data to identify stars with infrared excesses due to dusty disks and envelopes. Using the presence of an infrared excess as an indicator of youth, we show the distribution of young stars and protostars in the two molecular clouds. We find that roughly half of the stars are found in dense clusters surrounding the two regions of recent massive star formation in the Orion clouds, NGC 2024 and the Orion Nebula.

We present Very Large Array observations at 3.5 cm of the nearby young star TW Hya that show the emission is constant in time over weeks, months and years, and spatially resolved with peak brightness temperature ~10 K at ~0.25 (15 AU) resolution. These features are naturally explained if the emission mechanism at this wavelength is thermal emission from dust particles in the disk surrounding the star. To account quantitatively for the observations, we construct a self-consistent accretion disk model that incorporates a population of centimeter size particles that matches the long wavelength spectrum and spatial distribution. A substantial mass fraction of orbiting particles in the TW Hya disk must have agglomerated to centimeter size. These data provide the first clear indication that dust emission from protoplanetary disks may be observed at centimeter wavelengths, and that changes in the spectral slope of the dust emission may be detected, providing constraints on dust evolution and the planet formation process.

We carry out high resolution simulations of the inner regions of a wide angle wind driven bipolar outflow using an Adaptive Mesh Refinement code. Our code follows H-He gas with molecular, atomic and ionic components and the associated time dependent molecular chemistry and ionization dynamics with radiative cooling. Our simulations explore the nature of the outflow when a spherical wind expands into a rotating, collapsing envelope. We compare with key observational properties of the outflow system of Source I in the BN/KL region. Our calculations show that the wind evacuates a bipolar outflow cavity in the infalling envelope. We find the head of the outflow to be unstable and that it rapidly fragments into clumps. We resolve the dynamics of the strong shear layer which defines the side walls of the cavity. We conjecture that this layer is the likely site of maser emission and examine its morphology and rotational properties. The shell of swept up ambient gas that delineates the cavity edge retains its angular momentum. This rotation is roughly consistent with that observed in the Source I SiO maser spots. The observed proper motions and line-of-sight velocity are approximately reproduced by the model. The cavity shell at the base of the flow assumes an X-shaped morphology which is also consistent with Source I. We conclude that the wide opening angle of the outflow is evidence that a wide-angle wind drives the Source I outflow and not a collimated jet.

We present infrared photometry obtained with the IRAC camera on the Spitzer Space Telescope of a sample of 82 pre-main sequence stars and brown dwarfs in the Taurus star-forming region. We find a clear separation in some IRAC color-color diagrams between objects with and without disks. A few ``transition'' objects are noted, which correspond to systems in which the inner disk has been evacuated of small dust. Separating pure disk systems from objects with remnant protostellar envelopes is more difficult at IRAC wavelengths, especially for objects with infall at low rates and large angular momenta. Our results generally confirm the IRAC color classification scheme used in previous papers by Allen et al. and Megeath et al. to distinguish between protostars, T Tauri stars with disks, and young stars without (inner) disks. The observed IRAC colors are in good agreement with recent improved disk models, and in general accord with models for protostellar envelopes derived from analyzing a larger wavelength region. We also comment on a few Taurus objects of special interest. Our results should be useful for interpreting IRAC results in other, less well-studied star-forming regions.

We report the results of a sensitive K-band survey of Herbig Ae/Be disk sizes using the 85-m baseline Keck Interferometer. Targets were chosen to span the maximum range of stellar properties to probe the disk size dependence on luminosity and effective temperature. For most targets, the measured near-infrared sizes (ranging from 0.2 to 4 AU) support a simple disk model possessing a central optically-thin (dust-free) cavity, ringed by hot dust emitting at the expected sublimation temperatures (T_sub~1000-1500K). Furthermore, we find a tight correlation of disk size with source luminosity R propto L^(1/2) for Ae and late Be systems (valid over more than 2 decades in luminosity), confirming earlier suggestions based on lower-quality data. Interestingly, the inferred dust-free inner cavities of the highest luminosity sources (Herbig B0-B3 stars) are under-sized compared to predictions of the optically-thin cavity model, likely due to optically-thick gas within the inner AU.

Using the Infrared Array Camera (IRAC) aboard the Spitzer Space Telescope, we have obtained mid-infrared photometry of the least massive known brown dwarf in the Chamaeleon I star-forming region. For this young brown dwarf, OTS 44, we have constructed a spectral energy distribution (SED) from 0.8-8 um by combining the measurements at 3.6, 4.5, 5.8, and 8.0 um from IRAC with ground-based photometry at I, J, H, and K. The resulting SED for OTS 44 exhibits significant excess emission longward of 3 um relative to the SED expected from the photosphere of the brown dwarf. We have successfully modeled the source of this excess emission in terms of an irradiated viscous accretion disk with M'<=10^-10 M_sun/year. With a spectral type of M9.5 and a mass of ~15 M_Jup, OTS 44 is now the coolest and least massive brown dwarf observed to have a circumstellar disk. These measurements demonstrate that disks exist around brown dwarfs even down to the deuterium burning mass limit and the approximate upper mass limit of extrasolar planetary companions.

We present observations of Halpha emission line profiles taken at Magellan Observatory for a sample of 39 young low-mass stars and brown dwarfs in the Taurus and Chamaeleon I star forming regions. We have identified 11 new substellar accretors, more than tripling the number of known brown dwarfs with measurable accretion activity. These include the lowest-mass objects yet seen with accretion, with masses down to ~0.015 Msun. Using models of Halpha emission produced in magnetospheric accretion flows, the most widely applicable primary calibrator now available, we determine the first estimates of mass accretion rates for objects at such extremely low masses. For the six objects with masses < 0.03 Msun, we find accretion rates of ~ 5x10^{-12} Msun/yr, among the smallest yet measured. These new results continue the trend of decreasing mass accretion rate with decreasing (sub)stellar mass that we have noted previously for samples of more massive objects; the overall correlation is Mdot \~ M^{2.1}, and now extends over a mass range of over two orders of magnitude. Finally, the absence of a discontinuity in the distribution of accretion rates with mass tends to suggest that stars and brown dwarfs share similar formation histories.

We announce the release of a catalog of physical models of irradiated accretion disks around young stars based on the modelling techniques by D'Alessio et al. The WWW catalog includes ~ 3000 disk models for different central stars, disk sizes, inclinations, dust contents and mass accretion rates. For any of them, radial profiles of disk physical parameters and synthetic spectral energy distributions can be browsed and downloaded to compare with observations. It can be accessed at http://www-cfa.harvard.edu/youngstars/dalessio/ (US), http://www.astrosmo.unam.mx/~dalessio/ (Mexico), and at http://www.laeff.esa.es/models/dalessio/ (Spain).

We estimate the levels of turbulence in the envelopes of Class 0 and I protostars using a model based on measurements of the peak separation of double-peaked asymmetric line profiles. We use observations of 20 protostars of both Class 0 and I taken in the HCO+ (J=3-2) line that show the classic double-peaked profile. We find that some Class 0 sources show high levels of turbulence whilst others demonstrate much lower levels. In Class I protostars we find predominantly low levels of turbulence. The observations are consistent with a scenario in which Class 0 protostars form in a variety of environments and subsequently evolve into Class I protostars. The data do not appear to be consistent with a recently proposed scenario in which Class 0 protostars can only form in extreme environments.

We present a model of a dusty disk with an inner hole which accounts for the Spitzer Space Telescope Infrared Spectrograph observations of the low-mass pre-main sequence star CoKu Tau/4. We have modeled the mid-IR spectrum (between 8 and 25 mic) as arising from the inner wall of a disk. Our model disk has an evacuated inner zone of radius ~ 10 AU, with a dusty inner ``wall'', of half-height ~ 2 AU, that is illuminated at normal incidence by the central star. The radiative equilibrium temperature decreases from the inner disk edge outward through the optically-thick disk; this temperature gradient is responsible for the emission of the silicate bands at 10 and 20 mic. The observed spectrum is consistent with being produced by Fe-Mg amorphous glassy olivine and/or pyroxene, with no evidence of a crystalline component. The mid-infrared spectrum of CoKu Tau/4 is reminiscent of that of the much older star TW Hya, where it has been suggested that the significant clearing of its inner disk is due to planet formation. However, no inner disk remains in CoKu Tau/4, consistent with the star being a weak-emission (non-accreting) T Tauri star. The relative youth of CoKu Tau/4 (~ 1 Myr) may indicate much more rapid planet formation than typically assumed.

We present results of a large scale, multi-epoch optical survey of the Ori OB1 association, carried out with the QuEST camera at the Venezuela National Astronomical Observatory. We identify for the first time the widely spread low-mass, young population in the Orion OB1a and OB1b sub-associations. Candidate members were picked up by their variability in the V-band and position in color-magnitude diagrams. We obtained spectra to confirm membership. In a region spanning ~ 68 deg^2 we found 197 new young stars; of these, 56 are located in the Ori OB1a subassociation and 142 in Ori OB1b. Comparison with the spatial extent of molecular gas and extinction maps indicates that the subassociation Ori 1b is concentrated within a ring-like structure of radius ~2 deg (~15 pc at 440 pc), centered roughly on the star epsilon Ori in the Orion belt. The ring is apparent in 13CO and corresponds to a region with an extinction Av>=1. The stars exhibiting strong Ha emission, an indicator of active accretion, are found along this ring, while the center is populated with weak Ha emitting stars. In contrast, Ori OB1a is located in a region devoid of gas and dust. We identify a grouping of stars within a ~3 deg^2 area located in 1a, roughly clustered around the B2 star 25 Ori. The Herbig Ae/Be star V346 Ori is also associated with this grouping, which could be an older analog of sigma Ori. Using using several sets of evolutionary tracks we find an age of 7 - 10 Myr for Ori 1a and of ~4 - 6 Myr for Ori OB1b, consistent with previous estimates from OB stars. Indicators such as the equivalent width of Ha and near-IR excesses show that while a substantial fraction of accreting disks remain at ages ~5 Myr, inner disks are essentially dissipated by 10 Myr.

We analyze multi-band photometry of a subsample of low mass stars in the associations Ori OB1a and 1b discovered during the CIDA Orion Variability Survey, which have ages of 7 - 10 Myr and 3 - 5 Myr, respectively. We obtained UBVRcIc photometry at Mt. Hopkins for 6 Classical T Tauri stars (CTTS) and 26 Weak T Tauri stars (WTTS) in Ori OB1a, and for 21 CTTS and 2 WTTS in Ori OB1b. We also obtained L band photometry for 14 CTTS at Mt. Hopkins, and 10um and 18um photometry with OSCIR at Gemini for 6 CTTS; of these, all 6 were detected at 10um while only one was detected at 18um. We estimate mass accretion rates from the excess luminosity at U, and find that they are consistent with determinations for a number of other associations, with or without high mass star formation. The observed decrease of mass accretion rate with age is qualitatively consistent with predictions of viscous evolution of accretion disks. We find an overall decrease of disk emission from Taurus to Ori OB1b to Ori OB1a. This decrease implies that significant grain growth and settling towards the midplane has taken place in the inner disks of Ori OB1. We compare the SED of the star detected at both 10um and 18um with disk models for similar stellar and accretion parameters. We find that the low <= 18 um fluxes of this Ori OB1b star cannot be due to the smaller disk radius expected from viscous evolution in the presence of the FUV radiation fields from the OB stars in the association. Instead, we find that the disk of this star is essentially a flat disk, with little if any flaring, indicating a a significant degree of dust settling towards the midplane, as expected from dust evolution in protoplanetary disks.

We have carried out a study of the early type stars in nearby OB associations spanning an age range of $\sim$ 3 to 16 Myr, with the aim of determining the fraction of stars which belong to the Herbig Ae/Be class. We studied the B, A, and F stars in the nearby ($\le 500$ pc) OB associations Upper Scorpius, Perseus OB2, Lacerta OB1, and Orion OB1, with membership determined from Hipparcos data. We obtained spectra for 440 Hipparcos stars in these associations, from which we determined accurate spectral types, visual extinctions, effective temperatures, luminosities and masses, using Hipparcos photometry. Using colors corrected for reddening, we find that the Herbig Ae/Be stars and the Classical Be stars (CBe) occupy clearly different regions in the JHK diagram. Thus, we use the location on the JHK diagram, as well as the presence of emission lines and of strong 12 microns flux relative to the visual to identify the Herbig Ae/Be stars in the associations. We find that the Herbig Ae/Be stars constitute a small fraction of the early type stellar population even in the younger associations. Comparing the data from associations with different ages and assuming that the near-infrared excess in the Herbig Ae/Be stars arises from optically thick dusty inner disks, we determined the evolution of the inner disk frequency with age. We find that the inner disk frequency in the age range 3 - 10 Myr in intermediate mass stars is lower than that in the low mass stars (< 1msun). This indicates that the time-scales for disk evolution are much shorter in the intermediate mass stars, which could be a consequence of more efficient mechanisms of inner disk dispersal (viscous evolution, dust growth and settling toward the midplane).

We present two-dimensional simulations of finite, self-gravitating gaseous sheets. Unlike the case of infinite sheets, such configurations do not constitute equilibrium states but instead are subject to global collapse unless countered by pressure forces or rotation. The initial effect of finite geometry is to promote concentrations of material at the edges of the sheet. If the sheet is not perfectly circular, gravitational focussing results in enhanced concentrations of mass. In the second-most simple geometry, that of an elliptical outer boundary, the general result is collapse to a filamentary structure with the densest concentrations of mass at the ends of the filament. We suggest that these simple calculations have interesting implications for the gravitational evolution of overall molecular cloud structure, envisioning that such clouds might originate as roughly sheetlike sections of gas accumulated as a result of large-scale flows in the local interstellar medium. We show some examples of local clouds with overall filamentary shape and denser concentrations of mass and star clusters near the ends of the overall extended structure, suggestive of our simple ellipse collapse calculations. We suggest that cluster-forming gas is often concentrated as a result of gravity acting on irregular boundaries; this mechanism can result in very rapid infall of gas which may be of importance to the formation of massive stars. This picture suggests that much of the supersonic ``turbulence'' observed in molecular clouds might be gravitationally-generated. Our results may provide impetus for further theoretical explorations of global gravitational effects in molecular clouds and their implications for generating the substructure needed for fragmentation into stars and clusters.

We present evidence of magnetically-mediated disk accretion in Herbig Ae/Be stars. Magnetospheric accretion models of Balmer and sodium profiles calculated with appropriate stellar and rotational parameters are in qualitative agreement with the observed profiles of the Herbig Ae star UX Ori, and yield a mass accretion rate of ~ 10^{-8} Msun/yr. If more recent indications of an extremely large rotation rate for this object are correct, the magnetic field geometry must deviate from that of a standard dipole in order to produce line emission consistent with observed flux levels. Models of the associated accretion shock qualitatively explain the observed distribution of excess fluxes in the Balmer discontinuity for a large ensemble of Herbig Ae/Be stars, and imply typically small mass accretion rates, < 10^{-7} Msun/yr. In order for accretion to proceed onto the star, significant amounts of gas must exist inside the dust destruction radius, which is potentially problematic for recently advocated scenarios of "puffed" inner dust wall geometries. However, our models of the inner gas disk show that for the typical accretion rates we have derived, the gas should be generally optically thin, thus allowing direct stellar irradiation of the inner dust edge of the disk.

The optically-dark globule IC 1396A is revealed using Spitzer images at 3.6, 4.5, 5.8, 8, and 24 microns to be infrared-bright and to contain a set of previously unknown protostars. The mid-infrared colors of the 24 microns detected sources indicate several very young (Class I or 0) protostars and a dozen Class II stars. Three of the new sources (IC 1396A: gamma, delta, and epsilon) emit over 90% of their bolometric luminosities at wavelengths greater than 3 microns, and they are located within ~0.02 pc of the ionization front at the edge of the globule. Many of the sources have spectra that are still rising at 24 microns. The two previously-known young stars LkHa 349 a and c are both detected, with component c harboring a massive disk and component a being bare. Of order 5% of the mass of material in the globule is presently in the form of protostars in the 10^5 to 10^6 yr age range. This high star formation rate was likely triggered by radiation from a nearby O star.

We present two dimensional cylindrically symmetric hydrodynamic simulations and synthetic emission maps of a stellar wind propagating into an infalling, rotating environment. The resulting outflow morphology, collimation and stability observed in these simulations have relevance to the study of young stellar objects, Herbig-Haro jets and molecular outflows. Our code follows hydrogen gas with molecular, atomic and ionic components tracking the associated time dependent molecular chemistry and ionization dynamics with radiative cooling appropriate for a dense molecular gas. We present tests of the code as well as new simulations which indicate the presence of instabilities in the wind-blown bubble's swept-up shell.

We have used the Spitzer Infrared Spectrograph to observe seven members of the TW Hya association, the nearest stellar association whose age ($\sim$ 10 Myr) is similar to the timescales thought to apply to planet formation and disk dissipation. Only two of the seven targets display infrared excess emission, indicating that substantial amounts of dust still exist closer to the stars than is characteristic of debris disks; however, in both objects we confirm an abrupt short-wavelength edge to the excess, as is seen in disks with cleared-out central regions. The mid-infrared excesses in the spectra of Hen 3-600 and TW Hya include crystalline silicate emission features, indicating that the grains have undergone significant thermal processing. We offer a detailed comparison between the spectra of TW Hya and Hen 3-600, and a model that corroborates the spectral shape and our previous understanding of the radial structure of these protoplanetary disks.

We present observations of the star formation region NGC 7129 taken with the Multiband Imaging Photometer for Spitzer (MIPS). A significant population of sources, likely pre-main sequence members of the young stellar cluster, is revealed outside the central photoionization region. Combining with Infrared Array Camera (IRAC) and ground-based near-infrared images, we have obtained colors and spectral energy distributions for some 60 objects. The [3.6]-[4.5] vs. [8]-[24] color-color plane shows sources clustered at several different loci, which roughly correspond to the archetypal evolutionary sequence Class 0, I, II, and III. We obtain preliminary classifications for 36 objects, and find significant numbers of both Class I and II objects. Most of the pre-main sequence candidates are associated with the densest part of the molecular cloud surrounding the photoionization region, indicating active star formation over a broad area outside the central cluster. We discuss three Class II candidates that exhibit evidence of inner disk clearing, which would be some of the youngest known examples of a transition from accretion to optically thin quiescent disks.

We report initial results from IRAC observations of four young stellar clusters. These regions are part of a larger Spitzer survey of 31 young stellar groups and clusters within 1 kpc of the Sun. In each of the four clusters, there are between 39 and 85 objects with colors inconsistent with reddened stellar photospheres. We identify these objects as young stars with significant emission from circumstellar dust. Applying an analysis developed in a companion paper (Allen et al. 2004), we classify these objects as either pre-main sequence stars with disks (class II) or protostellar objects (class I). These show that the sites of recent star formation are distributed over multi-parsec size scales. In two clusters, Cepheus C and S140, we find protostars embedded in filamentary dark clouds seen against diffuse emission in the IRAC bands.

We report spectra obtained with the Spitzer Space Telescope in the wavelength range between 14 microns and 35 microns of 19 nearby main-sequence stars with infrared excesses. The six stars with strong dust emission show no recognizable spectral features, suggesting that the bulk of the emitting particles have diameters larger than 10 microns. If the observed dust results from collisional grinding of larger solids, we infer minimum masses of the parent body population between 0.004 of the Earth's mass and 0.06 of the Earth's mass. We estimate grain production rates of 10 Gg/s around lambda Boo and HR 1570; selective accretion of this matter may help explain their peculiar surface abundances. There appear to be inner truncations in the dust clouds at 48 AU, 11 AU, 52 AU and 54 AU around HR 333, HR 506, HR 1082 and HR 3927, respectively.

We compare the infrared colors predicted by theoretical models of protostellar envelopes and protoplanetary disks with initial observations of young stellar objects made with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope (Werner et al. 2004, Fazio et al. 2004). Disk and envelope models characterized by infall and/or accretion rates found in previous studies can quantitatively account for the range of IRAC colors found in four young embedded clusters: S 140, S 171, NGC 7129, and Cep C. The IRAC color-color diagram ([3.6]-[4.5] vs. [5.8]-[8.0]) can be used to help to distinguish between young stars with only disk emission and protostars with circumstellar envelopes.

We examine the formation of molecular gas behind shocks in atomic gas using a chemical/dynamical model, particular emphasis is given to constraints the chemistry places on the dynamical evolution. The most important result of this study is to stress the importance of shielding the molecular gas from the destructive effects of UV radiation. For shock ram pressures comparable to or exceeding typical local ISM pressures, self-shielding controls the formation time of H2 but CO formation requires shielding of the interstellar radiation field by dust grains. We find that the molecular hydrogen fractional abundance can become significant well before CO forms. The timescale for (CO) molecular cloud formation is not set by H2 formation, but rather by the timescale for accumulating a sufficient column density or extinction, A_V > 0.7. The local ratio of atomic to molecular gas (4:1), coupled with short estimates for cloud lifetimes (3-5 Myr), suggests that the timescales for accumulating molecular clouds from atomic material typically must be no longer than about 12-20 Myr. Based on the shielding requirement, this implies that the typical product of pre-shock density and velocity must be n*v > 20 cm^-3 km s^-1. Based on these results we find that flow-driven formation of molecular clouds in the local interstellar medium can occur sufficiently rapidly to account for observations. We also provide detailed predictions of atomic and molecular emission and absorption that track molecular cloud formation, with a view toward helping to verify cloud formation by shock waves. Finally, we provide an analytic solution for time-dependent H2 formation which may be of use in numerical hydrodynamic calculations.

We present a sequence of I-band images obtained at the Venezuela 1m Schmidt telescope during the outburst of the nebula recently discovered by J.W. McNeil in the Orion L1630 molecular cloud. We derive photometry spanning the pre-outburst state and the brightening itself, a unique record including 14 epochs and spanning a time scale of ~5 years. We constrain the beginning of the outburst at some time between Oct. 28 and Nov. 15, 2003. The light curve of the object at the vertex of the nebula, the likely exciting source of the outburst, reveals that it has brightened ~5 magnitudes in about 4 months. The time scale for the nebula to develop is consistent with the light travel time, indicating that we are observing light from the central source scattered by the ambient cloud into the line of sight. We also show recent FLWO optical spectroscopy of the exciting source and of the nearby HH 22. The spectrum of the source is highly reddened; in contrast, the spectrum of HH 22 shows a shock spectrum superimposed on a continuum, most likely due to reflected light from the exciting source reaching the HH object through a much less reddened path. The blue portion of this spectrum is consistent with an early B spectral type, similar to the early outburst spectrum of the FU Ori variable V1057 Cyg; we estimate a luminosity of L ~219 Lsun. The eruptive behavior of the McNeil nebula source, its spectroscopic characteristics and luminosity, suggest we may be witnessing an FU Ori event on its way to maximum. Further monitoring of this object will decide whether it qualifies as a member of this rare class of objects.

We present near-infrared spectra of the excess continuum emission from the innermost regions of classical T Tauri disks. In almost all cases, the shape of the excess is consistent with that of a single-temperature blackbody with T ~ 1400 K, similar to the expected dust sublimation temperature for typical dust compositions. The amount of excess flux roughly correlates with the accretion luminosity in objects with similar stellar properties. We compare our observations with the predictions of simple disk models having an inner rim located at the dust sublimation radius, including irradiation heating of the dust from both the stellar and accretion luminosities. The models yield inner rim radii in the range 0.07-0.54 AU, increasing with higher stellar and accretion luminosities. Using typical parameters which fit our observed sample, we predict a rim radius ~ 0.2 AU for the T Tauri star DG Tau, which agrees with recent Keck near-infrared interferometric measurements. For large mass accretion rates, the inner rim lies beyond the corotation radius at (or within) which magnetospheric accretion flows are launched, which implies that pure gaseous disks must extend inside the dust rim. Thus, for a significant fraction of young stars, dust cannot exist in the innermost disk, calling into question theories in which solid particles are ejected by a wind originating at the magnetospheric radius.

We discuss the importance of accretion in calculating disk models for young stellar objects. In particular, we show that a disk inner rim, irradiated by both the star and the accretion shocks at the stellar surface, can naturally explain recent observations of DG Tau with the Keck interferometer. We present models for two objects, with mass accretion rates differing by almost two orders of magnitude, to illustrate the effects of accretion on the overall disk structure and emission.

Protostellar core formation is probably much more dynamic, and magnetic fields are probably much less important, than has been previously assumed in the standard model of low-mass star formation. This revised picture has important consequences: it is easier to understand the observed rapidity of star formation in molecular clouds; cores are more likely to have structures favoring high infall rates at early times, helping to explain the differences between Class 0 and Class I protostars; and core structure and asymmetry will strongly favor post-collapse fragmentation into binary and multiple stellar systems.

The K dwarfs in the Pleiades fall nearly one half magnitude below a main sequence isochrone when plotted in a color-magnitude diagram utilizing V magnitude as the luminosity index and B-V as the color index. This peculiarity has been known for forty years but has gone unexplained and mostly ignored. When compared to Praesepe members, the Pleiades K dwarfs again are subluminous (or blue) in a color-magnitude diagram using B-V as the color index. However, using V-I as the color index, stars in the two clusters are coincident to M_V ~ 10; using V-K as the color index, Pleiades late K and M stars fall above the main sequence locus defined by Praesepe members. We believe that the anomalous spectral energy distributions for the Pleiades K dwarfs, as compared to older clusters, are a consequence of rapid stellar rotation and may be primarily due to spottedness. If so, the required areal filling factor for the cool component has to be very large (=> 50%). Weak-lined T Tauri stars have similar color anomalies, and we suspect this is a common feature of all very young K dwarfs (sp. type > K3). The peculiar spectral energy distribution needs to be considered in deriving accurate pre-main sequence isochrone-fitting ages for clusters like the Pleiades since the age derived will depend on the temperature index used.

New subarcsecond FUV observations of T Tau with HST/STIS show spatially resolved structures in the 2"x2" area around the star. The structures show in multiline emission of fluorescent H_2 pumped by Lyman alpha. One emission structure follows the cavity walls observed around T Tau N in scattered light in the optical. A temperature greater or equal to 1000K is required to have enough population in the H_2 to produce the observed fluorescent lines; in the cool environment of the T Tau system, shock heating is required to achieve this temperature at distances of a few tens of AU. Fluorescent H_2 along the cavity wall represents the best evidence to date for the action of low-density, wide-opening-angle outflows driving cavities into the molecular medium at scales smaller than 100 AU. A southern region of emission consists of two arcs, with shape and orientation similar to the arcs of H_2 2.12 microns and forbidden line emission crossing the outflow associated with the embedded system T Tau S. This region is located near the centroid of forbidden line emission at the blueshifted lobe of the N-S outflow.

Briceno et al. recently used optical imaging, data from the Two-Micron All-Sky Survey (2MASS), and follow-up spectroscopy to search for young low-mass stars and brown dwarfs in 8 square degrees of the Taurus star-forming region. By the end of that study, there remained candidate members of Taurus that lacked the spectroscopic observations needed to measure spectral types and determine membership. In this work, we have obtained spectroscopy of the 22 candidates that have A_V<=8, from which we find six new Taurus members with spectral types of M2.75 through M9. The new M9 source has the second latest spectral type of the known members of Taurus (~0.02 M_sun). Its spectrum contains extremely strong emission in H_alpha (W~950 A) as well as emission in He I 6678 A and the Ca II IR triplet. This is the least massive object known to exhibit emission in He I and Ca II, which together with the strong H_alpha are suggestive of intense accretion.

We present a study of accretion in a sample of 45 young, low mass objects in a variety of star forming regions and young associations, about half of which are likely substellar. Based primarily on the presence of broad, asymmetric Halpha emission, we have identified 13 objects (~30% of our sample) which are strong candidates for ongoing accretion. At least 3 of these are substellar. We do not detect significant continuum veiling in most of the accretors with late spectral types (M5-M7). Accretion shock models show that lack of measurable veiling allows us to place an upper limit to the mass accretion rates of <~ 10^{-10} Msun/yr. Using magnetospheric accretion models with appropriate (sub)stellar parameters, we can successfully explain the accretor Halpha emission line profiles, and derive quantitative estimates of accretion rates in the range 10^{-12} < Mdot < 10^{-9} Msun/yr. There is a clear trend of decreasing accretion rate with stellar mass, with mean accretion rates declining by 3-4 orders of magnitude over ~ 1 - 0.05 Msun.

We model the Class I source L1551 IRS 5, adopting a flattened infalling envelope surrounding a binary disk system and a circumbinary disk. With our composite model, we calculate self-consistently the spectral energy distribution of each component of the L1551 IRS 5 system, using additional constraints from recent observations by ISO, the water ice feature from observations with SpeX, the SCUBA extended spatial brightness distribution at sub-mm wavelengths, and the VLA spatial intensity distributions at 7 mm of the binary disks. We analyze the sensitivity of our results to the various parameters involved. Our results show that a flattened envelope collapse model is required to explain simultaneously the large scale fluxes and the water ice and silicate features. On the other hand, we find that the circumstellar disks are optically thick in the millimeter range and are inclined so that their outer parts hide the emission along the line of sight from their inner parts. We also find that these disks have lower mass accretion rates than the infall rate of the envelope.

Palla & Stahler have recently argued that star formation in Taurus and other nearby molecular clouds extends over a period of at least 10 Myr, implying quasi-static cloud evolution and star formation. Their conclusions contradict other recent results indicating that molecular clouds are transient objects and star formation proceeds rapidly. The Palla & Stahler picture implies that most molecular clouds should have extremely low rates of star formation, and that in such inactive stages the stellar initial mass function should be strongly skewed toward producing stars with masses $\gtrsim 1 \msun$; neither prediction is supported by observations. I show that the Palla & Stahler conclusions for Taurus depend almost entirely on a small number of stars with masses $\gtrsim 1 \msun$; the lower-mass stars show no evidence for such an extended period of star formation. I further show that most of the stars apparently older than 10 Myr in the direction of Taurus are probably foreground non-members. I also present birthline calculations which support the idea that the ages of the stars with masses $\gtrsim 1 \msun$ have been systematically overestimated because ``birthline'' age corrections have been underestimated; such birthlines would eliminate the need to postulate skewed initial mass functions. The simplest and most robust explanation of current observations characterizing the vast majority of young stars in molecular clouds is that cloud and star formation is rapid and dynamic.

The evolution of a stellar disk under the influence of viscous evolution, photoevaporation from the central source, and photoevaporation by external stars is studied. We take the typical parameters of TTSs and the Trapezium Cluster conditions. The photoionizing flux from the central source is assumed to arise both from the quiescent star and accretion shocks at the base of stellar magnetospheric columns, along which material from the disk accretes. The accretion flux is calculated self-consistently from the accretion mass loss rate. We find that the disk cannot be entirely removed using only viscous evolution and photoionization from the disk-star accretion shock. However, when FUV photoevaporation by external massive stars is included the disk is removed in 10^6 -10^7yr; and when EUV photoevaporation by external massive stars is included the disk is removed in 10^5 - 10^6yr. An intriguing feature of photoevaporation by the central star is the formation of a gap in the disk at late stages of the disk evolution. As the gap starts forming, viscous spreading and photoevaporation work in resonance. There is no gap formation for disks nearby external massive stars because the outer annuli are quickly removed by the dominant EUV flux. On the other hand, at larger, more typical distances (d>>0.03pc) from the external stars the flux is FUV dominated. As a consequence, the disk is efficiently evaporated at two different locations; forming a gap during the last stages of the disk evolution.

We present the results of a deep, optical/IR wide field imaging survey of selected fields in the nearby (d ~ 140 pc) Taurus star-forming region. We report the discovery of 9 new members with spectral types M5.75-M9.5. We derive an Initial Mass Function encompassing 54% of the known members in Taurus. Comparison with dense regions like the Trapezium Cluster in Orion shows that Taurus has produced x2 less brown dwarfs. We suggest that the lower frequency of brown dwarfs in Taurus may result from the low-density star-forming environment, leading to larger minimum Jeans masses.

The remarkably filamentary spatial distribution of young stars in the Taurus molecular cloud has significant implications for understanding low-mass star formation in relatively quiescent conditions. The large scale and regular spacing of the filaments suggests that small-scale turbulence is of limited importance, which could be consistent with driving on large scales by flows which produced the cloud. The small spatial dispersion of stars from gaseous filaments indicates that the low-mass stars are generally born with small velocity dispersions relative to their natal gas, of order the sound speed or less. The spatial distribution of the stars exhibits a mean separation of about 0.25 pc, comparable to the estimated Jeans length in the densest gaseous filaments, and is consistent with roughly uniform density along the filaments. The efficiency of star formation in filaments is much higher than elsewhere, with an associated higher frequency of protostars and accreting T Tauri stars. The protostellar cores generally are aligned with the filaments, suggesting that they are produced by gravitational fragmentation, resulting in initially quasi-prolate cores. Given the absence of massive stars which could strongly dominate cloud dynamics, Taurus provides important tests of theories of dispersed low-mass star formation and numerical simulations of molecular cloud structure and evolution.

We present recent results from our ongoing large scale variability survey of the Orion OB1 Association. In an area of ~ 25 square degrees we have unveiled new populations of low-mass young stars over a range of environments, and ages from 1-2 Myr in Ori OB 1b to roughly 10 Myr in Ori OB 1a. There is a lack of strong Ha emission and near-IR excesses in the young stars of Ori OB 1a, suggesting that accretion stops and disks dissipate for most solar type stars in a few Myr, probably caused by of the onset of planet formation. The absence of dust and gas in Ori OB 1a is consistent with the picture of star formation as a rapid process, in which molecular clouds dissipate in just a few Myr after the first stars are born.

We present two-dimensional, cylindrically symmetric simulations of hydrodynamic and magnetohydrodynamic (MHD) wide-angle winds interacting with a collapsing environment. These simulations have direct relevance to young stellar objects (YSOs). The results may also be of use in the study of collimated outflows from proto-planetary and planetary nebulae. We study a range of wind configurations consistent with asymptotic MHD wind collimation. The degree of collimation is parameterized by the ratio of the wind density at the pole to that of the equator. We find that a toroidal magnetic field can have a significant influence on the resulting outflow, giving rise to a very dense, jet-like flow in the post-shock region. The properties of the flow in this region are similar to the asymptotic state of a collimated MHD wind. We conclude that wide-angle MHD winds are quite likely capable of driving molecular outflows. Due to difficulty in treating MHD winds ab-initio in simulations we choose magnetic field strengths in the wind consistent slow magnetic rotators. While MHD launched winds will be in the fast rotator regime we discuss how our results, which rely on toroidal pinch effects, will hold for stronger field strengths.

We have developed a self-consistent model of the disk around the nearby 10 Myr old star TW Hya which matches the observed spectral energy distribution and 7mm images of the disk. The model requires a significant dust size evolution and a partially-evacuated inner disk region, as predicted by theories of planet formation. The outer disk, which extends to at least 140 AU in radius, is very optically thick at IR wavelengths and quite massive ~0.06 Msun for the relatively advanced age of this T Tauri star. This implies long viscous and dust evolution timescales, although dust must have grown to sizes of order ~1cm to explain the sub-mm and mm spectral slopes. In contrast, the negligible near-infrared excess emission of this system requires that the disk be optically thin inside ~4 AU.This inner region cannot be completely evacuated; we need ~0.5 lunar mass of ~1 micron particles remaining to produce the observed 10 micron silicate emission. Our model requires a distinct transition in disk properties at ~4 AU, separating the inner and outer disk. The inner edge of the optically-thick outer disk must be heated almost frontally by the star to account for the excess flux at mid-IR wavelengths. We speculate that this truncation of the outer disk may be the signpost of a developing gap due to the effects of a growing protoplanet; the gap is still presumably evolving because material still resides in it, as indicated by the silicate emission, the molecular hydrogen emission, and by the continued accretion onto the central star (albeit at a much lower rate than typical of younger T Tauri stars). TW Hya thus may become the Rosetta stone for our understanding of the evolution and dissipation of protoplanetary disks.

Observations of both star-forming regions and young, gas-free stellar associations indicate that most nearby molecular clouds form stars only over a short time span before dispersal; large-scale flows in the diffuse interstellar medium have the potential for forming clouds sufficiently rapidly, and for producing stellar populations with ages much less than the lateral crossing times of their host molecular clouds. We identify four important factors for understanding rapid star formation and short cloud lifetimes. First, much of the accumulation and dispersal of clouds near the solar circle might occur in the atomic phase; only the high-density portion of a cloud's lifecycle is spent in the molecular phase, thus helping to limit molecular cloud ``lifetimes''. Second, once a cloud achieves a high enough column density to form $\h2$ and CO, gravitational forces become larger than typical interstellar pressure forces; thus star formation can follow rapidly upon molecular gas formation and turbulent dissipation in limited areas of each cloud complex. Third, typical magnetic fields are not strong enough to prevent rapid cloud formation and gravitational collapse. Fourth, rapid dispersal of gas by newly-formed stars, and reduction of shielding by a small expansion of the cloud after the first events of star formation, might limit the length of the star formation epoch and the lifetime of a cloud in its molecular state. This picture emphasizes the importance of large-scale boundary conditions for understanding molecular cloud formation, and implies that star formation is a highly dynamic, rather than quasi-static, process.

We present an analysis of the accretion properties of the members of the TW Hydrae association. The Halpha emission line profile provides explicit evidence of an active accretion flow in 3 out of 23 observed members of the association; we use radiative transfer models to derive approximate accretion rates for these 3 objects. The resulting accretion rates for these 10 Myr-old stars (<10^-9 Msun/yr) are well over an order of magnitude lower than typical values in 1 Myr-old T Tauri stars. The small fraction of TW Hydrae association objects still accreting (13%, compared with ~70% in 1 Myr-old regions like Taurus), along with the very small accretion rates, points to significant disk evolution over 10 Myr, providing an important constraint on the timescales for planet formation.

We have conducted deep mid-infrared imaging of a relatively nearby sample of candidate Vega-like stars using the OSCIR instrument on the CTIO 4-meter and Keck II 10-meter telescopes. Our discovery of a spatially-resolved disk around HR 4796A has already been reported (Jayawardhana et al. 1998). Here we present imaging observations of the other members of the sample, including the discovery that only the primary in the HD 35187 binary system appears to harbor a substantial circumstellar disk and the possible detection of extended disk emission around 49 Ceti. We derive global properties of the dust disks, place constraints on their sizes, and discuss several interesting cases in detail. Although our targets are believed to be main sequence stars, we note that several have large infrared excesses compared to prototype Vega-like systems, and may therefore be somewhat younger. The disk size constraints we derive, in many cases, imply emission from relatively large ($\gtrsim$ 10$\mu$m) particles at mid-infrared wavelengths.

We present detailed models of irradiated T Tauri disks including dust grain growth with power-law size distributions. The models assume complete mixing between dust and gas and solve for the vertical disk structure self-consistentlyincluding the heating effects of stellar irradiation as well as local viscous heating. For a given total dust mass, grain growth is found to decrease the vertical height of the surface where the optical depth to the stellar radiation becomes unit and thus the local irradiation heating, while increasing the disk emission at mm and sub-mm wavelengths. The resulting disk models are less geometrically thick than our previous models assuming interstellar medium dust, and agree better with observed spectral energy distributions and images of edge-on disks, like HK Tau/c and HH 30. The implications of models with grain growth for determining disk masses from long-wavelength emission are considered.

We have carried out radiative transfer calculations of infalling, dusty envelopes surrounding embedded protostars to understand the observed properties of the recently identified ``Class 0'' sources. To match the far-infrared peaks in the spectral energy distributions of objects such as the prototype Class 0 source VLA 1623, pure collapse models require mass infall rates $\sim10^{-4}\msun$yr$^{-1}$. The radial intensity distributions predicted by such infall models are inconsistent with observations of VLA 1623 at sub-mm wavelengths, in agreement with the results of Andre et al. (1993) who found a density profile of $\rho \propto r^{-1/2}$ rather than the expected $\rho \propto r^{-3/2}$ gradient. To resolve this conflict, while still invoking infall to produce the outflow source at the center of VLA 1623, we suggest that the observed sub-mm intensity distribution is the sum of two components: an inner infall zone, plus an outer, more nearly constant-density region. This explanation of the observations requires that roughly half the total mass observed within 2000 AU radius of the source lies in a region external to the infall zone. The column densities for this external region are comparable to those found in the larger Oph A cloud within which VLA 1623 is embedded. The extreme environments of Class 0 sources lead us to suggest an alternative or additional interpretation of these objects: rather than simply concluding with Andre et al. that Class 0 objects only represent the earliest phases of protostellar collapse, and ultimately evolve into older ``Class I'' protostars, we suggest that many Class 0 sources could be the protostars of very dense regions. (Shortened)

We compare accretion shock models with optical and ultraviolet spectra of pre-main sequence stars to (1) make the first determinations of accretion rates in intermediate mass T Tauri stars from continuum emission, and (2) derive improved estimates of accretion rates and extinctions for continuum T Tauri stars. Our method extends the shock models developed by Calvet and Gullbring to enable comparisons with optical and archival International Ultraviolet Explorer ultraviolet spectra. We find good agreement between the observations and the model predictions, supporting the basic model of magnetospheric accretion shocks as well as previous determinations of accretion rates and interstellar reddening for the low mass T Tauri stars. The accretion rates determined for the intermediate-mass T Tauri stars agree well with values obtained through other methods that use near-infrared hydrogen line strengths (Muzerolle et al. 1998).

The presence of debris disks around $\sim$ 1-Gyr-old main sequence stars suggests that an appreciable amount of dust may persist even in mature planetary systems. Here we report the detection of dust emission from 55 Cancri, a star with one, or possibly two, planetary companions detected through radial velocity measurements. Our observations at 850$\mu$m and 450$\mu$m imply a dust mass of 0.0008-0.005 Earth masses, somewhat higher than that in the the Kuiper Belt of our solar system. The estimated temperature of the dust grains and a simple model fit both indicate a central disk hole of at least 10 AU in radius. Thus, the region where the planets are detected is likely to be significantly depleted of dust. Our results suggest that far-infrared and sub-millimeter observations are powerful tools for probing the outer regions of extrasolar planetary systems.

The Taurus molecular cloud complex is the paradigm for quiescent, low-density, isolated star formation. Yet the age distribution of its stellar population indicates that star formation is a rapid and dynamic process, inconsistent with the old picture of magnetically-controlled protostellar cloud collapse. Instead, Taurus seems to have formed stars in a manner qualitatively consistent with the rapid cloud formation and dispersal inferred for other, higher-density star-forming regions. I suggest that the Taurus clouds were formed rapidly by the collision of atomic gas streams in the interstellar medium, and that star formation ensued immediately afterward.

Using the 8k x 8k CCD Mosaic Camera on the 1m Schmidt telescope in Venezuela, we are conducting a large-scale, deep optical, multiepoch, photometric (BVRIHa) survey over 120 sq.deg. in the Orion OB association, aimed at identifying the low mass stellar populations with ages less than about 10 Myr. We present initial results for a 34 sq.deg. area spanning Orion 1b, 1a and the B Cloud. Using variability as our main selection criterion we derive much cleaner samples than with the usual single-epoch photometric selection, allowing us to attain a much higher efficiency in follow up spectroscopy and resulting in an preliminary list of 74 new low-mass (~ 0.4 Msun) pre-main sequence stars. Though preliminary, this list of new T Tauri stars already suggests that the fraction of accreting young stars in 1a is much lower than in 1b, which would be expected if 1a is indeed older than 1b. We are analyzing in detail the light curves of these new stars, and spectroscopy of further candidates is under way.

The interaction of a stellar or disk wind with a collapsing environment holds promise for explaining a variety of outflow phenomena observed around young stars. In this paper we present the first simulations of these interactions. The focus here is on exploring how ram pressure balance between wind and ambient gas and post-shock cooling affects the shape of the resulting outflows. In our models we explore the role of ram pressure and cooling by holding the wind speed constant and adjusting the ratio of the inflow mass flux to the wind mass flux (Mdot_a/Mdot_w) Assuming non-spherical cloud collapse, we find that relatively strong winds can carve out wide, conical outflow cavities and that relatively weak winds can be strongly collimated into jet-like structures. If the winds become weak enough, they can be cut off entirely by the infalling environment. We identify discrepancies between results from standard snowplow models and those presented here that have important implications for molecular outflows. We also present mass vs. velocity curves for comparison with observations.

We present new 10.8 and 18.2 micron images of HR 4796A, a young A0V star that was recently discovered to have a spectacular, nearly edge-on, circumstellar disk prominent at ~20 microns (Jayawardhana et al. 1998; Koerner et al. 1998). These new images, obtained with OSCIR at Keck II, show that the disk's size at 10 microns is comparable to its size at 18 microns. Therefore, the 18 micron-emitting dust may also emit some, or all, of the 10 micron radiation. Using these multi-wavelength images, we determine a "characteristic" diameter of 2-3 microns for the mid-infrared-emitting dust particles if they are spherical and composed of astronomical silicates. Particles this small are expected to be blown out of the system by radiation pressure in a few hundred years, and therefore these particles are unlikely to be primordial. Dynamical modeling of the disk (Wyatt et al. 2000) indicates that the disk surface density is relatively sharply peaked near 70 AU, which agrees with the mean annular radius deduced by Schneider et al. (1999) from their NICMOS images. We present evidence (~1.8 sigma significance) for a brightness asymmetry that may result from the presence of the hole and the gravitational perturbation of the disk particle orbits by the low-mass stellar companion or a planet. This "pericenter glow," which must still be confirmed, results from a very small (a few AU) shift of the disk's center of symmetry relative to the central star HR 4796A; one side of the inner boundary of the annulus is shifted towards HR 4796A, thereby becoming warmer and more infrared-emitting. The possible detection of pericenter glow implies that the detection of even complex dynamical effects of planets on disks is within reach.

We construct detailed vertical structure models of irradiated accretion disks around T Tauri stars with interstellar medium dust uniformly mixed with gas. The dependence of the structure and emission properties on mass accretion rate, viscosity parameter, and disk radius is explored using these models. The theoretical spectral energy distributions (SEDs) and images for all inclinations are compared with observations of the entire population of Classical T Tauri stars (CTTS) and Class I objects in Taurus. In particular, we find that the median near-infrared fluxes can be explained within the errors with the most recent values for the median accretion rates for CTTS. We further show that the majority of the Class I sources in Taurus cannot be Class II sources viewed edge-on because they are too luminous and their colors would be consistent with disks seen only in a narrow range of inclinations. Our models appear to be too geometrically thick at large radii, as suggested by: (a) larger far-infrared disk emission than in the typical SEDs of T Tauri stars; (b) wider dark dust lanes in the model images than in the images of HH30 and HK Tau/c; and (c) larger predicted number of stars extincted by edge-on disks than consistent with current surveys. The large thickness of the model is a consequence of the assumption that dust and gas are well-mixed, suggesting that some degree of dust settling may be required to explain the observations.

We suggest that molecular clouds can be formed on short time scales by compressions from large scale streams in the interstellar medium (ISM). In particular, we argue that the Taurus-Auriga complex, with filaments of 10-20 pc $\times$ 2-5 pc, most have been formed by H I flows in $\lesssim 3$Myr, explaining the absence of post-T Tauri stars in the region with ages $\gtrsim 3$ Myr. Observations in the 21 cm line of the H I `halos' around the Taurus molecular gas show many features (broad asymmetric profiles, velocity shifts of H I relative to $^{12}$CO) predicted by our MHD numerical simulations, in which large-scale H I streams collide to produce dense filamentary structures. This rapid evolution is possible because the H I flows producing and disrupting the cloud have much higher velocities (5-10 kms) than present in the molecular gas resulting from the colliding flows. The simulations suggest that such flows can occur from the global ISM turbulence without requiring a single triggering event such as a SN explosion.

A group of young, active stars in the vicinity of TW Hydrae has recently been identified as a possible physical association with a common origin. Given its proximity ($\sim$50 pc), age ($\sim$10 Myr) and abundance of binary systems, the TW Hya Association is ideally suited to studies of diversity and evolution of circumstellar disks. Here we present mid-infrared observations of 15 candidate members of the group, 11 of which have no previous flux measurements at wavelengths longer than 2$\mu$m. We report the discovery of a possible 10$\mu$m excess in CD -33$^{\circ}$7795, which may be due to a circumstellar disk or a faint, as yet undetected binary companion. Of the other stars, only TW Hya, HD 98800, Hen 3-600A, and HR 4796A -- all of which were detected by IRAS -- show excess thermal emission. Our 10$\mu$m flux measurements for the remaining members of the Association are consistent with photospheric emission, allowing us to rule out dusty inner disks. In light of these findings, we discuss the origin and age of the TW Hya Association as well as implications for disk evolution timescales.

We report Hubble Space Telescope near-infrared NICMOS observations of a remarkable low-luminosity Class I (protostellar) source in the Taurus Molecular Cloud. IRAS 04325+2402 exhibits a complex bipolar scattered light nebula. The central continuum source is resolved and may be multiple, or may be crossed by a small dust lane. Complex arcs seen in scattered light surround the central source; the physical nature of these structures is not clear, but they may reflect perturbations from multiple stellar sources or from time-dependent mass ejection. A second, resolved continuum source is found at a projected distance of approximately 1150 AU from the central region, near the edge of a nebular lobe probably produced by outflow. The images indicate that this second source is another low-luminosity young stellar object, seen nearly edge-on through a dusty disk and envelope system with disk diameter $\sim 60$ AU. We suggest that the scattered light ``streaks'' associated with this second source are limb-brightened outflow cavities in the dusty envelope, possibly perturbed by interaction with the outflow lobes of the main source. The nature of the companion is uncertain, since it is observed mostly in scattered light, but is most probably a very low mass star or brown dwarf, with a minimum luminosity of $\sim 10^{-2} \lsun$. Our results show that protostellar sources may have multiple centers of infall and non-aligned disks and outflows, even on relatively small scales.

We present high-resolution mid-infrared observations of the nearby late-type young binary system Hen 3-600. The binary, at a distance of $\sim$ 50 pc, could be a member of the TW Hydrae Association, the nearest known group of young stars, with an age of a few million years. Our images make it possible for the first time to determine which star in the pair, separated by 1.4'', harbors the mid-infrared excess detected by IRAS. In the near-infrared, where the radiation is primarily photospheric, Hen 3-600A (M3) and Hen 3-600B (M3.5) have a flux ratio of 1.6. At 4.8$\mu$m, 10.8$\mu$m, and 18.2$\mu$m, the primary becomes increasingly dominant over the secondary, suggesting that most of the circumstellar dust in the system resides around Hen 3-600A. Comparison of the spectral energy distribution (SED) of Hen 3-600A to the median SED of classical T Tauri stars suggests that its disk may be truncated by the secondary and provides tentative evidence for a central disk hole. The distribution of dust in the Hen 3-600 system may provide important clues to the formation and evolution of protoplanetary disks in close binaries.

We review the present knowledge of disk accretion in young low mass stars, and in particular, the mass accretion rate and its evolution with time. The methods used to obtain mass accretion rates from ultraviolet excesses and emission lines are described, and the current best estimates of mass accretion rate for Classical T Tauri stars and for objects still surrounded by infalling envelopes are given. We argue that the low mass accretion rates of the latter objects require episodes of high mass accretion rate to build the bulk of the star. Similarity solutions for viscous disk evolution suggest that the inner disk mass accretion rates can be self-consistently understood in terms of the disk mass and size if the viscosity parameter alpha is approximately 1e-2. Close companion stars may accelerate the disk accretion process, resulting in accretion onto the central star in less or equal than 1Myr; this may help explain the number of very young stars which are not currently surrounded by accretion disks (the weak emission T Tauri stars).

We report the codiscovery of the spatially-resolved dust disk of the Vega-like star HR 4796A. Images of the thermal dust emission at $\lambda = 18 \mu$m show an elongated structure approximately 200 AU in diameter surrounding the central A0V star. The position angle of the disk, $30^{\circ} \pm 10^{\circ}$, is consistent to the position angle of the M companion star, $225^{\circ}$, suggesting that the disk-binary system is being seen nearly along its orbital plane. The surface brightness distribution of the disk is consistent with the presence of an inner disk hole of approximately 50 AU radius, as was originally suggested by Jura et al. on the basis of the infrared spectrum. HR 4796 is a unique system among the Vega-like or $\beta$ Pictoris stars in that the M star companion (a weak-emission T Tauri star) shows that the system is relatively young, $\sim 8 \pm 3$ Myr. The inner disk hole may provide evidence for coagulation of dust into larger bodies on a timescale similar to that suggested for planet formation in the solar system.

We estimate here the age of one of the prototypes of the Beta Pic-like stars, Fomalhaut, based on the properties of its common proper motion companion Gl879. By combining constraints derived from the lithium abundance, rotational velocity, HR diagram position, and coronal activity we conclude that the age for Gl879, and hence the age for Fomalhaut, is 200$\pm$100 Myr. This age estimate agrees quite well with the completely independent age estimate derived directly from isochrone--fitting to Fomalhaut's position in an HR diagram, and thus confirms that circumstellar dust disks can persist in A stars for several hundred Myr.

We present solutions for the accretion disks and boundary layers in pre-main-sequence stars undergoing FU Orionis outbursts. These solutions differ from earlier disk solutions in that they include a self-consistent treatment of the boundary layer region. In a previous paper (Popham 1996), we showed that these stars should stop accreting angular momentum once they spin up to modest rotation rates. Here we show that for reasonable values of $\alpha$, these low angular momentum accretion rate solutions fit the spectra and line profiles observed in FU Orionis objects better than solutions with high rates of angular momentum accretion. We find solutions which fit the observations of FU Orionis and V1057 Cygni. These solutions have mass accretion rates of 2 and $1 \times 10^{-4} \msyr$, stellar masses of 0.7 and $0.5 \msun$, and stellar radii of 5.75 and $5.03 \rsun$, respectively. They also have modest stellar rotation rates $8-9 \xss$, comparable to the observed rotation rates of T Tauri stars, and angular momentum accretion rates of zero. This supports our earlier suggestion that FU Orionis outbursts may regulate the rotation rates of T Tauri stars.