Source author record

Philip C. Myers

Philip C. Myers 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

49works
4topics
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

49 published item(s)

preprint2022arXiv

The Magnetic Field in the Milky Way Filamentary Bone G47

Star formation primarily occurs in filaments where magnetic fields are expected to be dynamically important. The largest and densest filaments trace spiral structure within galaxies. Over a dozen of these dense ($\sim$10$^4$\,cm$^{-3}$) and long ($>$10\,pc) filaments have been found within the Milky Way, and they are often referred to as "bones." Until now, none of these bones have had their magnetic field resolved and mapped in their entirety. We introduce the SOFIA legacy project FIELDMAPS which has begun mapping $\sim$10 of these Milky Way bones using the HAWC+ instrument at 214\,$μ$m and 18$\farcs$2 resolution. Here we present a first result from this survey on the $\sim$60\,pc long bone G47. Contrary to some studies of dense filaments in the Galactic plane, we find that the magnetic field is often not perpendicular to the spine (i.e., the center-line of the bone). Fields tend to be perpendicular in the densest areas of active star formation and more parallel or random in other areas. The average field is neither parallel or perpendicular to the Galactic plane nor the bone. The magnetic field strengths along the spine typically vary from $\sim$20 to $\sim$100\,$μ$G. Magnetic fields tend to be strong enough to suppress collapse along much of the bone, but for areas that are most active in star formation, the fields are notably less able to resist gravitational collapse.

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.

preprint2022arXiv

Virial Clumps in Central Molecular Zone Clouds

CMZoom survey observations with the Submillimeter Array are analyzed to describe the virial equilibrium (VE) and star-forming potential of 755 clumps in 22 clouds in the Central Molecular Zone (CMZ) of the Milky Way. In each cloud, nearly all clumps follow the column-density-mass trend N~M^s, where s = 0.38 +- 0.03 is near the pressure-bounded limit s=1/3. This trend is expected when gravitationally unbound clumps in VE have similar velocity dispersion and external pressure. Nine of these clouds also harbor one or two distinctly more massive clumps. These properties allow a VE model of bound and unbound clumps in each cloud, where the most massive clump has the VE critical mass. These models indicate that 213 clumps have velocity dispersion 1-2 km s^(-1), mean external pressure 0.5-4 x 10^8 cm^(-3) K, bound clump fraction 0.06, and typical virial parameter alpha=4-15. These mostly unbound clumps may be in VE with their turbulent cloud pressure, possibly driven by inflow from the galactic bar. In contrast, most Sgr B2 clumps are bound according to their associated sources and N-M trends. When the CMZ clumps are combined into mass distributions, their typical power-law slope is analyzed with a stopped accretion model. It also indicates that most clumps are unbound and cannot grow significantly, due to their similar time scales of accretion and dispersal, ~0.2 Myr. Thus, virial and dynamical analyses of the most extensive clump census available indicate that star formation in the CMZ may be suppressed by a significant deficit of gravitationally bound clumps.

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

CS Depletion in Prestellar Cores

The CS molecule is known to be absorbed onto dust in the cold and dense conditions, causing it to get significantly depleted in the central region of cores. This study is aimed to investigate the depletion of the CS molecule using the optically thin C$^{34}$S molecular line observations. We mapped five prestellar cores, L1544, L1552, L1689B, L694-2, and L1197 using two molecular lines, C$^{34}$S $(J=2-1)$ and N$_2$H$^+$ $(J=1-0)$ with the NRO 45-m telescope, doubling the number of cores where the CS depletion was probed using C$^{34}$S. In most of our targets, the distribution of C$^{34}$S emission shows features that suggest that the CS molecule is generally depleted in the center of the prestellar cores. The radial profile of the CS abundance with respect to H$_2$ directly measured from the CS emission and the Herschel dust emission indicates that the CS molecule is depleted by a factor of $\sim$3 toward the central regions of the cores with respect to their outer regions. The degree of the depletion is found to be even more enhanced by an order of magnitude when the contaminating effect introduced by the presence of CS molecules in the surrounding envelope that lie along the line-of-sight is removed. Except for L1197 which is classified as relatively the least evolved core in our targets based on its observed physical parameters, we found that the remaining four prestellar cores are suffering from significant CS depletion at their central region regardless of the relative difference in their evolutionary status.

preprint2020arXiv

Efficient methanol production on the dark side of a prestellar core

We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particularly efficient. Methanol and sulphur monoxide are most abundant in a narrow zone that follows the eastern side of the core. This side is sheltered from the stronger external radiation field coming from the west. We show that photodissociation on the illuminated side can give rise to an asymmetric methanol distribution, but that the stark contrast observed in H-MM1 is hard to explain without assuming enhanced desorption on the shaded side. The region of the brightest emission has a wavy structure that rolls up at one end. This is the signature of Kelvin-Helmholtz instability occurring in sheared flows. We suggest that in this zone, methanol and sulphur are released as a result of grain-grain collisions induced by shear vorticity.

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

Relative Alignment between Dense Molecular Cores and Ambient Magnetic Field: The Synergy of Numerical Models and Observations

The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular clouds. We find that the results on relative alignment between cores and background magnetic field in synthetic observations slightly disagree with those measured in fully 3D simulation data, which is partly because cores identified in projected 2D maps tend to coexist within filamentary structures, while 3D cores are generally more rounded. In addition, we examine the progression of magnetic field from pc- to core-scale in the simulation, which is consistent with the anisotropic core formation model that gas preferably flow along the magnetic field toward dense cores. When comparing the observed cores identified from the GBT Ammonia Survey (GAS) and Planck polarization-inferred magnetic field orientations, we find that the relative core-field alignment has a regional dependence among different clouds. More specifically, we find that dense cores in the Taurus molecular cloud tend to align perpendicular to the background magnetic field, while those in Perseus and Ophiuchus tend to have random (Perseus) or slightly parallel (Ophiuchus) orientations with respect to the field. We argue that this feature of relative core-field orientation could be used to probe the relative significance of the magnetic field within the cloud.

preprint2020arXiv

Star-Gas Surface Density Correlations in Twelve Nearby Molecular Clouds I: Data Collection and Star-Sampled Analysis

We explore the relation between the stellar mass surface density and the mass surface density of molecular hydrogen gas in twelve nearby molecular clouds that are located at $<$1.5 kpc distance. The sample clouds span an order of magnitude range in mass, size, and star formation rates. We use thermal dust emission from $Herschel$ maps to probe the gas surface density and the young stellar objects from the most recent $Spitzer$ Extended Solar Neighborhood Archive (SESNA) catalog to probe the stellar surface density. Using a star-sampled nearest neighbor technique to probe the star-gas surface density correlations at the scale of a few parsecs, we find that the stellar mass surface density varies as a power-law of the gas mass surface density, with a power-law index of $\sim$2 in all the clouds. The consistent power-law index implies that star formation efficiency is directly correlated with gas column density, and no gas column density threshold for star formation is observed. We compare the observed correlations with the predictions from an analytical model of thermal fragmentation, and with the synthetic observations of a recent hydrodynamic simulation of a turbulent star-forming molecular cloud. We find that the observed correlations are consistent for some clouds with the thermal fragmentation model and can be reproduced using the hydrodynamic simulations.

preprint2020arXiv

Ubiquitous $\rm NH_3$ supersonic component in L1688 coherent cores

Context : Star formation takes place in cold dense cores in molecular clouds. Earlier observations have found that dense cores exhibit subsonic non-thermal velocity dispersions. In contrast, CO observations show that the ambient large-scale cloud is warmer and has supersonic velocity dispersions. Aims : We aim to study the ammonia ($\rm NH_3$) molecular line profiles with exquisite sensitivity towards the coherent cores in L1688 in order to study their kinematical properties in unprecedented detail. Methods : We used $\rm NH_3$ (1,1) and (2,2) data from the first data release (DR1) in the Green Bank Ammonia Survey (GAS). We first smoothed the data to a larger beam of 1' to obtain substantially more extended maps of velocity dispersion and kinetic temperature, compared to the DR1 maps. We then identified the coherent cores in the cloud and analysed the averaged line profiles towards the cores. Results : For the first time, we detected a faint (mean $\rm NH_3$(1,1) peak brightness $<$0.25 K in $T_{MB}$), supersonic component towards all the coherent cores in L1688. We fitted two components, one broad and one narrow, and derived the kinetic temperature and velocity dispersion of each component. The broad components towards all cores have supersonic linewidths ($\mathcal{M}_S \ge 1$). This component biases the estimate of the narrow dense core component's velocity dispersion by $\approx$28% and the kinetic temperature by $\approx$10%, on average, as compared to the results from single-component fits. Conclusions : Neglecting this ubiquitous presence of a broad component towards all coherent cores causes the typical single-component fit to overestimate the temperature and velocity dispersion. This affects the derived detailed physical structure and stability of the cores estimated from $\rm NH_3$ observations.

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

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

Characteristic structure of star-forming clouds

This paper gives a new way to diagnose the star-forming potential of a molecular cloud region from the probability density function of its column density (N-pdf). It gives expressions for the column density and mass profiles of a symmetric filament having the same N-pdf as a filamentary region. The central concentration of this characteristic filament can distinguish regions and can quantify their fertility for star formation. Profiles are calculated for N-pdfs which are pure lognormal, pure power law, or a combination. In relation to models of singular polytropic cylinders, characteristic filaments can be unbound, bound, or collapsing depending on their central concentration. Such filamentary models of the dynamical state of N-pdf gas are more relevant to star-forming regions than are models of spherical collapse. The star formation fertility of a bound or collapsing filament is quantified by its mean mass accretion rate when in radial free fall. For a given mass per length, the fertility increases with the filament mean column density and with its initial concentration. In selected regions the fertility of their characteristic filaments increases with the level of star formation.

preprint2015arXiv

Detection of Infall in the Protostar B335 with ALMA

Observations of the isolated globule B335 with ALMA have yielded absorption features against the continuum that are redshifted from the systemic velocity in both HCN and HCO$^+$ lines. These features provide unambiguous evidence for infall toward a central luminosity source. Previously developed models of inside-out collapse can match the observed line profiles of HCN and HCO$^+$ averaged over the central 50 AU. At the new distance of 100 pc, the inferred infall radius is 0.012 pc, the mass infall rate is $3 \times 10^{-6}$ solar masses per year, the age is 50,000 years, and the accumulated mass in the central zone is 0.15 solar masses, most of which must be in the star or in parts of a disk that are opaque at 0.8 mm. The continuum detection indicates an optically thin mass (gas and dust) of only $7.5\times 10^{-4}$ solar masses in the central region, consistent with only a very small disk mass.

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

A Spitzer View of Mon OB1 East/NGC 2264

We present Spitzer 3.6, 4.5, 5.8, 8.0, and 24 micron images of the Mon OB1 East giant molecular cloud, which contains the young star forming region NGC 2264, as well as more extended star formation. With Spitzer data and 2MASS photometry, we identify and classify young stellar objects (YSOs) with dusty circumstellar disks and/or envelopes in Mon OB1 East by their infrared-excess emission and study their distribution with respect to cloud material. We find a correlation between the local surface density of YSOs and column density of molecular gas as traced by dust extinction that is roughly described as a power law in these quantities. NGC 2264 follows a power law index of ~2.7, exhibiting a large YSO surface density for a given gas column density. Outside of NGC 2264 where the surface density of YSOs is lower, the power law is shallower and the region exhibits a larger gas column density for a YSO surface density, suggesting the star formation is more recent. In order to measure the fraction of cloud members with circumstellar disks/envelopes, we estimate the number of diskless pre-main sequence stars by statistical removal of background star detections. We find that the disk fraction of the NGC 2264 region is 45%, while the surrounding more distributed regions show a disk fraction of 19%. This may be explained by the presence an older, more dispersed population of stars. In total, the Spitzer observations provide evidence for heterogenous, non-coeval star formation throughout the Mon OB1 cloud.

preprint2014arXiv

IC348-SMM2E: a Class 0 proto-brown dwarf candidate forming as a scaled-down version of low-mass stars

We report on Submillimeter Array observations of the 870 micron continuum and CO(3-2), 13CO(2-1) and C18O(2-1) line emission of a faint object, SMM2E, near the driving source of the HH797 outflow in the IC348 cluster. The continuum emission shows an unresolved source for which we estimate a mass of gas and dust of 30 Mjup, and the CO(3-2) line reveals a compact bipolar outflow centred on SMM2E, and barely seen also in 13CO(2-1). In addition, C18O(2-1) emission reveals hints of a possible rotating envelope/disk perpendicular to the outflow, for which we infer a dynamical mass of ~16 Mjup. In order to further constrain the accreted mass of the object, we gathered data from Spitzer, Herschel, and new and archive submillimetre observations, and built the Spectral Energy Distribution (SED). The SED can be fitted with one single modified black-body from 70 micron down to 2.1 cm, using a dust temperature of ~24 K, a dust emissivity index of 0.8, and an envelope mass of ~35 Mjup. The bolometric luminosity is 0.10 Lsun, and the bolometric temperature is 35 K. Thus, SMM2E is comparable to the known Class 0 objects in the stellar domain. An estimate of the final mass indicates that SMM2E will most likely remain substellar, and the SMM2E outflow force matches the trend with luminosity known for young stellar objects. Thus, SMM2E constitutes an excellent example of a Class 0 proto-brown dwarf candidate which forms as a scaled-down version of low-mass stars. Finally, SMM2E seems to be part of a wide (~2400 AU) multiple system of Class 0 sources.

preprint2014arXiv

The Evolution of Protostars: Insights from Ten Years of Infrared Surveys with Spitzer and Herschel

Stars form from the gravitational collapse of dense molecular cloud cores. In the protostellar phase, mass accretes from the core onto a protostar, likely through an accretion disk, and it is during this phase that the initial masses of stars and the initial conditions for planet formation are set. Over the past decade, new observational capabilities provided by the Spitzer Space Telescope and Herschel Space Observatory have enabled wide-field surveys of entire star-forming clouds with unprecedented sensitivity, resolution, and infrared wavelength coverage. We review resulting advances in the field, focusing both on the observations themselves and the constraints they place on theoretical models of star formation and protostellar evolution. We also emphasize open questions and outline new directions needed to further advance the field.

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

Origin of the dense core mass function in contracting filaments

Mass functions of starless dense cores (CMFs) may arise from contraction and dispersal of core-forming filaments. In an illustrative model, a filament contracts radially by self-gravity, increasing the mass of its cores. During this contraction, FUV photoevaporation and ablation by shocks and winds disperse filament gas and limit core growth. The stopping times of core growth are described by a waiting-time distribution. The initial filament column density profile and the resulting CMF each match recent Herschel observations in detail. Then low-mass cores have short growth ages and arise from the innermost filament gas, while massive cores have long growth ages and draw from more extended filament gas. The model fits the initial density profile and CMF best for mean core density 2 10^4 cm^-3 and filament dispersal time scale 0.5 Myr. Then the typical core mass, radius, mean column density, and contraction speed are respectively 0.8 solar masses, 0.06 pc, 6 10^21 cm^-2, and 0.07 km s^-1, also in accord with observed values.

preprint2013arXiv

Protostar mass functions in young clusters

In an improved model of protostar mass functions (PMFs), protostars gain mass from isothermal cores in turbulent clumps. Their mass accretion rate is similar to Shu accretion at low mass, and to reduced Bondi accretion at high mass. Accretion durations follow a simple expression in which higher-mass protostars accrete for longer times. These times are set by ejections, stellar feedback, and gravitational competition, which terminate accretion and reduce its efficiency. The mass scale is the mass of a critically stable isothermal core. In steady state, the PMF approaches a power law at high mass due to competition between clump accretion and accretion stopping. The power law exponent is the ratio of the time scales of accretion and accretion stopping. The luminosity function (PLF) peaks near 1 L_Sun, due to inefficient accretion of core gas. Models fit observed PLFs in four large embedded clusters. These indicate that their underlying PMFs may be top-heavy compared to the IMF, depending on the model of protostar radius.

preprint2012arXiv

Mass and luminosity evolution of young stellar objects

A model of protostar mass and luminosity evolution in clusters gives new estimates of cluster age, protostar birthrate, accretion rate and mean accretion time. The model assumes constant protostar birthrate, core-clump accretion, and equally likely accretion stopping. Its parameters are set to reproduce the initial mass function, and to match protostar luminosity distributions in nearby star-forming regions. It obtains cluster ages and birthrates from the observed numbers of protostars and pre-main sequence (PMS) stars, and from the modal value of the protostar luminosity. In 31 embedded clusters and complexes the global cluster age is 1-3 Myr, matching available estimates based on optical spectroscopy and evolutionary tracks. This method of age estimation is simpler than optical spectroscopy, and is more useful for young embedded clusters where optical spectrocopy is not possible. In the youngest clusters, the protostar fraction decreases outward from the densest gas, indicating that the local star-forming age increases outward from a few 0.1 Myr in small protostar-dominated zones to a few Myr in large PMS-dominated zones.

preprint2011arXiv

Filamentary condensations in a young cluster

New models are presented for star-forming condensations in clusters. In each model, the condensation mass increases linearly with radius on small scales, and more rapidly on large scales, as in "thermal-nonthermal" models. Spherical condensations with this structure form protostars which match the IMF if their infall is subject to equally likely stopping. However such spherical models do not match the filamentary nature of cluster gas, and they are too extended to form protostars having high mass and short spacing. Two hybrid models are presented, which are spherical on small scales and filamentary on large scales. In and around clusters, cores embedded in linear filaments match the elongation of cluster gas, and the central concentration of low-mass stars. In cluster centers, condensations require a low volume filling factor to produce massive stars with short spacing. These may have stellate shape, where cores are nodes of filamentary networks, as seen in some simulations of colliding flows and of collapsing turbulent clumps. A dense configuration of such stellate condensations may be indistinguishable from a clump forming multiple protostars via filamentary flow paths.

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

Internal Motions in Starless Dense Cores

This paper discusses the statistics of internal motions in starless dense cores and the relation of these motions to core density and evolution. Four spectral lines from three molecular species are analyzed from single-pointing and mapped observations of several tens of starless cores. Blue asymmetric profiles are dominant, indicating that inward motions are prevalent in sufficiently dense starless cores. These blue profiles are found to be more abundant, and their asymmetry is bluer, at core positions with stronger $\rm N_2H^+$ line emission or higher column density. Thirty three starless cores are classified into four types according to the blue and red shifts of the lines in their molecular line maps. Among these cores, contracting motions dominate: 19 are classified as contracting, 3 as oscillating, 3 as expanding, and 8 as static. Contracting cores have inward motions all over the core with predominance of those motions near the region of peak density. Cores with the bluest asymmetry tend to have greater column density than other cores and all five cores with peak column density $> \rm 6\times 10^{21} cm^{-2}$ are found to be contracting. This suggests that starless cores are likely to have contracting motions if they are sufficiently condensed. Our classification of the starless cores may indicate a sequence of core evolution in the sense that column density increases from static to contracting cores: the static cores in the earliest stage, the expanding and/or the oscillating cores in the next, and the contracting cores in the latest stage.

preprint2011arXiv

Mass distributions of stars and cores in young groups and clusters

We investigate the relation of the stellar initial mass function (IMF) and the dense core mass function (CMF), using stellar masses and positions in 14 well-studied young groups. Initial column density maps are computed by replacing each star with a model initial core having the same star formation efficiency (SFE). For each group the SFE, core model, and observational resolution are varied to produce a realistic range of initial maps. A clumpfinding algorithm parses each initial map into derived cores, derived core masses, and a derived CMF. The main result is that projected blending of initial cores causes derived cores to be too few and too massive. The number of derived cores is fewer than the number of initial cores by a mean factor 1.4 in sparse groups and 5 in crowded groups. The mass at the peak of the derived CMF exceeds the mass at the peak of the initial CMF by a mean factor 1.0 in sparse groups and 12.1 in crowded groups. These results imply that in crowded young groups and clusters, the mass distribution of observed cores may not reliably predict the mass distribution of protostars which will form in those cores.

preprint2011arXiv

Star formation in dense clusters

A model of core-clump accretion with equally likely stopping describes star formation in the dense parts of clusters, where models of isolated collapsing cores may not apply. Each core accretes at a constant rate onto its protostar, while the surrounding clump gas accretes as a power of protostar mass. Short accretion flows resemble Shu accretion, and make low-mass stars. Long flows resemble reduced Bondi accretion and make massive stars. Accretion stops due to environmental processes of dynamical ejection, gravitational competition, and gas dispersal by stellar feedback, independent of initial core structure. The model matches the field star IMF from 0.01 to more than 10 solar masses. The core accretion rate and the mean accretion duration set the peak of the IMF, independent of the local Jeans mass. Massive protostars require the longest accretion durations, up to 0.5 Myr. The maximum protostar luminosity in a cluster indicates the mass and age of its oldest protostar. The distribution of protostar luminosities matches those in active star-forming regions if protostars have a constant birthrate but not if their births are coeval. For constant birthrate, the ratio of YSOs to protostars indicates the star-forming age of a cluster, typically ~1 Myr. The protostar accretion luminosity is typically less than its steady spherical value by a factor of ~2, consistent with models of episodic disk accretion.

preprint2011arXiv

Variations in the Mass Functions of Clustered and Isolated Young Stellar Objects

We analyze high quality, complete stellar catalogs for four young (roughly 1 Myr) and nearby (within ~300 pc) star-forming regions: Taurus, Lupus3, ChaI, and IC348, which have been previously shown to have stellar groups whose properties are similar to those of larger clusters such as the ONC. We find that stars at higher stellar surface densities within a region or belonging to groups tend to have a relative excess of more massive stars, over a wide range of masses. We find statistically significant evidence for this result in Taurus and IC348 as well as the ONC. These differences correspond to having typically a ~10 - 20% higher mean mass in the more clustered environment. Stars in ChaI show no evidence for a trend with either surface density or grouped status, and there are too few stars in Lupus3 to make any definitive interpretation. Models of clustered star formation do not typically extend to sufficiently low masses or small group sizes in order for their predictions to be tested but our results suggest that this regime is important to consider.

preprint2010arXiv

Direct observation of a sharp transition to coherence in Dense Cores

We present NH3 observations of the B5 region in Perseus obtained with the Green Bank Telescope (GBT). The map covers a region large enough (~11'x14') that it contains the entire dense core observed in previous dust continuum surveys. The dense gas traced by NH3(1,1) covers a much larger area than the dust continuum features found in bolometer observations. The velocity dispersion in the central region of the core is small, presenting subsonic non-thermal motions which are independent of scale. However, it is thanks to the coverage and high sensitivity of the observations that we present the detection, for the first time, of the transition between the coherent core and the dense but more turbulent gas surrounding it. This transition is sharp, increasing the velocity dispersion by a factor of 2 in less than 0.04 pc (the 31" beam size at the distance of Perseus, ~250 pc). The change in velocity dispersion at the transition is ~3 km/s/pc. The existence of the transition provides a natural definition of dense core: the region with nearly-constant subsonic non-thermal velocity dispersion. From the analysis presented here we can not confirm nor rule out a corresponding sharp density transition.

preprint2010arXiv

Star-forming gas in young clusters

Initial conditions for star formation in clusters are estimated for protostars whose masses follow the initial mass function (IMF) from 0.05 to 10 solar masses. Star-forming infall is assumed equally likely to stop at any moment, due to gas dispersal dominated by stellar feedback. For spherical infall, the typical initial condensation must have a steep density gradient, as in low-mass cores, surrounded by a shallower gradient, as in the clumps around cores. These properties match observed column densities in cluster-forming regions when the mean infall stopping time is 0.05 Myr and the accretion efficiency is 0.5. The infall duration increases with final protostar mass, from 0.01 to 0.3 Myr, and the mass accretion rate increases from 3 to 300 x 10^(-6) solar masses/yr. The typical spherical accretion luminosity is ~5 solar luminosities, reducing the luminosity problem to a factor ~3. The initial condensation density gradient changes from steep to shallow at radius 0.04 pc, enclosing 0.9 solar masses, with mean column density 2 x 10^(22) cm^(-2), and with effective central temperature 16 K. These initial conditions are denser and warmer than those for isolated star formation.

preprint2010arXiv

The Mass-Size Relation from Clouds to Cores. II. Solar Neighborhood Clouds

We measure the mass and size of cloud fragments in several molecular clouds continuously over a wide range of spatial scales (0.05 < r / pc < 3). Based on the recently developed "dendrogram-technique", this characterizes dense cores as well as the enveloping clouds. "Larson's 3rd Law" of constant column density, m(r) = C*r^2, is not well suited to describe the derived mass-size data. Solar neighborhood clouds not forming massive stars (< 10 M_sun; Pipe Nebula, Taurus, Perseus, and Ophiuchus) obey m(r) < 870 M_sun (r / pc)^1.33 . In contrast to this, clouds forming massive stars (Orion A, G10.15$-$0.34, G11.11$-$0.12) do exceed the aforementioned relation. Thus, this limiting mass-size relation may approximate a threshold for the formation of massive stars. Across all clouds, cluster-forming cloud fragments are found to be---at given radius---more massive than fragments devoid of clusters. The cluster-bearing fragments are found to roughly obey a mass-size law m = C*r^1.27 (where the exponent is highly uncertain in any given cloud, but is certainly smaller than 1.5).

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.

preprint2010arXiv

The Spitzer c2d Survey of Weak-Line T Tauri Stars. III. The Transition from Primordial Disks to Debris Disks

We present 3.6 to 70 μm Spitzer photometry of 154 weak-line T Tauri stars (WTTS) in the Chamaeleon, Lupus, Ophiuchus and Taurus star formation regions, all of which are within 200 pc of the Sun. For a comparative study, we also include 33 classical T Tauri stars (CTTS) which are located in the same star forming regions. Spitzer sensitivities allow us to robustly detect the photosphere in the IRAC bands (3.6 to 8 μm) and the 24 μm MIPS band. In the 70 μm MIPS band, we are able to detect dust emission brighter than roughly 40 times the photosphere. These observations represent the most sensitive WTTS survey in the mid to far infrared to date, and reveal the frequency of outer disks (r = 3-50 AU) around WTTS. The 70 μm photometry for half the c2d WTTS sample (the on-cloud objects), which were not included in the earlier papers in this series, Padgett et al. (2006) and Cieza et al. (2007), are presented here for the first time. We find a disk frequency of 19% for on-cloud WTTS, but just 5% for off- cloud WTTS, similar to the value reported in the earlier works. WTTS exhibit spectral energy distributions (SEDs) that are quite diverse, spanning the range from optically thick to optically thin disks. Most disks become more tenuous than Ldisk/L* = 2 x 10^-3 in 2 Myr, and more tenuous than Ldisk/L* = 5 x 10^-4 in 4 Myr.

preprint2010arXiv

Young Stellar Groups and Their Most Massive Stars

We analyze the masses and spatial distributions of fourteen young stellar groups in Taurus, Lupus3, ChaI, and IC348. These nearby groups, which typically contain 20 to 40 members, have membership catalogs complete to ~0.02 M_sun, and are sufficiently young that their locations should be similar to where they formed. These groups show five properties seen in clusters having many more stars and much greater surface density of stars: (1) a broad range of masses, (2) a concentration of the most massive star towards the centre of the group, (3) an association of the most massive star with a high surface density of lower-mass stars, (4) a correlation of the mass of the most massive star with the total mass of the group, and (5) the distribution of a large fraction of the mass in a small fraction of the stars.

preprint2009arXiv

On the distribution of protostar masses

The distribution of protostar masses is studied for core-environment systems whose duration of infall follows a waiting-time distribution. Each core-environment system has a continuous density profile with no barrier to mass flow. The core is an isothermal sphere and the environment is a filament, a layer, or a uniform medium. The infall is terminated by gas dispersal due to outflows and turbulence. The distribution of infall durations is a declining exponential, the simplest waiting-time distribution. The resulting distribution of protostar masses closely resembles the initial mass function, provided the environment density is sufficiently high, and the distribution of initial core masses is sufficiently narrow. The high-mass tail of the mass function increases strongly with environment density and weakly with environment dimension. Isolated regions of low environment density form protostars of low mass from within the parent core. In contrast, clustered regions of high environment density form protostars of low mass from core gas, and protostars of high mass from core and environment gas.

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.

preprint2008arXiv

The Spitzer c2d Legacy Results: Star Formation Rates and Efficiencies; Evolution and Lifetimes

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

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.

preprint2007arXiv

The Spitzer c2d Survey of Large, Nearby, Insterstellar Clouds. IX. The Serpens YSO Population As Observed With IRAC and MIPS

We discuss the results from the combined IRAC and MIPS c2d Spitzer Legacy observations of the Serpens star-forming region. In particular we present a set of criteria for isolating bona fide young stellar objects, YSO's, from the extensive background contamination by extra-galactic objects. We then discuss the properties of the resulting high confidence set of YSO's. We find 235 such objects in the 0.85 deg^2 field that was covered with both IRAC and MIPS. An additional set of 51 lower confidence YSO's outside this area is identified from the MIPS data combined with 2MASS photometry. We describe two sets of results, color-color diagrams to compare our observed source properties with those of theoretical models for star/disk/envelope systems and our own modeling of the subset of our objects that appear to be star+disks. These objects exhibit a very wide range of disk properties, from many that can be fit with actively accreting disks to some with both passive disks and even possibly debris disks. We find that the luminosity function of YSO's in Serpens extends down to at least a few x .001 Lsun or lower for an assumed distance of 260 pc. The lower limit may be set by our inability to distinguish YSO's from extra-galactic sources more than by the lack of YSO's at very low luminosities. A spatial clustering analysis shows that the nominally less-evolved YSO's are more highly clustered than the later stages and that the background extra-galactic population can be fit by the same two-point correlation function as seen in other extra-galactic studies. We also present a table of matches between several previous infrared and X-ray studies of the Serpens YSO population and our Spitzer data set.

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

preprint2002arXiv

Dense Cores in Dark Clouds. XIV. N2H+(1-0) maps of dense cloud cores

We present results of an extensive mapping survey of N2H+(1-0) in about 60 low mass cloud cores already mapped in the NH3(1,1) inversion transition line. The survey has been carried out at the FCRAO antenna with an angular resolution about 1.5 times finer than the previous ammonia observations. Cores with stars typically have map sizes about a factor of two smaller for N2H+ than for NH3, indicating the presence of denser and more centrally concentrated gas compared to starless cores. Significant correlations are found between NH3 and N2H+ column densities and excitation temperatures in starless cores, but not in cores with stars, suggesting a different chemical evolution of the two species. Velocity gradients range between 0.5 and 6 km/s/pc, similar to what has been found with NH3 data. ``Local'' velocity gradients show significant variation in both magnitude and direction, suggesting the presence of complexmotions not interpretable as simple solid body rotation. Integrated intensity profiles of starless cores present a ``central flattening'' and are consistent with a spherically symmetric density law n ~ r^{-1.2} for r < ~0.03 pc and n ~ r^{-2} at larger r. Cores with stars are better modelled with single density power laws with n ~ r^{-2}. Line widths change across the core but we did not find a general trend. The deviation in line width correlates with the mean line width, suggesting that the line of sight contains ~ 10 coherence lengths. The corresponding value of the coherence length, ~ 0.01 pc, is similar to the expected cutoff wavelength for MHD waves. This similarity may account for the increased ``coherence'' of line widths on small scales. Despite of the finer angular resolution, the majority of N2H+ and NH3 maps show a similar ``simple'' structure, with single peaks and no elongation.

preprint1999arXiv

A Catalogue of Optically Selected Cores

We present a new catalogue of 406 dense cores optically selected by using the STScI Digitized Sky Survey (DSS). In this catalogue 306 cores have neither an Embedded YSO (EYSO) nor a Pre-Main-Sequence (PMS) star, 94 cores have EYSOs (1 core has both an EYSO and a PMS star), and 6 cores have PMS star only. Our sample of dense cores in the catalogue is fairly complete within a category of northern Lynds class 5, 6 clouds, and southern Hartley et al. (1986)'s class A clouds, providing a database useful for the systematic study of dense cores. Most of the cores listed in the catalogue have diameters between $0.05 - 0.36$ pc with a mean of $\sim 0.24$ pc. The sizes ($\sim 0.33$ pc in the mean) of cores with EYSOs are found to be usually larger than the sizes ($\sim 0.22 $ pc in the mean) of starless cores. The typical mean gas density of the cores is $\sim7\times 10^3 cm^{-3}$. Most of the cores are more likely elongated than spherical (mean aspect ratio: $\sim 2.4$). The ratio of the number of cores with EYSOs to the number of starless cores for our sample is about 0.3, suggesting that the typical lifetime of starless cores is $0.3-1.6$ Myr, about 3 times longer than the duration of the Class 0 and Class I phases. This lifetime is shorter than expected from models of ambipolar diffusion, by factors of 2-44.