Source author record

Shih-Ping Lai

Shih-Ping Lai 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

41works
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

41 published item(s)

preprint2022arXiv

A VLA View of the Flared, Asymmetric Disk Around the Class 0 Protostar L1527 IRS

We present high resolution Karl G. Jansky Very Large Array (VLA) observations of the protostar L1527 IRS at 7 mm, 1.3 cm, and 2 cm wavelengths. We detect the edge-on dust disk at all three wavelengths and find that it is asymmetric, with the southern side of the disk brighter than the northern side. We confirm this asymmetry through analytic modeling and also find that the disk is flared at 7 mm. We test the data against models including gap features in the intensity profile, and though we cannot rule such models out, they do not provide a statistically significant improvement in the quality of fit to the data. From these fits, we can however place constraints on allowed properties of any gaps that could be present in the true, underlying intensity profile. The physical nature of the asymmetry is difficult to associate with physical features due to the edge-on nature of the disk, but could be related to spiral arms or asymmetries seen in other imaging of more face-on disks.

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

Magnetic fields and outflows in the large Bok globule CB 54

We have observed the large Bok globule CB 54 in 850$μ$m polarised light using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT). We find that the magnetic field in the periphery of the globule shows significant, ordered deviation from the mean field direction in the globule centre. This deviation appears to correspond with the extended but relatively weak $^{12}$CO outflow emanating from the Class 0 sources at the centre of the globule. Energetics analysis suggests that if the outflow is reshaping the magnetic field in the globule's periphery, then we can place an upper limit of $<27\,μ$G on the magnetic field strength in the globule's periphery. Comparison with archival Planck and CARMA measurements shows that the field in the centre of the globule is consistent over several orders of magnitude in size scale, and oriented parallel to the density structure in the region in projection. We thus hypothesise that while non-thermal motions in the region may be sub-Alfvénic, the magnetic field is subdominant to gravity over a wide range of size scales. Our results suggest that even a relatively weak outflow may be able to significantly reshape magnetic fields in star-forming regions on scales $> 0.1$ pc.

preprint2022arXiv

Revealing the dust grain polarization properties as a function of extinction and distance towards NGC 1893

Dust polarization observations at optical wavelengths help to understand the dust grain properties and trace the plane-of-the-sky component of the magnetic field. In this study, we make use of the $I$-band polarization data acquired from AIMPOL along with the distances ($d$) and extinction ($A_{V}$) data to study the variation of polarization fraction ($P$) as a function of $A_{V}$ and $d$ towards the star-forming region, NGC 1893. We employ a broken power-law fit and Bayesian analysis on extinction ($A_{V}$) versus polarization efficiency ($P$/$A_{V}$) and distance ($d$) versus rate of polarization ($P$/$d$). We find that $P$/$A_{V}$ shows a break at an extinction of $\sim$0.9 mag, whereas $P/d$ exhibits a break at a distance of $\sim$1.5 kpc. Based on these, we categorize the dust towards NGC 1893 into two populations: (i) foreground dust confined to $A_{V}$ $<$ $\sim$1 mag and distance up to $\sim$2 kpc and (ii) Perseus spiral arm dust towards NGC 1893 characterized with $A_{V}$ $>$ $\sim$1 mag and distance beyond $\sim$2 kpc. Foreground dust exhibits higher polarization efficiency but a lower polarization rate, whereas Perseus dust shows a lower polarization efficiency but a slightly higher polarization rate. Hence, we suggest that while polarization efficiency reveals the dust grain alignment, the rate of polarization infers about the distribution of dust grains towards NGC 1893. Further, we also shed a light on the spatial variation of intrinsic polarization and magnetic field orientation, and other parameters within the intra-cluster medium of NGC 1893.

preprint2021arXiv

Accretion Flows or Outflow Cavities? Uncovering the Gas Dynamics around Lupus 3-MMS

Understanding how material accretes onto the rotationally supported disk from the surrounding envelope of gas and dust in the youngest protostellar systems is important for describing how disks are formed. Magnetohydrodynamic simulations of magnetized, turbulent disk formation usually show spiral-like streams of material (accretion flows) connecting the envelope to the disk. However, accretion flows in these early stages of protostellar formation still remain poorly characterized due to their low intensity and possibly some extended structures are disregarded as being part of the outflow cavity. We use ALMA archival data of a young Class 0 protostar, Lupus 3-MMS, to uncover four extended accretion flow-like structures in C$^{18}$O that follow the edges of the outflows. We make various types of position-velocity cuts to compare with the outflows and find the extended structures are not consistent with the outflow emission, but rather more consistent with a simple infall model. We then use a dendrogram algorithm to isolate five sub-structures in position-position-velocity space. Four out of the five sub-structures fit well ($>$95%) with our simple infall model, with specific angular momenta between $2.7-6.9\times10^{-4}\,$km$\,$s$^{-1}\,$pc and mass-infall rates of $0.5-1.1\times10^{-6}\,M_{\odot}\,$yr$^{-1}$. Better characterization of the physical structure in the supposed "outflow-cavities" is important to disentangle the true outflow cavities and accretion flows.

preprint2021arXiv

An observational correlation between magnetic field, angular momentum and fragmentation in the envelopes of Class 0 protostars?

To assess the potential role of magnetic fields in regulating the envelope rotation and the fragmentation of Class 0 protostars, we carried out observations of the dust polarized emission at 0.87 mm with the SMA, in the envelopes of a large sample of 20 Class 0 protostars. We estimate the mean magnetic field orientation over the central 1000 au envelope scales and compared it to that of the protostellar outflow in order to study the relation between their misalignment and the kinematics of the circumstellar gas. We discover a strong relationship between the misalignment of the magnetic field orientation with the outflow and the amount of angular momentum observed at similar scales in the protostellar envelope, revealing a potential link between the kinetic and the magnetic energy at envelope scales. The relation could be driven by favored B misalignments in more dynamical envelopes or a dependence of the envelope dynamics with the large-scale B initial configuration. Comparing the trend with the presence of fragmentation, we observe that single sources are mostly associated with conditions of low angular momentum in the inner envelope and good alignment of the magnetic field with protostellar outflows, at intermediate scales. Our results suggest that the properties of the magnetic field in protostellar envelopes bear a tight relationship with the rotating-infalling gas directly involved in the star and disk formation: we find that it may not only influence the fragmentation of protostellar cores into multiple stellar systems, but also set the conditions establishing the pristine properties of planet-forming disks.

preprint2020arXiv

Determining the physical conditions of extremely young Class 0 circumbinary disk around VLA1623A

We present detailed analysis of high-resolution C18O (2-1), SO (88-77), CO (3-2) and DCO+ (3-2) data obtained by the Atacama Large Millimeter/sub-millimeter Array (ALMA) towards a Class 0 Keplerian circumbinary disk around VLA1623A, which represents one of the most complete analysis towards a Class 0 source. From the dendrogram analysis, we identified several accretion flows feeding the circumbinary disk in a highly anisotropic manner. Stream-like SO emission around the circumbinary disk reveals the complicated shocks caused by the interactions between the disk, accretion flows and outflows. A wall-like structure is discovered south of VLA1623B. The discovery of two outflow cavity walls at the same position traveling at different velocities suggests the two outflows from both VLA1623A and VLA1623B overlays on top of each other in the plane of sky. Our detailed flat and flared disk modeling shows that Cycle 2 C18O J = 2-1 data is inconsistent with the combined binary mass of 0.2 Msun as suggested by early Cycle 0 studies. The combined binary mass for VLA1623A should be modified to 0.3 ~ 0.5 Msun.

preprint2020arXiv

Formation of the Hub-Filament System G33.92+0.11: Local Interplay between Gravity, Velocity, and Magnetic Field

The formation of filaments in molecular clouds is an important process in star formation. Hub-filament systems (HFSs) are a transition stage connecting parsec-scale filaments and proto-clusters. Understanding the origin of HFSs is crucial to reveal how star formation proceeds from clouds to cores. Here, we report JCMT POL-2 850 $μ$m polarization and IRAM 30-m C$^{18}$O (2-1) line observations toward the massive HFS G33.92+0.11. The 850 $μ$m continuum map reveals four major filaments converging to the center of G33.92+0.11 with numerous short filaments connecting to the major filaments at local intensity peaks. We estimate the local orientations of filaments, magnetic field, gravity, and velocity gradients from observations, and we examine their correlations based on their local properties. In the high-density areas, our analysis shows that the filaments tend to align with the magnetic field and local gravity. In the low-density areas, we find that the local velocity gradients tend to be perpendicular to both the magnetic field and local gravity, although the filaments still tend to align with local gravity. A global virial analysis suggests that the gravitational energy overall dominates the magnetic and kinematic energy. Combining local and global aspects, we conclude that the formation of G33.92+0.11 is predominantly driven by gravity, dragging and aligning the major filaments and magnetic field on the way to the inner dense center. Traced by local velocity gradients in the outer diffuse areas, ambient gas might be accreted onto the major filaments directly or via the short filaments.

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

Physical and chemical modeling of the starless core L1512

The deuterium fractionation in starless cores gives us a clue to estimate their lifetime scales, thus allowing us to distinguish between different dynamical theories of core formation. Cores also seem to be subject to a differential N2 and CO depletion which was not expected from models. We aim to make a survey of 10 cores to estimate their lifetime scales and depletion profiles in detail. After L183, in Serpens, we present the second cloud of the series, L1512 in Auriga. To constrain the lifetime scale, we perform chemical modeling of the deuteration profiles across L1512 based on dust extinction measurements from near-infrared observations and non-local thermal equilibrium radiative transfer with multiple line observations of N2H+, N2D+, DCO+, C18O, and 13CO, plus H2D+ (1$_{10}$--1$_{11}$). We find a peak density of 1.1$\times$10$^5$ cm$^{-3}$ and a central temperature of 7.5$\pm$1 K, which are respectively higher and lower compared with previous dust emission studies. The depletion factors of N2H+ and N2D+ are 27$^{+17}_{-13}$ and 4$^{+2}_{-1}$ in L1512, intermediate between the two other more advanced and denser starless core cases, L183 and L1544. These factors also indicate a similar freeze-out of N2 in L1512, compared to the two others despite a peak density one to two orders of magnitude lower. Retrieving CO and N2 abundance profiles with the chemical model, we find that CO has a depletion factor of $\sim$430-870 and the N2 profile is similar to that of CO unlike towards L183. Therefore, L1512 has probably been living long enough so that N2 chemistry has reached steady state. N2H+ modeling remains compulsory to assess the precise physical conditions in the center of cold starless cores, rather than dust emission. L1512 is presumably older than 1.4 Myr. Therefore, the dominating core formation mechanism should be ambipolar diffusion for this source.

preprint2020arXiv

Probing the Temperature Structure of Optically Thick Disks Using Polarized Emission of Aligned Grains

Polarized continuum emission from aligned grains in disks around young stellar objects can be used to probe the magnetic field, radiation anisotropy, or drift between dust and gas, depending on whether the non-spherical grains are aligned magnetically, radiatively or mechanically. We show that it can also be used to probe another key disk property -- the temperature gradient -- along sight lines that are optically thick, independent of the grain alignment mechanism. We first illustrate the technique analytically using a simple 1D slab model, which yields an approximate formula that relates the polarization fraction to the temperature gradient with respect to the optical depth tau at the tau=1 surface. The formula is then validated using models of stellar irradiated disks with and without accretion heating. The promises and challenges of the technique are illustrated with a number of Class 0 and I disks with ALMA dust polarization data, including NGC 1333 IRAS4A1, IRAS 16293B, BHB 07-11, L1527, HH 212 and HH 111. We find, in particular, that the sight lines passing through the near-side of a highly inclined disk trace different temperature gradient directions than those through the far-side, which can lead to a polarization orientation on the near-side that is orthogonal to that on the far-side, and that the HH 111 disk may be such a case. Our technique for probing the disk temperature gradient through dust polarization can complement other methods, particularly those using molecular lines.

preprint2019arXiv

Resolving Linear Polarization due to Emission and Extinction of Aligned Dust Grains on NGC1333 IRAS4A with JVLA and ALMA

We report high angular resolution observations of linearly polarized dust emission towards the Class 0 young stellar object (YSO) NGC1333 IRAS4A (hereafter, IRAS4A) using the Karl G. Jansky Very Large Array (JVLA) at K (11.5-16.7 mm), Ka (8.1-10.3 mm), and Q bands (6.3-7.9 mm), and using the Atacama Large Millimeter Array (ALMA) at Band 6 (1.2 mm) and Band 7 (0.85-0.89 mm). On 100-1000 AU scales, all of these observations consistently trace the hourglass shaped magnetic field topology as shown in the previous studies. In the innermost 100 AU region of IRAS4A1, the polarization position angles (E-field) detected at 6.3-16.7 mm are consistent, however, are nearly 90 degrees offset from those detected at 1.2 mm and 0.85-0.89 mm. Such a 90 degree offset may be explained by that the inner $\sim$100 AU area is optically thick at wavelengths shorter than $\sim$1.5 mm, whereby the observations probe the absorption of aligned dust against the weakly or unpolarized warm dust emission from the innermost region. This can also consistently explain why the highest angular resolution ALMA images at Band 7 show that the polarization percentage increases with dust brightness temperature in the inner $\sim$100 AU region of IRAS4A1. Following this interpretation and assuming that the dust grains are aligned with the magnetic fields, the inferred magnetic field position angle based on the 90$^{\circ}$ rotated at 6.3-7.9 mm in the central peak of IRAS4A1 is $\sim-22^{\circ}$, which is approximately consistent with the outflow direction $\sim-9^{\circ}$.

preprint2016arXiv

Detection of Linearly Polarized 6.9 mm Continuum Emission from the Class 0 Young Stellar Object NGC1333 IRAS4A

We report new Karl G. Jansky Very Large Array (JVLA), 0$"$.5 angular resolution observations of linearly polarized continuum emission at 6.9 mm, towards the Class 0 young stellar object (YSO) NGC1333 IRAS4A. This target source is a collapsing dense molecular core, which was resolved at short wavelengths to have hourglass shaped B-field configuration. We compare these 6.9 mm observations with previous polarization Submillimeter Array (SMA) observations at 0.88 mm, which have comparable angular resolution ($\sim$0$"$7). We found that at the same resolution, the observed polarization position angles at 6.9 mm are slightly deviated from those observed at 0.88 mm. Due to the lower optical depth of the emission at 6.9 mm, and the potential effect of dust grain growth, the new JVLA observations are likely probing B-field alignments in regions interior to those sampled by the previous polarization observations at higher frequencies. Our understanding can be improved by more sensitive observations, and observations for the more extended spatial scales.

preprint2016arXiv

Helical Magnetic Fields in the NGC1333 IRAS 4A Protostellar Outflows

We present Submillimeter Array (SMA) polarization observations of the CO $J$ = 3--2 line toward the NGC1333 IRAS 4A. The CO Stokes $I$ maps at an angular resolution of $\sim$1$\arcsec$ reveal two bipolar outflows from the binary sources of the NGC 1333 IRAS 4A. The kinematic features of the CO emission can be modeled by wind-driven outflows at $\sim$ 20$\arcdeg$ inclined from the plane of the sky. Close to the protostars the CO polarization, at an angular resolution of $\sim$$2\farcs3$, has a position angle approximately parallel to the magnetic field direction inferred from the dust polarizations. The CO polarization direction appears to vary smoothly from an hourglass field around the core to an arc-like morphology wrapping around the outflow, suggesting a helical structure of magnetic fields that inherits the poloidal fields at the launching point and consists of toroidal fields at a farther distance of outflow. The helical magnetic field is consistent with the theoretical expectations for launching and collimating outflows from a magnetized rotating disk. Considering that the CO polarized emission is mainly contributed from the low-velocity and low-resolution data, the helical magnetic field is likely a product of the wind-envelope interaction in the wind-driven outflows. The CO data reveal a PA of $\sim$ 30$\arcdeg$ deflection in the outflows. The variation in the CO polarization angle seems to correlate with the deflections. We speculate that the helical magnetic field contributes to $\sim$ 10$\arcdeg$ deflection of the outflows by means of Lorenz force.

preprint2016arXiv

Interferometric mapping of magnetic fields: The ALMA view of the massive star forming clump W43-MM1

Here we present the first results from ALMA observations of 1 mm polarized dust emission towards the W43-MM1 high mass star forming clump. We have detected a highly fragmented filament with source masses ranging from 14Msun to 312Msun, where the largest fragment, source A, is believed to be one of the most massive in our Galaxy. We found a smooth, ordered, and detailed polarization pattern throughout the filament which we used to derived magnetic field morphologies and strengths for 12 out of the 15 fragments detected ranging from 0.2 to 9 mG. The dynamical equilibrium of each fragment was evaluated finding that all the fragments are in a super-critical state which is consistent with previously detected infalling motions towards W43-MM1. Moreover, there are indications suggesting that the field is being dragged by gravity as the whole filament is collapsing.

preprint2016arXiv

IRAS 16253-2429: the First Proto-Brown Dwarf Binary Candidate Identified through Dynamics of Jets

The formation mechanism of brown dwarfs (BDs) is one of the long-standing problems in star formation because the typical Jeans mass in molecular clouds is too large to form these substellar objects. To answer this question, it is crucial to study a BD at the embedded phase. IRAS 16253-2429 is classified as a very low luminosity object (VeLLO) with internal luminosity 0.1 Lsun. VeLLOs are believed to be very low-mass protostars or even proto-BDs. We observed the jet/outflow driven by IRAS 16253-2429 in CO (2-1), (6-5), and (7-6) using the IRAM 30 m and APEX telescopes and the SMA in order to study its dynamical features and physical properties. Our SMA map reveals two protostellar jets, indicating the existence of a proto-binary system as implied by the precessing jet detected in H2 emission. We detect a wiggling pattern in the position-velocity diagrams along the jet axes, which is likely due to the binary orbital motion. Based on this, we derive the current mass of the binary as ~0.032 Msun. Given the low envelope mass, IRAS 16253-2429 will form a binary that probably consist of one or two BDs. Furthermore, we found that the outflow force as well as the mass accretion rate are very low based on the multi-transition CO observations, which suggests that the final masses of the binary components are at the stellar/substellar boundary. Since IRAS 16253 is located in an isolated environment, we suggest that BDs can form through fragmentation and collapse like low-mass stars.

preprint2016arXiv

The intrinsic abundance ratio and X-factor of CO isotopologues in L1551 shielded from FUV photodissociation

L1551 is chosen because it is relatively isolated in the Taurus molecular cloud shielded from FUV photons, providing an ideal environment for studying the target properties. Our observations cover ~40'x40' with resolution ~30", which are the maps with highest spatial dynamical range to date. We derive the X(13CO)/X(C18O) value on the sub-parsec scales in the range of ~3-27 with a mean value of 8.0+-2.8. Comparing to the visual extinction map derived from the Herschel observations, we found that the abundance ratio reaches its maximum at low Av (i.e., Av ~ 1-4mag), and decreases to the typical solar system value of 5.5 inside L1551 MC. The high X(13CO)/X(C18O) value at the boundary of the cloud is most likely due to the selective FUV photodissociation of C18O. This is in contrast with Orion-A where its internal OB stars keep the abundance ratio at a high level greater than ~10. In addition, we explore the variation of the X-factor, because it is an uncertain but widely used quantity in extragalactic studies. We found that X-factor is proportional to N(H2)^1.0 which is consistent with previous simulations. Excluding the high density region, the average X-factor is similar to the Milky Way average value.

preprint2015arXiv

A low-mass protostar's disk-envelope interface: disk-shadowing evidence from ALMA DCO+ observations of VLA1623

Due to instrumental limitations and a lack of disk detections, the structure between the envelope and the rotationally supported disk has been poorly studied. This is now possible with ALMA through observations of CO isotopologs and tracers of freezeout. Class 0 sources are ideal for such studies given their almost intact envelope and young disk. The structure of the disk-envelope interface of the prototypical Class 0 source, VLA1623A which has a confirmed Keplerian disk, is constrained from ALMA observations of DCO+ 3-2 and C18O 2-1. The physical structure of VLA1623 is obtained from the large-scale SED and continuum radiative transfer. An analytic model using a simple network coupled with radial density and temperature profiles is used as input for a 2D line radiative transfer calculation for comparison with the ALMA Cycle 0 12m array and Cycle 2 ACA observations of VLA1623. DCO+ emission shows a clumpy structure bordering VLA1623A's Keplerian disk, suggesting a cold ring-like structure at the disk-envelope interface. The radial position of the observed DCO+ peak is reproduced in our model only if the region's temperature is between 11-16K, lower than expected from models constrained by continuum and SED. Altering the density has little effect on the DCO+ position, but increased density is needed to reproduce the disk traced in C18O. The DCO+ emission around VLA1623A is the product of shadowing of the envelope by the disk. Disk-shadowing causes a drop in the gas temperature outside of the disk on >200AU scales, encouraging deuterated molecule production. This indicates that the physical structure of the disk-envelope interface differs from the rest of the envelope, highlighting the drastic impact that the disk has on the envelope and temperature structure. The results presented here show that DCO+ is an excellent cold temperature tracer.

preprint2015arXiv

Planck cold clumps in the $λ$ Orionis complex: I. Discovery of an extremely young Class 0 protostellar object and a proto-brown dwarf candidate in a bright rimmed clump PGCC G192.32-11.88

We are performing a series of observations with ground-based telescopes toward Planck Galactic cold clumps (PGCCs) in the $λ$ Orionis complex in order to systematically investigate the effects of stellar feedback. In the particular case of PGCC G192.32-11.88, we discovered an extremely young Class 0 protostellar object (G192N) and a proto-brown dwarf candidate (G192S). G192N and G192S are located in a gravitationally bound bright-rimmed clump. The velocity and temperature gradients seen in line emission of CO isotopologues indicate that PGCC G192.32-11.88 is externally heated and compressed. G192N probably has the lowest bolometric luminosity ($\sim0.8$ L$_{\sun}$) and accretion rate (6.3$\times10^{-7}$ M$_{\sun}$~yr$^{-1}$) when compared with other young Class 0 sources (e.g. PACS Bright Red sources (PBRs)) in the Orion complex. It has slightly larger internal luminosity ($0.21\pm0.01$ L$_{\sun}$) and outflow velocity ($\sim$14 km~s$^{-1}$) than the predictions of first hydrostatic cores (FHSCs). G192N might be among the youngest Class 0 sources, which are slightly more evolved than a FHSC. Considering its low internal luminosity ($0.08\pm0.01$ L$_{\odot}$) and accretion rate (2.8$\times10^{-8}$ M$_{\sun}$~yr$^{-1}$), G192S is an ideal proto-brown dwarf candidate. The star formation efficiency ($\sim$0.3\%-0.4\%) and core formation efficiency ($\sim$1\%) in PGCC G192.32-11.88 are significantly smaller than in other giant molecular clouds or filaments, indicating that the star formation therein is greatly suppressed due to stellar feedback.

preprint2015arXiv

Properties of the Molecular Cores of Low Luminosity Objects

We present a survey toward 16 Low Luminosity Objects (LLOs with an internal luminosity, Lint, lower than 0.2 Lsun) with N2H+ (1-0), N2H+ (3-2), N2D+ (3-2), HCO+ (3-2) and HCN (3-2) using the Arizona Radio Observatory Kitt Peak 12m Telescope and Submillimeter Telescope. Our goal is to probe the nature of these faint protostars which are believed to be either very low mass or extremely young protostars. We find that the N2D+/N2H+ column density ratios of LLOs are similar to those of typical starless cores and Class 0 objects. The N2D+/N2H+ column density ratios are relatively high (> 0.05) for LLOs with kinetic temperatures less than 10 K in our sample. The distribution of N2H+ (1-0) line widths spreads between that of starless cores and young Class 0 objects. If we use the line width as a dynamic evolutionary indicator, LLOs are likely young Class 0 protostellar sources. We further use the optically thick tracers, HCO+ (3-2) and HCN (3-2), to probe the infall signatures of our targets. We derive the asymmetry parameters from both lines and estimate the infall velocities by fitting the HCO+ (3-2) spectra with two-layer models. As a result, we identify eight infall candidates based on the infall velocities and seven candidates have infall signatures supported by asymmetry parameters from at least one of HCO+ (3-2) and HCN (3-2).

preprint2015arXiv

The Inner Debris Structure in the Fomalhaut Planetary System

Fomalhaut plays an important role in the study of debris disks and small bodies in other planetary systems. The proximity and luminosity of the star make key features of its debris, like the water ice-line, accessible. Here we present ALMA cycle 1, 870 μm (345 GHz) observations targeted at the inner part of the Fomalhaut system with a synthesized beam of 0.45"x0.37" (~3 AU linear resolution at the distance of Fomalhaut) and a rms of 26 μJy/beam. The high angular resolution and sensitivity of the ALMA data enable us to place strong constraints on the nature of the warm excess revealed by Spitzer and Herschel observations. We detect a point source at the star position with a total flux consistent with thermal emission from the stellar photosphere. No structures that are brighter than 3σ are detected in the central 15 AU x 15 AU region. Modeling the spectral energy distribution using parameters expected for a dust-producing planetesimal belt indicates a radial location in the range ~8-15 AU. This is consistent with the location where ice sublimates in Fomalhaut, i.e., an asteroid-belt analog. The 3σ upper limit for such a belt is <1.3 mJy at 870 μm. We also interpret the 2 and 8-13 μm interferometric measurements to reveal the structure in the inner 10 AU region as dust naturally connected to this proposed asteroid belt by Poynting-Robertson drag, dust sublimation, and magnetically trapped nano grains.

preprint2015arXiv

The Spitzer c2d Survey of Large, Nearby, Interstellar Clouds. XII. The Perseus YSO Population as Observed with IRAC and MIPS

The Spitzer Space Telescope mapped the Perseus molecular cloud complex with IRAC and MIPS as part of the c2d Spitzer Legacy project. This paper combines the observations from both instruments giving an overview of low-mass star formation across Perseus from 3.6 to 70 micron. We provide an updated list of young stellar objects with new classifications and source fluxes from previous works, identifying 369 YSOs in Perseus with the Spitzer dataset. By synthesizing the IRAC and MIPS maps of Perseus and building on the work of previous papers in this series (Jorgensen et al. 2006, Rebull et al. 2007), we present a current census of star formation across the cloud and within smaller regions. 67% of the YSOs are associated with the young clusters NGC 1333 and IC 348. The majority of the star formation activity in Perseus occurs in the regions around the clusters, to the eastern and western ends of the cloud complex. The middle of the cloud is nearly empty of YSOs despite containing regions of high visual extinction. The western half of Perseus contains three-quarters of the total number of embedded YSOs (Class 0+I and Flat SED sources) in the cloud and nearly as many embedded YSOs as Class II and III sources. Class II and III greatly outnumber Class 0+I objects in eastern Perseus and IC 348. These results are consistent with previous age estimates for the clusters. Across the cloud, 56% of YSOs and 91% of the Class 0+I and Flat sources are in areas where Av > 5 mag, indicating a possible extinction threshold for star formation.

preprint2014arXiv

Magnetic field structure in the Flattened Envelope and Jet in the young protostellar system HH 211

HH 211 is a young Class 0 protostellar system, with a flattened envelope, a possible rotating disk, and a collimated jet. We have mapped it with the Submillimeter Array in 341.6 GHz continuum and SiO J=8-7 at ~ 0.6 resolution. The continuum traces the thermal dust emission in the flattened envelope and the possible disk. Linear polarization is detected in the continuum in the flattened envelope. The field lines implied from the polarization have different orientations, but they are not incompatible with current gravitational collapse models, which predict different orientation depending on the region/distance. Also, we might have detected for the first time polarized SiO line emission in the jet due to the Goldreich-Kylafis effect. Observations at higher sensitivity are needed to determine the field morphology in the jet.

preprint2014arXiv

Magnetic Fields and Massive Star Formation

Massive stars ($M > 8$ \msun) typically form in parsec-scale molecular clumps that collapse and fragment, leading to the birth of a cluster of stellar objects. We investigate the role of magnetic fields in this process through dust polarization at 870 $μ$m obtained with the Submillimeter Array (SMA). The SMA observations reveal polarization at scales of $\lsim$ 0.1 pc. The polarization pattern in these objects ranges from ordered hour-glass configurations to more chaotic distributions. By comparing the SMA data with the single dish data at parsec scales, we found that magnetic fields at dense core scales are either aligned within $40^\circ$ of or perpendicular to the parsec-scale magnetic fields. This finding indicates that magnetic fields play an important role during the collapse and fragmentation of massive molecular clumps and the formation of dense cores. We further compare magnetic fields in dense cores with the major axis of molecular outflows. Despite a limited number of outflows, we found that the outflow axis appears to be randomly oriented with respect to the magnetic field in the core. This result suggests that at the scale of accretion disks ($\lsim 10^3$ AU), angular momentum and dynamic interactions possibly due to close binary or multiple systems dominate over magnetic fields. With this unprecedentedly large sample massive clumps, we argue on a statistical basis that magnetic fields play an important role during the formation of dense cores at spatial scale of 0.01 - 0.1 pc in the context of massive star and cluster star formation.

preprint2014arXiv

The Importance of the Magnetic Field from an SMA-CSO-Combined Sample of Star-Forming Regions

Submillimeter dust polarization measurements of a sample of 50 star-forming regions, observed with the SMA and the CSO covering pc-scale clouds to mpc-scale cores, are analyzed in order to quantify the magnetic field importance. The magnetic field misalignment $δ$ -- the local angle between magnetic field and dust emission gradient -- is found to be a prime observable, revealing distinct distributions for sources where the magnetic field is preferentially aligned with or perpendicular to the source minor axis. Source-averaged misalignment angles $\langle|δ|\rangle$ fall into systematically different ranges, reflecting the different source-magnetic field configurations. Possible bimodal $\langle|δ|\rangle$-distributions are found for the separate SMA and CSO samples. Combining both samples broadens the distribution with a wide maximum peak at small $\langle|δ|\rangle$-values. Assuming the 50 sources to be representative, the prevailing source-magnetic field configuration is one that statistically prefers small magnetic field misalignments $|δ|$. When interpreting $|δ|$ together with an MHD force equation, as developed in the framework of the polarization-intensity gradient method, a sample-based scaling fits the magnetic field tension-to-gravity force ratio $\langleΣ_B\rangle$ versus $\langle|δ|\rangle$ with $\langleΣ_B\rangle = 0.116 \cdot \exp(0.047\cdot \langle|δ|\rangle)\pm 0.20$ (mean error), providing a way to estimate the relative importance of the magnetic field, only based on measurable field misalignments $|δ|$. The force ratio $Σ_B$ discriminates systems that are collapsible on average ($\langle Σ_B\rangle <1$) from other molecular clouds where the magnetic field still provides enough resistance against gravitational collapse ($\langle Σ_B\rangle >1$) (abridged).

preprint2013arXiv

A Keplerian disk around a Class 0 source: ALMA observations of VLA1623A

Context: Rotationally supported disks are critical in the star formation process. The questions of when do they form and what factors influence or hinder their formation have been studied but are largely unanswered. Observations of early stage YSOs are needed to probe disk formation. Aims: VLA1623 is a triple non-coeval protostellar system, with a weak magnetic field perpendicular to the outflow, whose Class 0 component, VLA1623A, shows a disk-like structure in continuum with signatures of rotation in line emission. We aim to determine whether this structure is in part or in whole a rotationally supported disk, i.e. a Keplerian disk, and what are its characteristics. Methods: ALMA Cycle 0 Early Science 1.3 mm continuum and C$^{18}$O (2-1) observations in the extended configuration are presented here and used to perform an analysis of the disk-like structure using PV diagrams and thin disk modelling with the addition of foreground absorption. Results: The PV diagrams of the C$^{18}$O line emission suggest the presence of a rotationally supported component with a radius of at least 50 AU. Kinematical modelling of the line emission shows that the disk out to 180 AU is actually rotationally supported, with the rotation being well described by Keplerian rotation out to at least 150 AU, and the central source mass to be $\sim$0.2 M$_{sun}$ for an inclination of 55$^{\circ}$. Pure infall and conserved angular momentum rotation models are excluded. Conclusions: VLA1623A, a very young Class 0 source, presents a disk with an outer radius $R_{\rm out}$ = 180 AU with a Keplerian velocity structure out to at least 150 AU. The weak magnetic fields and recent fragmentation in this region of rho Ophiuchus may have played a lead role in the formation of the disk.

preprint2013arXiv

Detection of a Magnetized Disk around a Very Young Protostar

We present subarcsecond resolution polarimetric observations of the 878 mum thermal dust continuum emission obtained with the Submillimeter Array (SMA) towards the IRAS 16293-2422 protostellar binary system. We report the detection of linearly polarized dust emission arising from the circumstellar disk associated with the IRAS 16293-2422 B protostar. The fractional polarization of ~1.4$% is only slightly lower than that expected from theoretical calculations in such disks. The magnetic field structure on the plane of the sky derived from the dust polarization suggests a complex magnetic field geometry in the disk, possibly associated with a rotating disk that is wrapping the field lines as expected from the simulations. The polarization around IRAS 16293-2422 A at subarcsecond angular resolution is only marginally detected.

preprint2013arXiv

Disentangling the entangled: Observations and analysis of the triple non-coeval protostellar system VLA1623

Commonplace at every evolutionary stage, Multiple Protostellar Systems (MPSs) are thought to be formed through fragmentation, but it is unclear when and how. The youngest MPSs, which have not yet undergone much evolution, provide important constraints to this question. It is then of interest to disentangle early stage MPSs. In this letter we present the results of our work on VLA1623 using our observations and archival data from the Submillimeter Array (SMA). Our continuum and line observations trace VLA1623's components, outflow and envelope, revealing unexpected characteristics. We construct the SED for each component using the results of our work and data from literature, as well as derive physical parameters from continuum and perform a simple kinematical analysis of the circumstellar material. Our results show VLA1623 to be a triple non-coeval system composed of VLA1623A, B & W, with each source driving its own outflow and unevenly distributed circumstellar material. From the SED, physical parameters and IR emission we conclude that VLA1623A & W are Class 0 and I protostars, respectively, and together drive the bulk of the observed outflow. Furthermore, we find two surprising results, first the presence of a rotating disk-like structure about VLA1623A with indications of pure Keplerian rotation, which, if real, would make it one of the first evidence of Keplerian disk structures around Class 0 protostars. Second, we find VLA1623B to be a bonafide extremely young protostellar object between the starless core and Class 0 stages.

preprint2013arXiv

Populations of Young Stellar Objects in Nearby Molecular Clouds

We develop a new method to identify YSOs from star-forming regions using the photometry data from Spitzer's c2d Legacy Project. The aim is to obtain YSO lists as complete as possible for studying the statistical properties, such as Star Formation Rate (SFR) and lifetimes of YSOs in different evolutionary stages. The largest obstacle for identifying YSOs comes from background galaxies with similar SEDs to YSOs. Traditionally, selected color-color and color-magnitude criteria are used to separate YSOs and galaxies. However, since there is no obvious boundary between YSOs and galaxies in Color-Color Diagrams (CCDs) and Color-Magnitude Diagrams (CMDs), those criteria may exclude faint YSOs near the boundary. In this paper, we separate the YSOs and galaxies in multi-dimensional (Multi-D) magnitude space, which is equivalent to using all variations of CMDs simultaneously. Comparing sources from molecular clouds to Spitzer's SWIRE data, which have negligible amount of YSOs, we can naturally identify YSO candidates locating outside of the galaxy populated regions in the Multi-D space. In the five c2d-surveyed clouds, we select 322 new YSO candidates (YSOc), miss/exclude 33 YSOc compared to Evans et al. 2009 and result in 1313 YSOc in total. As a result, SFR increases 28% correspondingly, but the lifetimes of YSOs in different evolutionary stages remain unchanged. Comparing to theories Krumholz & McKee 2005, our derived SFR suggests that star formation in large scale is dominated by supersonic turbulence rather than magnetic fields. Furthermore, we identify 7 new Very Low Luminosity Objects.

preprint2012arXiv

Emission Mechanism of "Green Fuzzies" in High-mass Star Forming Regions

The Infrared Array Camera (IRAC) on the Spitzer Space Telescope has revealed that a number of high-mass protostars are associated with extended mid-infrared emission, particularly prominent at 4.5-micron. These are called "Green Fuzzy" emission or "Extended Green Objects". We present color analysis of this emission toward six nearby (d=2-3 kpc) well-studied high-mass protostars and three candidate high-mass protostars identified with the Spitzer GLIMPSE survey. In our color-color diagrams most of the sources show a positive correlation between the [3.6]-[4.5] and [3.5]-[5.8] colors along the extinction vector in all or part of the region. We compare the colors with those of scattered continuum associated with the low-mass protostar L 1527, modeled scattered continuum in cavities, shocked emission associated with low-mass protostars, modeled H2 emission for thermal and fluorescent cases, and modeled PAH emission. Of the emission mechanisms discussed above, scattered continuum provides the simplest explanation for the observed linear correlation. In this case, the color variation within each object is attributed to different foreground extinctions at different positions. Alternative possible emission mechanisms to explain this correlation may be a combination of thermal and fluorescent H2 emission in shocks, and a combination of scattered continuum and thermal H2 emission, but detailed models or spectroscopic follow-up are required to further investigate this possibility. Our color-color diagrams also show possible contributions from PAHs in two objects. However, none of our sample show clear evidence for PAH emission directly associated with the high-mass protostars, several of which should be associated with ionizing radiation. This suggests that those protostars are heavily embedded even at mid-infrared wavelengths.

preprint2010arXiv

The magnetic field in the NGC 2024 FIR 5 dense core

We used the Submillimeter Array (SMA) to observe the thermal polarized dust emission from the protostellar source NGC 2024 FIR 5. The polarized emission outlines a partial hourglass morphology for the plane-of-sky component of the core magnetic field. Our data are consistent with previous BIMA maps, and the overall magnetic field geometries obtained with both instruments are similar. We resolve the main core into two components, FIR 5A and FIR 5B. A possible explanation for the asymmetrical field lies in depolarization effects due to the lack of internal heating from FIR 5B source, which may be in a prestellar evolutionary state. The field strength was estimated to be 2.2 mG, in agreement with previous BIMA data. We discuss the influence of a nearby H{\sc ii} region over the field lines at scales of $\sim 0.01$ pc. Although the hot component is probably compressing the molecular gas where the dust core is embedded, it is unlikely that the radiation pressure exceeds the magnetic tension. Finally, a complex outflow morphology is observed in CO (3 $\rightarrow$ 2) maps. Unlike previous maps, several features associated with dust condensations other than FIR 5 are detected.

preprint2010arXiv

Tracing Turbulent Ambipolar Diffusion in Molecular Clouds

Though flux freezing is a good approximation frequently assumed for molecular clouds, ambipolar diffusion (AD) is inevitable at certain scales. The scale at which AD sets in can be a crucial parameter for turbulence and the star formation process. However, both observation and simulation of AD are very challenging and our knowledge of it is very limited. We proposed earlier (Li and Houde 2008) that the difference between ion and neutral velocity spectra is a signature of turbulent AD and can be used to estimate the AD scales and magnetic field strength. Here we present observational evidence showing that this difference between the velocity dispersions from coexistent ions and neutrals is indeed correlated with magnetic field strength.

preprint2009arXiv

IRAS 16293: A "Magnetic" Tale of Two Cores

We present polarization observations of the dust continuum emission from IRAS 16293 which is a Class 0 protostar and is known to have at least two components, source A and B. These measurements were conducted by the Submillimeter Array (SMA) at a frequency of ~341 GHz and with high angular resolution (~2-3 arcseconds). We find that the large scale global direction of the field, which is perpendicular to the observed polarization, appears to be along the dust ridge where the emission peaks. On smaller scales we find that the field structure is significantly different for the two components of the binary. The first component, source A, shows a magnetic field structure which is "hourglass" shaped as predicted from theoretical models of low mass star formation in the presence of strong magnetic fields. However, the other component, source B, shows a relatively ordered magnetic field with no evidence of any deformation. There is an observed decrease in polarization towards the center. Our calculations show that in IRAS 16293 the magnetic energy is stronger than the turbulent energy but is approximately similar to the centrifugal energy. Our results provide additional evidence to show that the two protostars appear to be in different stages during their evolution.

preprint2009arXiv

The Evolution of Density Structure of Starless and Protostellar Cores

We present a near-infrared extinction study of nine dense cores at evolutionary stages between starless to Class I. Our results show that the density structure of all but one observed cores can be modeled with a single power law rho \propto r^p between ~ 0.2R-R of the cores. The starless cores in our sample show two different types of density structures, one follows p ~ -1.0 and the other follows p ~ -2.5, while the protostellar cores all have p ~ -2.5. The similarity between the prestellar cores with p ~ -2.5 and protostellar cores implies that those prestellar cores could be evolving towards the protostellar stage. The slope of p ~ -2.5 is steeper than that of an singular isothermal sphere, which may be interpreted with the evolutionary model of cores with finite mass.

preprint2008arXiv

Evolution of Magnetic Fields in High Mass Star Formation: SMA dust polarization image of the UCHII region G5.89-0.39

We report high angular resolution (3") Submillimeter Array (SMA) observations of the molecular cloud associated with the Ultra-Compact HII region G5.89-0.39. Imaged dust continuum emission at 870 micron reveals significant linear polarization. The position angles (PAs) of the polarization vary enormously but smoothly in a region of 2x10^4 AU. Based on the distribution of the PAs and the associated structures, the polarized emission can be separated roughly into two components. The component "x" is associated with a well defined dust ridge at 870 micron, and is likely tracing a compressed B field. The component "o" is located at the periphery of the dust ridge and is probably from the original B field associated with a pre-existing extended structure. The global B field morphology in G5.89, as inferred from the PAs, is clearly disturbed by the expansion of the HII region and the molecular outflows. Using the Chandrasekhar-Fermi method, we estimate from the smoothness of the field structures that the B field strength in the plane of sky can be no more than 2-3 mG. We then compare the energy densities in the radiation, the B field, and the mechanical motions as deduced from the C^17O 3-2 line emission. We conclude that the B field structures are already overwhelmed and dominated by the radiation, outflows, and turbulence from the newly formed massive stars.

preprint2007arXiv

The Spitzer c2d Survey of Large, Nearby, Interstellar Clouds VIII. Serpens Observed with MIPS

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.

preprint2007arXiv

The Spitzer c2d Survey of Weak-line T Tauri Stars II: New Constraints on the Timescale for Planet Building

One of the central goals of the Spitzer Legacy Project ``From Molecular Cores to Planet-forming Disks'' (c2d) is to determine the frequency of remnant circumstellar disks around weak-line T Tauri stars (wTTs) and to study the properties and evolutionary status of these disks. Here we present a census of disks for a sample of over 230 spectroscopically identified wTTs located in the c2d IRAC (3.6, 4.5, 4.8, and 8.0 um) and MIPS (24 um) maps of the Ophiuchus, Lupus, and Perseus Molecular Clouds. We find that ~20% of the wTTs in a magnitude limited subsample have noticeable IR-excesses at IRAC wavelengths indicating the presence of a circumstellar disk. The disk frequencies we find in these 3 regions are ~3-6 times larger than that recently found for a sample of 83 relatively isolated wTTs located, for the most part, outside the highest extinction regions covered by the c2d IRAC and MIPS maps. The disk fractions we find are more consistent with those obtained in recent Spitzer studies of wTTs in young clusters such as IC 348 and Tr 37. From their location in the H-R diagram, we find that, in our sample, the wTTs with excesses are among the younger part of the age distribution. Still, up to ~50% of the apparently youngest stars in the sample show no evidence of IR excess, suggesting that the circumstellar disks of a sizable fraction of pre-main-sequence stars dissipate in a timescale of ~1 Myr. We also find that none of the stars in our sample apparently older than ~10 Myrs have detectable circumstellar disks at wavelengths < 24 um. Also, we find that the wTTs disks in our sample exhibit a wide range of properties (SED morphology, inner radius, L_DISK/L*, etc) which bridge the gaps observed between the cTTs and the debris disk regimes.

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 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).