Source author record

Tyler L. Bourke

Tyler L. Bourke appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

51works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

51 published item(s)

preprint2026arXiv

The Structure and Kinematics of Three Class 0 Protostellar Jets from JWST

We present observations of jets within 2000 au of three deeply embedded protostars using 2.9-27 micron observations with JWST. These observations show the morphologies and kinematics of the collimated jets from three protostars, the low-mass Class 0 protostars B335 and HOPS 153, and the intermediate-mass protostar HOPS 370. These jets are traced by shock-ionized fine-structure line emission observed with the JWST NIRSpec and MIRI IFUs. We find that [Fe II] emission traces the full extent of the inner 1000 to 2000 au of the jets, depending on distance to the protostar, while other ions mostly trace isolated shocked knots. The jets show evidence of wiggling motion in the plane of the sky as well as asymmetries between blue and red-shifted lobes. The widths of the jets increase non-monotonically with distance from the central protostar, with opening angles ranging from 2.1 degrees to < 10.1 degrees for the three protostars in the sample. The jets have total velocities ranging from 147 to 184 km/s after correcting for disk inclination. For B335, an 8-month gap between NIRSpec and MIRI MRS observations enabled measurement of the tangential velocity of a shocked knot; in combination with the radial velocity, this shows that the jet has a different inclination than the outflow cavity. We find multiple knots before and during a recent outburst in B335, although the knots were more frequent during the burst. The asymmetries between blue- and red-shifted lobes strongly suggest complex interactions between the circumstellar disks and magnetic fields.

preprint2022arXiv

B-fields in Star-Forming Region Observations (BISTRO): Magnetic Fields in the Filamentary Structures of Serpens Main

We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~μ$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities.

preprint2022arXiv

Effects of magnetic field orientations in dense cores on gas kinematics in protostellar envelopes

Theoretically, misalignment between the magnetic field and rotational axis in a dense core is considered to be dynamically important in the star formation process, however, extent of this influence remains observationally unclear. For a sample of 32 Class 0 and I protostars in the Perseus Molecular Cloud, we analyzed gas motions using C$^{18}$O data from the SMA MASSES survey and the magnetic field structures using 850 $μ$m polarimetric data from the JCMT BISTRO-1 survey and archive. We do not find any significant correlation between the velocity gradients in the C$^{18}$O emission in the protostellar envelopes at a 1,000 au scale and the misalignment between the outflows and magnetic field orientations in the dense cores at a 4,000 au scale, and there is also no correlation between the velocity gradients and the angular dispersions of the magnetic fields. However, a significant dependence on the misalignment angles emerges after we normalize the rotational motion by the infalling motion, where the ratios increase from $\lesssim1$ to $\gtrsim1$ with increasing misalignment angles. This suggests that the misalignment could prompt angular momentum transportation to the envelope scale but is not a dominant factor in determining the envelope rotation, and other parameters, like mass accretion in protostellar sources, also play an important role. These results remain valid after taking into account projection effects. The comparison between our estimated angular momentum in the protostellar envelopes and the sizes of the known protostellar disks suggests that significant angular momentum is likely lost between radii of $\sim$1,000-100 au in protostellar envelopes.

preprint2022arXiv

The Central 1000 au of a Pre-stellar Core Revealed with ALMA. II. Almost Complete Freeze-out

Pre-stellar cores represent the initial conditions in the process of star and planet formation. Their low temperatures ($<$10 K) allow the formation of thick icy dust mantles, which will be partially preserved in the future protoplanetary disks, ultimately affecting the chemical composition of planetary systems. Previous observations have shown that carbon- and oxygen-bearing species, in particular CO, are heavily depleted in pre-stellar cores due to the efficient molecular freeze-out onto the surface of cold dust grains. However, N-bearing species such as NH$_3$ and, in particular, its deuterated isotopologues, appear to maintain high abundances where CO molecules are mainly in solid phase. Thanks to ALMA, we present here the first clear observational evidence of NH$_2$D freeze-out toward the L1544 pre-stellar core, suggestive of the presence of a"complete-depletion zone" within a $\simeq$1800 au radius, in agreement with astrochemical pre-stellar core model predictions. Our state-of-the-art chemical model coupled with a non-LTE radiative transfer code demonstrates that NH$_2$D becomes mainly incorporated in icy mantles in the central 2000 au and starts freezing-out already at $\simeq$7000 au. Radiative transfer effects within the pre-stellar core cause the NH$_2$D(1$_{11}$-1$_{01}$) emission to appear centrally concentrated, with a flattened distribution within the central $\simeq$3000 au, unlike the 1.3 mm dust continuum emission which shows a clear peak within the central $\simeq$1800 au. This prevented NH$_2$D freeze-out to be detected in previous observations, where the central 1000 au cannot be spatially resolved.

preprint2022arXiv

The SKA as a prebiotic molecule detector

One of the theories for the origin of life proposes that a significant fraction of prebiotic material could have arrived to Earth from outer space between 4.1 and 3.8 billion years ago. This suggests that those prebiotic compounds could have originated in interstellar space, to be later on incorporated to small Solar-system bodies and planetesimals. The recent discovery of prebiotic molecules such as hydroxylamine and ethanolamine in the interstellar medium, strongly supports this hypothesis. However, some species such as sugars, key for the synthesis of ribonucleotides and for metabolic processes, remain to be discovered in space. The unmatched sensitivity of the Square Kilometer Array (SKA) at centimeter wavelengths will be able to detect even more complex and heavier prebiotic molecules than existing instrumentation. In this contribution, we illustrate the potential of the SKA to detect simple sugars with three and four carbon atoms, using a moderate investment of observing time.

preprint2022arXiv

The Twisted Magnetic Field of the Protobinary L483

We present H-band (1.65 $μ$m) and SOFIA HAWC+ 154 $μ$m polarization observations of the low-mass core L483. Our H-band observations reveal a magnetic field that is overwhelmingly in the E-W direction, which is approximately parallel to the bipolar outflow that is observed in scattered IR light and in single-dish $^{12}$CO observations. From our 154 $μ$m data, we infer a $\sim$ 45$^{\circ}$ twist in the magnetic field within the inner 5" (1000 au) of L483. We compare these new observations with published single-dish 350 $μ$m polarimetry and find that the 10,000 au scale H-band data match the smaller scale 350 $μ$m data, indicating that the collapse of L483 is magnetically regulated on these larger scales. We also present high-resolution 1.3 mm ALMA data of L483 which reveals it is a close binary star with a separation of 34 au. The plane of the binary of L483 is observed to be approximately parallel to the twisted field in the inner 1000 au. Comparing this result to the $\sim$ 1000 au protostellar envelope, we find that the envelope is roughly perpendicular to the 1000 au HAWC+ field. Using the data presented, we speculate that L483 initially formed as a wide binary and the companion star migrated to its current position, causing an extreme shift in angular momentum thereby producing the twisted magnetic field morphology observed. More observations are needed to further test this scenario.

preprint2021arXiv

Evolution and Kinematics of Protostellar Envelopes in the Perseus Molecular Cloud

We present a comprehensive analysis on the evolution of envelopes surrounding protostellar systems in the Perseus molecular cloud using data from the MASSES survey. We focus our attention to the C$^{18}$O(2--1) spectral line, and we characterize the shape, size, and orientation of 54 envelopes and measure their fluxes, velocity gradients, and line widths. To look for evolutionary trends, we compare these parameters to the bolometric temperature Tbol, a tracer of protostellar age. We find evidence that the angular difference between the elongation angle of the C$^{18}$O envelope and the outflow axis direction generally becomes increasingly perpendicular with increasing Tbol, suggesting the envelope evolution is directly affected by the outflow evolution. We show that this angular difference changes at Tbol = $53 \pm 20$ K, which includes the conventional delineation between Class 0 and I protostars of 70K. We compare the C$^{18}$O envelopes with larger gaseous structures in other molecular clouds and show that the velocity gradient increases with decreasing radius ($|\mathcal{G}| \sim R^{-0.72 \pm 0.06}$). From the velocity gradients we show that the specific angular momentum follows a power law fit $J/M \propto R^{1.83 \pm 0.05}$ for scales from 1pc down to $\sim$500 au, and we cannot rule out a possible flattening out at radii smaller than $\sim$1000 au.

preprint2020arXiv

JCMT POL-2 and BISTRO Survey observations of magnetic fields in the L1689 molecular cloud

We present 850$μ$m polarization observations of the L1689 molecular cloud, part of the nearby Ophiuchus molecular cloud complex, taken with the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT). We observe three regions of L1689: the clump L1689N which houses the IRAS 16293-2422 protostellar system, the starless clump SMM-16, and the starless core L1689B. We use the Davis-Chandrasekhar-Fermi method to estimate plane-of-sky field strengths of $366\pm 55$ $μ$G in L1689N, $284\pm 34$ $μ$G in SMM-16, and $72\pm 33$ $μ$G in L1689B, for our fiducial value of dust opacity. These values indicate that all three regions are likely to be magnetically trans-critical with sub-Alfvénic turbulence. In all three regions, the inferred mean magnetic field direction is approximately perpendicular to the local filament direction identified in $Herschel$ Space Telescope observations. The core-scale field morphologies for L1689N and L1689B are consistent with the cloud-scale field morphology measured by the $Planck$ Space Observatory, suggesting that material can flow freely from large to small scales for these sources. Based on these magnetic field measurements, we posit that accretion from the cloud onto L1689N and L1689B may be magnetically regulated. However, in SMM-16, the clump-scale field is nearly perpendicular to the field seen on cloud scales by $Planck$, suggesting that it may be unable to efficiently accrete further material from its surroundings.

preprint2020arXiv

Magnetic Field Structure in Spheroidal Star-Forming Clouds. II. Estimating Field Structure from Observed Maps

This paper presents models to estimate the structure of density and magnetic field strength in spheroidal condensations, from maps of their column density and their polarization of magnetically aligned dust grains. The density model is obtained by fitting a column density map with an embedded p = 2 Plummer spheroid of any aspect ratio and inclination. The magnetic properties are based on the density model, on the Davis-Chandrasekhar-Fermi (DCF) model of Alfvénic fluctuations, and on the Spheroid Flux Freezing (SFF) model of mass and flux conservation in Paper I. The field strength model has the resolution of the column density map, which is finer than the resolution of the DCF estimate of field strength. The models are applied to ALMA observations of the envelope of the protostar BHR71 IRS1. Column density fits give the density model, from (2.0 +- 0.4) x 10^5 cm^-3 to (7 +- 1) x 10^7 cm^-3 . The density model predicts the field directions map, which fits the polarization map best within 1100 au, with standard deviation of angle differences 17°. In this region the DCF mean field strength is 0.7 +- 0.2 mG and the envelope mass is supercritical, with ratio of mass to magnetic critical mass 1.5 +- 0.4. The SFF field strength profile scales with the DCF field strength, from 60 x 10μG to 4+-1 mG. The spatial resolution of the SFF field strength estimate is finer than the DCF resolution by a factor ~7, and the peak SFF field strength exceeds the DCF field strength by a factor ~5.

preprint2020arXiv

Magnetism Science with the Square Kilometre Array

The Square Kilometre Array (SKA) will answer fundamental questions about the origin, evolution, properties, and influence of magnetic fields throughout the Universe. Magnetic fields can illuminate and influence phenomena as diverse as star formation, galactic dynamics, fast radio bursts, active galactic nuclei, large-scale structure, and Dark Matter annihilation. Preparations for the SKA are swiftly continuing worldwide, and the community is making tremendous observational progress in the field of cosmic magnetism using data from a powerful international suite of SKA pathfinder and precursor telescopes. In this contribution, we revisit community plans for magnetism research using the SKA, in the light of these recent rapid developments. We focus in particular on the impact that new radio telescope instrumentation is generating, thus advancing our understanding of key SKA magnetism science areas, as well as the new techniques that are required for processing and interpreting the data. We discuss these recent developments in the context of the ultimate scientific goals for the SKA era.

preprint2019arXiv

Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) -- Full Data Release

We present and release the full dataset for the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. This survey used the Submillimeter Array (SMA) to image the 74 known protostars within the Perseus molecular cloud. The SMA was used in two array configurations to capture outflows for scales $>$30$^{\prime\prime}$ ($>$9000 au) and to probe scales down to $\sim$1$^{\prime\prime}$ ($\sim$300 au). The protostars were observed with the 1.3 mm and 850 $μ$m receivers simultaneously to detect continuum at both wavelengths and molecular line emission from CO(2-1), $^{13}$CO(2-1), C$^{18}$O(2-1), N$_2$D$^+$(3-2), CO(3-2), HCO$^+$(4-3), and H$^{13}$CO$^+$(4-3). Some of the observations also used the SMA's recently upgraded correlator, SWARM, whose broader bandwidth allowed for several more spectral lines to be observed (e.g., SO, H$_2$CO, DCO$^+$, DCN, CS, CN). Of the main continuum and spectral tracers observed, 84% of the images and cubes had emission detected. The median C$^{18}$O(2-1) linewidth is $\sim$1.0 km s$^{-1}$, which is slightly higher than those measured with single-dish telescopes at scales of 3000-20000 au. Of the 74 targets, six are suggested to be first hydrostatic core candidates, and we suggest that L1451-mm is the best candidate. We question a previous continuum detection toward L1448 IRS2E. In the SVS13 system, SVS13A certainly appears to be the most evolved source, while SVS13C appears to be hotter and more evolved than SVS13B. The MASSES survey is the largest publicly available interferometric continuum and spectral line protostellar survey to date, and is largely unbiased as it only targets protostars in Perseus. All visibility ($uv$) data and imaged data are publicly available at https://dataverse.harvard.edu/dataverse/full_MASSES/.

preprint2019arXiv

Understanding the origin of the magnetic field morphology in the wide-binary protostellar system BHR 71

We present 1.3 mm ALMA observations of polarized dust emission toward the wide-binary protostellar system BHR 71 IRS1 and IRS2. IRS1 features what appears to be a natal, hourglass-shaped magnetic field. In contrast, IRS2 exhibits a magnetic field that has been affected by its bipolar outflow. Toward IRS2, the polarization is confined mainly to the outflow cavity walls. Along the northern edge of the redshifted outflow cavity of IRS2, the polarized emission is sandwiched between the outflow and a filament of cold, dense gas traced by N$_2$D$^+$, toward which no dust polarization is detected. This suggests that the origin of the enhanced polarization in IRS2 is the irradiation of the outflow cavity walls, which enables the alignment of dust grains with respect to the magnetic field -- but only to a depth of ~300 au, beyond which the dust is cold and unpolarized. However, in order to align grains deep enough in the cavity walls, and to produce the high polarization fraction seen in IRS2, the aligning photons are likely to be in the mid- to far-infrared range, which suggests a degree of grain growth beyond what is typically expected in very young, Class 0 sources. Finally, toward IRS1 we see a narrow, linear feature with a high (10-20%) polarization fraction and a well ordered magnetic field that is not associated with the bipolar outflow cavity. We speculate that this feature may be a magnetized accretion streamer; however, this has yet to be confirmed by kinematic observations of dense-gas tracers.

preprint2016arXiv

A Catalog of Low-Mass Star-Forming Cores Observed with SHARC-II at 350 microns

We present a catalog of low-mass dense cores observed with the SHARC-II instrument at 350 microns. Our observations have an effective angular resolution of 10", approximately 2.5 times higher than observations at the same wavelength obtained with the Herschel Space Observatory, albeit with lower sensitivity, especially to extended emission. The catalog includes 81 maps covering a total of 164 detected sources. For each detected source, we tabulate basic source properties including position, peak intensity, flux density in fixed apertures, and radius. We examine the uncertainties in the pointing model applied to all SHARC-II data and conservatively find that the model corrections are good to within ~3", approximately 1/3 of the SHARC-II beam. We examine the differences between two array scan modes and find that the instrument calibration, beam size, and beam shape are similar between the two modes. We also show that the same flux densities are measured when sources are observed in the two different modes, indicating that there are no systematic effects introduced into our catalog by utilizing two different scan patterns during the course of taking observations. We find a detection rate of 95% for protostellar cores but only 45% for starless cores, and demonstrate the existence of a SHARC-II detection bias against all but the most massive and compact starless cores. Finally, we discuss the improvements in protostellar classification enabled by these 350 micron observations.

preprint2016arXiv

An ALMA Search for Substructure, Fragmentation, and Hidden Protostars in Starless Cores in Chamaeleon I

We present an Atacama Large Millimeter/submillimeter Array (ALMA) 106 GHz (Band 3) continuum survey of the complete population of dense cores in the Chamaeleon I molecular cloud. We detect a total of 24 continuum sources in 19 different target fields. All previously known Class 0 and Class I protostars in Chamaeleon I are detected, whereas all of the 56 starless cores in our sample are undetected. We show that the Spitzer+Herschel census of protostars in Chamaeleon I is complete, with the rate at which protostellar cores have been misclassified as starless cores calculated as <1/56, or < 2%. We use synthetic observations to show that starless cores collapsing following the turbulent fragmentation scenario are detectable by our ALMA observations when their central densities exceed ~10^8 cm^-3, with the exact density dependent on the viewing geometry. Bonnor-Ebert spheres, on the other hand, remain undetected to central densities at least as high as 10^10 cm^-3. Our starless core non-detections are used to infer that either the star formation rate is declining in Chamaeleon I and most of the starless cores are not collapsing, matching the findings of previous studies, or that the evolution of starless cores are more accurately described by models that develop less substructure than predicted by the turbulent fragmentation scenario, such as Bonnor-Ebert spheres. We outline future work necessary to distinguish between these two possibilities.

preprint2016arXiv

Misalignment of Outflow Axes in the Proto-Multiple Systems in Perseus

We investigate the alignment between outflow axes in nine of the youngest binary/multiple systems in the Perseus Molecular Cloud. These systems have typical member spacing larger than 1000 AU. For outflow identification, we use 12CO(2-1) and 12CO(3-2) data from a large survey with the Sub- millimeter Array: Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES). The distribution of outflow orientations in the binary pairs is consistent with random or preferentially anti-aligned distributions, demonstrating that these outflows are misaligned. This result suggests that these systems are possibly formed in environments where the distribution of angular momentum is complex and disordered, and these systems do not come from the same co-rotating structures or from an initial cloud with aligned vectors of angular momentum.

preprint2015arXiv

A Sub-arcsecond Survey Toward Class 0 Protostars in Perseus: Searching for Signatures of Protostellar Disks

We present a CARMA 1.3 mm continuum survey toward 9 Class 0 protostars in the Perseus molecular cloud at $\sim$0.3$^{\prime\prime}$ (70 AU) resolution. This study approximately doubles the number of Class 0 protostars observed with spatial resolutions $<$ 100 AU at millimeter wavelengths, enabling the presence of protostellar disks and proto-binary systems to be probed. We detect flattened structures with radii $>$ 100 AU around 2 sources (L1448 IRS2 and Per-emb-14) and these sources may be strong disk candidates. Marginally-resolved structures with position angles within 30$^{\circ}$ of perpendicular to the outflow are found toward 3 protostars (L1448 IRS3C, IRAS 03282+3035, and L1448C) and are considered disk candidates. Two others (L1448 IRS3B and IRAS 03292+3039) have resolved structure, possibly indicative of massive inner envelopes or disks; L1448 IRS3B also has a companion separated by 0.9$^{\prime\prime}$ ($\sim$210 AU). IC348-MMS does not have well-resolved structure and the candidate first hydrostatic core L1451-MMS is marginally resolved on 1$^{\prime\prime}$ scales. The strong disk candidate sources were followed-up with C$^{18}$O ($J=2\rightarrow1$) observations, detecting velocity gradients consistent with rotation, but it is unclear if the rotation is Keplerian. We compare the observed visibility amplitudes to radiative transfer models, finding that visibility amplitude ratios suggest a compact component (possibly a disk) is necessary for 5 of 9 Class 0 sources; envelopes alone may explain the other 4 systems. We conclude that there is evidence for the formation of large disks in the Class 0 phase with a range of radii and masses dependent upon their initial formation conditions.

preprint2015arXiv

Infrared and Radio observations of a small group of protostellar objects in the molecular core, L1251-C

We present a multi-wavelength observational study of a low-mass star-forming region, L1251-C, with observational results at wavelengths from the near-infrared to the millimeter. Spitzer Space Telescope observations confirmed that IRAS 22343+7501 is a small group of protostellar objects. The extended emission to east-west direction with its intensity peak at the center of L1251A has been detected at 350 and 850 micron with the CSO and JCMT telescopes, tracing dense envelope materials around L1251A. The single-dish data from the KVN and TRAO telescopes show inconsistencies between the intensity peaks of several molecular line emission and that of the continuum emission, suggesting complex distributions of molecular abundances around L1251A. The SMA interferometer data, however, show intensity peaks of CO 2-1 and 13CO 2-1 located at the position of IRS 1, which is both the brightest source in IRAC image and the weakest source in the 1.3 mm dust continuum map. IRS 1 is the strongest candidate for the driving source of the newly detected compact CO 2-1 outflow. Over the whole region (14' by 14') of L125l-C, 3 Class I and 16 Class II sources have been detected, including three YSOs in L1251A. A comparison with the average projected distance among 19 YSOs in L1251-C and that among 3 YSOs in L1251A suggests L1251-C is an example of low-mass cluster formation, where protostellar objects are forming in a small group.

preprint2015arXiv

Mass Assembly of Stellar Systems and Their Evolution with the SMA (MASSES). Multiplicity and the Physical Environment in L1448N

We present continuum and molecular line observations at 230 GHz and 345 GHz from the Sub-millimeter Array (SMA) toward three protostars in the Perseus L1448N region. The data are from the large project "Mass Assembly of Stellar Systems and their Evolution with the SMA" (MASSES). Three dust continuum sources, Source B, Source NW, and Source A, are detected at both frequencies. These sources have corresponding emission peaks in C18O (J=2-1), 13CO (J=2-1), and HCO+ (J=4-3), and have offsets with N2D+ (J=3-2) peaks. High angular resolution data from a complimentary continuum survey with the Karl G. Jansky Very Large Array show that Source B is associated with three 8 mm continuum objects, Source NW with two, and Source A remains single. These results suggest that multiplicity in L1448N exists at different spatial scales from a few thousand AU to < 100 AU. Velocity gradients in each source obtained from two-dimensional fits to the SMA C18O emission are found to be perpendicular to within 20 degrees of the outflow directions as revealed by 12CO (J=2-1). We have observed that Sources B and NW with multiplicity have higher densities than Source A without multiplicity. This suggests that thermal Jeans fragmentation can be relevant in the fragmentation process. However, we have not observed a difference in the ratio between rotational and gravitational energy between sources with and without multiplicity. We also have not observed a trend between non-thermal velocity dispersions and the level of fragmentation. Our study has provided the first direct and comprehensive comparison between multiplicity and core properties in low-mass protostars, although based on small number statistics.

preprint2014arXiv

Dynamical structure of the inner 100 AU of the deeply embedded protostar IRAS 16293-2422

A fundamental question about the early evolution of low-mass protostars is when circumstellar disks may form. High angular resolution observations of molecular transitions in the (sub)millimeter wavelength windows make it possible to investigate the kinematics of the gas around newly-formed stars, for example to identify the presence of rotation and infall. IRAS 16293-2422 was observed with the extended Submillimeter Array (eSMA) resulting in subarcsecond resolution (0.46" x 0.29", i.e. $\sim$ 55 $\times$ 35~AU) images of compact emission from the C$^{17}$O (3-2) and C$^{34}$S (7-6) transitions at 337~GHz (0.89~mm). To recover the more extended emission we have combined the eSMA data with SMA observations of the same molecules. The emission of C$^{17}$O (3-2) and C$^{34}$S (7-6) both show a velocity gradient oriented along a northeast-southwest direction with respect to the continuum marking the location of one of the components of the binary, IRAS16293A. Our combined eSMA and SMA observations show that the velocity field on the 50--400~AU scales is consistent with a rotating structure. It cannot be explained by simple Keplerian rotation around a single point mass but rather needs to take into account the enclosed envelope mass at the radii where the observed lines are excited. We suggest that IRAS 16293-2422 could be among the best candidates to observe a pseudo-disk with future high angular resolution observations.

preprint2014arXiv

The Spitzer Survey of Interstellar Clouds in the Gould Belt. VI. The Auriga-California Molecular Cloud observed with IRAC and MIPS

We present observations of the Auriga-California Molecular Cloud (AMC) at 3.6, 4.5, 5.8, 8.0, 24, 70 and 160 micron observed with the IRAC and MIPS detectors as part of the Spitzer Gould Belt Legacy Survey. The total mapped areas are 2.5 sq-deg with IRAC and 10.47 sq-deg with MIPS. This giant molecular cloud is one of two in the nearby Gould Belt of star-forming regions, the other being the Orion A Molecular Cloud (OMC). We compare source counts, colors and magnitudes in our observed region to a subset of the SWIRE data that was processed through our pipeline. Using color-magnitude and color-color diagrams, we find evidence for a substantial population of 166 young stellar objects (YSOs) in the cloud, many of which were previously unknown. Most of this population is concentrated around the LkHalpha 101 cluster and the filament extending from it. We present a quantitative description of the degree of clustering and discuss the fraction of YSOs in the region with disks relative to an estimate of the diskless YSO population. Although the AMC is similar in mass, size and distance to the OMC, it is forming about 15 - 20 times fewer stars.

preprint2014arXiv

Tomography of Galactic star-forming regions and spiral arms with the Square Kilometer Array

Very Long Baseline Interferometry (VLBI) at radio wavelengths can provide astrometry accurate to 10 micro-arcseconds or better (i.e. better than the target GAIA accuracy) without being limited by dust obscuration. This means that unlike GAIA, VLBI can be applied to star-forming regions independently of their internal and line-of-sight extinction. Low-mass young stellar objects (particularly T Tauri stars) are often non-thermal compact radio emitters, ideal for astrometric VLBI radio continuum experiments. Existing observations for nearby regions (e.g. Taurus, Ophiuchus, or Orion) demonstrate that VLBI astrometry of such active T Tauri stars enables the reconstruction of both the regions' 3D structure (through parallax measurements) and their internal kinematics (through proper motions, combined with radial velocities). The extraordinary sensitivity of the SKA telescope will enable similar "tomographic mappings" to be extended to regions located several kpc from Earth, in particular to nearby spiral arm segments. This will have important implications for Galactic science, galactic dynamics and spiral structure theories.

preprint2013arXiv

Filamentary Accretion Flows in the Embedded Serpens South Protocluster

One puzzle in understanding how stars form in clusters is the source of mass -- is all of the mass in place before the first stars are born, or is there an extended period when the cluster accretes material which can continuously fuel the star formation process? We use a multi-line spectral survey of the southern filament associated with the Serpens South embedded cluster-forming region in order to determine if mass is accreting from the filament onto the cluster, and whether the accretion rate is significant. Our analysis suggests that material is flowing along the filament's long axis at a rate of ~30Msol/Myr (inferred from the N2H+ velocity gradient along the filament), and radially contracting onto the filament at ~130Msol/Myr (inferred from HNC self-absorption). These accretion rates are sufficient to supply mass to the central cluster at a similar rate to the current star formation rate in the cluster. Filamentary accretion flows may therefore be very important in the ongoing evolution of this cluster.

preprint2013arXiv

SMA Observations of Class 0 Protostars: A High-Angular Resolution Survey of Protostellar Binary Systems

We present high angular resolution 1.3 mm and 850 um dust continuum data obtained with the Submillimeter Array toward 33 Class 0 protostars in nearby clouds (distance < 500 pc), which represents so far the largest survey toward protostellar binary/multiple systems. The median angular resolution in the survey is 2.5 arcsec, while the median linear resolution is approximately 600 AU. Compact dust continuum emission is observed from all sources in the sample. Twenty-one sources in the sample show signatures of binarity/multiplicity, with separations ranging from 50 to 5000 AU. The numbers of singles, binaries, triples, and quadruples in the sample are 12, 14, 5, and 2, respectively. The derived multiplicity frequency (MF) and companion star fraction (CSF) for Class 0 protostars are 0.64+/-0.08 and 0.91+/-0.05, respectively, with no correction for completeness. The derived MF and CSF in this survey are approximately two times higher than the values found in the binary surveys toward Class I YSOs, and approximately three (for MF) and four (for CSF) times larger than the values found among MS stars, with a similar range of separations. Furthermore, the observed fraction of high order multiple systems to binary systems in Class 0 protostars (0.50+/-0.09) is also larger than the fractions found in Class I YSOs (0.31+/-0.07) and MS stars (< 0.2). These results suggest that binary properties evolve as protostars evolve, as predicted by numerical simulations. The distribution of separations for Class 0 protostellar binary systems shows a general trend in which companion star fraction increases with decreasing companion separation. We find that 67%+/-8% of the protobinary systems have circumstellar mass ratios below 0.5, implying that unequal-mass systems are preferred in the process of binary star formation. We suggest an empirical sequential fragmentation picture for binary star formation.

preprint2013arXiv

VLA and CARMA Observations of Protostars in the Cepheus Clouds: Sub-arcsecond Proto-Binaries Formed via Disk Fragmentation

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 $λ= 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 $λ= 1.3$ mm with a best resolution of $\sim0\farcs3$ (100 AU), and the CB230 IRS1 observations were taken at $λ= 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.

preprint2012arXiv

A First Look at the Auriga-California Giant Molecular Cloud With Herschel and the CSO: Census of the Young Stellar Objects and the Dense Gas

We have mapped the Auriga/California molecular cloud with the Herschel PACS and SPIRE cameras and the Bolocam 1.1 mm camera on the Caltech Submillimeter Observatory (CSO) with the eventual goal of quantifying the star formation and cloud structure in this Giant Molecular Cloud (GMC) that is comparable in size and mass to the Orion GMC, but which appears to be forming far fewer stars. We have tabulated 60 compact 70/160um sources that are likely pre-main-sequence objects and correlated those with Spitzer and WISE mid-IR sources. At 1.1 mm we find 18 cold, compact sources and discuss their properties. The most important result from this part of our study is that we find a modest number of additional compact young objects beyond those identified at shorter wavelengths with Spitzer. We also describe the dust column density and temperature structure derived from our photometric maps. The column density peaks at a few x 10^22 cm^-2 (N_H2) and is distributed in a clear filamentary structure along which nearly all the pre-main-sequence objects are found. We compare the YSO surface density to the gas column density and find a strong non-linear correlation between them. The dust temperature in the densest parts of the filaments drops to ~10K from values ~ 14--15K in the low density parts of the cloud. We also derive the cumulative mass fraction and probability density function of material in the cloud which we compare with similar data on other star-forming clouds.

preprint2012arXiv

Detection of the simplest sugar, glycolaldehyde, in a solar-type protostar with ALMA

Glycolaldehyde (HCOCH2OH) is the simplest sugar and an important intermediate in the path toward forming more complex biologically relevant molecules. In this paper we present the first detection of 13 transitions of glycolaldehyde around a solar-type young star, through Atacama Large Millimeter Array (ALMA) observations of the Class 0 protostellar binary IRAS 16293-2422 at 220 GHz (6 transitions) and 690 GHz (7 transitions). The glycolaldehyde lines have their origin in warm (200-300 K) gas close to the individual components of the binary. Glycolaldehyde co-exists with its isomer, methyl formate (HCOOCH3), which is a factor 10-15 more abundant toward the two sources. The data also show a tentative detection of ethylene glycol, the reduced alcohol of glycolaldehyde. In the 690 GHz data, the seven transitions predicted to have the highest optical depths based on modeling of the 220 GHz lines all show red-shifted absorption profiles toward one of the components in the binary (IRAS16293B) indicative of infall and emission at the systemic velocity offset from this by about 0.2" (25 AU). We discuss the constraints on the chemical formation of glycolaldehyde and other organic species - in particular, in the context of laboratory experiments of photochemistry of methanol-containing ices. The relative abundances appear to be consistent with UV photochemistry of a CH3OH-CO mixed ice that has undergone mild heating. The order of magnitude increase in line density in these early ALMA data illustrate its huge potential to reveal the full chemical complexity associated with the formation of solar system analogs.

preprint2012arXiv

Revealing The Millimeter Environment of the New FU Orionis Candidate HBC722 with the Submillimeter Array

We present 230 GHz Submillimeter Array continuum and molecular line observations of the newly discovered FUor candidate HBC722. We report the detection of seven 1.3 mm continuum sources in the vicinity of HBC722, none of which correspond to HBC722 itself. We compile infrared and submillimeter continuum photometry of each source from previous studies and conclude that three are Class 0 embedded protostars, one is a Class I embedded protostar, one is a Class I/II transition object, and two are either starless cores or very young, very low luminosity protostars or first hydrostatic cores. We detect a northwest-southeast outflow, consistent with the previous detection of such an outflow in low-resolution, single-dish observations, and note that its axis may be precessing. We show that this outflow is centered on and driven by one of the nearby Class 0 sources rather than HBC722, and find no conclusive evidence that HBC722 itself is driving an outflow. The non-detection of HBC722 in the 1.3 mm continuum observations suggests an upper limit of 0.02 solar masses for the mass of the circumstellar disk. This limit is consistent with typical T Tauri disks and with a disk that provides sufficient mass to power the burst.

preprint2012arXiv

SMA and Spitzer Observations of Bok Glouble CB17: A Candidate First Hydrostatic Core?

We present high angular resolution SMA and Spitzer observations toward the Bok globule CB17. SMA 1.3mm dust continuum images reveal within CB17 two sources with an angular separation of about 21" (about 5250 AU at a distance of 250 pc). The northwestern continuum source, referred to as CB17 IRS, dominates the infrared emission in the Spitzer images, drives a bipolar outflow extending in the northwest-southeast direction, and is classified as a low luminosity Class0/I transition object (L_bol ~ 0.5 L_sun). The southeastern continuum source, referred to as CB17 MMS, has faint dust continuum emission in the SMA 1.3mm observations (about 6 sigma detection; ~3.8 mJy), but is not detected in the deep Spitzer infrared images at wavelengths from 3.6 to 70 micron. Its bolometric luminosity and temperature, estimated from its spectral energy distribution, are less than 0.04 L_sun and 16 K, respectively. The SMA CO(2-1) observations suggest that CB17 MMS may drive a low-velocity molecular outflow (about 2.5 km/s), extending in the east-west direction. Comparisons with prestellar cores and Class0 protostars suggest that CB17 MMS is more evolved than prestellar cores but less evolved than Class0 protostars. The observed characteristics of CB17 MMS are consistent with the theoretical predictions from radiative/magneto hydrodynamical simulations of a first hydrostatic core, but there is also the possibility that CB17 MMS is an extremely low luminosity protostar deeply embedded in an edge-on circumstellar disk. Further observations are needed to study the properties of CB17 MMS and to address more precisely its evolutionary stage.

preprint2011arXiv

Arcsecond resolution images of the chemical structure of the low-mass protostar IRAS 16293-2422

It remains a key challenge to establish the molecular content of different components of low-mass protostars, like their envelopes and disks, and how this depends on the evolutionary stage and/or environment of the young stars. Observations at submillimeter wavelengths provide a direct possibility to study the chemical composition of low-mass protostars through transitions probing temperatures up to a few hundred K in the gas surrounding these sources. This paper presents a large molecular line survey of the deeply embedded protostellar binary IRAS 16293-2422 from the Submillimeter Array (SMA) - including images of individual lines down to approximately 1.5-3" (190-380 AU) resolution. More than 500 individual transitions are identified related to 54 molecular species (including isotopologues) probing temperatures up to about 550 K. Strong chemical differences are found between the two components in the protostellar system with a separation between, in particular, the sulfur- and nitrogen-bearing species and oxygen-bearing complex organics. The action of protostellar outflow on the ambient envelope material is seen in images of CO and SiO and appear to influence a number of other species, including (deuterated) water, HDO. The effects of cold gas-phase chemistry is directly imaged through maps of CO, N2D+ and DCO+, showing enhancements of first DCO+ and subsequently N2D+ in the outer envelope where CO freezes-out on dust grains.

preprint2011arXiv

Detection of a Bipolar Molecular Outflow Driven by a Candidate First Hydrostatic Core

We present new 230 GHz Submillimeter Array observations of the candidate first hydrostatic core Per-Bolo 58. We report the detection of a 1.3 mm continuum source and a bipolar molecular outflow, both centered on the position of the candidate first hydrostatic core. The continuum detection has a total flux density of 26.6 +/- 4.0 mJy, from which we calculate a total (gas and dust) mass of 0.11 +/- 0.05 Msun and a mean number density of 2.0 +/- 1.6 X 10^7 cm-3. There is some evidence for the existence of an unresolved component in the continuum detection, but longer-baseline observations are required in order to confirm the presence of this component and determine whether its origin lies in a circumstellar disk or in the dense inner envelope. The bipolar molecular outflow is observed along a nearly due east-west axis. The outflow is slow (characteristic velocity of 2.9 km/s), shows a jet-like morphology (opening semi-angles ~8 degrees for both lobes), and extends to the edges of the primary beam. We calculate the kinematic and dynamic properties of the outflow in the standard manner and compare them to several other protostars and candidate first hydrostatic cores with similarly low luminosities. We discuss the evidence both in support of and against the possibility that Per-Bolo 58 is a first hydrostatic core, and we outline future work needed to further evaluate the evolutionary status of this object.

preprint2011arXiv

Initial conditions for star formation in clusters: physical and kinematical structure of the starless core Oph A-N6

We present high spatial (<300 AU) and spectral (0.07 km/s) resolution Submillimeter Array observations of the dense starless cluster core Oph A-N6, in the 1 mm dust continuum and the 3-2 line of N2H+ and N2D+. The dust continuum observations reveal a compact source not seen in single-dish observations, of size ~1000 AU and mass 0.005-0.01 M\odot. The combined line and single-dish observations reveal a core of size 3000 \times 1400 AU elongated in a NW-SE direction, with almost no variation in either line width or line center velocity across the map, and very small non-thermal motions. The deuterium fraction has a peak value of ~0.15 and is >0.05 over much of the core. The N2H+ column density profile across the major axis of Oph A-N6 is well represented by an isothermal cylinder, with temperature 20 K, peak density 7.1 \times 10^6 cm^{-3}, and N2H+ abundance 2.7 \times 10^{-10}. The mass of Oph A-N6 is estimated to be 0.29 M\odot, compared to a value of 0.18 M\odot from the isothermal cylinder analysis, and 0.63 M\odot for the critical mass for fragmentation of an isothermal cylinder. Compared to isolated low-mass cores, Oph A-N6 shows similar narrow line widths and small velocity variation, with a deuterium fraction similar to "evolved" dense cores. It is significantly smaller than isolated cores, with larger peak column and volume density. The available evidence suggests Oph A-N6 has formed through the fragmentation of the Oph A filament and is the precursor to a low-mass star. The dust continuum emission suggests it may already have begun to form a star.

preprint2011arXiv

The Spitzer c2d Survey of Nearby Dense Cores XI : Infrared and Submillimeter Observations of CB130

We present new observations of the CB130 region, composed of three separate cores. Using the \textit{Spitzer Space Telescope} we detected a Class 0 and a Class II object in one of these, CB130-1. The observed photometric data from \textit{Spitzer} and ground-based telescopes are used to establish the physical parameters of the Class 0 object. SED fitting with a radiative transfer model shows that the luminosity of the Class 0 object is 0.14 $-$ 0.16 L$_{\odot}$, which is a low luminosity for a protostellar object. In order to constrain the chemical characteristics of the core having the low luminosity object, we compare our molecular line observations to models of lines including abundance variations. We tested both ad hoc step function abundance models and a series of self-consistent chemical evolution models. In the chemical evolution models, we consider a continuous accretion model and an episodic accretion model to explore how variable luminosity affects the chemistry. The step function abundance models can match observed lines reasonably well. The best fitting chemical evolution model requires episodic accretion and the formation of CO$_2$ ice from CO ice during the low luminosity periods. This process removes C from the gas phase, providing a much improved fit to the observed gas-phase molecular lines and the CO$_2$ ice absorption feature. Based on the chemical model result, the low luminosity of CB130-1 is explained better as a quiescent stage between episodic accretion bursts rather than being at the first hydrostatic core stage.

preprint2011arXiv

The Spitzer Survey of Interstellar Clouds in the Gould Belt. III. A Multi-Wavelength View of Corona Australis

We present Spitzer Space Telescope IRAC and MIPS observations of a 0.85 deg^2 field including the Corona Australis (CrA) star-forming region. At a distance of 130 pc, CrA is one of the closest regions known to be actively forming stars, particularly within its embedded association, the Coronet. Using the Spitzer data, we identify 51 young stellar objects (YSOs) in CrA which include sources in the well-studied Coronet cluster as well as distributed throughout the molecular cloud. Twelve of the YSOs discussed are new candidates, one of which is located in the Coronet. Known YSOs retrieved from the literature are also added to the list, and a total of 116 candidate YSOs in CrA are compiled. Based on these YSO candidates, the star formation rate is computed to be 12 M_o Myr^-1, similar to that of the Lupus clouds. A clustering analysis was also performed, finding that the main cluster core, consisting of 68 members, is elongated (having an aspect ratio of 2.36), with a circular radius of 0.59 pc and mean surface density of 150 pc^-2. In addition, we analyze outflows and jets in CrA by means of new CO and H_2 data. We present 1.3 mm interferometric continuum observations made with the Submillimeter Array (SMA) covering R CrA, IRS 5, IRS 7, and IRAS 18595-3712 (IRAS 32). We also present multi-epoch H_2 maps and detect jets and outflows, study their proper motions, and identify exciting sources. The Spitzer and ISAAC/VLT observations of IRAS 32 show a bipolar precessing jet, which drives a CO (2-1) outflow detected in the SMA observations. There is also clear evidence for a parsec-scale precessing outflow, E-W oriented, and originating in the SMA 2 region, likely driven by SMA 2 or IRS 7A.

preprint2011arXiv

The Spitzer Survey of Interstellar Clouds in the Gould Belt. IV. Lupus V and VI Observed with IRAC and MIPS

We present Gould's Belt (GB) Spitzer IRAC and MIPS observations of the Lupus V and VI clouds and discuss them in combination with near-infrared (2MASS) data. Our observations complement those obtained for other Lupus clouds within the frame of the Spitzer "Core to Disk" (c2d) Legacy Survey. We found 43 Young Stellar Object (YSO) candidates in Lupus V and 45 in Lupus VI, including 2 transition disks, using the standard c2d/GB selection method. None of these sources was classified as a pre-main sequence star from previous optical, near-IR and X-ray surveys. A large majority of these YSO candidates appear to be surrounded by thin disks (Class III; ~79% in Lupus V and ~87% in Lupus VI). These Class III abundances differ significantly from those observed for the other Lupus clouds and c2d/GB surveyed star-forming regions, where objects with optically thick disks (Class II) dominate the young population. We investigate various scenarios that can explain this discrepancy. In particular, we show that disk photo-evaporation due to nearby OB stars is not responsible for the high fraction of Class III objects. The gas surface densities measured for Lupus V and VI lies below the star-formation threshold (AV {\eqsim}8.6 mag), while this is not the case for other Lupus clouds. Thus, few Myrs older age for the YSOs in Lupus V and VI with respect to other Lupus clouds is the most likely explanation of the high fraction of Class III objects in these clouds, while a higher characteristic stellar mass might be a contributing factor. Better constraints on the age and binary fraction of

preprint2010arXiv

Extreme Active Molecular Jets in L1448C

The protostellar jet driven by L1448C was observed in the SiO J=8-7 and CO J=3-2 lines and 350 GHz dust continuum at ~1" resolution with the Submillimeter Array (SMA). A narrow jet from the northern source L1448C(N) was observed in the SiO and the high-velocity CO. The jet consists of a chain of emission knots with an inter-knot spacing of ~2" (500 AU) and a semi-periodic velocity variation. The innermost pair of knots, which are significant in the SiO map but barely seen in the CO, are located at ~1" (250 AU) from the central source, L1448C(N). Since the dynamical time scale for the innermost pair is only ~10 yr, SiO may have been formed in the protostellar wind through the gas-phase reaction, or been formed on the dust grain and directly released into the gas phase by means of shocks. It is found that the jet is extremely active with a mechanical luminosity of ~7 L_sun, which is comparable to the bolometric luminosity of the central source (7.5 L_sun). The mass accretion rate onto the protostar derived from the mass-loss rate is ~10^{-5} M_sun/yr. Such a high mass accretion rate suggests that the mass and the age of the central star are 0.03-0.09 M_sun and (4-12)x10^3 yr, respectively, implying that the central star is in the very early stage of protostellar evolution. The low-velocity CO emission delineates two V-shaped shells with a common apex at L1448C(N). The kinematics of these shells are reproduced by the model of a wide opening angle wind. The co-existence of the highly-collimated jets and the wide-opening angle shells can be explained by the unified X-wind model" in which highly-collimated jet components correspond to the on-axis density enhancement of the wide-opening angle wind. The CO $J$=3--2 map also revealed the second outflow driven by the southern source L1448C(S) located at ~8.3" (2000 AU) from L1448C(N).

preprint2010arXiv

First Interferometric Images of the 36 GHz Methanol Masers in the DR21 Complex

Class I methanol masers are believed to be produced in the shock-excited environment around star-forming regions. Many authors have argued that the appearance of various subsets of class I masers may be indicative of specific evolutionary stages of star formation or excitation conditions. Until recently, however, no major interferometer was capable of imaging the important 36 GHz transition. We report on Expanded Very Large Array observations of the 36 GHz methanol masers and Submillimeter Array observations of the 229 GHz methanol masers in DR21(OH), DR21N, and DR21W. The distribution of 36 GHz masers in the outflow of DR21(OH) is similar to that of the other class I methanol transitions, with numerous multitransition spatial overlaps. At the site of the main continuum source in DR21(OH), class I masers at 36 and 229 GHz are found in virtual overlap with class II 6.7 GHz masers. To the south of the outflow, the 36 GHz masers are scattered over a large region but usually do not appear coincident with 44 GHz masers. In DR21W we detect an "S-curve" signature in Stokes V that implies a large value of the magnetic field strength if interpreted as due to Zeeman splitting, suggesting either that class I masers may exist at higher densities than previously believed or that the direct Zeeman interpretation of S-curve Stokes V profiles in class I masers may be incorrect. We find a diverse variety of different maser phenomena in these sources, suggestive of differing physical conditions among them.

preprint2010arXiv

Grain growth across protoplanetary discs: 10-micron silicate feature versus millimetre slope

Young stars are formed within dusty discs. The grains in the disc are originally of the same size as interstellar dust. Models predict that these grains will grow in size through coagulation. Observations of the silicate features at micron wavelengths are consistent with growth to micron sizes whereas the slope of the SED at longer wavelengths traces growth up to mm sizes. We here look for a correlation between these two grain growth indicators. A large sample of T-Tauri and Herbig-Ae/Be stars was observed with the Spitzer Space Telescope at 5-13 micron; a subsample was observed at mm wavelengths. We complement this subsample with data from the literature to maximise the overlap between micron and mm observations and search for correlations. Synthetic spectra are produced to determine which processes may produce the dust evolution. Dust disc masses in the range <1 to 7 x 10^-4 MSun are obtained. Most sources have a mm spectral slope consistent with grain growth. There is a tentative correlation between the 10-micron silicate feature and the mm slope of the SED. The observed sources seem to be grouped per star-forming region in the micron-vs-mm diagram. The modelling results show that the 10-micron feature becomes flatter and subsequently the mm slope becomes shallower. Grain size distributions shallower than that of the ISM and/or bright central stars are required to explain specific features. Settling of larger grains towards the disc midplane affects the 10-micron feature, but hardly the mm slope. The tentative correlation between the strength of the 10-micron feature and the mm slope suggests that the inner and outer disc evolve simultaneously. Dust with a mass dominated by mm-sized grains is required to explain the shallowest mm slopes. Other processes besides grain growth may also be responsible for the removal of small grains.

preprint2010arXiv

L1448 IRS2E: A candidate first hydrostatic core

Intermediate between the prestellar and Class 0 protostellar phases, the first core is a quasi-equilibrium hydrostatic object with a short lifetime and an extremely low luminosity. Recent MHD simulations suggest that the first core can even drive a molecular outflow before the formation of the second core (i.e., protostar). Using the Submillimeter Array and the Spitzer Space Telescope, we present high angular resolution observations towards the embedded dense core IRS2E in L1448. We find that source L1448 IRS2E is not visible in the sensitive Spitzer infrared images (at wavelengths from 3.6 to 70 um), and has weak (sub-)millimeter dust continuum emission. Consequently, this source has an extremely low bolometric luminosity (< 0.1 L_sun). Infrared and (sub-)millimeter observations clearly show an outflow emanating from this source; L1448 IRS2E represents thus far the lowest luminosity source known to be driving a molecular outflow. Comparisons with prestellar cores and Class 0 protostars suggest that L1448 IRS2E is more evolved than prestellar cores but less evolved than Class 0 protostars, i.e., at a stage intermediate between prestellar cores and Class 0 protostars. All these results are consistent with the theoretical predictions of the radiative/magneto hydrodynamical simulations, making L1448 IRS2E the most promising candidate of the first hydrostatic core revealed so far.

preprint2010arXiv

The Spitzer c2d Survey of Nearby Dense Cores. IX. Discovery of a Very Low Luminosity Object Driving a Molecular Outflow in the Dense Core L673-7

We present new infrared, submillimeter, and millimeter observations of the dense core L673-7 and report the discovery of a low-luminosity, embedded Class 0 protostar driving a molecular outflow. L673-7 is seen in absorption against the mid-infrared background in 5.8, 8, and 24 micron Spitzer images, allowing for a derivation of the column density profile and total enclosed mass of L673-7, independent of dust temperature assumptions. Estimates of the core mass from these absorption profiles range from 0.2-4.5 solar masses. Millimeter continuum emission indicates a mass of about 2 solar masses, both from a direct calculation assuming isothermal dust and from dust radiative transfer models constrained by the millimeter observations. We use dust radiative transfer models to constrain the internal luminosity of L673-7, defined to be the luminosity of the central source and excluding the luminosity from external heating, to be 0.01-0.045 solar luminosities, with 0.04 solar luminosities the most likely value. L673-7 is thus classified as a very low luminosity object (VeLLO), and is among the lowest luminosity VeLLOs yet studied. We calculate the kinematic and dynamic properties of the molecular outflow in the standard manner, and we show that the expected accretion luminosity based on these outflow properties is greater than or equal to 0.36 solar luminosities. The discrepancy between this expected accretion luminosity and the internal luminosity derived from dust radiative transfer models indicates that the current accretion rate is much lower than the average rate over the lifetime of the outflow. Although the protostar embedded within L673-7 is consistent with currently being substellar, it is unlikely to remain as such given the substantial mass reservoir remaining in the core.

preprint2009arXiv

PROSAC: A Submillimeter Array survey of low-mass protostars. II. The mass evolution of envelopes, disks, and stars from the Class 0 through I stages

The key question about early protostellar evolution is how matter is accreted from the large-scale molecular cloud, through the circumstellar disk onto the central star. A sample of 20 Class 0 and I protostars has been observed in continuum at (sub)millimeter wavelengths at high angular resolution with the Submillimeter Array. Using detailed dust radiative transfer models, we have developed a framework for disentangling the continuum emission from the envelopes and disks, and from that estimated their masses. For the Class I sources in the sample, HCO+ 3-2 line emission has furthermore been observed with the Submillimeter Array. Four of these sources show signs of Keplerian rotation, constraining the masses of the central stars. Both Class 0 and I protostars are surrounded by disks with typical masses of about 0.05 M_sun. No evidence is found for a correlation between the disk mass and evolutionary stage of the young stellar objects. This contrasts the envelope mass, which decreases sharply from 1 M_sun in the Class 0 stage to <0.1 M_sun in the Class I stage. Typically, the disks have masses that are 1-10% of the corresponding envelope masses in the Class 0 stage and 20-60% in the Class I stage. For the Class I sources for which Keplerian rotation is seen, the central stars contain 70-98% of the total mass in the star-disk-envelope system, confirming that these objects are late in their evolution through the embedded protostellar stages. Theoretical models tend to overestimate the disk masses relative to the stellar masses in the late Class I stage. The results argue in favor of a picture in which circumstellar disks are formed early during the protostellar evolution and rapidly process material accreted from the larger scale envelope onto the central star.

preprint2009arXiv

Rotation and Outflow motions in the very low-mass Class 0 protostellar system HH 211 at subarcsecond resolution

HH 211 is a nearby young protostellar system with a highly collimated jet. We have mapped it in 352 GHz continuum, SiO (J=8-7), and HCO+ (J=4-3) emission at up to ~ 0.2" resolution with the Submillimeter Array (SMA). The continuum source is now resolved into two sources, SMM1 and SMM2, with a separation of ~ 84 AU. SMM1 is seen at the center of the jet, probably tracing a (inner) dusty disk around the protostar driving the jet. SMM2 is seen to the southwest of SMM1 and may trace an envelope-disk around a small binary companion. A flattened envelope-disk is seen in HCO+ around SMM1 with a radius of ~ 80 AU perpendicular to the jet axis. Its velocity structure is consistent with a rotation motion and can be fitted with a Keplerian law that yields a mass of ~ 50+-15 Jupiter mass (a mass of a brown dwarf) for the protostar. Thus, the protostar could be the lowest mass source known to have a collimated jet and a rotating flattened envelope-disk. A small-scale (~ 200 AU) low-speed (~ 2 km/s) outflow is seen in HCOP+ around the jet axis extending from the envelope-disk. It seems to rotate in the same direction as the envelope-disk and may carry away part of the angular momentum from the envelope-disk. The jet is seen in SiO close to ~ 100 AU from SMM1. It is seen with a "C-shaped" bending. It has a transverse width of <~ 40 AU and a velocity of ~ 170+-60 km/s. A possible velocity gradient is seen consistently across its innermost pair of knots, with ~ 0.5 km/s at ~ 10 AU, consistent with the sense of rotation of the envelope-disk. If this gradient is an upper limit of the true rotational gradient of the jet, then the jet carries away a very small amount of angular momentum of ~ 5 AU km/s and thus must be launched from the very inner edge of the disk near the corotation radius.

preprint2009arXiv

The Spitzer c2d Survey of Nearby Dense Cores VII: Chemistry and Dynamics in L43

We present results from the Spitzer Space Telescope and molecular line observations of 9 species toward the dark cloud L43. The Spitzer images and molecular line maps suggest it has a starless core and a Class I protostar evolving in the same environment. CO depletion is seen in both sources, and DCO+ lines are stronger toward the starless core. With a goal of testing the chemical characteristics from pre- to protostellar stages, we adopt an evolutionary chemical model to calculate the molecular abundances and compare with our observations. Among the different model parameters we tested, the best-fit model suggests a longer total timescale at the pre-protostellar stage, but with faster evolution at the later steps with higher densities.

preprint2009arXiv

The Spitzer c2d Survey of Nearby Dense Cores: Jet and Molecular Outflow Associated with a YSO in core A of L1251

A long infrared jet has been discovered by the Spitzer c2d legacy program in core A of L1251. It is associated with a very embedded Class 0 object with an accretion luminosity of about 0.9 Lsun derived by radiative transfer model fitting to the observed SED. Comparing the observed IRAC colors along the infrared jet with those calculated from a model of an admixture of gas with a power-law temperature distribution indicates that the jet is possibly created by a paraboloidal bow shock propagating into the ambient medium of n(H_2)=10^5 cm^{-3}. In addition, the variation of the power-law index along the jet suggests that the portion of hot gas decreases with distance from the jet engine. The molecular outflow in this region has been mapped for the first time using CO data. From the calculated outflow momentum flux, a very strong lower limit to the average accretion luminosity is 3.6 sin i/cos^3 i Lsun, indicative of a decrease in the accretion rate with time.

preprint2008arXiv

The Spitzer c2d Survey of Nearby Dense Cores. V. Discovery of a VeLLO in the "Starless" Dense Core L328

This paper reports the discovery of a Very Low Luminosity Object (VeLLO) in the "starless" dense core L328, using the Spitzer Space Telescope and ground based observations from near-infrared to millimeter wavelengths. The Spitzer 8 micron image indicates that L328 consists of three subcores of which the smallest one may harbor a source, L328-IRS while two other subcores remain starless. L328-IRS is a Class 0 protostar according to its bolometric temperature (44 K) and the high fraction ~72 % of its luminosity emitted at sub-millimeter wavelengths. Its inferred "internal luminosity" (0.04 - 0.06 Lsun) using a radiative transfer model under the most plausible assumption of its distance as 200 pc is much fainter than for a typical protostar, and even fainter than other VeLLOs studied previously. Note, however, that its inferred luminosity may be uncertain by a factor of 2-3 if we consider two extreme values of the distance of L328-IRS (125 or 310 pc). Low angular resolution observations of CO do not show any clear evidence of a molecular outflow activity. But broad line widths toward L328, and Spitzer and near-infrared images showing nebulosity possibly tracing an outflow cavity, strongly suggest the existence of outflow activity. Provided that an envelope of at most ~0.1 Msunis the only mass accretion reservoir for L328-IRS, and the star formation efficiency is close to the canonical value ~30%, L328-IRS has not yet accreted more than 0.05 Msun. At the assumed distance of 200 pc, L328-IRS is destined to be a brown dwarf.

preprint2007arXiv

PROSAC: A Submillimeter Array Survey of Low-Mass Protostars. I. Overview of Program: Envelopes, Disks, Outflows and Hot Cores

This paper presents a large spectral line and continuum survey of 8 deeply embedded, low-mass protostellar cores using the Submillimeter Array. Each source was observed in high excitation lines of some of the most common molecular species, CO, HCO+, CS, SO, H2CO, CH3OH and SiO. Line emission from 11 species originating from warm and dense gas have been imaged at high angular resolution (1-3"; typically 200-600 AU) together with continuum emission at 230 GHz (1.3 mm) and 345 GHz (0.8 mm). Compact continuum emission is observed for all sources which likely originates in marginally optically thick circumstellar disks, with typical lower limits to their masses of 0.1 M_sun (1-10% of the masses of their envelopes) and having a dust opacity law with beta approximately 1. Prominent outflows are present in CO 2-1 observations in all sources: the most diffuse outflows are found in the sources with the lowest ratios of disk-to-envelope mass, and it is suggested that these sources are in a phase where accretion of matter from the envelope has almost finished and the remainder of the envelope material is being dispersed by the outflows. Other characteristic dynamical signatures are found with inverse P Cygni profiles indicative of infalling motions seen in the 13CO 2-1 lines toward NGC1333-IRAS4A and -IRAS4B. Outflow-induced shocks are present on all scales in the protostellar environments and are most clearly traced by the emission of CH3OH in NGC1333-IRAS4A and -IRAS4B. These observations suggest that the emission of CH3OH and H2CO from these proposed "hot corinos" are related to the shocks caused by the protostellar outflows. Only one source, NGC1333-IRAS2A, has evidence for hot, compact CH3OH emission coincident with the embedded protostar.

preprint2006arXiv

SHARC-II Mapping of Spitzer c2d Small Clouds and Cores

We present the results of a submillimeter survey of 53 low-mass dense cores with the Submillimeter High Angular Resolution Camera II (SHARC-II). The survey is a follow-up project to the Spitzer Legacy Program ``From Molecular Cores to Planet-Forming Disks'', with the purpose being to create a complete data set of nearby low-mass dense cores from the infrared to the millimeter. We present maps of 52 cores at 350 microns and three cores at 450 microns, two of which were observed at both wavelengths. Of these 52 cores, 41 were detected by SHARC-II: 32 contained one submillimeter source while 9 contained multiple sources. For each submillimeter source detected, we report various source properties including source position, fluxes in various apertures, size, aspect ratio, and position angle. For the 12 cores that were not detected we present upper limits. The sources detected by SHARC-II have, on average, smaller sizes at the 2sigma contours than those derived from longer-wavelength bolometer observations. We conclude that this is not caused by a failure to integrate long enough to detect the full extent of the core; instead it arises primarily from the fact that the observations presented in this survey are insensitive to smoothly varying extended emission. We find that SHARC-II observations of low-mass cores are much better suited to distinguishing between starless and protostellar cores than observations at longer wavelengths. Very Low Luminosity Objects, a new class of objects being discovered by the Spitzer Space Telescope in cores previously classified as starless, look very similar at 350 microns to other cores with more luminous protostars.

preprint2006arXiv

The Spitzer c2d Survey of Large, Nearby, Interstellar Clouds. III. Perseus Observed with IRAC

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.

preprint2006arXiv

The Spitzer c2d survey of Large, Nearby, Interstellar Clouds. V. Chamaeleon II Observed with IRAC

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.

preprint2006arXiv

The Spitzer c2d Survey of Nearby Dense Cores: II: Discovery of a Low Luminosity Object in the "Evolved Starless Core" L1521F

We present Spitzer Space Telescope observations of the "evolved starless core" L1521F which reveal the presence of a very low luminosity object (L < 0.07 Lsun). The object, L1521F-IRS, is directly detected at mid-infrared wavelengths (>5 micron) but only in scattered light at shorter infrared wavelengths, showing a bipolar nebula oriented east-west which is probably tracing an outflow cavity. The nebula strongly suggests that L1521F-IRS is embedded in the L1521F core. Thus L1521F-IRS is similar to the recently discovered L1014-IRS and the previously known IRAM 04191 in its substellar luminosity and dense core environment. However these objects differ significantly in their core density, core chemistry, and outflow properties, and some may be destined to be brown dwarfs rather than stars.

preprint2006arXiv

The Spitzer c2d Survey of Nearby Dense Cores: III: Low Mass Star Formation in a Small Group, L1251B

We present a comprehensive study of a low-mass star-forming region,L1251B, at wavelengths from the near-infrared to the millimeter. L1251B, where only one protostar, IRAS 22376+7455, was known previously, is confirmed to be a small group of protostars based on observations with the Spitzer Space Telescope. The most luminous source of L1251B is located 5" north of the IRAS position. A near-infrared bipolar nebula, which is not associated with the brightest object and is located at the southeast corner of L1251B, has been detected in the IRAC bands. OVRO and SMA interferometric observations indicate that the brightest source and the bipolar nebula source in the IRAC bands are deeply embedded disk sources.Submillimeter continuum observations with single-dish telescopes and the SMA interferometric observations suggest two possible prestellar objects with very high column densities. Outside of the small group, many young stellar object candidates have been detected over a larger region of 12' x 12'. Extended emission to the east of L1251B has been detected at 850 micron; this "east core" may be a site for future star formation since no point source has been detected with IRAC or MIPS. This region is therefore a possible example of low-mass cluster formation, where a small group of pre- and protostellar objects (L1251B) is currently forming, alongside a large starless core (the east core).