Search

We present Atacama Large Millimeter/sub-millimeter Array (ALMA) Cycle-5 observations of HBC 494, as well as calculations of the kinematic and dynamic variables which represent the object's wide-angle bipolar outflows. HBC 494 is a binary FU Orionis type object located in the Orion A molecular cloud. We take advantage of combining the ALMA main array, Atacama Compact Array (ACA), and Total Power (TP) array in order to map HBC 494's outflows and thus, estimate their kinematic parameters with higher accuracy in comparison to prior publications. We use $^{12}$CO, $^{13}$CO, C$^{18}$O and SO observations to describe the object's outflows, envelope, and disc, as well as estimate the mass, momentum, and kinetic energy of the outflows. After correcting for optical opacity near systemic velocities, we estimate a mass of $3.0\times10^{-2}$ M$_{\odot}$ for the southern outflow and $2.8\times10^{-2}$ M$_{\odot}$ for the northern outflow. We report the first detection of a secondary outflow cavity located approximately $15$" north of the central binary system, which could be a remnant of a previous large-scale accretion outburst. Furthermore, we find CO spatial features in HBC 494's outflows corresponding to position angles of $\sim35^{\circ}$ and $\sim145^{\circ}$. This suggests that HBC 494's outflows are most likely a composite of overlapping outflows from two different sources, i.e., HBC 494a and HBC 494b, the two objects in the binary system., Comment: 26 pages, 24 figures, accepted to MNRAS

The internal velocity structure within dense gaseous cores plays a crucial role in providing the initial conditions for star formation in molecular clouds. However, the kinematic properties of dense gas at core scales (~0.01 - 0.1 pc) has not been extensively characterized because of instrument limitations until the unique capabilities of GBT-Argus became available. The ongoing GBT-Argus Large Program, Dynamics in Star-forming Cores (DiSCo) thus aims to investigate the origin and distribution of angular momenta of star-forming cores. DiSCo will survey all starless cores and Class 0 protostellar cores in the Perseus molecular complex down to ~0.01 pc scales with < 0.05 km/s velocity resolution using the dense gas tracer N$_2$H$^+$. Here, we present the first datasets from DiSCo toward the B1 and NGC 1333 regions in Perseus. Our results suggest that a dense core's internal velocity structure has little correlation with other core-scale properties, indicating these gas motions may be originated externally from cloud-scale turbulence. These first datasets also reaffirm the ability of GBT-Argus for studying dense core velocity structure and provided an empirical basis for future studies that address the angular momentum problem with a statistically broad sample., Comment: 17 pages, 12 figures, accepted by MNRAS

We investigate the physical structure and conditions of photodissociation regions (PDRs) and molecular gas within the Pillars of Creation in the Eagle Nebula using SOFIA FEEDBACK observations of the [C II] 158 micron line. These observations are velocity resolved to 0.5 km s$^{-1}$ and are analyzed alongside a collection of complimentary data with similar spatial and spectral resolution: the [O I] 63 micron line, also observed with SOFIA, and rotational lines of CO, HCN, HCO$^{+}$, CS, and N$_2$H$^{+}$. Using the superb spectral resolution of SOFIA, APEX, CARMA, and BIMA, we reveal the relationships between the warm PDR and cool molecular gas layers in context of the Pillars' kinematic structure. We assemble a geometric picture of the Pillars and their surroundings informed by illumination patterns and kinematic relationships and derive physical conditions in the PDRs associated with the Pillars. We estimate an average molecular gas density $n_{{\rm H}_2} \sim 1.3 \times 10^5$ cm$^{-3}$ and an average atomic gas density $n_{\rm H} \sim 1.8 \times 10^4$ cm$^{-3}$ and infer that the ionized, atomic, and molecular phases are in pressure equilibrium if the atomic gas is magnetically supported. We find pillar masses of 103, 78, 103, and 18 solar masses for P1a, P1b, P2, and P3 respectively, and evaporation times of $\sim$1-2 Myr. The dense clumps at the tops of the pillars are currently supported by the magnetic field. Our analysis suggests that ambipolar diffusion is rapid and these clumps are likely to collapse within their photoevaporation timescales., Comment: 42 pages, 16 figures. Accepted for publication in The Astronomical Journal

A space-based far-infrared interferometer could work synergistically with the James Webb Space Telescope (JWST) and the Atacama Large Millimeter Array (ALMA) to revolutionize our understanding of the astrophysical processes leading to the formation of habitable planets and the co-evolution of galaxies and their central supermassive black holes. Key to these advances are measurements of water in its frozen and gaseous states, observations of astronomical objects in the spectral range where most of their light is emitted, and access to critical diagnostic spectral lines, all of which point to the need for a far-infrared observatory in space. The objects of interest - circumstellar disks and distant galaxies - typically appear in the sky at sub-arcsecond scales, which rendered all but a few of them unresolvable with the successful and now-defunct 3.5-m Herschel Space Observatory, the largest far-infrared telescope flown to date. A far-infrared interferometer with maximum baseline length in the tens of meters would match the angular resolution of JWST at 10x longer wavelengths and observe water ice and water-vapor emission, which ALMA can barely do through the Earth's atmosphere. Such a facility was conceived and studied two decades ago. Here we revisit the science case for a space-based far-infrared interferometer in the era of JWST and ALMA and summarize the measurement capabilities that will enable the interferometer to achieve a set of compelling scientific objectives. Common to all the science themes we consider is a need for sub-arcsecond image resolution., Comment: 11 pages, 12 figures. To be published in Proc. SPIE Conf. 12686 "Instruments: Scientific Returns and Conceptual Designs"

Annular structures, or rings and gaps, in disks around pre-main sequence stars have been detected in abundance towards Class II objects ~1,000,000 years in age. These structures are often interpreted as evidence of planet formation, with planet-mass bodies carving rings and gaps in the disk. This implies that planet formation may already be underway in even younger disks in the Class I phase, when the protostar is still embedded in a larger scale dense envelope of gas and dust. While younger disks likely play an important role in the onset of planet formation, only within the past decade have detailed properties of disks in the youngest star-forming phases begun to be observed. Here we present 1.3 mm dust emission observations with 5 au resolution that show four annular substructures in the disk of the young (<500,000 years) protostar IRS 63. IRS 63, a single Class I source located in the nearby Ophiuchus molecular cloud (at a distance of 144 pc), is one of the brightest Class I protostars at (sub)millimeter wavelengths that also has a relatively large disk (>50 au). Multiple annular substructures observed towards disks at young times can act as an early foothold for dust grain growth, a prerequisite of planet formation. Whether planets already exist or not in the disk of IRS 63, it is clear that the planet formation process begins in the young protostellar phases, earlier than predicted by current planet formation theories., Comment: Accepted for publication in Nature, appearing in the 8 Oct 2020 issue, DOI: 10.1038/s41586-020-2779-6

In typical environments of star-forming clouds, converging supersonic turbulence generates shock-compressed regions, and can create strongly-magnetized sheet-like layers. Numerical MHD simulations show that within these post-shock layers, dense filaments and embedded self-gravitating cores form via gathering material along the magnetic field lines. As a result of the preferred-direction mass collection, a velocity gradient perpendicular to the filament major axis is a common feature seen in simulations. We show that this prediction is in good agreement with recent observations from the CARMA Large Area Star Formation Survey (CLASSy), from which we identified several filaments with prominent velocity gradients perpendicular to their major axes. Highlighting a filament from the northwest part of Serpens South, we provide both qualitative and quantitative comparisons between simulation results and observational data. In particular, we show that the dimensionless ratio $C_v \equiv {\Delta v_h}^2/(GM/L)$, where $\Delta v_h$ is half of the observed perpendicular velocity difference across a filament, and $M/L$ is the filament's mass per unit length, can distinguish between filaments formed purely due to turbulent compression and those formed due to gravity-induced accretion. We conclude that the perpendicular velocity gradient observed in the Serpens South northwest filament can be caused by gravity-induced anisotropic accretion of material from a flattened layer. Using synthetic observations of our simulated filaments, we also propose that a density-selection effect may explain observed subfilaments (one filament breaking into two components in velocity space) as reported in Dhabal et al. (2018)., Comment: 11 pages, 7 figures, accepted for publication in MNRAS

We present the first results of high-spectral resolution (0.023 km/s) N$_2$H$^+$ observations of dense gas dynamics at core scales (~0.01 pc) using the recently commissioned Argus instrument on the Green Bank Telescope (GBT). While the fitted linear velocity gradients across the cores measured in our targets nicely agree with the well-known power-law correlation between the specific angular momentum and core size, it is unclear if the observed gradients represent core-scale rotation. In addition, our Argus data reveal detailed and intriguing gas structures in position-velocity (PV) space for all 5 targets studied in this project, which could suggest that the velocity gradients previously observed in many dense cores actually originate from large-scale turbulence or convergent flow compression instead of rigid-body rotation. We also note that there are targets in this study with their star-forming disks nearly perpendicular to the local velocity gradients, which, assuming the velocity gradient represents the direction of rotation, is opposite to what is described by the classical theory of star formation. This provides important insight on the transport of angular momentum within star-forming cores, which is a critical topic on studying protostellar disk formation., Comment: 14 pages, 10 figures, accepted by MNRAS

We use NH3 inversion transitions to trace the dense gas in the NGC 1333 region of the Perseus molecular cloud. NH3(1,1) and NH3(2,2) maps covering an area of 102 square arcminutes at an angular resolution of ~3.7" are produced by combining VLA interferometric observations with GBT single dish maps. The combined maps have a spectral resolution of 0.14 km/s and a sensitivity of 4 mJy/beam. We produce integrated intensity maps, peak intensity maps and dispersion maps of NH3(1,1) and NH3(2,2) and a line-of-sight velocity map of NH3(1,1). These are used to derive the optical depth for the NH3(1,1) main component, the excitation temperature of NH3(1,1), and the rotational temperature, kinetic temperature and column density of NH3 over the mapped area. We compare these observations with the CARMA J=1-0 observations of N2H+ and H13CO+ and conclude that they all trace the same material in these dense star forming regions. From the NH3(1,1) velocity map, we find that a velocity gradient ridge extends in an arc across the entire southern part of NGC 1333. We propose that a large scale turbulent cell is colliding with the cloud, which could result in the formation of a layer of compressed gas. This region along the velocity gradient ridge is dotted with Class 0/I YSOs, that could have formed from local overdensities in the compressed gas leading to gravitational instabilities. The NH3(1,1) velocity dispersion map also has relatively high values along this region, thereby substantiating the shock layer argument.

We present H13CO+ (J=1-0) and HNC (J=1-0) maps of regions in Serpens South, Serpens Main and NGC 1333 containing filaments. We also observe the Serpens regions using H13CN (J=1-0). These dense gas tracer molecular line observations carried out with CARMA have an angular resolution of ~7", a spectral resolution of ~0.16 km/s and a sensitivity of 50-100 mJy/beam. Although the large scale structure compares well with the Herschel dust continuum maps, we resolve finer structure within the filaments identified by Herschel. The H13CO+ emission distribution agrees with the existing CARMA N2H+ (J=1-0) maps; so they trace the same morphology and kinematics of the filaments. The H13CO+ maps additionally reveal that many regions have multiple structures partially overlapping in the line-of-sight. In two regions, the velocity differences are as high as 1.4 m/s. We identify 8 filamentary structures having typical widths of 0.03-0.08 pc in these tracers. At least 50% of the filamentary structures have distinct velocity gradients perpendicular to their major axis with average values in the range 4-10 km/s/pc. These findings are in support of the theoretical models of filament formation by 2-D inflow in the shock layer created by colliding turbulent cells. We also find evidence of velocity gradients along the length of two filamentary structures; the gradients suggest that these filaments are inflowing towards the cloud core., Comment: 30 pages, 16 figures

Gravitational forces are expected to excite spiral density waves in protoplanetary disks, disks of gas and dust orbiting young stars. However, previous observations that showed spiral structure were not able to probe disk midplanes, where most of the mass is concentrated and where planet formation takes place. Using the Atacama Large Millimeter/submillimeter Array we detected a pair of trailing symmetric spiral arms in the protoplanetary disk surrounding the young star Elias 2-27. The arms extend to the disk outer regions and can be traced down to the midplane. These millimeter-wave observations also reveal an emission gap closer to the star than the spiral arms. We argue that the observed spirals trace shocks of spiral density waves in the midplane of this young disk., Comment: This is our own version of the manuscript, the definitive version was published in Science (DOI: 10.1126/science.aaf8296) on September 30, 2016. Posted to the arxiv for non-commercial use

We present a 3 mm spectral line and continuum survey of L1451 in the Perseus Molecular Cloud. These observations are from the CARMA Large Area Star Formation Survey (CLASSy), which also imaged Barnard 1, NGC 1333, Serpens Main and Serpens South. L1451 is the survey region with the lowest level of star formation activity---it contains no confirmed protostars. HCO+, HCN, and N2H+ (J=1-0) are all detected throughout the region, with HCO+ the most spatially widespread, and molecular emission seen toward 90% of the area above N(H_2) column densities of 1.9x10^21 cm^-2. HCO+ has the broadest velocity dispersion, near 0.3 km/s on average, compared to ~0.15 km/s for the other molecules, thus representing a range from supersonic to subsonic gas motions. Our non-binary dendrogram analysis reveals that the dense gas traced by each molecule has similar hierarchical structure, and that gas surrounding the candidate first hydrostatic core (FHSC), L1451-mm, and other previously detected single-dish continuum clumps have similar hierarchical structure; this suggests that different sub-regions of L1451 are fragmenting on the pathway to forming young stars. We determined the three-dimensional morphology of the largest detectable dense gas structures to be relatively ellipsoidal compared to other CLASSy regions, which appeared more flattened at largest scales. A virial analysis shows the most centrally condensed dust structures are likely unstable against collapse. Additionally, we identify a new spherical, centrally condensed N2H+ feature that could be a new FHSC candidate. The overall results suggest L1451 is a young region starting to form its generation of stars within turbulent, hierarchical structures., Comment: Accepted to The Astrophysical Journal (ApJ), 45 pages, 24 figures (some with reduced resolution in this preprint); Project website is at http://carma.astro.umd.edu/classy

We present new results from the Disks@EVLA program for two young stars: CY Tau and DoAr 25. We trace continuum emission arising from their circusmtellar disks from spatially resolved observations, down to tens of AU scales, at {\lambda} = 0.9, 2.8, 8.0, and 9.8 mm for DoAr25 and at {\lambda} = 1.3, 2.8, and 7.1 mm for CY Tau. Additionally, we constrain the amount of emission whose origin is different from thermal dust emission from 5 cm observations. Directly from interferometric data, we find that observations at 7 mm and 1 cm trace emission from a compact disk while millimeter-wave observations trace an extended disk structure. From a physical disk model, where we characterize the disk structure of CY Tau and DoAr 25 at wavelengths shorter than 5 cm, we find that (1) dust continuum emission is optically thin at the observed wavelengths and over the spatial scales studied, (2) a constant value of the dust opacity is not warranted by our observations, and (3) a high-significance radial gradient of the dust opacity spectral index, {\beta}, is consistent with the observed dust emission in both disks, with low-{\beta} in the inner disk and high-{\beta} in the outer disk. Assuming that changes in dust properties arise solely due to changes in the maximum particle size (amax), we constrain radial variations of amax in both disks, from cm-sized particles in the inner disk (R < 40 AU) to millimeter sizes in the outer disk (R > 80 AU). These observational constraints agree with theoretical predictions of the radial-drift barrier, however, fragmentation of dust grains could explain our amax(R) constraints if these disks have lower turbulence and/or if dust can survive high-velocity collisions., Comment: 18 pages, 15 figures, accepted for publication at ApJ

Double-Fourier interferometry is the most viable path to sub-arcsecond spatial resolution for future astronomical instruments that will observe the universe at far-infrared wavelengths. The double transform spatio-spectral interferometry couples pupil plane beam combination with detector arrays to enable imaging spectroscopy of wide fields, that will be key to accomplishing top-level science goals. The wide field of view and the necessity for these instruments to fly above the opaque atmosphere create unique characteristics and requirements compared to instruments on ground-based telescopes. In this paper, we discuss some characteristics of single-baseline spatio-spectral interferometers. We investigate the impact of intensity and optical path difference noise on the interferogram and the spectral signal-to-noise ratio. We apply our findings to the special case of the Balloon Experimental Twin Telescope for Infrared Interferometry (BETTII), a balloon payload that will be a first application of this technique at far-infrared wavelengths on a flying platform.

We present continuum observations at 1.3 and 2.7 mm using the Combined Array for Research in Millimeter-wave Astronomy (CARMA) toward six protoplanetary disks in the Taurus molecular cloud: CI Tau, DL Tau, DO Tau, FT Tau, Haro 6-13, and HL Tau. We constrain physical properties of the disks with Bayesian inference using two disk models; flared power-law disk model and flared accretion disk model. Comparing the physical properties, we find that the more extended disks are less flared and that the dust opacity spectral index (beta) is smaller in the less massive disks. In addition, disks with a steeper mid-plane density gradient have a smaller beta, which suggests that grains grow and radially move. Furthermore, we compare the two disk models quantitatively and find that the accretion disk model provides a better fit overall. We also discuss the possibilities of substructures on three extended protoplanetary disks., Comment: 45 pages, 10 figures, 7 tables, to be published in ApJ

Phase fluctuations introduced by the atmosphere are the main limiting factor in attaining diffraction limited performance in extended interferometric arrays at millimeter and submillimeter wavelengths. We report the results of C-PACS, the Combined Array for Research in Millimeter-Wave Astronomy Paired Antenna Calibration System. We present a systematic study of several hundred test observations taken during the 2009-2010 winter observing season where we utilize CARMA's eight 3.5-m antennas to monitor an atmospheric calibrator while simultaneously acquiring science observations with 6.1-m and 10.4-m antennas on baselines ranging from a few hundred meters to ~2 km. We find that C-PACS is systematically successful at improving coherence on long baselines under a variety of atmospheric conditions. We find that the angular separation between the atmospheric calibrator and target source is the most important consideration, with consistently successful phase correction at CARMA requiring a suitable calibrator located $\lesssim$6$^\circ$ away from the science target. We show that cloud cover does not affect the success of C-PACS. We demonstrate C-PACS in typical use by applying it to the observations of the nearby very luminous infrared galaxy Arp 193 in $^{12}$CO(2-1) at a linear resolution of ~70 pc (0.12" x 0.18"), 3 times better than previously published molecular maps of this galaxy. We resolve the molecular disk rotation kinematics and the molecular gas distribution and measure the gas surface densities and masses on 90 pc scales. We find that molecular gas constitutes $\sim30\%$ of the dynamical mass in the inner 700 pc of this object with a surface density $\sim10^4 M_\odot$ pc$^{-2}$; we compare these properties to those of the starburst region of NGC 253., Comment: 19 pages, 18 figures

We present observations of N2H+(1-0), HCO+(1-0), and HCN(1-0) toward the Serpens Main molecular cloud from the CARMA Large Area Star Formation Survey (CLASSy). We mapped 150 square arcminutes of Serpens Main with an angular resolution of 7 arcsecs. The gas emission is concentrated in two subclusters (the NW and SE subclusters). The SE subcluster has more prominent filamentary structures and more complicated kinematics compared to the NW subcluster. The majority of gas in the two subclusters has subsonic to sonic velocity dispersions. We applied a dendrogram technique with N2H+(1-0) to study the gas structures; the SE subcluster has a higher degree of hierarchy than the NW subcluster. Combining the dendrogram and line fitting analyses reveals two distinct relations: a flat relation between nonthermal velocity dispersion and size, and a positive correlation between variation in velocity centroids and size. The two relations imply a characteristic depth of 0.15 pc for the cloud. Furthermore, we have identified six filaments in the SE subcluster. These filaments have lengths of 0.2 pc and widths of 0.03 pc, which is smaller than a characteristic width of 0.1 pc suggested by Herschel observations. The filaments can be classified into two types based on their properties. The first type, located in the northeast of the SE subcluster, has larger velocity gradients, smaller masses, and nearly critical mass-per-unit-length ratios. The other type, located in the southwest of the SE subcluster, has the opposite properties. Several YSOs are formed along two filaments which have supercritical mass per unit length ratios, while filaments with nearly critical mass-per-unit-length ratios are not associated with YSOs, suggesting that stars are formed on gravitationally unstable filaments., Comment: Accepted to ApJ. 38 pages, 16 figures, 5 tables

Magnetic fields in accretion disks play a dominant role during the star formation process but have hitherto been observationally poorly constrained. Field strengths have been inferred on T Tauri stars themselves and possibly in the innermost part of the accretion disk, but the strength and morphology of the field in the bulk of the disk have not been observed. Unresolved measurements of polarized emission (arising from elongated dust grains aligned perpendicular to the field) imply average fields aligned with the disks. Theoretically, the fields are expected to be largely toroidal, poloidal, or a mixture of the two, which imply different mechanisms for transporting angular momentum in the disks of actively accreting young stars such as HL Tau. Here we report resolved measurements of the polarized 1.25 mm continuum emission from HL Tau's disk. The magnetic field on a scale of 80 AU is coincident with the major axis (~210 AU diameter) of the disk. From this we conclude that the magnetic field inside the disk at this scale cannot be dominated by a vertical component, though a purely toroidal field does not fit the data well either. The unexpected morphology suggests that the magnetic field's role for the accretion of a T Tauri star is more complex than the current theoretical understanding., Comment: Accepted for publication in Nature

We present sub-arcsecond resolution observations of continuum emission associated with the GG Tau quadruple star system at wavelengths of 1.3, 2.8, 7.3, and 50 mm. These data confirm that the GG Tau A binary is encircled by a circumbinary ring at a radius of 235 AU with a FWHM width of ~60 AU. We find no clear evidence for a radial gradient in the spectral shape of the ring, suggesting that the particle size distribution is spatially homogeneous on angular scales of ~0.1". A central point source, likely associated with the primary component (GG Tau Aa), exhibits a composite spectrum from dust and free-free emission. Faint emission at 7.3 mm is observed toward the low-mass star GG Tau Ba, although its origin remains uncertain. Using these measurements of the resolved, multifrequency emission structure of the GG Tau A system, models of the far-infrared to radio spectrum are developed to place constraints on the grain size distribution and dust mass in the circumbinary ring. The non-negligible curvature present in the ring spectrum implies a maximum particle size of 1-10 mm, although we are unable to place strong constraints on the distribution shape. The corresponding dust mass is 30-300 M_earth, at a temperature of 20-30 K. We discuss how this significant concentration of relatively large particles in a narrow ring at a large radius might be produced in a local region of higher gas pressures (i.e., a particle "trap") located near the inner edge of the circumbinary disk., Comment: ApJ, in press

We present dust continuum observations of the protoplanetary disk surrounding the pre-main sequence star AS 209, spanning more than an order of magnitude in wavelength from 0.88 to 9.8 mm. The disk was observed with sub-arcsecond angular resolution (0.2"-0.5") to investigate radial variations in its dust properties. At longer wavelengths, the disk emission structure is notably more compact, providing model-independent evidence for changes in the grain properties across the disk. We find that physical models which reproduce the disk emission require a radial dependence of the dust opacity \kappa_{\nu}. Assuming that the observed wavelength-dependent structure can be attributed to radial variations in the dust opacity spectral index (\beta), we find that \beta(R) increases from \beta<0.5 at \sim20 AU to \beta>1.5 for R>80 AU, inconsistent with a constant value of \beta\ across the disk (at the 10\sigma\ level). Furthermore, if radial variations of \kappa_{\nu} are caused by particle growth, we find that the maximum size of the particle-size distribution (a_{max}) increases from sub-millimeter-sized grains in the outer disk (R>70 AU) to millimeter and centimeter-sized grains in the inner disk regions (R< 70 AU). We compare our observational constraint on a_{max}(R) with predictions from physical models of dust evolution in proto-planetary disks. For the dust composition and particle-size distribution investigated here, our observational constraints on a_{max}(R) are consistent with models where the maximum grain size is limited by radial drift., Comment: 7 pages, 5 figures, Accepted to ApJL

We present results of high-resolution imaging toward HL Tau by the Combined Array for Research in Millimeter-wave Astronomy (CARMA). We have obtained 1.3 and 2.7 mm dust continua with an angular resolution down to 0.13 arc second. Through model fitting to the two wavelength data simultaneously in Bayesian inference using a flared viscous accretion disk model, we estimate the physical properties of HL Tau, such as density distribution, dust opacity spectral index, disk mass, disk size, inclination angle, position angle, and disk thickness. HL Tau has a circumstellar disk mass of 0.13 solar mass, a characteristic radius of 79 AU, an inclination of 40 degree, and a position angle of 136 degree. Although a thin disk model is preferred by our two wavelength data, a thick disk model is needed to explain the high mid- and far-infrared emission of the HL Tau spectral energy distribution. This could imply large dust grains settled down on the mid plane with fine dust grains mixed with gas. The HL Tau disk is likely gravitationally unstable and can be fragmented between 50 and 100 AU of radius. However, we did not detect dust thermal continuum supporting the protoplanet candidate claimed by a previous study using observations of the Very Large Array at 1.3 cm., Comment: 27 pages, 9 figures, 1 table, to be published in ApJ

We present observations in eight wavebands from 1.25-24 microns of four dense cores: L204C-2, L1152, L1155C-2, and L1228. Our goals are to study the YSO population of these cores and to measure the mid-infrared extinction law. With our combined near-infrared and Spitzer photometry, we classify each source in the cores as, among other things, background stars, galaxies, or embedded young stellar objects (YSOs). L1152 contains three YSOs and L1228 has seven, but neither L204C-2 nor L1155C-2 appear to contain any YSOs. We estimate an upper limit of 7x10^-5 to 5x10^-4 solar luminosities for any undiscovered YSOs in our cores. We also compute the line-of-sight extinction law towards each background star. These measurements are averaged spatially, to create chi-squared maps of the changes in the mid-infrared extinction law throughout our cores, and also in different ranges of extinction. From the chi-squared maps we identify two small regions in L1152 and L1228 where the outflows in those cores appear to be destroying the larger dust grains, thus altering the extinction law in those regions. On average, however, our extinction law is relatively flat from 3.6 to 24 microns for all ranges of extinction and in all four cores. From 3.6 to 8 microns this law is consistent with a dust model that includes larger dust grains than the diffuse interstellar medium, which suggests grain growth has occurred in our cores. At 24 microns, our extinction law is 2-4 times higher than predicted by dust models. However, it is similar to other empirical measurements., Comment: Accepted by ApJ. 46 pages, 20 figures, 6 tables in preprint format. For a version with full-resolution figures, see http://peggysue.as.utexas.edu/SIRTF/PAPERS/

We present dust opacity spectral indexes (beta) of the youngest protostellar systems (so-called Class 0 sources), L1448 IRS 2, L1448 IRS 3, and L1157, obtained between 1.3 mm and 2.7 mm continua, using the Combined Array for Research in Millimeter-wave Astronomy (CARMA). The unprecedented compact configuration and image fidelity of CARMA allow a better detection of the dust continuum emission from Class 0 sources, with a less serious missing flux problem normally associated with interferometry. Through visibility-modeling at both 1.3 mm and 2.7 mm simultaneously, as well as image- and visibility-comparison, we show that beta of the three Class 0 sources are around or smaller than 1, indicating that dust grains have already significantly grown at the Class 0 stage. In addition, we find a radial dependence of beta, which implies faster grain growth in the denser central regions and/or dust segregation. Density distributions of the Class 0 sources are also addressed by visibility-modeling., Comment: 16 pages; to be published in ApJ

(Abridged) The c2d Spitzer Legacy project obtained images and photometry with both IRAC and MIPS instruments for five large, nearby molecular clouds. This paper combines information drawn from studies of individual clouds into a combined and updated statistical analysis of star formation rates and efficiencies, numbers and lifetimes for SED classes, and clustering properties. Current star formation efficiencies range from 3% to 6%. Taken together, the five clouds are producing about 260 solar masses of stars per Myr. The star formation surface density is more than an order of magnitude larger than would be predicted from the Kennicutt relation used in extragalactic studies. Measured against the dense gas probed by the maps of dust continuum emission, the efficiencies are much higher, and the current stock of dense cores would be exhausted in 1.8 Myr on average. The derived lifetime for the Class I phase is 0.44 to 0.54 Myr, considerably longer than some estimates. Similarly, the lifetime for the Class 0 SED class, 0.10 to 0.16 Myr, is longer than early estimates. The great majority (90%) of young stars lie within loose clusters with at least 35 members and a stellar density of 1 solar mass per cubic pc. Accretion at the sound speed from an isothermal sphere over the lifetime derived for the Class I phase could build a star of about 0.25 solar masses, given an efficiency of 0.3. Our data confirm and aggravate the "luminosity problem" for protostars. Our results strongly suggest that accretion is time variable, with prolonged periods of very low accretion. Based on a very simple model and this sample of sources, half the mass of a star would be accreted during only 7% of the Class I lifetime, as represented by the eight most luminous objects., Comment: 41 pages, 14 figures, ApJSupp, in press. Sample tables included for tables 8-12. Full tables and higher quality figures available from http://peggysue.as.utexas.edu/SIRTF/

We compute the mid-infrared extinction law from 3.6-24 microns in three molecular clouds: Ophiuchus, Perseus, and Serpens, by combining data from the "Cores to Disks" Spitzer Legacy Science program with deep JHKs imaging. Using a new technique, we are able to calculate the line-of-sight extinction law towards each background star in our fields. With these line-of-sight measurements, we create, for the first time, maps of the chi-squared deviation of the data from two extinction law models. Because our chi-squared maps have the same spatial resolution as our extinction maps, we can directly observe the changing extinction law as a function of the total column density. In the Spitzer IRAC bands, 3.6-8 microns, we see evidence for grain growth. Below $A_{K_s} = 0.5$, our extinction law is well-fit by the Weingartner & Draine (2001) $R_V = 3.1$ diffuse interstellar medium dust model. As the extinction increases, our law gradually flattens, and for $A_{K_s} >= 1$, the data are more consistent with the Weingartner & Draine $R_V = 5.5$ model that uses larger maximum dust grain sizes. At 24 microns, our extinction law is 2-4 times higher than the values predicted by theoretical dust models, but is more consistent with the observational results of Flaherty et al. (2007). Lastly, from our chi-squared maps we identify a region in Perseus where the IRAC extinction law is anomalously high considering its column density. A steeper near-infrared extinction law than the one we have assumed may partially explain the IRAC extinction law in this region., Comment: 38 pages, 19 figures in pre-print format. Accepted for publication in ApJ. A version with full-resolution figures can be found here: http://peggysue.as.utexas.edu/SIRTF/

Using interferometric data from BIMA observations, combined with detailed modeling in Fourier space of the physical structures predicted by models, we constrain the circumstellar envelope parameters for four Class 0 young stellar objects, as well as their embedded circumstellar disks. The envelopes of these objects are still undergoing collapse, and theoretical collapse models can be compared to the observations. Since it has been suggested in a previous study that both the Larson-Penston and Shu similarity solutions underestimate the age of the system, we adopt Tassis & Mouschovias' model of the collapse process, which includes all relevant magnetic fields effects. The results of the model fitting show a good consistency between theory and data; furthermore, no age problem exists since the Tassis & Mouschovias' model is age independent for the first 255 kyr. Although the majority of the continuum dust emission arises from the circumstellar envelopes, these objects have well known outflows, which suggest the presence of circumstellar disks. At the highest resolution, most of the large-scale envelope emission is resolved out by interferometry, but the small-scale residual emission remains, making it difficult to observe only the compact disk component. By modeling the emission of the envelope and subtracting it from the total emission, we constrain the disk masses in our four systems to be comparable to or smaller than the typical disk masses for T Tauri systems., Comment: 23 pages, 4 figures, accepted for publication in the Astrophysical Journal

We present maps of 14.4 deg^2 of the Ophiuchus dark clouds observed by the Spitzer Space Telescope Multiband Imaging Photometer for Spitzer (MIPS). These high quality maps depict both numerous point sources as well as extended dust emission within the star-forming and non-star-forming portions of these clouds. Using PSF-fitting photometry, we detect 5779 sources at 24 um and 81 sources at 70 um at the 10 sigma level of significance. Three hundred twenty-three candidate young stellar objects (YSOs) were identified according to their positions on the MIPS/2MASS K versus K$-$[24] color-magnitude diagrams as compared to 24 um detections in the SWIRE extragalactic survey. We find that more than half of the YSO candidates, and almost all the ones with protostellar Class I spectral energy distributions, are confined to the known cluster and aggregates., Comment: 55 pages; 22 figures; accepted by Astrophysical Journal

By providing sensitive sub-arcsecond images and integral field spectroscopy in the 25 - 400 micron wavelength range, a far-IR interferometer will revolutionize our understanding of planetary system formation, reveal otherwise-undetectable planets through the disk perturbations they induce, and spectroscopically probe the atmospheres of extrasolar giant planets in orbits typical of most of the planets in our solar system. The technical challenges associated with interferometry in the far-IR are greatly relaxed relative to those encountered at shorter wavelengths or when starlight nulling is required. A structurally connected far-IR interferometer with a maximum baseline length of 36 m can resolve the interesting spatial structures in nascent and developed exoplanetary systems and measure exozodiacal emission at a sensitivity level critical to TPF-I mission planning. The Space Infrared Interferometric Telescope was recommended in the Community Plan for Far-IR/Submillimeter Space Astronomy, studied as a Probe-class mission, and estimated to cost 800M dollars. The scientific communities in Europe, Japan, and Canada have also demonstrated a keen interest in far-IR interferometry through mission planning workshops and technology research, suggesting the possibility of an international collaborative effort., Comment: 7 pages, 4 figures, submitted to the Exoplanet Task Force (AAAC), 30 March 2007

We present maps of 7.78 square degrees of the Lupus molecular cloud complex at 24, 70, and $160\:\mu$m. They were made with the Spitzer Space Telescope's Multiband Imaging Photometer for Spitzer (MIPS) instrument as part of the Spitzer Legacy Program, ``From Molecular Cores to Planet-Forming Disks'' (c2d). The maps cover three separate regions in Lupus, denoted I, III, and IV. We discuss the c2d pipeline and how our data processing differs from it. We compare source counts in the three regions with two other data sets and predicted star counts from the Wainscoat model. This comparison shows the contribution from background galaxies in Lupus I. We also create two color magnitude diagrams using the 2MASS and MIPS data. From these results, we can identify background galaxies and distinguish them from probable young stellar objects. The sources in our catalogs are classified based on their spectral energy distribution (SED) from 2MASS and Spitzer wavelengths to create a sample of young stellar object candidates. From 2MASS data, we create extinction maps for each region and note a strong corresponence between the extinction and the $160\:\mu$m emission. The masses we derived in each Lupus cloud from our extinction maps are compared to masses estimated from $^{13}$CO and C$^{18}$O and found to be similar to our extinction masses in some regions, but significantly different in others. Finally, based on our color-magnitude diagrams, we selected 12 of our reddest candidate young stellar objects for individual discussion. Five of the 12 appear to be newly-discovered YSOs., Comment: 15 pages, 17 figures, uses emulateapj.cls. Accepted for publication in ApJ. A version with high-quality figures can be found at http://peggysue.as.utexas.edu/SIRTF/

We present maps of 1.5 square degrees of the Serpens dark cloud at 24, 70, and 160\micron observed with the Spitzer Space Telescope MIPS Camera. More than 2400 compact sources have been extracted at 24um, nearly 100 at 70um, and 4 at 160um. We estimate completeness limits for our 24um survey from Monte Carlo tests with artificial sources inserted into the Spitzer maps. We compare source counts, colors, and magnitudes in the Serpens cloud to two reference data sets, a 0.50 deg^2 set on a low-extinction region near the dark cloud, and a 5.3 deg^2 subset of the SWIRE ELAIS N1 data that was processed through our pipeline. These results show that there is an easily identifiable population of young stellar object candidates in the Serpens Cloud that is not present in either of the reference data sets. We also show a comparison of visual extinction and cool dust emission illustrating a close correlation between the two, and find that the most embedded YSO candidates are located in the areas of highest visual extinction.

We present deep IRAC images that highlight the scattered light emission around many of the youngest protostars, the so-called Class 0 sources, in L1448. By comparison of the data with a Monte Carlo radiative transfer code \citep{whitney2003a}, we demonstrate for the first time that the observed infrared light from these objects is consistent with scattered light from the central protostar. The scattered light escapes out the cavity, carved by molecular outflows, in the circumstellar envelope. In particular, we observe prominent scattered light nebulae associated with the Class 0 sources: L1448-mm, L1448 IRS 2, and 3B, as well as a Class I source: IRS 3A. We use a grid of models with probable protostellar properties to generate model spectral energy distributions (SEDs) and images for bands sensitive to this scattered light: J, H, Ks, and \textit{Spitzer} IRAC bands. By simultaneously fitting SEDs and images of the outflow cavities, we are able to model geometric parameters, i.e. inclination angle and opening angle, and loosely constrain physical parameters. The opening angle may be an important indicator of the evolutionary state of a source. We compare our results for Class 0 sources to similar studies of Class I sources. There may be a transition phase from Class 0 to Class I when a source has an opening angle between 20$^o$ to 30$^o$. It is important to note that while the best fit model image and SED do not fully describe the sources, the fits generally describe the circumstellar structure of Class 0 sources in L1448., Comment: Accepted ApJ, 42 pages, 13 figures

We present IRAC (3.6, 4.5, 5.8, and 8.0 micron) observations of the Chamaeleon II molecular cloud. The observed area covers about 1 square degree defined by $A_V >2$. Analysis of the data in the 2005 c2d catalogs reveals a small number of sources (40) with properties similar to those of young stellaror substellar objects (YSOs). The surface density of these YSO candidates is low, and contamination by background galaxies appears to be substantial, especially for sources classified as Class I or flat SED. We discuss this problem in some detail and conclude that very few of the candidate YSOs in early evolutionary stages are actually in the Cha II cloud. Using a refined set of criteria, we define a smaller, but more reliable, set of 24 YSO candidates., Comment: 19 pages, 10 figures, in press ApJ

Infrared images of the dark cloud core B59 were obtained with the Spitzer Space Telescope as part of the "Cores to Disks" Legacy Science project. Photometry from 3.6-70 microns indicates at least 20 candidate low-mass young stars near the core, more than doubling the previously known population. Out of this group, 13 are located within about 0.1 pc in projection of the molecular gas peak, where a new embedded source is detected. Spectral energy distributions span the range from small excesses above photospheric levels to rising in the mid-infrared. One other embedded object, probably associated with the millimeter source B59-MMS1, with a bolometric luminosity L(bol) roughly 2 L(sun), has extended structure at 3.6 and 4.5 microns, possibly tracing the edges of an outflow cavity. The measured extinction through the central part of the core is A(V) greater than of order 45 mag. The B59 core is producing young stars with a high efficiency.

We present high resolution 0.45" x 0.32" observations from the BIMA array toward the Herbig Ae system MWC 480 in the lambda = 1.4 mm dust continuum. We resolve a circumstellar disk of radius ~170 AU and constrain the disk parameters by comparing the observations to flat disk models. These results show that the typical fit parameters of the disk, such as the mass, Md ~ 0.04-0.18 Mo, and the surface density power law index, p=0.5 or 1, are comparable to those of the lower mass T Tauri stars. The dust in the MWC 480 disk can be modeled as processed dust material (beta ~ 0.8), similar to the Herbig Ae star CQ Tau disk; the fitted disk parameters are also consistent with less-evolved dust (beta ~ 1.2). The possibility of grain growth in the MWC 480 circumstellar disk is supported by the acceptable fits with beta ~ 0.8. The surface density power-law profiles of p=0.5 and p=1 can be easily fit to the MWC 480 disk; however, a surface density power-law profile similiar to the minimum mass solar nebula model p=1.5 is ruled out at an 80% confidence level., Comment: ApJ Accepted. 17 pages

We present observations of 3.86 sq. deg. of the Perseus molecular cloud complex with the Spitzer Space Telescope Infrared Array Camera (IRAC). The maps show strong extended emission arising from shocked H2 in outflows in the region and from polycyclic aromatic hydrocarbon features. More than 120,000 sources are extracted toward the cloud. Based on their IRAC colors and comparison to off-cloud and extragalactic fields, we identify 400 candidate young stellar objects. About two thirds of these are associated with the young clusters IC348 and NGC1333, while the last third is distributed over the remaining cloud. We classify the young stellar objects according to the traditional scheme based on the slope of their spectral energy distributions. Significant differences are found for the numbers of embedded Class I objects relative to more evolved Class II objects in IC348, NGC1333 and the remaining cloud with the embedded Class I and "flat spectrum" YSOs constituting 14%, 36% and 47% of the total number of YSOs identified in each of these regions. The high number of Class I objects in the extended cloud (61% of the Class I objects in the entire cloud) suggests that a significant fraction of the current star formation is occuring outside the two main clusters. Finally we discuss a number of outflows and identify their driving sources, including the known deeply embedded Class 0 sources outside the two major clusters. The Class 0 objects are found to be detected by Spitzer and have very red [3.6]-[4.5] colors but do not show similarly red [5.8]-[8.0] colors. The Class 0 objects are easily identifiable in color-color diagrams plotting these two colors but are problematic to extract automatically due to the extended emission from shocked gas or scattered light in cavities related to the associated outflows., Comment: Accepted for publication in ApJ (49 pages, 18 figures). Abstract abridged. Version with full resolution figures available at http://peggysue.as.utexas.edu/SIRTF/PAPERS/pap50.pub.pdf

The Wide-Field Imaging Interferometry Testbed (WIIT) will provide valuable information for the development of space-based interferometers. This laboratory instrument operates at optical wavelengths and provides the ability to test operational algorithms and techniques for data reduction of interferometric data. Here we present some details of the system design and implementation, discuss the overall performance of the system to date, and present our plans for future development of WIIT. In order to make best use of the interferometric data obtained with this system, it is critical to limit uncertainties within the system and to accurately understand possible sources of error. The WIIT design addresses these criteria through a number of ancillary systems. The use of redundant metrology systems is one of the most important features of WIIT, and provides knowledge of the delay line position to better than 10 nm. A light power detector is used to monitor the brightness of our light sources to ensure that small fluctuations in brightness do not affect overall performance. We have placed temperature sensors on critical components of the instrument, and on the optical table, in order to assess environmental effects on the system. The use of these systems provides us with estimates of the overall system uncertainty, and allows an overall characterization of the results to date. These estimates allow us to proceed forward with WIIT, adding rotation stages for 2-D interferometry. In addition, they suggest possible avenues for system improvement. Funding for WIIT is provided by NASA Headquarters through the ROSS/SARA Program and by the Goddard Space Flight Center through the IR&D Program., Comment: 11 pages, 9 figures

We present two-dimensional hydrodynamical simulations of both jet- and wind-driven models for protostellar outflows in order to make detailed comparisons to the kinematics of observed molecular outflows. Comparing the different simulations with observations, we find that some outflows, e.g., HH 212, show features consistent with the jet-driven model, while others, e.g., VLA 05487, are consistent with the wind-driven model., Comment: 38 pages, 14 figures, accepted for publication in ApJ

We analyze the dynamics of the shell produced when a bow shock from a collimated jet propagates into the surrounding medium. Under interstellar conditions, the shock is radiative, and a ballistic approximation for the shell flow is appropriate beyond the working surface where the jet impacts its surroundings. The solution is then determined by the ambient and jet densities and velocities and by the momentum impulse applied in the working surface. Using estimates for these impulses (confirmed by separate numerical simulations), we obtain solutions for the shell structure, and for the range of velocities in the shell at any point. We provide predictions for the position-velocity and mass-velocity relations expected for plane-of-sky bow-shock shells, and for the bulk shell properties. In a companion paper, we show that these analytic solutions are in excellent agreement with the results of direct numerical simulations. We argue that classical molecular (CO) outflows cannot be purely jet-driven, because the bow-shock shell solutions are much too elongated compared with observations. Finally, we suggest that the ``spur'' structures seen in position-velocity diagrams of observed molecular outflows are the manifestation of internal bow shocks which may be fit with our simple dynamical models., Comment: 20 pages,4 figures, accepted for publications in ApJ

We have modeled the emission from dust in pre-protostellar cores, including a self-consistent calculation of the temperature distribution for each input density distribution. Model density distributions include Bonnor-Ebert spheres and power laws. The Bonnor-Ebert spheres fit the data well for all three cores we have modeled. The dust temperatures decline to very low values ($\Td \sim 7$ K) in the centers of these cores, strongly affecting the dust emission. Compared to earlier models that assume constant dust temperatures, our models indicate higher central densities and smaller regions of relatively constant density. Indeed, for L1544, a power-law density distribution, similar to that of a singular, isothermal sphere, cannot be ruled out. For the three sources modeled herein, there seems to be a sequence of increasing central condensation, from L1512 to L1689B to L1544. The two denser cores, L1689B and L1544, have spectroscopic evidence for contraction, suggesting an evolutionary sequence for pre-protostellar cores., Comment: 22 pages, 9 figures, Ap. J. accepted, uses emulateapj5.sty

We present the results of a 2.7 mm continuum interferometric survey of 24 young stellar objects in 11 fields. The target objects range from deeply embedded Class 0 sources to optical T Tauri sources. This is the first sub-arcsecond survey of the 2.7 mm dust continuum emission from young, embedded stellar systems. The images show a diversity of structure and complexity. The optically visible T Tauri stars (DG Tauri, HL Tauri, GG Tauri,and GM Aurigae) have continuum emission dominated by compact, less than 1", circumstellar disks. The more embedded near-infrared sources (SVS13 and L1551 IRS5) have continuum emission that is extended and compact. The embedded sources (L1448 IRS3, NGC1333 IRAS2, NGC1333 IRAS4, VLA1623, and IRAS 16293-2422) have continuum emission dominated by the extended envelope, typically more than 85%. In fact, in many of the deeply embedded systems it is difficult to uniquely isolate the disk emission component from the envelope extending inward to AU size scales. All of the target embedded objects are in multiple systems with separations on scales of 30" or less. Based on the system separation, we place the objects into three categories: separate envelope (separation > 6500 AU), common envelope (separation 150-3000 AU), and common disk (separation < 100 AU). These three groups can be linked with fragmentation events during the star formation process: separate envelopes from prompt initial fragmentation and the separate collapse of a loosely condensed cloud, common envelopes from fragmentation of a moderately centrally condensed spherical system, and common disk from fragmentation of a high angular momentum circumstellar disk., Comment: 47 Pages, 18 Figures, ApJ accepted

A major goal of modern astrophysics is to understand the processes by which the universe evolved from its initial simplicity, as seen in measurements of the Cosmic Microwave Background, to the universe we see today, with complexity on all scales. While the diffuse background measurements of COBE reveal the importance of the far infrared and submillimeter in early galaxy and star formation, the understanding of the development of complex structure requires high resolution imaging and spectroscopy. We present a concept for a space mission called SPECS, the Submillimeter Probe of the Evolution of Cosmic Structure, which provides these capabilities. SPECS is a cold, spatial and spectral Michelson interferometer with adjustable baselines ranging up to 1 km. It has Hubble sensitivity and angular resolution in the far IR and submillimeter, spectral resolution up to 10,000, and a 14' field of view. SPECS will be able to image thermal dust continuum and infrared cooling and diagnostic line emission over a wide range of redshifts, providing extinction-free astrophysical probes of young galaxies and early cosmic structures, and measures of the luminosity and heavy element formation history of the universe. SPECS would also have the potential to improve vastly our knowledge of protostars, protoplanetary systems, Active Galactic Nuclei and other objects in the local universe. We recommend that a concerted effort be made during the next decade to develop and test certain critical technologies (photon-counting far IR detectors, formation flying spacecraft, cold, lightweight mirrors, and active coolers), so that SPECS can be deployed in about 2015., Comment: 29 pages, 8 figures; see http://www.gsfc.nasa.gov/astro/specs

We present sub-arcsecond resolution imaging of the $\lambda$ = 2.7 mm continuum emission from the young, embedded system L1551 IRS5 using the nine-element, high-resolution configuration of the BIMA array. The observed emission arises from two compact sources separated by $0\farcs$35, coinciding with the two sources seen at $\lambda$ = 2 cm and $\lambda$ = 1.3 cm. When the high resolution data is combined with data from two compact configurations, L1551 IRS5 is argued to consist of a protobinary system separated by $\sim$50 AU with individual circumstellar disks, a circumbinary structure, and a large-scale envelope. The characteristic masses of the components are: 0.024 M$_{\sun}$ for the northern circumstellar disk, 0.009 M$_{\sun}$ for the southern circumstellar disk, 0.04 M$_{\sun}$ for the circumbinary material, and 0.28 M$_{\sun}$ for the envelope., Comment: 11 pages, 2 figures, latex