Far-infrared observations of young clusters embedded in the R Coronae Austrinae and RHO Ophiuchi dark clouds

Multicolor far infrared maps in two nearby dark clouds, R Coronae Austrinae and rho Ophiuchi, were made in order to investigate the individual contribution of low mass stars to the energetics and dynamics of the surrounding gas and dust. Emission from cool dust associated with five low mass stars in Cr A and four in rho Oph was detected; their far infrared luminosities range from 2 far infrared luminosities L. up to 40 far infrared luminosities. When an estimate of the bolometric luminosity was possible, it was found that typically more than 50% of the star's energy was radiated longward of 20 micrometers. meaningful limits to the far infrared luminosities of an additional eleven association members in Cr A and two in rho Oph were also obtained. The dust optical depth surrounding the star R Cr A appears to be asymmetric and may control the dynamics of the surrounding molecular gas. The implications of the results for the cloud energetics and star formation efficiency in these two clouds are discussed

Low-Mass Star Formation and the Initial Mass Function in the Rho Ophiuchi Cloud Core

We have obtained moderate-resolution (R=800-1200) K-band spectra for ~100 stars within and surrounding the cloud core of rho Oph. We have measured spectral types and continuum veilings and have combined this information with results from new deep imaging. The IMF peaks at about 0.4 M_sun and slowly declines to the hydrogen burning limit with a slope of ~-0.5 in logarithmic units (Salpeter is +1.35). Our lower limits on the numbers of substellar objects demonstrate that the IMF probably does not fall more steeply below the hydrogen burning limit, at least down to ~0.02 M_sun. We then make the first comparison of mass functions of stars and pre-stellar clumps (Motte, Andre, & Neri) measured in the same region. The similar behavior of the two mass functions in rho Oph supports the suggestion of Motte et al. and Testi & Sargent that the stellar mass function in young clusters is a direct product of the process of cloud fragmentation. After considering the effect of extinction on the SED classifications of the sample, we find that ~17% of the rho Oph stars are Class I, implying ~0.1 Myr for the lifetime of this stage. In spectra separated by two years, we observe simultaneous variability in the Br gamma emission and K-band continuum veiling for two stars, where the hydrogen emission is brighter in the more heavily veiled data. This behavior indicates that the disk may contribute significantly to continuous K-band emission, in contrast to the proposal that the infalling envelope always dominates. Our detection of strong 2 micron veiling (r_K=1-4) in several Class II and III stars, which should have disks but little envelope material, further supports this proposition.Comment: 35 pages, 14 figures, accepted to Ap

A Mid-Infrared Imaging Survey of Embedded Young Stellar Objects in the Rho Ophiuchi Cloud Core

Results of a comprehensive, new, ground-based mid-infrared imaging survey of the young stellar population of the Rho Ophiuchi cloud are presented. Data were acquired at the Palomar 5-m and at the Keck 10-m telescopes with the MIRLIN and LWS instruments, at 0.25 arcsec and 0.25 arcsec resolutions, respectively. Of 172 survey objects, 85 were detected. Among the 22 multiple systems observed, 15 were resolved and their individual component fluxes determined. A plot of the frequency distribution of the detected objects with SED spectral slope shows that YSOs spend ~400,000 yr in the Flat Spectrum phase, clearing out their remnant infall envelopes. Mid-infrared variability is found among a significant fraction of the surveyed objects, and is found to occur for all SED classes with optically thick disks. Large-amplitude near-infrared variability, also found for all SED classes with optically thick disks, seems to occur with somewhat higher frequency at the earlier evolutionary stages. Although a general trend of mid-infrared excess and NIR veiling exists proceeding through SED classes, with Class I objects generally exhibiting K-veilings > 1, Flat Spectrum objects with K-veilings > 0.58, and Class III objects with K-veilings =0, Class II objects exhibit the widest range of K-band veiling values, 0-4.5. However, the highly variable value of veiling that a single source can exhibit in any of the SED classes in which active disk accretion can take place is striking, and is direct observational evidence for highly time-variable accretion activity in disks. Finally, by comparing mid-infrared vs. near-infrared excesses in a subsample with well-determined effective temperatures and extinction values, disk clearing mechanisms are explored. The results are consistent with disk clearing proceeding from the inside-out.Comment: 18 pages + 5 tables + 7 figure

Optical Spectroscopy of the Surface Population of the rho Ophiuchi Molecular Cloud: The First Wave of Star Formation

We present the results of optical spectroscopy of 139 stars obtained with the Hydra multi-object spectrograph. The objects extend over a 1.3 square degree area surrounding the main cloud of the rho Oph complex. The objects were selected from narrowband images to have H alpha in emission. Using the presence of strong H alpha emission, lithium absorption, location in the Hertzsprung-Russell diagram, or previously reported x-ray emission, we were able to identify 88 objects as young stars associated with the cloud. Strong H alpha emission was confirmed in 39 objects with line widths consistent with their origin in magnetospheric accretion columns. Two of the strongest emission-line objects are young, x-ray emitting brown dwarf candidates with M8 spectral types. Comparisons of the bolometric luminosities and effective temperatures with theoretical models suggest a medianage for this population of 2.1 Myr which is signifcantly older than the ages derived for objects in the cloud core. It appears that these stars formed contemporaneously with low mass stars in the Upper Scorpius subgroup, likely triggered by massive stars in the Upper-Centaurus subgroup.Comment: 35 pages of postscript which includes seven figures (some of which are multi-panel) and four postscript tables. Astronomical Journal (in press

Proper Motion of H2O Masers in IRAS 20050+2720 MMS1: An AU Scale Jet Associated with An Intermediate-Mass Class 0 Source

We conducted a 4 epoch 3 month VLBA proper motion study of H$_2$O masers toward an intermediate-mass class 0 source IRAS 20050+2720 MMS1 (d=700 pc). From milli-arcsecond (mas) resolution VLBA images, we found two groups of H2O maser spots at the center of the submillimeter core of MMS1. One group consists of more than $\sim 50$ intense maser spots; the other group consisting of several weaker maser spots is located at 18 AU south-west of the intense group. Distribution of the maser spots in the intense group shows an arc-shaped structure which includes the maser spots that showed a clear velocity gradient. The spatial and velocity structures of the maser spots in the arc-shape did not significantly change through the 4 epochs. Furthermore, we found a relative proper motion between the two groups. Their projected separation increased by 1.13+/-0.11 mas over the 4 epochs along a line connecting them. The spatial and velocity structures of the intense group and the relative proper motions strongly suggest that the maser emission is associated with a protostellar jet. Comparing the observed LSR velocities with calculated radial velocities from a simple biconical jet model, we conclude that the most of the maser emission are likely to be associated with an accelerating biconical jet which has large opening angle. The large opening angle of the jet traced by the masers would support the hypothesis that poor jet collimation is an inherent property of luminous (proto)stars.Comment: 14 pages, 10 figures, Fig.3 was downsized significantly. accepted for publication in A&

On the fidelity of the core mass functions derived from dust column density data

Aims: We examine the recoverability and completeness limits of the dense core mass functions (CMFs) derived for a molecular cloud using extinction data and a core identification scheme based on two-dimensional thresholding. Methods: We performed simulations where a population of artificial cores was embedded into the variable background extinction field of the Pipe nebula. We extracted the cores from the simulated extinction maps, constructed the CMFs, and compared them to the input CMFs. The simulations were repeated using a variety of extraction parameters and several core populations with differing input mass functions and differing degrees of crowding. Results: The fidelity of the observed CMF depends on the parameters selected for the core extraction algorithm for our background. More importantly, it depends on how crowded the core population is. We find that the observed CMF recovers the true CMF reliably when the mean separation of cores is larger than their mean diameter (f>1). If this condition holds, the derived CMF is accurate and complete above M > 0.8-1.5 Msun, depending on the parameters used for the core extraction. In the simulations, the best fidelity was achieved with the detection threshold of 1 or 2 times the rms-noise of the extinction data, and with the contour level spacings of 3 times the rms-noise. Choosing larger threshold and wider level spacings increases the limiting mass. The simulations show that when f>1.5, the masses of individual cores are recovered with a typical uncertainty of 25-30 %. When f=1 the uncertainty is ~60 %. In very crowded cases where f<1 the core identification algorithm is unable to recover the masses of the cores adequately. For the cores of the Pipe nebula f~2.0 and therefore the use of the method in that region is justified.Comment: 9 pages, 6 figures, accepted for publication in A&

VLT/NACO adaptive optics imaging of the TY CrA system - A fourth stellar component candidate detected

We report the detection of a possible subsolar mass companion to the triple young system TY CrA using the NACO instrument at the VLT UT4 during its commissioning. Assuming for TY CrA a distance similar to that of the close binary system HD 176386, the photometric spectral type of this fourth stellar component candidate is consistent with an ~M4 star. We discuss the dynamical stability of this possible quadruple system as well as the possible location of dusty particles inside or outside the system.Comment: 4 pages, 2 figures postscrip

The Structure of the DoAr 25 Circumstellar Disk

We present high spatial resolution (< 0.3" = 40$AU) Submillimeter Array observations of the 865 micron continuum emission from the circumstellar disk around the young star DoAr 25. Despite its bright millimeter emission, this source exhibits only a comparatively small infrared excess and low accretion rate, suggesting that the material and structural properties of the inner disk may be in an advanced state of evolution. A simple model of the physical conditions in the disk is derived from the submillimeter visibilities and the complete spectral energy distribution using a Monte Carlo radiative transfer code. For the standard assumption of a homogeneous grain size distribution at all disk radii, the results indicate a shallow surface density profile,$\Sigma \propto r^{-p}$ with p = 0.34, significantly less steep than a steady-state accretion disk (p = 1) or the often adopted minimum mass solar nebula (p = 1.5). Even though the total mass of material is large (M_d = 0.10 M_sun), the densities inferred in the inner disk for such a model may be too low to facilitate any mode of planet formation. However, alternative models with steeper density gradients (p = 1) can explain the observations equally well if substantial grain growth in the planet formation region (r < 40 AU) has occurred. We discuss these data in the context of such models with dust properties that vary with radius and highlight their implications for understanding disk evolution and the early stages of planet formation.Comment: ApJL in pres