Source author record

Neal J. Evans II

Neal J. Evans II appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.SR astro-ph astro-ph.CO astro-ph.EP astro-ph.IM astro-ph.HE

Catalog footprint

What is connected

100works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2026arXiv

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

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

preprint2022arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- XI. From inflow to infall in hub-filament systems

We investigate the presence of hub-filament systems in a large sample of 146 active proto-clusters, using H$^{13}$CO$^{+}$ J=1-0 molecular line data obtained from the ATOMS survey. We find that filaments are ubiquitous in proto-clusters, and hub-filament systems are very common from dense core scales ($\sim$0.1 pc) to clump/cloud scales ($\sim$1-10 pc). The proportion of proto-clusters containing hub-filament systems decreases with increasing dust temperature ($T_d$) and luminosity-to-mass ratios ($L/M$) of clumps, indicating that stellar feedback from H{\sc ii} regions gradually destroys the hub-filament systems as proto-clusters evolve. Clear velocity gradients are seen along the longest filaments with a mean velocity gradient of 8.71 km s$^{-1}$pc$^{-1}$ and a median velocity gradient of 5.54 km s$^{-1}$pc$^{-1}$. We find that velocity gradients are small for filament lengths larger than $\sim$1~pc, probably hinting at the existence of inertial inflows, although we cannot determine whether the latter are driven by large-scale turbulence or large-scale gravitational contraction. In contrast, velocity gradients below $\sim$1~pc dramatically increase as filament lengths decrease, indicating that the gravity of the hubs or cores starts to dominate gas infall at small scales. We suggest that self-similar hub-filament systems and filamentary accretion at all scales may play a key role in high-mass star formation.

preprint2022arXiv

Nobeyama Survey of Inward Motions toward Cores in Orion Identified by SCUBA-2

In this study, 36 cores (30 starless and 6 protostellar) identified in Orion were surveyed to search for inward motions. We used the Nobeyama 45 m radio telescope, and mapped the cores in the $J = 1\rightarrow0$ transitions of HCO$^+$, H$^{13}$CO$^+$, N$_2$H$^+$, HNC, and HN$^{13}$C. The asymmetry parameter $δV$, which was the ratio of the difference between the HCO$^+$ and H$^{13}$CO$^+$ peak velocities to the H$^{13}$CO$^+$ line width, was biased toward negative values, suggesting that inward motions were more dominant than outward motions. Three starless cores (10% of all starless cores surveyed) were identified as cores with blue-skewed line profiles (asymmetric profiles with more intense blue-shifted emission), and another two starless cores (7%) were identified as candidate blue-skewed line profiles. The peak velocity difference between HCO$^+$ and H$^{13}$CO$^+$ of them was up to 0.9 km s$^{-1}$, suggesting that some inward motions exceeded the speed of sound for the quiescent gas ($\sim10-17$ K). The mean of $δV$ of the five aforementioned starless cores was derived to be $-$0.5$\pm$0.3. One core, G211.16$-$19.33North3, observed using the ALMA ACA in DCO$^+$ $J = 3\rightarrow2$ exhibited blue-skewed features. Velocity offset in the blue-skewed line profile with a dip in the DCO$^+$ $J = 3\rightarrow2$ line was larger ($\sim 0.5$ km s$^{-1}$) than that in HCO$^+$ $J = 1\rightarrow0$ ($\sim 0.2$ km s$^{-1}$), which may represent gravitational acceleration of inward motions. It seems that this core is at the last stage in the starless phase, judging from the chemical evolution factor version 2.0 (CEF2.0).

preprint2022arXiv

Slow Star Formation in the Milky Way: Theory Meets Observations

The observed star formation rate of the Milky Way can be explained by applying a metallicity-dependent factor to convert CO luminosity to molecular gas mass and a star formation efficiency per free-fall time that depends on the virial parameter of a molecular cloud. These procedures also predict the trend of star formation rate surface density with Galactocentric radius. The efficiency per free-fall time variation with virial parameter plays the major role in bringing theory into agreement with observations for the total star formation rate, while the metallicity dependence of the CO luminosity to mass conversion is most notable in the variation with Galactocentric radius. Application of these changes resolves a factor of over 100 discrepancy between observed and theoretical star formation rates that has been known for nearly 50 years.

preprint2022arXiv

Tracers of Dense Gas in the Outer Galaxy

We have mapped HCN and HCO$^{+}$ (J = 1 $\to$ 0) line emission toward a sample of seven star-forming regions (with 12 + log[O/H] range from 8.34 to 8.69) in the outer Milky Way (Galactocentric distance > 9.5 kpc), using the 14-meter radio telescope of the Taeduk Radio Astronomy Observatory (TRAO). We compare these two molecular lines with other conventional tracers of dense gas, millimeter-wave continuum emission from dust and extinction thresholds ($A_{V} \geq 8$ mag), inferred from the $^{13}$CO line data. HCN and HCO$^{+}$ correlate better with the millimeter emission than with the extinction criterion. A significant amount of luminosity comes from regions below the extinction criterion and outside the millimeter clump for all the clouds. The average fraction of HCN luminosity from within the regions with $A_{V} \geq 8$ mag is $0.343\pm0.225$; for the regions of millimeter emission, it is $0.478\pm0.149$. Based on a comparison with column density maps from Herschel, HCN and HCO$^{+}$ trace dense gas in high column density regions better than does $^{13}$CO. HCO$^{+}$ is less concentrated than HCN for outer Galaxy targets, in contrast with the inner Galaxy sample, suggesting that metallicity may affect the interpretation of tracers of dense gas. The conversion factor between the dense gas mass ($M_{dense}$) and line luminosities of HCN and HCO$^{+}$, when integrated over the whole cloud, is comparable with factors used in extragalactic studies.

preprint2021arXiv

Atomic Shocks in the Outflow of L1551 IRS 5 Identified with SOFIA-upGREAT Observations of [OI]

We present velocity resolved SOFIA/upGREAT observations of [OI] and [CII] lines toward a Class I protostar, L1551 IRS 5, and its outflows. The SOFIA observations detect [OI] emission toward only the protostar and [CII] emission toward the protostar and the red-shifted outflow. The [OI] emission has a width of $\sim$100 km s$^{-1}$ only in the blue-shifted velocity, suggesting an origin in shocked gas. The [CII] lines are narrow, consistent with an origin in a photodissociation region. Differential dust extinction from the envelope due to the inclination of the outflows is the most likely cause of the missing red-shifted [OI] emission. Fitting the [OI] line profile with two Gaussian components, we find one component at the source velocity with a width of $\sim$20 km s$^{-1}$ and another extremely broad component at -30 km s$^{-1}$ with a width of 87.5 km s$^{-1}$, the latter of which has not been seen in L1551 IRS 5. The kinematics of these two components resemble cavity shocks in molecular outflows and spot shocks in jets. Radiative transfer calculations of the [OI], high-J CO, and H$_2$O lines in the cavity shocks indicate that [OI] dominates the oxygen budget, making up more than 70% of the total gaseous oxygen abundance and suggesting [O]/[H] of $\sim$1.5$\times$10$^{-4}$. Attributing the extremely broad [OI] component to atomic winds, we estimate the intrinsic mass loss rate of (1.3$\pm$0.8)$\times$ 10$^{-6}$ M$_{\odot}$ yr$^{-1}$. The intrinsic mass loss rates derived from low-J CO, [OI], and HI are similar, supporting the model of momentum-conserving outflows, where the atomic wind carries most momentum and drives the molecular outflows.

preprint2021arXiv

Molecular Cloud Cores with High Deuterium Fractions: Nobeyama Mapping Survey

We present the results of on-the-fly mapping observations of 44 fields containing 107 SCUBA-2 cores in the emission lines of molecules, N$_2$H$^+$, HC$_3$N, and CCS at 82$-$94 GHz using the Nobeyama 45-m telescope. This study aimed at investigating the physical properties of cores that show high deuterium fractions and might be close to the onset of star formation. We found that the distributions of the N$_2$H$^+$ and HC$_3$N line emissions are approximately similar to that of 850-$μ$m dust continuum emission, whereas the CCS line emission is often undetected or is distributed in a clumpy structure surrounding the peak position of the 850-$μ$m dust continuum emission. Occasionally (12%), we observe the CCS emission which is an early-type gas tracer toward the young stellar object, probably due to local high excitation. Evolution toward star formation does not immediately affect nonthermal velocity dispersion.

preprint2020arXiv

ALMA ACA and Nobeyama observations of two Orion cores in deuterated molecular lines

We mapped two molecular cloud cores in the Orion A cloud with the ALMA ACA 7-m Array and with the Nobeyama 45-m radio telescope. These cores have bright N$_2$D$^+$ emission in single-pointing observations with the Nobeyama 45-m radio telescope, have relatively high deuterium fraction, and are thought to be close to the onset of star formation. One is a star-forming core, and the other is starless. These cores are located along filaments observed in N$_2$H$^+$, and show narrow linewidths of 0.41 km s$^{-1}$ and 0.45 km s$^{-1}$ in N$_2$D$^+$, respectively, with the Nobeyama 45-m telescope. Both cores were detected with the ALMA ACA 7m Array in the continuum and molecular lines at Band 6. The starless core G211 shows clumpy structure with several sub-cores, which in turn show chemical differences. Also, the sub-cores in G211 have internal motions that are almost purely thermal. The starless sub-core G211D, in particular, shows a hint of the inverse P Cygni profile, suggesting infall motion. The star-forming core G210 shows an interesting spatial feature of two N$_2$D$^+$ peaks of similar intensity and radial velocity located symmetrically with respect to the single dust continuum peak. One interpretation is that the two N$_2$D$^+$ peaks represent an edge-on pseudo-disk. The CO outflow lobes, however, are not directed perpendicular to the line connecting both N$_2$D$^+$ peaks.

preprint2020arXiv

Large-scale molecular gas distribution in the M17 cloud complex: dense gas conditions of massive star formation?

The non-uniform distribution of gas and protostars in molecular clouds is caused by combinations of various physical processes that are difficult to separate. We explore this non-uniform distribution in the M17 molecular cloud complex that hosts massive star formation activity using the $^{12}$CO ($J=1-0$) and $^{13}$CO ($J=1-0$) emission lines obtained with the Nobeyama 45m telescope. Differences in clump properties such as mass, size, and gravitational boundedness reflect the different evolutionary stages of the M17-H{\scriptsize II} and M17-IRDC clouds. Clumps in the M17-H{\scriptsize II} cloud are denser, more compact, and more gravitationally bound than those in M17-IRDC. While M17-H{\scriptsize II} hosts a large fraction of very dense gas (27\%) that has column density larger than the threshold of $\sim$ 1 g cm$^{-2}$ theoretically predicted for massive star formation, this very dense gas is deficient in M17-IRDC (0.46\%). Our HCO$^+$ ($J=1-0$) and HCN ($J=1-0$) observations with the TRAO 14m telescope, {\nlqb trace all gas with column density higher than $3\times 10^{22}$ cm$^{-2}$}, confirm the deficiency of high density ($\gtrsim 10^5$ cm$^{-3}$) gas in M17-IRDC. Although M17-IRDC is massive enough to potentially form massive stars, its deficiency of very dense gas and gravitationally bound clumps can explain the current lack of massive star formation.

preprint2020arXiv

Star Formation Occurs in Dense Gas, but What Does "Dense" Mean?

We report results of a project to map HCN and HCO+ J = 1-0 emission toward a sample of molecular clouds in the inner Galaxy, all containing dense clumps that are actively engaged in star formation. We compare these two molecular line tracers with millimeter continuum emission and extinction, as inferred from 13CO, as tracers of dense gas in molecular clouds. The fraction of the line luminosity from each tracer that comes from the dense gas, as measured by AV > 8 mag, varies substantially from cloud to cloud. In all cases, a substantial fraction (in most cases, the majority) of the total luminosity arises in gas below the AV > 8 mag threshold and outside the region of strong mm continuum emission. Measurements of the luminosity of HCN toward other galaxies will likely be dominated by such gas at lower surface density. Substantial, even dominant, contributions to the total line luminosity can arise in gas with densities typical of the cloud as a whole (densities about 100 per cubic cm). Defining the dense clump from the HCN or HCO+ emission itself, similarly to previous studies, leads to a wide range of clump properties, with some being considerably larger and less dense than in previous studies. HCN and HCO+ have similar ability to trace dense gas for the clouds in this sample. For the two clouds with low virial parameters, the 13CO is definitely a worse tracer of the dense gas, but for the other four, it is equally good (or bad) at tracing dense gas.

preprint2016arXiv

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

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

preprint2016arXiv

Probing Planet Forming Zones with Rare CO Isotopologues

The gas near the midplanes of planet-forming protostellar disks remains largely unprobed by observations due to the high optical depth of commonly observed molecules such as CO and H$_2$O. However, rotational emission lines from rare molecules may have optical depths near unity in the vertical direction, so that the lines are strong enough to be detected, yet remain transparent enough to trace the disk midplane. Here we present a chemical model of an evolving T-Tauri disk and predict the optical depths of rotational transitions of $^{12}$C$^{16}$O, $^{13}$C$^{16}$O, $^{12}$C$^{17}$O and $^{12}$C$^{18}$O. The MRI-active disk is primarily heated by the central star due to the formation of the dead zone. CO does not freeze out in our modeled region within $70$AU around a sunlike star. However, the abundance of CO decreases because of the formation of complex organic molecules (COM), producing an effect that can be misinterpreted as the "snow line". These results are robust to variations in our assumptions about the evolution of the gas to dust ratio. The optical depths of low-order rotational lines of C$^{17}$O are around unity, making it possible to see into the disk midplane using C$^{17}$O. Combining observations with modeled C$^{17}$O$/$H$_2$ ratios, like those we provide, can yield estimates of protoplanetary disks' gas masses.

preprint2016arXiv

The abundance and thermal history of water ice in the disk surrounding HD142527 from the DIGIT Herschel Key Program

The presence or absence of ice in protoplanetary disks is of great importance for the formation of planets. By enhancing the solid surface density and increasing the sticking efficiency, ice catalyzes the rapid formation of planetesimals and decreases the time scale for giant planet core accretion. Aims: In this paper we analyse the composition of the outer disk around the Herbig star HD~142527. We focus on the composition of the water ice, but also analyse the abundances of previously proposed minerals. Methods: We present new Herschel far infrared spectra and a re-reduction of archival data from the Infrared Space Observatory (ISO). We model the disk using full 3D radiative transfer to obtain the disk structure. Also, we use an optically thin analysis of the outer disk spectrum to obtain firm constraints on the composition of the dust component. Results: The water ice in the disk around HD~142527 contains a large reservoir of crystalline water ice. We determine the local abundance of water ice in the outer disk (i.e. beyond 130\,AU). The re-reduced ISO spectrum differs significantly from that previously published, but matches the new Herschel spectrum at their common wavelength range. In particular, we do not detect any significant contribution from carbonates or hydrous silicates, in contrast to earlier claims. Conclusions: The amount of water ice detected in the outer disk requires $\sim80\,$\% of the oxygen atoms. This is comparable to the water ice abundance in the outer solar system, in comets and in dense interstellar clouds. The water ice is highly crystalline while the temperatures where we detect it are too low to crystallize the water on relevant time scales. We discuss the implications of this finding.

preprint2016arXiv

The CDF Archive: Herschel PACS and SPIRE Spectroscopic Data Pipeline and Products for Protostars and Young Stellar Objects

We present the COPS-DIGIT-FOOSH (CDF) Herschel spectroscopy data product archive, and related ancillary data products, along with data fidelity assessments, and a user-created archive in collaboration with the Herschel-PACS and SPIRE ICC groups. Our products include datacubes, contour maps, automated line fitting results, and best 1-D spectra products for all protostellar and disk sources observed with PACS in RangeScan mode for two observing programs: the DIGIT Open Time Key Program (KPOT_nevans_1 and SDP_nevans_1; PI: N. Evans), and the FOOSH Open Time Program (OT1_jgreen02_2; PI: J. Green). In addition, we provide our best SPIRE-FTS spectroscopic products for the COPS Open Time Program (OT2_jgreen02_6; PI: J. Green) and FOOSH sources. We include details of data processing, descriptions of output products, and tests of their reliability for user applications. We identify the parts of the dataset to be used with caution. The resulting absolute flux calibration has improved in almost all cases. Compared to previous reductions, the resulting rotational temperatures and numbers of CO molecules have changed substantially in some sources. On average, however, the rotational temperatures have not changed substantially (< 2%), but the number of warm (Trot ~ 300 K) CO molecules has increased by about 18%.

preprint2016arXiv

The inflow signature toward different evolutionary phases of massive star formation

We analyzed both HCN J=1-0 and HNC J=1-0 line profiles to study the inflow motions in different evolutionary stages of massive star formation: 54 infrared dark clouds (IRDCs), 69 high-mass protostellar object (HMPOs), and 54 ultra-compact HII regions (UCHIIs). The inflow asymmetry in HCN spectra seems to be prevalent throughout all the three evolutionary phases, with IRDCs showing the largest excess in blue profile. In the case of HNC spectra, the prevalence of blue sources does not appear, excepting for IRDCs. We suggest that this line is not appropriate to trace inflow motion in evolved stages of massive star formation because the abundance of HNC decreases at high temperatures. This result spotlights the importance of considering chemistry in the dynamics study of massive star-forming regions. The fact that the IRDCs show the highest blue excess in both transitions indicates that the most active inflow occurs in the early phase of star formation, i.e., the IRDC phase rather than in the later phases. However, mass is still inflowing onto some UCHIIs. We also found that the absorption dips of the HNC spectra in 6 out of 7 blue sources are red-shifted relative to their systemic velocities. These red-shifted absorption dips may indicate global collapse candidates, although mapping observations with better resolution are needed to examine this feature in more detail.

preprint2016arXiv

The Mid-Infrared Evolution of the FU Orionis Disk

We present new SOFIA-FORCAST observations obtained in Feburary 2016 of the archetypal outbursting low mass young stellar object FU Orionis, and compare the continuum, solid state, and gas properties with mid-IR data obtained at the same wavelengths in 2004 with Spitzer-IRS. In this study, we conduct the first mid-IR spectroscopic comparison of an FUor over a long time period. Over a 12 year period, UBVR monitoring indicates that FU Orionis has continued its steady decrease in overall brightness by ~ 14%. We find that this decrease in luminosity occurs only at wavelengths < 20 microns. In particular, the continuum short ward of the silicate emission complex at 10 microns exhibits a ~ 12% (~ 3 sigma) drop in flux density, but no apparent change in slope; both the Spitzer and SOFIA spectra are consistent with a 7200 K blackbody. Additionally, the detection of water absorption is consistent with the Spitzer spectrum. The silicate emission feature at 10 microns continues to be consistent with unprocessed grains, unchanged over 12 years. We conclude that either the accretion rate in FU Orionis has decreased by ~ 12-14% over this time baseline, or that the inner disk has cooled, but the accretion disk remains in a superheated state outside of the innermost region.

preprint2015arXiv

"Dust, Ice, and Gas In Time" (DIGIT) Herschel Observations of GSS30-IRS1 in Ophiuchus

As a part of the "Dust, Ice, and Gas In Time" (DIGIT) key program on Herschel, we observed GSS30-IRS1, a Class I protostar located in Ophiuchus (d = 120 pc), with Herschel/Photodetector Array Camera and Spectrometer (PACS). More than 70 lines were detected within a wavelength range from 50 micron to 200 micron, including CO, H2O, OH, and two atomic [O I] lines at 63 and 145 micron. The [C II] line, known as a tracer of externally heated gas by the interstellar radiation field, is also detected at 158 micron. All lines, except [O I] and [C II], are detected only at the central spaxel of 9.4" X 9.4". The [O I] emissions are extended along a NE-SW orientation, and the [C II] line is detected over all spaxels, indicative of external PDR. The total [C II] intensity around GSS30 reveals that the far-ultraviolet radiation field is in the range of 3 to 20 G0, where G0 is in units of the Habing Field, 1.6 X 10^{-3} erg cm^{-2} s^{-1}. This enhanced external radiation field heats the envelope of GSS30-IRS1, causing the continuum emission to be extended, unlike the molecular emission. The best-fit continuum model of GSS30-IRS1 with the physical structure including flared disk, envelope, and outflow shows that the internal luminosity is 10 Lsun, and the region is externally heated by a radiation field enhanced by a factor of 130 compared to the standard local interstellar radiation field.

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

Evidence for Decay of Turbulence by MHD Shocks in Molecular Clouds via CO Emission

We utilize observations of sub-millimeter rotational transitions of CO from a Herschel Cycle 2 open time program ("COPS", PI: J. Green) to identify previously predicted turbulent dissipation by magnetohydrodynamic (MHD) shocks in molecular clouds. We find evidence of the shocks expected for dissipation of MHD turbulence in material not associated with any protostar. Two models fit about equally well: model 1 has a density of 10$^{3}$ cm$^{-3}$, a shock velocity of 3 km s$^{-1}$, and a magnetic field strength of 4 $μ$G; model 2 has a density of 10$^{3.5}$ cm$^{-3}$, a shock velocity of $2$ km s$^{-1}$, and a magnetic field strength of 8 $μ$G. Timescales for decay of turbulence in this region are comparable to crossing times. Transitions of CO up to $J$ of 8, observed close to active sites of star formation, but not within outflows, can trace turbulent dissipation of shocks stirred by formation processes. Although the transitions are difficult to detect at individual positions, our Herschel-SPIRE survey of protostars provides a grid of spatially-distributed spectra within molecular clouds. We averaged all spatial positions away from known outflows near seven protostars. We find significant agreement with predictions of models of turbulent dissipation in slightly denser (10$^{3.5}$ cm$^{-3}$) material with a stronger magnetic field (24 $μ$G) than in the general molecular cloud.

preprint2015arXiv

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

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

preprint2015arXiv

Infrared emission from tidal disruption events --- probing the pc-scale dust content around galactic nuclei

Recent UV-optical surveys have been successful in finding tidal disruption events (TDEs), in which a star is tidally disrupted by a supermassive black hole (BH). These TDEs release a huge amount of radiation energy ~ 10^51-52 erg into the circum-nuclear medium. If the medium is dusty, most of the radiation energy will be absorbed by dust grains within ~ 1 pc from the BH and re-radiated in the infrared. We calculate the dust emission lightcurve from a 1-D radiative transfer model, taking into account the time-dependent heating, cooling and sublimation of dust grains. We show that the dust emission peaks at 3-10 microns and has typical luminosities ~ 10^42-43 erg/s (with sky covering factor of dusty clouds ranging from 0.1-1). This is detectable by current generation of telescopes. In the near future, James Webb Space Telescope will be able to perform photometric and spectroscopic measurements, in which silicate or polycyclic aromatic hydrocarbon (PAH) features may be found. Observations at rest-frame wavelength > 2 microns have only been reported from two TDE candidates, SDSS J0952+2143 and Swift J1644+57. Although consistent with the dust emission from TDEs, the mid-infrared fluxes of the two events may be from other sources. Long-term monitoring is needed to draw a firm conclusion. We also point out two nearby TDE candidates (ASSASN-14ae and -14li) where the dust emission may be currently detectable. The dust infrared emission can give a snapshot of the pc-scale dust content around weakly- or non-active galactic nuclei, which is hard to probe otherwise.

preprint2015arXiv

The Bolocam Galactic Plane Survey. XI. Temperatures and Substructure of Galactic Clumps Based on 350 micron Observations

We present 107 maps of continuum emission at 350 microns from Galactic molecular clumps. Observed sources were mainly selected from the Bolocam Galactic Plane Survey (BGPS) catalog, with 3 additional maps covering star forming regions in the outer Galaxy. The higher resolution of the SHARC-II images (8.5" beam) compared with the 1.1 mm images from BGPS (33" beam) allowed us to identify a large population of smaller substructures within the clumps. A catalog is presented for the 1386 sources extracted from the 350 micron maps. The color temperature distribution of clumps based on the two wavelengths has a median of 13.3 K and mean of 16.3 +- 0.4 K, assuming an opacity law index of 1.7. For the structures with good determination of color temperatures, the mean ratio of gas temperature, determined from NH3 observations, to dust color temperature is 0.88 and the median ratio is 0.76. About half the clumps have more than two substructures and 22 clumps have more than 10. The fraction of the mass in dense substructures seen at 350 microns compared to the mass of their parental clump is ~0.19, and the surface densities of these substructures are, on average, 2.2 times those seen in the clumps identified at 1.1 mm. For a well-characterized sample, 88 structures (31%) exceed a surface density of 0.2 g cm^(-2), and 18 (6%) exceed 1.0 g cm^(-2), thresholds for massive star formation suggested by theorists.

preprint2015arXiv

The Bolocam Galactic Plane Survey. XIII. Physical Properties and Mass Functions of Dense Molecular Cloud Structures

We use the distance probability density function (DPDF) formalism of Ellsworth-Bowers et al. (2013, 2015) to derive physical properties for the collection of 1,710 Bolocam Galactic Plane Survey (BGPS) version 2 sources with well-constrained distance estimates. To account for Malmquist bias, we estimate that the present sample of BGPS sources is 90% complete above 400 $M_\odot$ and 50% complete above 70 $M_\odot$. The mass distributions for the entire sample and astrophysically motivated subsets are generally fitted well by a lognormal function, with approximately power-law distributions at high mass. Power-law behavior emerges more clearly when the sample population is narrowed in heliocentric distance (power-law index $α= 2.0\pm0.1$ for sources nearer than 6.5 kpc and $α= 1.9\pm0.1$ for objects between 2 kpc and 10 kpc). The high-mass power-law indices are generally $1.85 \leq α\leq 2.05$ for various subsamples of sources, intermediate between that of giant molecular clouds and the stellar initial mass function. The fit to the entire sample yields a high-mass power-law $\hatα = 1.94_{-0.10}^{+0.34}$. Physical properties of BGPS sources are consistent with large molecular cloud clumps or small molecular clouds, but the fractal nature of the dense interstellar medium makes difficult the mapping of observational categories to the dominant physical processes driving the observed structure. The face-on map of the Galactic disk's mass surface density based on BGPS dense molecular cloud structures reveals the high-mass star-forming regions W43, W49, and W51 as prominent mass concentrations in the first quadrant. Furthermore, we present a 0.25-kpc resolution map of the dense gas mass fraction across the Galactic disk that peaks around 5%.

preprint2015arXiv

The Gould Belt 'MISFITS' Survey: The Real Solar Neighborhood Protostars

We present an HCO$^{+}$ J=3-2 survey of Class 0+I and Flat SED young stellar objects (YSOs) found in the Gould Belt clouds by surveys with Spitzer. Our goal is to provide a uniform Stage 0+I source indicator for these embedded protostar candidates. We made single point HCO$^{+}$ J = 3-2 measurements toward the source positions at the CSO and APEX of 546 YSOs (89% of the Class 0+I + Flat SED sample). Using the criteria from van Kempen et al. (2009), we classify sources as Stage 0+I or bona fide protostars and find that 84% of detected sources meet the criteria. We recommend a timescale for the evolution of Stage 0+I (embedded protostars) of 0.54 Myr. We find significant correlations of HCO$^{+}$ integrated intensity with $α$ and $T_{bol}$ but not with $L_{bol}$. The detection fraction increases smoothly as a function of $α$ and $L_{bol}$, while decreasing smoothly with $T_{bol}$. Using the Stage 0+I sources tightens the relation between protostars and high extinction regions of the cloud; 89% of Stage I sources lie in regions with $A_{V}$ >8 mag. Class 0+I and Flat SED YSOs that are not detected in HCO$^{+}$ have, on average, a factor of ~2 higher $T_{bol}$ and a factor of ~5 lower $L_{bol}$ than YSOs with HCO$^{+}$ detections. We find less YSO contamination, defined as the number of undetected YSOs divided by the total number surveyed, for sources with $T_{bol}\lesssim$ 600 K and $L_{bol} \gtrsim$ 1 $L_{\odot}$. The contamination percentage is >90% at $A_{V}$< 4 mag and decreases as $A_{V}$ increases.

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.

preprint2015arXiv

Young Stellar Objects in the Gould Belt

We present the full catalog of Young Stellar Objects (YSOs) identified in the 18 molecular clouds surveyed by the Spitzer Space Telescope "cores to disks" (c2d) and "Gould Belt" (GB) Legacy surveys. Using standard techniques developed by the c2d project, we identify 3239 candidate YSOs in the 18 clouds, 2966 of which survive visual inspection and form our final catalog of YSOs in the Gould Belt. We compile extinction corrected SEDs for all 2966 YSOs and calculate and tabulate the infrared spectral index, bolometric luminosity, and bolometric temperature for each object. We find that 326 (11%), 210 (7%), 1248 (42%), and 1182 (40%) are classified as Class 0+I, Flat-spectrum, Class II, and Class III, respectively, and show that the Class III sample suffers from an overall contamination rate by background AGB stars between 25% and 90%. Adopting standard assumptions, we derive durations of 0.40-0.78 Myr for Class 0+I YSOs and 0.26-0.50 Myr for Flat-spectrum YSOs, where the ranges encompass uncertainties in the adopted assumptions. Including information from (sub)millimeter wavelengths, one-third of the Class 0+I sample is classified as Class 0, leading to durations of 0.13-0.26 Myr (Class 0) and 0.27-0.52 Myr (Class I). We revisit infrared color-color diagrams used in the literature to classify YSOs and propose minor revisions to classification boundaries in these diagrams. Finally, we show that the bolometric temperature is a poor discriminator between Class II and Class III YSOs.

preprint2014arXiv

Herschel Key Program, "Dust, Ice, and Gas In Time" (DIGIT): the origin of molecular and atomic emission in low-mass protostars in Taurus

Six low-mass embedded sources (L1489, L1551-IRS5, TMR1, TMC1-A, L1527, and TMC1) in Taurus have been observed with Herschel-PACS to cover the full spectrum from 50 to 210 $μ$m as part of the Herschel key program, "Dust, Ice, and Gas In Time (DIGIT)". The relatively low intensity of the interstellar radiation field surrounding Taurus minimizes contamination of the [C II] emission associated with the sources by diffuse emission from the cloud surface, allowing study of the [C II] emission from the source. In several sources, the [C II] emission is distributed along the outflow, as is the [O I] emission. The atomic line luminosities correlate well with each other, as do the molecular lines, but the atomic and molecular lines correlate poorly. The relative contribution of CO to the total gas cooling is constant at $\sim$30 %, while the cooling fraction by H$_2$O varies from source to source, suggesting different shock properties resulting in different photodissociation levels of H$_2$O. The gas with a power-law temperature distribution with a moderately high density can reproduce the observed CO fluxes, indicative of CO close to LTE. However, H$_2$O is mostly subthermally excited. L1551-IRS5 is the most luminous source ($L_{\rm bol}=24.5$ L$_\odot$) and the [O I] 63.1 $μ$m line accounts for more than 70 % of its FIR line luminosity, suggesting complete photodissociation of H$_2$O by a J-shock. In L1551-IRS5, the central velocity shifts of the [O I] line, which exceed the wavelength calibration uncertainty ($\sim$70 km s$^{-1}$) of PACS, are consistent with the known red- and blue-shifted outflow direction.

preprint2014arXiv

Star Formation Relations in Nearby Molecular Clouds

We test some ideas for star formation relations against data on local molecular clouds. On a cloud by cloud basis, the relation between the surface density of star formation rate and surface density of gas divided by a free-fall time, calculated from the mean cloud density, shows no significant correlation. If a crossing time is substituted for the free-fall time, there is even less correlation. Within a cloud, the star formation rate volume and surface densities increase rapidly with the corresponding gas densities, faster than predicted by models using the free-fall time defined from the local density. A model in which the star formation rate depends linearly on the mass of gas above a visual extinction of 8 mag describes the data on these clouds, with very low dispersion. The data on regions of very massive star formation, with improved star formation rates based on free-free emission from ionized gas, also agree with this linear relation.

preprint2014arXiv

The Bolocam Galactic Plane Survey. XII. Distance Catalog Expansion Using Kinematic Isolation of Dense Molecular Cloud Structures With 13CO(1-0)

We present an expanded distance catalog for 1,710 molecular cloud structures identified in the Bolocam Galactic Plane Survey (BGPS) version 2, representing a nearly threefold increase over the previous BGPS distance catalog. We additionally present a new method for incorporating extant data sets into our Bayesian distance probability density function (DPDF) methodology. To augment the dense-gas tracers (e.g., HCO+(3-2), NH3(1,1)) used to derive line-of-sight velocities for kinematic distances, we utilize the Galactic Ring Survey 13CO(1-0) data to morphologically extract velocities for BGPS sources. The outline of a BGPS source is used to select a region of the GRS 13CO data, along with a reference region to subtract enveloping diffuse emission, to produce a line profile of 13CO matched to the BGPS source. For objects with a HCO+(3-2) velocity, \approx 95% of the new 13CO(1-0) velocities agree with that of the dense gas. A new prior DPDF for kinematic distance ambiguity (KDA) resolution, based on a validated formalism for associating molecular cloud structures with known objects from the literature, is presented. We demonstrate this prior using catalogs of masers with trigonometric parallaxes and HII regions with robust KDA resolutions. The distance catalog presented here contains well-constrained distance estimates for 20% of BGPS V2 sources, with typical distance uncertainties \lesssim 0.5 kpc. Approximately 75% of the well-constrained sources lie within 6 kpc of the Sun, concentrated in the Scutum-Centarus arm. Galactocentric positions of objects additionally trace out portions of the Sagittarius, Perseus, and Outer arms in the first and second Galactic quadrants, and we also find evidence for significant regions of interarm dense gas.

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.

preprint2014arXiv

The Gould Belt Very Large Array Survey IV: The Taurus-Auriga complex

We present a multi-epoch radio study of the Taurus-Auriga star-forming complex made with the Karl G. Jansky Very Large Array at frequencies of 4.5 GHz and 7.5 GHz. We detect a total of 610 sources, 59 of which are related to young stellar objects and 18 to field stars. The properties of 56\% of the young stars are compatible with non-thermal radio emission. We also show that the radio emission of more evolved young stellar objects tends to be more non-thermal in origin and, in general, that their radio properties are compatible with those found in other star forming regions. By comparing our results with previously reported X-ray observations, we notice that young stellar objects in Taurus-Auriga follow a Güdel-Benz relation with $κ$=0.03, as we previously suggested for other regions of star formation. In general, young stellar objects in Taurus-Auriga and in all the previous studied regions seem to follow this relation with a dispersion of $\sim1$ dex. Finally, we propose that most of the remaining sources are related with extragalactic objects but provide a list of 46 unidentified radio sources whose radio properties are compatible with a YSO nature.

preprint2014arXiv

The Star Formation Relation for Regions in the Galactic Plane: The Effect of Spatial Resolution

We examined the relations between molecular gas surface density and star formation rate surface density in a 11 square degree region of the Galactic Plane. Dust continuum at 1.1 mm from the Bolocam Galactic Plane Survey and 22 micron emission from the WISE All-sky survey were used as tracers of molecular gas and star formation rate, respectively, across Galactic longitude of 31.5 > l > 20.5 and Galactic latitude of 0.5 > b > -0.5. The relation was studied over a range of resolutions from 33 arcsecond to 20 arcminute by convolving images to larger scales. The pixel-by-pixel correlation between 1.1 mm and 22 micron increases rapidly at small scales and levels off at the scale of 5 - 8 arcminute. We studied the star formation relation based on pixel-by-pixel analysis and 1.1 mm and 22 micron peaks analysis. The star formation relation was found to be nearly linear with no significant changes in the form of the relation across all spatial scales and lie above the extragalactic relation from Kennicutt (1998). The average gas depletion time is approximately 200 Myr and does not change significantly at different scales, but the scatter in the depletion time decreases as the scale increases.

preprint2013arXiv

ALMA observations of the massive molecular outflow G331.512-0.103

The object of this study is one of the most energetic and luminous molecular outflows known in the Galaxy, G331.512-0.103. Observations with ALMA Band 7 (350 GHz; 0.86 mm) reveal a very compact, extremely young bipolar outflow and a more symmetric outflowing shocked shell surrounding a very small region of ionized gas. The velocities of the bipolar outflow are about 70 km s^{-1} on either side of the systemic velocity. The expansion velocity of the shocked shell is ~24 km s^{-1}, implying a crossing time of about 2000 yrs. Along the symmetry axis of the outflow, there is a velocity feature, which could be a molecular "bullet" of high-velocity dense material. The source is one of the youngest examples of massive molecular outflow found associated with a high-mass star.

preprint2013arXiv

An Analysis of the Environments of FU Orionis Objects with Herschel

We present Herschel-HIFI, SPIRE, and PACS 50-670 μm imaging and spectroscopy of six FU Orionis-type objects and candidates (FU Orionis, V1735 Cyg, V1515 Cyg, V1057 Cyg, V1331 Cyg, and HBC 722), ranging in outburst date from 1936-2010, from the "FOOSH" (FU Orionis Objects Surveyed with Herschel) program, as well as ancillary results from Spitzer-IRS and the Caltech Submillimeter Observatory. In their system properties (Lbol, Tbol, line emission), we find that FUors are in a variety of evolutionary states. Additionally, some FUors have features of both Class I and II sources: warm continuum consistent with Class II sources, but rotational line emission typical of Class I, far higher than Class II sources of similar mass/luminosity. Combining several classification techniques, we find an evolutionary sequence consistent with previous mid-IR indicators. We detect [O I] in every source at luminosities consistent with Class 0/I protostars, much greater than in Class II disks. We detect transitions of 13CO (J_up of 5 to 8) around two sources (V1735 Cyg and HBC 722) but attribute them to nearby protostars. Of the remaining sources, three (FU Ori, V1515 Cyg, and V1331 Cyg) exhibit only low-lying CO, but one (V1057 Cyg) shows CO up to J = 23 - 22 and evidence for H2O and OH emission, at strengths typical of protostars rather than T Tauri stars. Rotational temperatures for "cool" CO components range from 20-81 K, for ~ 10^50 total CO molecules. We detect [C I] and [N II] primarily as diffuse emission.

preprint2013arXiv

DIGIT survey of far-infrared lines from protoplanetary discs: II. CO

CO is an important component of a protoplanetary disc as it is one of the most abundant gas phase species. Furthermore, observations of the CO transitions can be used as a diagnostic of the gas, tracing conditions in both the inner and outer disc. We present Herschel/PACS spectroscopy of a sample of 22 Herbig Ae/Be (HAEBEs) and 8 T Tauri stars (TTS), covering the pure rotational CO transitions from J = 14-13 up to J = 49-48. CO is detected in only 5 HAEBEs, AB Aur, HD 36112, HD 97048, HD 100546 and IRS 48, and in 4 TTS, AS 205, S CrA, RU Lup, and DG Tau. The highest transition seen is J = 36-35, with Eup of 3669 K, detected in HD 100546 and DG Tau. We construct rotational diagrams for the discs with at least 3 CO detections to derive Trot, and find average temperatures of 270 K for the HAEBEs and 485 K for the TTS. HD 100546 requires an extra temperature component with Trot ~ 900-1000 K, suggesting a range of temperatures in its disc atmosphere, consistent with thermo-chemical disc models. In HAEBEs, the objects with CO detections all have flared discs in which the gas and dust are thermally decoupled. We use a small model grid to analyse our observations and find that an increased amount of flaring means higher line flux, as it increases the mass in warm gas. CO is not detected in our flat discs as the emission is below the detection limit. We find that HAEBE sources with CO detections have high Luv and strong PAH emission, again connected to the heating of the gas. In TTS, the objects with CO detections are all sources with evidence for a disc wind or outflow. For both groups of objects, sources with CO detections generally have high UV luminosity (either stellar in HAEBEs or due to accretion in TTS), but this is not a sufficient condition for the detection of the far-IR CO lines.

preprint2013arXiv

DIGIT survey of far-infrared lines from protoplanetary disks I

[abridged] We present far-infrared spectroscopic observations of PMS stars taken with Herschel/PACS as part of the DIGIT key project. The sample includes 22 Herbig AeBe and 8 T Tauri sources. Multiple atomic fine structure and molecular lines are detected at the source position: [OI], [CII], CO, OH, H_2O, CH^+. The most common feature is the [OI] 63micron line detected in almost all of the sources followed by OH. In contrast with CO, OH is detected toward both Herbig AeBe groups (flared and non-flared sources). An isothermal LTE slab model fit to the OH lines indicates column densities of 10^13 < N_OH < 10^16 cm^-2, emitting radii 15 < r < 100 AU and excitation temperatures 100 < T_ex < 400 K. The OH emission thus comes from a warm layer in the disk at intermediate stellar distances. Warm H_2O emission is detected through multiple lines toward the T Tauri systems AS 205, DG Tau, S CrA and RNO 90 and three Herbig AeBe systems HD 104237, HD 142527, HD 163296 (through line stacking). Overall, Herbig AeBe sources have higher OH/H_2O abundance ratios across the disk than do T Tauri disks, from near- to far-infrared wavelengths. Far-infrared CH^+ emission is detected toward HD 100546 and HD 97048. The slab model suggests moderate excitation (T_ex ~ 100 K) and compact (r ~ 60 AU) emission in the case of HD 100546. The [CII] emission is spatially extended in all sources where the line is detected. This suggests that not all [CII] emission is associated with the disk and that there is a substantial contribution from diffuse material around the young stars. The flux ratios of the atomic fine structure lines are consistent with a disk origin for the oxygen lines for most of the sources.

preprint2013arXiv

Embedded Protostars in the Dust, Ice, and Gas In Time (DIGIT) Key Program: Continuum SEDs, and an Inventory of Characteristic Far-Infrared Lines from PACS Spectroscopy

We present 50-210 um spectral scans of 30 Class 0/I protostellar sources, obtained with Herschel-PACS, and 0.5-1000 um SEDs, as part of the Dust, Ice, and Gas in Time (DIGIT) Key Program. Some sources exhibit up to 75 H2O lines ranging in excitation energy from 100-2000 K, 12 transitions of OH, and CO rotational lines ranging from J=14-13 up to J=40-39. [O I] is detected in all but one source in the entire sample; among the sources with detectable [O I] are two Very Low Luminosity Objects (VeLLOs). The mean 63/145 um [O I] flux ratio is 17.2 +/- 9.2. The [O I] 63 um line correlates with Lbol, but not with the time-averaged outflow rate derived from low-J CO maps. [C II] emission is in general not local to the source. The sample Lbol increased by 1.25 (1.06) and Tbol decreased to 0.96 (0.96) of mean (median) values with the inclusion of the Herschel data. Most CO rotational diagrams are characterized by two optically thin components (<N> = (0.70 +/- 1.12) x 10^49 total particles). N_CO correlates strongly with Lbol, but neither Trot nor N_CO(warm)/N_CO(hot) correlates with Lbol, suggesting that the total excited gas is related to the current source luminosity, but that the excitation is primarily determined by the physics of the interaction (e.g., UV- heating/shocks). Rotational temperatures for H2O (<Trot> = 194 +/- 85 K) and OH (<Trot> = 183 +/- 117 K) are generally lower than for CO, and much of the scatter in the observations about the best fit is attributed to differences in excitation conditions and optical depths amongst the detected lines.

preprint2013arXiv

L1448-MM observations by the Herschel Key program, "Dust, Ice, and Gas In Time" (DIGIT)

We present Herschel/PACS observations of L1448-MM, a Class 0 protostar with a prominent outflow. Numerous emission lines are detected at 55 < lambda < 210 micrometer including CO, OH, H2O, and [O I]. We investigate the spatial distribution of each transition to find that lines from low energy levels tend to distribute along the outflow direction while lines from high energy levels peak at the central spatial pixel. Spatial maps reveal that OH emission lines are formed in a relatively small area, while [O I] emission is extended. According to the rotational diagram analysis, the CO emission can be fitted by two (warm and hot) temperature components. For H2O, the ortho-to-para ratio is close to 3. The non-LTE LVG calculations suggest that CO and H2O lines could instead be formed in a high kinetic temperature (T > 1000 K) environment, indicative of a shock origin. For OH, IR-pumping processes play an important role in the level population. The molecular emission in L1448-MM is better explained with a C-shock model, but the atomic emission of PACS [O I] and Spitzer/IRS [Si II] emission is not consistent with C-shocks, suggesting multiple shocks in this region. Water is the major line coolant of L1448-MM in the PACS wavelength range, and the best-fit LVG models predict that H2O and CO emit (50-80)% of their line luminosity in the PACS wavelength range.

preprint2013arXiv

Physical characteristics of G331.5-0.1: The luminous central region of a Giant Molecular Cloud

We report molecular line and dust continuum observations toward the high-mass star forming region G331.5-0.1, one of the most luminous regions of massive star-formation in the Milky Way, located at the tangent region of the Norma spiral arm, at a distance of 7.5 kpc. Molecular emission was mapped toward the G331.5-0.1 GMC in the CO (J=1-0) and C18O (J=1-0) lines with NANTEN, while its central region was mapped in CS (J=2-1 and J=5-4) with SEST, and in CS (J=7-6) and 13CO (J=3-2) with ASTE. Continuum emission mapped at 1.2 mm with SIMBA and at 0.87 mm with LABOCA reveal the presence of six compact and luminous dust clumps, making this source one of the most densely populated central regions of a GMC in the Galaxy. The dust clumps are associated with molecular gas and they have the following average properties: size of 1.6 pc, mass of 3.2x10^3 Msun, molecular hydrogen density of 3.7x10^4 cm^{-3}, dust temperature of 32 K, and integrated luminosity of 5.7x10^5 Lsun, consistent with values found toward other massive star forming dust clumps. The CS and 13CO spectra show the presence of two velocity components: a high-velocity component at ~ -89 km s^{-1}, seen toward four of the clumps, and a low-velocity component at ~ -101 km s^{-1} seen toward the other two clumps. Radio continuum emission is present toward four of the molecular clumps, with spectral index estimated for two of them of 0.8+-0.2 and 1.2+-0.2. A high-velocity molecular outflow is found at the center of the brightest clump, with a line width of 26 km s^{-1} (FWHM) in CS (J=7-6). Observations of SiO (J=7-6 and J=8-7), and SO (J_K=8_8-7_7 and J_K=8_7-7_6) lines provide estimates of the gas rotational temperature toward this outflow >120 K and >75 K, respectively.

preprint2013arXiv

Testing 24 micron and Infrared Luminosity as Star Formation Tracers for Galactic Star Forming Regions

We have tested some relations for star formation rates used in extra-galactic studies for regions within the Galaxy. In nearby molecular clouds, where the IMF is not fully-sampled, the dust emission at 24 micron greatly underestimates star formation rates (by a factor of 100 on average) when compared to star formation rates determined from counting YSOs. The total infrared emission does no better. In contrast, the total far-infrared method agrees within a factor of 2 on average with star formation rates based on radio continuum emission for massive, dense clumps that are forming enough massive stars to have the total infrared luminosity exceed 10^4.5 Lsun. The total infrared and 24 micron also agree well with each other for both nearby, low-mass star forming regions and the massive, dense clumps regions.

preprint2013arXiv

The Bolocam Galactic Plane Survey IX: Data Release 2 and Outer Galaxy Extension

We present a re-reduction and expansion of the Bolocam Galactic Plane Survey, first presented by Aguirre et al. (2011) and Rosolowsky et al. (2010). The BGPS is a 1.1 mm survey of dust emission in the Northern galactic plane, covering longitudes -10 < \ell < 90 and latitudes |b| < 0.5 with a typical 1-σRMS sensitivity of 30-100 mJy in a 33" beam. Version 2 of the survey includes an additional 20 square degrees of coverage in the 3rd and 4th quadrants and 2 square degrees in the 1st quadrant. The new data release has improved angular recovery, with complete recovery out to 80" and partial recovery to 300", and reduced negative bowls around bright sources resulting from the atmospheric subtraction process. We resolve the factor of 1.5 flux calibration offset between the v1.0 data release and other data sets and determine that there is no offset between v2.0 and other data sets. The v2.0 pointing accuracy is tested against other surveys and demonstrated to be accurate and an improvement over v1.0. We present simulations and tests of the pipeline and its properties, including measurements of the pipeline's angular transfer function. The Bolocat cataloging tool was used to extract a new catalog, which includes 8594 sources, with 591 in the expanded regions. We have demonstrated that the Bolocat 40" and 80" apertures are accurate even in the presence of strong extended background emission. The number of sources is lower than in v1.0, but the amount of flux and area included in identified sources is larger.

preprint2013arXiv

The Bolocam Galactic Plane Survey. VIII. A Mid-Infrared Kinematic Distance Discrimination Method

We present a new distance estimation method for dust-continuum-identified molecular cloud clumps. Recent (sub-)millimeter Galactic plane surveys have cataloged tens of thousands of these objects, but detailed study of their physical properties requires robust distance determinations. We derive Bayesian distance probability density functions (DPDFs) for 770 objects from the Bolocam Galactic Plane Survey in the longitude range 7.5 < l < 65. The DPDF formalism is based on kinematic distances, and uses external data sets to place prior distance probabilities to resolve the kinematic distance ambiguity (KDA) for objects in the inner Galaxy. We present priors related to the mid-infrared absorption of dust in dense molecular regions and the Galactic distribution of molecular gas. By assuming a numerical model of Galactic mid-infrared emission and simple radiative transfer, we match the morphology of millimeter thermal dust emission with mid-infrared absorption to compute a prior DPDF for distance discrimination. Selecting objects first from (sub-)millimeter source catalogs avoids a bias towards the darkest infrared dark clouds (IRDCs) and extends the range of heliocentric distance and contrast range probed by mid-infrared extinction. We derive well-constrained KDA resolutions for 618 molecular cloud clumps, with approximately 15% placed at or beyond the tangent distance. Objects sufficiently dark to be cataloged as IRDCs are generally placed at the near kinematic distance. Distance comparisons with Galactic Ring Survey KDA resolutions yield a 92% agreement, and a Galactic face-on view reveals sections of the Sagittarius and Scutum-Centaurus Arms. This KDA-resolution method for large catalogs of sources through the combination of (sub-)millimeter and mid-infrared observations of molecular cloud clumps is generally applicable to other dust-continuum Galactic plane surveys.

preprint2013arXiv

The Gould's Belt Very Large Array Survey I: The Ophiuchus complex

We present large-scale ($\sim$ 2000 square arcminutes), deep ($\sim$ 20 $μ$Jy), high-resolution ($\sim$ 1$''$) radio observations of the Ophiuchus star-forming complex obtained with the Karl G.\ Jansky Very Large Array at $λ$ = 4 and 6 cm. In total, 189 sources were detected, 56 of them associated with known young stellar sources, and 4 with known extragalactic objects; the other 129 remain unclassified, but most of them are most probably background quasars. The vast majority of the young stars detected at radio wavelengths have spectral types K or M, although we also detect 4 objects of A/F/B types and 2 brown dwarf candidates. At least half of these young stars are non-thermal (gyrosynchrotron) sources, with active coronas characterized by high levels of variability, negative spectral indices, and (in some cases) significant circular polarization. As expected, there is a clear tendency for the fraction of non-thermal sources to increase from the younger (Class 0/I or flat spectrum) to the more evolved (Class III or weak line T Tauri) stars. The young stars detected both in X-rays and at radio wavelengths broadly follow a Güdel-Benz relation, but with a different normalization than the most radio-active types of stars. Finally, we detect a $\sim$ 70 mJy compact extragalactic source near the center of the Ophiuchus core, which should be used as gain calibrator for any future radio observations of this region.

preprint2013arXiv

The Luminosities of Protostars in the Spitzer c2d and Gould Belt Legacy Clouds

Motivated by the long-standing "luminosity problem" in low-mass star formation whereby protostars are underluminous compared to theoretical expectations, we identify 230 protostars in 18 molecular clouds observed by two Spitzer Space Telescope Legacy surveys of nearby star-forming regions. We compile complete spectral energy distributions, calculate Lbol for each source, and study the protostellar luminosity distribution. This distribution extends over three orders of magnitude, from 0.01 Lsun - 69 Lsun, and has a mean and median of 4.3 Lsun and 1.3 Lsun, respectively. The distributions are very similar for Class 0 and Class I sources except for an excess of low luminosity (Lbol < 0.5 Lsun) Class I sources compared to Class 0. 100 out of the 230 protostars (43%) lack any available data in the far-infrared and submillimeter (70 um < wavelength < 850 um) and have Lbol underestimated by factors of 2.5 on average, and up to factors of 8-10 in extreme cases. Correcting these underestimates for each source individually once additional data becomes available will likely increase both the mean and median of the sample by 35% - 40%. We discuss and compare our results to several recent theoretical studies of protostellar luminosities and show that our new results do not invalidate the conclusions of any of these studies. As these studies demonstrate that there is more than one plausible accretion scenario that can match observations, future attention is clearly needed. The better statistics provided by our increased dataset should aid such future work.

preprint2013arXiv

The molecular ISM in luminous infrared galaxies: a 3mm line survey of Arp 157

In this work we present a search for molecular line emission in the interacting luminous infrared galaxy Arp 157 (NGC 520). We report the detection of 19 different lines, from 12 different molecular species in our 3mm line survey and supporting 2mm and 1mm observations. We investigated the profiles of the detected lines, and find all species show asymmetric emission lines, with profiles that differ systematically as a function of critical density. We derive the physical conditions within the gas reservoir of this system using isotope ratios, an excitation analysis, and rotation diagrams. We suggest that in this source we are observing a warm (Tex > 23K) edge-on, asymmetric starburst ring, plus a less dense central molecular component, which is more optically thin. We compare our molecular abundances with several well studied galaxies, and with Galactic regions, in order to determine the primary driver of the chemistry in Arp 157. We highlight some significant differences between the abundances of molecules in this source and in the more extreme ultra-luminous infrared galaxy Arp 220. This suggests that different physical mechanisms may become important in driving gas chemistry at different points in the evolution of a merger induced starburst. Of the nearby galaxies we considered, NGC 253 (whose nuclear chemistry is dominated by large-scale shocks) and NGC 1068 (whose nuclear chemistry is X-ray dominated) have the most similar molecular abundances to Arp 157. The Galactic region Sgr B2(OH) (where shock desorption of ices is thought to be important) may also have similar abundances to Arp 157. We hence postulate that shocks and/or the generation of hard radiation fields by the ongoing merger (or a starburst wind), are driving the molecular chemistry in Arp 157.

preprint2013arXiv

The Star Formation Rate of Molecular Clouds

We review recent advances in the analytical and numerical modeling of the star formation rate in molecular clouds and discuss the available observational constraints. We focus on molecular clouds as the fundamental star formation sites, rather than on the larger-scale processes that form the clouds and set their properties. Molecular clouds are shaped into a complex filamentary structure by supersonic turbulence, with only a small fraction of the cloud mass channeled into collapsing protostars over a free-fall time of the system. In recent years, the physics of supersonic turbulence has been widely explored with computer simulations, leading to statistical models of this fragmentation process, and to the prediction of the star formation rate as a function of fundamental physical parameters of molecular clouds, such as the virial parameter, the rms Mach number, the compressive fraction of the turbulence driver, and the ratio of gas to magnetic pressure. Infrared space telescopes, as well as ground-based observatories have provided unprecedented probes of the filamentary structure of molecular clouds and the location of forming stars within them.

preprint2013arXiv

The VIRUS-P Exploration of Nearby Galaxies (VENGA): Survey Design, Data Processing, and Spectral Analysis Methods

We present the survey design, data reduction, and spectral fitting pipeline for the VIRUS-P Exploration of Nearby Galaxies (VENGA). VENGA is an integral field spectroscopic survey, which maps the disks of 30 nearby spiral galaxies. Targets span a wide range in Hubble type, star formation activity, morphology, and inclination. The VENGA data-cubes have 5.6'' FWHM spatial resolution, ~5A FWHM spectral resolution, sample the 3600A-6800A range, and cover large areas typically sampling galaxies out to ~0.7 R_25. These data-cubes can be used to produce 2D maps of the star formation rate, dust extinction, electron density, stellar population parameters, the kinematics and chemical abundances of both stars and ionized gas, and other physical quantities derived from the fitting of the stellar spectrum and the measurement of nebular emission lines. To exemplify our methods and the quality of the data, we present the VENGA data-cube on the face-on Sc galaxy NGC 628 (a.k.a. M 74). The VENGA observations of NGC 628 are described, as well as the construction of the data-cube, our spectral fitting method, and the fitting of the stellar and ionized gas velocity fields. We also propose a new method to measure the inclination of nearly face-on systems based on the matching of the stellar and gas rotation curves using asymmetric drift corrections. VENGA will measure relevant physical parameters across different environments within these galaxies, allowing a series of studies on star formation, structure assembly, stellar populations, chemical evolution, galactic feedback, nuclear activity, and the properties of the interstellar medium in massive disk galaxies.

preprint2013arXiv

Warm gas towards young stellar objects in Corona Australis - Herschel/PACS observations from the DIGIT key programme

The effects of external irradiation on the chemistry and physics in the protostellar envelope around low-mass young stellar objects are poorly understood. The Corona Australis star-forming region contains the R CrA dark cloud, comprising several low-mass protostellar cores irradiated by an intermediate-mass young star. We study the effects on the warm gas and dust in a group of low-mass young stellar objects from the irradiation by the young luminous Herbig Be star R CrA. Herschel/PACS far-infrared datacubes of two low-mass star-forming regions in the R CrA dark cloud are presented. The distribution of CO, OH, H2O, [C II], [O I], and continuum emission is investigated. We have developed a deconvolution algorithm which we use to deconvolve the maps, separating the point-source emission from the extended emission. We also construct rotational diagrams of the molecular species. By deconvolution of the Herschel data, we find large-scale (several thousand AU) dust continuum and spectral line emission not associated with the point sources. Similar rotational temperatures are found for the warm CO ($282\pm4$ K), hot CO ($890\pm84$ K), OH ($79\pm4$ K), and H2O ($197\pm7$ K) emission, respectively, in the point sources and the extended emission. The rotational temperatures are also similar to what is found in other more isolated cores. The extended dust continuum emission is found in two ridges similar in extent and temperature to molecular mm emission, indicative of external heating from the Herbig Be star R CrA. Our results show that a nearby luminous star does not increase the molecular excitation temperatures in the warm gas around a young stellar object (YSO). However, the emission from photodissociation products of H2O, such as OH and O, is enhanced in the warm gas associated with these protostars and their surroundings compared to similar objects not suffering from external irradiation.

preprint2012arXiv

A Herschel Survey of Cold Dust in Disks Around Brown Dwarfs and Low-Mass Stars

We report the complete photometric results from our Herschel study which is the first comprehensive program to search for far-infrared emission from cold dust around young brown dwarfs. We surveyed 50 fields containing 51 known or suspected brown dwarfs and very low mass stars that have evidence of circumstellar disks based on Spitzer photometry and/or spectroscopy. The objects with known spectral types range from M3 to M9.5. Four of the candidates were subsequently identified as extragalactic objects. Of the remaining 47 we have successfully detected 36 at 70micron and 14 at 160micron with S/N greater than 3, as well as several additional possible detections with low S/N. The objects exhibit a range of [24]--[70] micron colors suggesting a range in mass and/or structure of the outer disk. We present modeling of the spectral energy distributions of the sample and discuss trends visible in the data. Using two Monte Carlo radiative transfer codes we investigate disk masses and geometry. We find a very wide range in modeled total disk masses from less than 1e-6 solar masses up to 1e-3 solar masses with a median disk mass of order 3e-5 solar masses, suggesting that the median ratio of disk mass to central object mass may be lower than for T Tauri stars. The disk scale heights and flaring angles, however, cover a range consistent with those seen around T Tauri stars. The host clouds in which the young brown dwarfs and low-mass stars are located span a range in estimated age from ~1-3 Myr to ~10 Myr and represent a variety of star-forming environments. No obvious dependence on cloud location or age is seen in the disk properties, though the statistical significance of this conclusion is not strong.

preprint2012arXiv

CO2 Ice toward Low-luminosity, Embedded Protostars: Evidence for Episodic Mass Accretion via Chemical History

We present Spitzer IRS spectroscopy of CO2 ice bending mode spectra at 15.2 micrometer toward 19 young stellar objects with luminosity lower than 1 Lsun (3 with luminosity lower than 0.1 Lsun). Ice on dust grain surfaces can encode the history of heating because pure CO2 ice forms only at elevated temperature, T > 20 K, and thus around protostars of higher luminosity. Current internal luminosities of YSOs with L < 1 Lsun do not provide the conditions needed to produce pure CO2 ice at radii where typical envelopes begin. The presence of detectable amounts of pure CO2 ice would signify a higher past luminosity. Many of the spectra require a contribution from a pure, crystalline CO2 component, traced by the presence of a characteristic band splitting in the 15.2 micrometer bending mode. About half of the sources (9 out of 19) in the low luminosity sample have evidence for pure CO2 ice, and six of these have significant double-peaked features, which are very strong evidence of pure CO2 ice. The presence of the pure CO2 ice component indicates that the dust temperature, and hence luminosity of the central star/accretion disk system, must have been higher in the past. An episodic accretion scenario, in which mixed CO-CO2 ice is converted to pure CO2 ice during each high luminosity phase, explains the presence of pure CO2 ice, the total amount of CO2 ice, and the observed residual C18O gas.

preprint2012arXiv

Fragmentation and Evolution of Molecular Clouds. III: The Effect of Dust and Gas Energetics

Dust and gas energetics are incorporated into a cluster-scale simulation of star formation in order to study the effect of heating and cooling on the star formation process. We build on our previous work by calculating separately the dust and gas temperatures. The dust temperature is set by radiative equilibrium between heating by embedded stars and radiation from dust. The gas temperature is determined using an energy-rate balance algorithm which includes molecular cooling, dust-gas collisional energy transfer, and cosmic-ray ionization. The fragmentation proceeds roughly similarly to simulations in which the gas temperature is set to the dust temperature, but there are differences. The structure of regions around sink particles have properties similar to those of Class 0 objects, but the infall speeds and mass accretion rates were, on average, higher than those seen for regions forming only low-mass stars. The gas and dust temperature have complex distributions not well modeled by approximations that ignore the detailed thermal physics. There is no simple relationship between density and kinetic temperature. In particular, high density regions have a large range of temperatures, determined by their location relative to heating sources. The total luminosity underestimates the star formation rate at these early stages, before ionizing sources are included, by an order of magnitude. As predicted in our previous work, a larger number of intermediate mass objects form when improved thermal physics is included, but the resulting IMF still has too few low mass stars. However, if we consider recent evidence on core-to-star efficiencies, the match to the IMF is improved.

preprint2012arXiv

Star Formation in the Milky Way and Nearby Galaxies

We review progress over the past decade in observations of large-scale star formation, with a focus on the interface between extragalactic and Galactic studies. Methods of measuring gas contents and star formation rates are discussed, and updated prescriptions for calculating star formation rates are provided. We review relations between star formation and gas on scales ranging from entire galaxies to individual molecular clouds.

preprint2012arXiv

Submillimeter Follow-up of WISE-Selected Hyperluminous Galaxies

We have used the Caltech Submillimeter Observatory (CSO) to follow-up a sample of WISE-selected, hyperluminous galaxies, so called W1W2-dropout galaxies. This is a rare (~ 1000 all-sky) population of galaxies at high redshift (peaks at z=2-3), that are faint or undetected by WISE at 3.4 and 4.6 um, yet are clearly detected at 12 and 22 um. The optical spectra of most of these galaxies show significant AGN activity. We observed 14 high-redshift (z > 1.7) W1W2-dropout galaxies with SHARC-II at 350 to 850 um, with 9 detections; and observed 18 with Bolocam at 1.1 mm, with five detections. Warm Spitzer follow-up of 25 targets at 3.6 and 4.5 um, as well as optical spectra of 12 targets are also presented in the paper. Combining WISE data with observations from warm Spitzer and CSO, we constructed their mid-IR to millimeter spectral energy distributions (SEDs). These SEDs have a consistent shape, showing significantly higher mid-IR to submm ratios than other galaxy templates, suggesting a hotter dust temperature. We estimate their dust temperatures to be 60-120 K using a single-temperature model. Their infrared luminosities are well over 10^{13} Lsun. These SEDs are not well fitted with existing galaxy templates, suggesting they are a new population with very high luminosity and hot dust. They are likely among the most luminous galaxies in the Universe. We argue that they are extreme cases of luminous, hot dust-obscured galaxies (DOGs), possibly representing a short evolutionary phase during galaxy merging and evolution. A better understanding of their long-wavelength properties needs ALMA as well as Herschel data.

preprint2012arXiv

The First Hyper-Luminous Infrared Galaxy Discovered by WISE

We report the discovery by the Wide-field Infrared Survey Explorer of the z = 2.452 source WISE J181417.29+341224.9, the first hyperluminous source found in the WISE survey. WISE 1814+3412 is also the prototype for an all-sky sample of ~1000 extremely luminous "W1W2-dropouts" (sources faint or undetected by WISE at 3.4 and 4.6 microns and well detected at 12 or 22 microns). The WISE data and a 350 micron detection give a minimum bolometric luminosity of 3.7 x 10^13 Lsun, with ~10^14 Lsun plausible. Followup images reveal four nearby sources: a QSO and two Lyman Break Galaxies (LBGs) at z = 2.45, and an M dwarf star. The brighter LBG dominates the bolometric emission. Gravitational lensing is unlikely given the source locations and their different spectra and colors. The dominant LBG spectrum indicates a star formation rate ~300 Msun/yr, accounting for < 10% of the bolometric luminosity. Strong 22 micron emission relative to 350 microns implies that warm dust contributes significantly to the luminosity, while cooler dust normally associated with starbursts is constrained by an upper limit at 1.1 mm. Radio emission is ~10x above the far-infrared/radio correlation, indicating an active galactic nucleus is present. An obscured AGN combined with starburst and evolved stellar components can account for the observations. If the black hole mass follows the local M_BH-bulge mass relation, the implied Eddington ratio is >~4. WISE 1814+3412 may be a heavily obscured object where the peak AGN activity occurred prior to the peak era of star formation.

preprint2012arXiv

The Herschel DIGIT Survey of Weak-line T Tauri Stars: implications for disk evolution and dissipation

As part of the "Dust, Ice, and Gas In Time (DIGIT)" Herschel Open Time Key Program, we present Herschel photometry (at 70, 160, 250, 350 and 500 micron) of 31 Weak-Line T Tauri star (WTTS) candidates in order to investigate the evolutionary status of their circumstellar disks. Thirteen stars in our sample had circumstellar disks previously known from infrared observations at shorter wavelengths, while eighteen of them had no previous evidence for a disk. We detect a total of 15 disks as all previously known disks are detected at one or more Herschel wavelengths and two additional disks are identified for the first time. The spectral energy distributions (SEDs) of our targets seem to trace the dissipation of the primordial disk and the transition to the debris disk regime. Seven of the 15 disks appear to be optically thick primordial disks, including two objects with SEDs indistinguishable from those of typical Classical T Tauri stars, four objects that have significant deficit of excess emission at all IR wavelengths, and one "pre-transitional" object with a known gap in the disk. Despite their previous WTTS classification, we find that the seven targets in our sample with optically thick disks show evidence for accretion. The remaining eight disks have weaker IR excesses similar to those of optically thin debris disks. Six of them are warm and show significant 24 micron Spitzer excesses, while the last two are newly identified cold debris-like disks with photospheric 24 micron fluxes, but significant excess emission at longer wavelengths. The Herschel photometry also places strong constraints on the non-detections, where systems with F70/F70,star > 5 - 15 and L,disk/L,star > 1xE-3 to 1xE-4 can be ruled out. We present preliminary models for both the optically thick and optically thin disks and discuss our results in the context of the evolution and dissipation of circumstellar disks.

preprint2012arXiv

The VIRUS-P Exploration of Nearby Galaxies (VENGA): The Xco Gradient in NGC 628

We measure the radial profile of the 12CO(1-0) to H_2 conversion factor (Xco) in NGC 628. The Hαemission from the VENGA integral field spectroscopy is used to map the star formation rate surface density (Σ_{SFR}). We estimate the molecular gas surface density (Σ_{H2}) from Σ_{SFR} by inverting the molecular star formation law (SFL), and compare it to the CO intensity to measure Xco. We study the impact of systematic uncertainties by changing the slope of the SFL, using different SFR tracers (Hαvs. far-UV plus 24μm), and CO maps from different telescopes (single-dish and interferometers). The observed Xco profile is robust against these systematics, drops by a factor of 2 from R~7 kpc to the center of the galaxy, and is well fit by a gradient Δlog(Xco)=0.06\pm0.02 dex kpc^-1. We study how changes in Xco follow changes in metallicity, gas density, and ionization parameter. Theoretical models show that the gradient in Xco can be explained by a combination of decreasing metallicity, and decreasing Σ_{H2} with radius. Photoelectric heating from the local UV radiation field appears to contribute to the decrease of Xco in higher density regions. Our results show that galactic environment plays an important role at setting the physical conditions in star forming regions, in particular the chemistry of carbon in molecular complexes, and the radiative transfer of CO emission. We caution against adopting a single Xco value when large changes in gas surface density or metallicity are present.

preprint2011arXiv

A Herschel Search For Cold Dust in Brown Dwarf Disks: First Results

We report initial results from a {\it Herschel} program to search for far-infrared emission from cold dust around a statistically significant sample of young brown dwarfs. The first three objects in our survey are all detected at 70\micron, and we report the first detection of a brown dwarf at 160\micron. The flux densities are consistent with the presence of substantial amounts of cold dust in the outer disks around these objects. We modeled the SED's with two different radiative transfer codes. We find that a broad range of model parameters provides a reasonable fit to the SED's, but that the addition of our 70\micron, and especially the 160\micron\ detection enables strong lower limits to be placed on the disk masses since most of the mass is in the outer disk. We find likely disk masses in the range of a few $\times 10^{-6}$ to $10^{-4}$ \msun. Our models provide a good fit to the SED's and do not require dust settling.

preprint2011arXiv

A Mid-Infrared Census of Star Formation Activity in Bolocam Galactic Plane Survey Sources

We present the results of a search for mid-infrared signs of star formation activity in the 1.1 mm sources in the Bolocam Galactic Plane Survey (BGPS). We have correlated the BGPS catalog with available mid-IR Galactic plane catalogs based on the Spitzer Space Telescope GLIMPSE legacy survey and the Midcourse Space Experiment (MSX) Galactic plane survey. We find that 44% (3,712 of 8,358) of the BGPS sources contain at least one mid-IR source, including 2,457 of 5,067 (49%) within the area where all surveys overlap (10 deg < l < 65 deg). Accounting for chance alignments between the BGPS and mid-IR sources, we conservatively estimate that 20% of the BPGS sources within the area where all surveys overlap show signs of active star formation. We separate the BGPS sources into four groups based on their probability of star formation activity. Extended Green Objects (EGOs) and Red MSX Sources (RMS) make up the highest probability group, while the lowest probability group is comprised of "starless" BGPS sources which were not matched to any mid-IR sources. The mean 1.1 mm flux of each group increases with increasing probability of active star formation. We also find that the "starless" BGPS sources are the most compact, while the sources with the highest probability of star formation activity are on average more extended with large skirts of emission. A subsample of 280 BGPS sources with known distances demonstrates that mass and mean H_2 column density also increase with probability of star formation activity.

preprint2011arXiv

Disentangling the Environment of the FU Orionis Candidate HBC 722 with Herschel

We analyze the submillimeter emission surrounding the new FU Orionis-type object, HBC 722. We present the first epoch of observations of the active environs of HBC 722, with imaging and spectroscopy from PACS, SPIRE, and HIFI aboard the Herschel Space Observatory, as well as CO J= 2-1 and 350 um imaging (SHARC-II) with the Caltech Submillimeter Observatory. The primary source of submillimeter continuum emission in the region -- 2MASS 20581767+4353310 -- is located 16$\arcsec$ south-southeast of the optical flaring source while the optical and near-IR emission is dominated by HBC 722. A bipolar outflow extends over HBC 722; the most likely driver is the submillimeter source. We detect warm (100 K) and hot (246 K) CO emission in the surrounding region, evidence of outflow-driven heating in the vicinity. The region around HBC 722 itself shows little evidence of heating driven by the new outbursting source itself.

preprint2011arXiv

Herschel Observations of the T Cha Transition Disk: Constraining the Outer Disk Properties

T Cha is a nearby (d = 100 pc) transition disk known to have an optically thin gap separating optically thick inner and outer disk components. Huelamo et al. (2011) recently reported the presence of a low-mass object candidate within the gap of the T Cha disk, giving credence to the suspected planetary origin of this gap. Here we present the Herschel photometry (70, 160, 250, 350, and 500 micron) of T Cha from the "Dust, Ice, and Gas in Time" (DIGIT) Key Program, which bridges the wavelength range between existing Spitzer and millimeter data and provide important constraints on the outer disk properties of this extraordinary system. We model the entire optical to millimeter wavelength spectral energy distribution (SED) of T Cha (19 data points between 0.36 and 3300 micron without any major gaps in wavelength coverage). T Cha shows a steep spectral slope in the far-IR, which we find clearly favors models with outer disks containing little or no dust beyond 40 AU. The full SED can be modeled equally well with either an outer disk that is very compact (only a few AU wide) or a much larger one that has a very steep surface density profile. That is, T Cha's outer disk seems to be either very small or very tenuous. Both scenarios suggest a highly unusual outer disk and have important but different implications for the nature of T Cha. Spatially resolved images are needed to distinguish between the two scenarios.

preprint2011arXiv

Ices in starless and starforming cores

Icy grain mantles are commonly observed through infrared spectroscopy toward dense clouds, cloud cores, protostellar envelopes and protoplanetary disks. Up to 80% of the available oxygen, carbon and nitrogen are found in such ices; the most common ice constituents - H2O, CO2 and CO - are second in abundance only to H2 in many star forming regions. In addition to being a molecular reservoir, ice chemistry is responsible for much of the chemical evolution from H2O to complex, prebiotic molecules. Combining the existing ISO, Spitzer, VLT and Keck ice data results in a large sample of ice sources (\sime80) that span all stages of star formation and a large range of protostellar luminosities (<0.1-105 L\odot). Here we summarize the different techniques that have been applied to mine this ice data set on information on typical ice compositions in different environments and what this implies about how ices form and evolve during star and planet formation. The focus is on how to maximize the use of empirical constraints from ice observations, followed by the application of information from experiments and models. This strategy is used to identify ice bands and to constrain which ices form early during cloud formation, which form later in the prestellar core and which require protostellar heat and/or UV radiation to form. The utility of statistical tests, survival analysis and ice maps is highlighted; the latter directly reveals that the prestellar ice formation takes place in two phases, associated with H2O and CO ice formation, respectively, and that most protostellar ice variation can be explained by differences in the prestellar CO ice formation stage. Finally, special attention is paid to the difficulty of observing complex ices directly and how gas observations, experiments and models help in constraining this ice chemistry stage.

preprint2011arXiv

The Bolocam Galactic Plane Survey V: HCO+ and N2H+ Spectroscopy of 1.1 mm Dust Continuum Sources

We present the results of observations of 1882 sources in the Bolocam Galactic Plane Survey (BGPS) at 1.1 mm with the 10m Heinrich Hertz Telescope simultaneously in HCO+ J=3-2 and N2H+ J=3-2. We detect 77% of these sources in HCO^+ and 51% in N2H+ at greater than 3$σ$. We find a strong correlation between the integrated intensity of both dense gas tracers and the 1.1 mm dust emission of BGPS sources. We determine kinematic distances for 529 sources (440 in the first quadrant breaking the distance ambiguity and 89 in the second quadrant) We derive the size, mass, and average density for this subset of clumps. The median size of BGPS clumps is 0.75 pc with a median mass of 330 M$_{\odot}$ (assuming T_{Dust}=20 K). The median HCO+ linewidth is 2.9 km s$^{-1}$ indicating that BGPS clumps are dominated by supersonic turbulence or unresolved kinematic motions. We find no evidence for a size-linewidth relationship for BGPS clumps. We analyze the effects of the assumed dust temperature on the derived clump properties with a Monte Carlo simulation and we find that changing the temperature distribution will change the median source properties (mass, volume-averaged number density, surface density) by factors of a few. The observed differential mass distribution has a power-law slope that is intermediate between that observed for diffuse CO clouds and the stellar IMF. BGPS clumps represent a wide range of objects (from dense cores to more diffuse clumps) and are typically characterized by larger sizes and lower densities than previously published surveys of high-mass star-forming regions. This collection of objects is a less-biased sample of star-forming regions in the Milky Way that likely span a wide range of evolutionary states.

preprint2011arXiv

The Bolocam Galactic Plane Survey. VII. Characterizing the Properties of Massive Star-Forming Regions

We present the results of a GBT survey of NH3(1,1), (2,2), (3,3) lines towards 631 Bolocam Galactic Plane Survey (BGPS) sources at a range of Galactic longitudes in the inner Galaxy. We have detected the NH3(1,1) line towards 72% of our targets (456), demonstrating that the high column density features identified in the BGPS and other continuum surveys accurately predict the presence of dense gas. We have determined kinematic distances and resolved the distance ambiguity for all BGPS sources detected in NH3. The BGPS sources trace the locations of the Scutum and Sagittarius spiral arms, with the number of sources peaking between Galactocentric radii of 4-5 kpc. We measure the physical properties of each source and find that depending on the distance, BGPS sources are primarily clumps, with some cores and clouds. We have examined the physical properties as a function of Galactocentric distance, and find a mean gas kinetic temperature of 15.6 K, and that the NH3 column density and abundance decrease by nearly an order of magnitude between Galactocentric radii of 3 - 11 kpc. Comparing sources at similar distances demonstrates that the physical properties are indistinguishable, which suggests a similarity in clump structure across the Galactic disk. We have also compared the BGPS sources to criteria for efficient star formation presented independently by Heiderman et al. and Lada et al., and massive star formation presented by Kauffmann et al. 48% of our sample should be forming stars (including massive stars) with high efficiency, and 87% contain subregions that should be efficiently forming stars. Indeed, we find that 67% of the sample exhibit signs of star formation activity based on an association with a mid-IR source.

preprint2011arXiv

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

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

preprint2011arXiv

The Spitzer ice legacy: Ice evolution from cores to protostars

Ices regulate much of the chemistry during star formation and account for up to 80% of the available oxygen and carbon. In this paper, we use the Spitzer c2d ice survey, complimented with data sets on ices in cloud cores and high-mass protostars, to determine standard ice abundances and to present a coherent picture of the evolution of ices during low- and high-mass star formation. The median ice composition H2O:CO:CO2:CH3OH:NH3:CH4:XCN is 100:29:29:3:5:5:0.3 and 100:13:13:4:5:2:0.6 toward low- and high-mass protostars, respectively, and 100:31:38:4:-:-:- in cloud cores. In the low-mass sample, the ice abundances with respect to H2O of CH4, NH3, and the component of CO2 mixed with H2O typically vary by <25%, indicative of co-formation with H2O. In contrast, some CO and CO2 ice components, XCN and CH3OH vary by factors 2-10 between the lower and upper quartile. The XCN band correlates with CO, consistent with its OCN- identification. The origin(s) of the different levels of ice abundance variations are constrained by comparing ice inventories toward different types of protostars and background stars, through ice mapping, analysis of cloud-to-cloud variations, and ice (anti-)correlations. Based on the analysis, the first ice formation phase is driven by hydrogenation of atoms, which results in a H2O-dominated ice. At later prestellar times, CO freezes out and variations in CO freeze-out levels and the subsequent CO-based chemistry can explain most of the observed ice abundance variations. The last important ice evolution stage is thermal and UV processing around protostars, resulting in CO desorption, ice segregation and formation of complex organic molecules. The distribution of cometary ice abundances are consistent with with the idea that most cometary ices have a protostellar origin.

preprint2011arXiv

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

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

preprint2011arXiv

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

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

preprint2010arXiv

A Spitzer c2d Legacy Survey to Identify and Characterize Disks with Inner Dust Holes

Understanding how disks dissipate is essential to studies of planet formation. However, identifying exactly how dust and gas dissipates is complicated due to difficulty in finding objects clearly in the transition of losing their surrounding material. We use Spitzer IRS spectra to examine 35 photometrically-selected candidate cold disks (disks with large inner dust holes). The infrared spectra are supplemented with optical spectra to determine stellar and accretion properties and 1.3mm photometry to measure disk masses. Based on detailed SED modeling, we identify 15 new cold disks. The remaining 20 objects have IRS spectra that are consistent with disks without holes, disks that are observed close to edge-on, or stars with background emission. Based on these results, we determine reliable criteria for identifying disks with inner holes from Spitzer photometry and examine criteria already in the literature. Applying these criteria to the c2d surveyed star-forming regions gives a frequency of such objects of at least 4% and most likely of order 12% of the YSO population identified by Spitzer. We also examine the properties of these new cold disks in combination with cold disks from the literature. Hole sizes in this sample are generally smaller than for previously discovered disks and reflect a distribution in better agreement with exoplanet orbit radii. We find correlations between hole size and both disk and stellar masses. Silicate features, including crystalline features, are present in the overwhelming majority of the sample although 10 micron feature strength above the continuum declines for holes with radii larger than ~7 AU. In contrast, PAHs are only detected in 2 out of 15 sources. Only a quarter of the cold disk sample shows no signs of accretion, making it unlikely that photoevaporation is the dominant hole forming process in most cases.

preprint2010arXiv

c2d Spitzer IRS spectra of embedded low-mass young stars: gas-phase emission lines

A survey of mid-IR gas-phase emission lines of H2, H2O and various atoms toward a sample of 43 embedded low-mass young stars in nearby star-forming regions is presented. The sources are selected from the Spitzer "Cores to Disks" (c2d) legacy program. The environment of embedded protostars is complex both in its physical structure (envelopes, outflows, jets, protostellar disks) and the physical processes (accretion, irradiation by UV and/or X-rays, excitation through slow and fast shocks) which take place. A key point is to spatially resolve the emission in the Spitzer-IRS spectra. An optimal extraction method is used to separate both spatially unresolved (compact, up to a few 100 AU) and spatially resolved (extended, 1000 AU or more) emission from the IRS spectra. The results are compared with the c2d disk sample and literature PDR and shock models to address the physical nature of the sources. Both compact and extended emission features are observed. Warm (Tex few 100 K) H2, observed through the pure rotational H2 S(0), S(1) and S(2) lines, and [S I] 25 mu emission is observed primarily in the extended component. [S I] is observed uniquely toward truly embedded sources and not toward disks. On the other hand hot (Tex>700 K) H2, observed primarily through the S(4) line, and [Ne II] emission is seen mostly in the compact component. [Fe II] and [Si II] lines are observed in both spatial components. Hot H2O emission is found in the compact component of some sources. The observed emission on >=1000 AU scales is characteristic of PDR emission and likely originates in the outflow cavities in the remnant envelope created by the stellar wind and jets from the embedded young stars. Weak shocks along the outflow wall can also contribute. The compact emission is likely of mixed origin, comprised of optically thick circumstellar disk and/or jet/outflow emission from the protostellar object.

preprint2010arXiv

Low-mass Star Formation: Observations

I briefly review recent observations of regions forming low mass stars. The discussion is cast in the form of seven questions that have been partially answered, or at least illuminated, by new data. These are the following: where do stars form in molecular clouds; what determines the IMF; how long do the steps of the process take; how efficient is star formation; do any theories explain the data; how are the star and disk built over time; and what chemical changes accompany star and planet formation. I close with a summary and list of open questions.

preprint2010arXiv

Observational Constraints on Submillimeter Dust Opacity

Infrared extinction maps and submillimeter dust continuum maps are powerful probes of the density structure in the envelope of star-forming cores. We make a direct comparison between infrared and submillimeter dust continuum observations of the low-mass Class 0 core, B335, to constrain the ratio of submillimeter to infrared opacity (\kaprat) and the submillimeter opacity power-law index ($κ\propto λ^{-β}$). Using the average value of theoretical dust opacity models at 2.2 \micron, we constrain the dust opacity at 850 and 450 \micron . Using new dust continuum models based upon the broken power-law density structure derived from interferometric observations of B335 and the infall model derived from molecular line observations of B335, we find that the opacity ratios are $\frac{κ_{850}}{κ_{2.2}} = (3.21 - 4.80)^{+0.44}_{-0.30} \times 10^{-4}$ and $\frac{κ_{450}}{κ_{2.2}} = (12.8 - 24.8)^{+2.4}_{-1.3} \times 10^{-4}$ with a submillimeter opacity power-law index of $β_{smm} = (2.18 - 2.58)^{+0.30}_{-0.30}$. The range of quoted values are determined from the uncertainty in the physical model for B335. For an average 2.2 \micron\ opacity of $3800 \pm 700$ cm$^2$g$^{-1}$, we find a dust opacity at 850 and 450 \micron\ of $κ_{850} = (1.18 - 1.77)^{+0.36}_{-0.24}$ and $κ_{450} = (4.72 - 9.13)^{+1.9}_{-0.98}$ cm$^2$g$^{-1}$ of dust. These opacities are from $(0.65 - 0.97) κ^{\rm{OH}5}_{850}$ of the widely used theoretical opacities of Ossenkopf and Henning for coagulated ice grains with thin mantles at 850\micron.

preprint2010arXiv

The Bolocam Galactic Plane Survey -- III. Characterizing Physical Properties of Massive Star-Forming Regions in the Gemini OB1 Molecular Cloud

We present the 1.1 millimeter Bolocam Galactic Plane Survey (BGPS) observations of the Gemini OB1 molecular cloud complex, and targeted ammonia observations of the BGPS sources. When paired with molecular spectroscopy of a dense gas tracer, millimeter observations yield physical properties such as masses, radii, mean densities, kinetic temperatures and line widths. We detect 34 distinct BGPS sources above 5-sigma=0.37 Jy/beam with corresponding 5-sigma detections in the ammonia (1,1) transition. Eight of the objects show water maser emission (20%). We find a mean millimeter source FWHM of 1.12 pc, and a mean kinetic temperature of 20 K for the sample of 34 BGPS sources. The observed ammonia line widths are dominated by non-thermal motions, typically found to be a few times the thermal sound speed expected for the derived kinetic temperature. We calculate the mass for each source from the millimeter flux assuming the sources are isothermal and find a mean isothermal mass within a 120" aperture of 230 +/- 180 solar masses. We find a total mass of 8,400 solar masses for all BGPS sources in the Gemini OB1 molecular cloud, representing 6.5% of the cloud mass. By comparing the millimeter isothermal mass to the virial mass within a radius equal to the mm source size calculated from the ammonia line widths, we find a mean virial parameter (M_vir/M_iso) of 1.0 +/- 0.9 for the sample. We find mean values for the distributions of column densities of 10^22 cm^-2 for H_2, and 3.0x10^14 cm^-2 for ammonia, giving a mean ammonia abundance of 3.0x10^-8 relative to H_2. We find volume-averaged densities on the order of 10^3-10^4 cm^-3. The sizes and densities suggest that in the Gem OB1 region the BGPS is detecting the clumps from which stellar clusters form, rather than smaller, higher density cores where single stars or small multiple systems form.

preprint2010arXiv

The Bolocam Galactic Plane Survey IV: 1.1 and 0.35 mm Dust Continuum Emission in the Galactic Center Region

The Bolocam Galactic Plane Survey (BGPS) data for a six square degree region of the Galactic plane containing the Galactic center is analyzed and compared to infrared and radio continuum data. The BGPS 1.1 mm emission consists of clumps interconnected by a network of fainter filaments surrounding cavities, a few of which are filled with diffuse near-IR emission indicating the presence of warm dust or with radio continuum characteristic of HII regions or supernova remnants. New 350 μm images of the environments of the two brightest regions, Sgr A and B, are presented. Sgr B2 is the brightest mm-emitting clump in the Central Molecular Zone and may be forming the closest analog to a super star cluster in the Galaxy. The Central Molecular Zone (CMZ) contains the highest concentration of mm and sub-mm emitting dense clumps in the Galaxy. Most 1.1 mm features at positive longitudes are seen in silhouette against the 3.6 to 24 μm background observed by the Spitzer Space Telescope. However, only a few clumps at negative longitudes are seen in absorption, confirming the hypothesis that positive longitude clumps in the CMZ tend to be on the near-side of the Galactic center, consistent with the suspected orientation of the central bar in our Galaxy. Some 1.1 mm cloud surfaces are seen in emission at 8 μm, presumably due to polycyclic aromatic hydrocarbons (PAHs). A ~0.2\degree (~30 pc) diameter cavity and infrared bubble between l \approx 0.0\degree and 0.2\degree surrounds the Arches and Quintuplet clusters and Sgr A. The bubble contains several clumpy dust filaments that point toward Sgr A\ast; its potential role in their formation is explored. [abstract truncated]

preprint2010arXiv

The Bolocam Galactic Plane Survey: Survey Description and Data Reduction

We present the Bolocam Galactic Plane Survey (BGPS), a 1.1 mm continuum survey at 33" effective resolution of 170 square degrees of the Galactic Plane visible from the northern hemisphere. The survey is contiguous over the range -10.5 < l < 90.5, |b| < 0.5 and encompasses 133 square degrees, including some extended regions |b| < 1.5. In addition to the contiguous region, four targeted regions in the outer Galaxy were observed: IC1396, a region towards the Perseus Arm, W3/4/5, and Gem OB1. The BGPS has detected approximately 8400 clumps over the entire area to a limiting non-uniform 1-sigma noise level in the range 11 to 53 mJy/beam in the inner Galaxy. The BGPS source catalog is presented in a companion paper (Rosolowsky et al. 2010). This paper details the survey observations and data reduction methods for the images. We discuss in detail the determination of astrometric and flux density calibration uncertainties and compare our results to the literature. Data processing algorithms that separate astronomical signals from time-variable atmospheric fluctuations in the data time-stream are presented. These algorithms reproduce the structure of the astronomical sky over a limited range of angular scales and produce artifacts in the vicinity of bright sources. Based on simulations, we find that extended emission on scales larger than about 5.9' is nearly completely attenuated (> 90%) and the linear scale at which the attenuation reaches 50% is 3.8'. Comparison with other millimeter-wave data sets implies a possible systematic offset in flux calibration, for which no cause has been discovered. This presentation serves as a companion and guide to the public data release through NASA's Infrared Processing and Analysis Center (IPAC) Infrared Science Archive (IRSA). New data releases will be provided through IPAC IRSA with any future improvements in the reduction.

preprint2010arXiv

The c2d Spitzer Spectroscopic Survey of Ices Around Low-Mass Young Stellar Objects. IV. NH3 and CH3OH

NH3 and CH3OH are key molecules in astrochemical networks leading to the formation of more complex N- and O-bearing molecules, such as CH3CN and HCOOCH3. Despite a number of recent studies, little is known about their abundances in the solid state. (...) In this work, we investigate the ~ 8-10 micron region in the Spitzer IRS (InfraRed Spectrograph) spectra of 41 low-mass young stellar objects (YSOs). These data are part of a survey of interstellar ices in a sample of low-mass YSOs studied in earlier papers in this series. We used both an empirical and a local continuum method to correct for the contribution from the 10 micron silicate absorption in the recorded spectra. In addition, we conducted a systematic laboratory study of NH3- and CH3OH-containing ices to help interpret the astronomical spectra. We clearly detect a feature at ~9 micron in 24 low-mass YSOs. Within the uncertainty in continuum determination, we identify this feature with the NH3 nu_2 umbrella mode, and derive abundances with respect to water between ~2 and 15%. Simultaneously, we also revisited the case of CH3OH ice by studying the nu_4 C-O stretch mode of this molecule at ~9.7 micron in 16 objects, yielding abundances consistent with those derived by Boogert et al. 2008 (hereafter paper I) based on a simultaneous 9.75 and 3.53 micron data analysis. Our study indicates that NH3 is present primarily in H2O-rich ices, but that in some cases, such ices are insufficient to explain the observed narrow FWHM. The laboratory data point to CH3OH being in an almost pure methanol ice, or mixed mainly with CO or CO2, consistent with its formation through hydrogenation on grains. Finally, we use our derived NH3 abundances in combination with previously published abundances of other solid N-bearing species to find that up to 10-20 % of nitrogen is locked up in known ices.

preprint2010arXiv

The Nature of Transition Circumstellar Disks I. The Ophiuchus Molecular Cloud

We have obtained millimeter wavelength photometry, high-resolution optical spectroscopy and adaptive optics near-infrared imaging for a sample of 26 Spitzer-selected transition circumstellar disks. All of our targets are located in the Ophiuchus molecular cloud (d ~125 pc) and have Spectral Energy Distributions (SEDs) suggesting the presence of inner opacity holes. We use these ground-based data to estimate the disk mass, multiplicity, and accretion rate for each object in our sample in order to investigate the mechanisms potentially responsible for their inner holes. We find that transition disks are a heterogeneous group of objects, with disk masses ranging from < 0.6 to 40 M_JUP and accretion rates ranging from <1E-11 1E-7 M_solar/yr, but most tend to have much lower masses and accretion rates than "full disks" (i.e., disks without opacity holes). Eight of our targets have stellar companions: 6 of them are binaries and the other 2 are triple systems. In four cases, the stellar companions are close enough to suspect they are responsible for the inferred inner holes. We find that 9 of our 26 targets have low disk mass (< 2.5 M_JUP) and negligible accretion (< 1E-11) M_solar/yr, and are thus consistent with photoevaporating (or photoevaporated) disks. Four of these 9 non-accreting objects have fractional disk luminosities < 1E-3 and could already be in a debris disk stage. Seventeen of our transition disks are accreting. Thirteen of these accreting objects are consistent with grain growth. The remaining 4 accreting objects have SEDs suggesting the presence of sharp inner holes, and thus are excellent candidates for harboring giant planets.

preprint2010arXiv

The Spatial Distribution of Star Formation in the Solar Neighbourhood: Do all stars form in clusters?

We present a global study of low mass, young stellar object (YSO) surface densities in nearby (< 500 pc) star forming regions based on a comprehensive collection of Spitzer Space Telescope surveys. We show that the distribution of YSO surface densities in the solar neighbourhood is a smooth distribution, being adequately described by a lognormal function from a few to 10^3 YSOs per pc^2, with a peak at 22 stars/pc^2 and a dispersion of 0.85. We do not find evidence for multiple discrete modes of star-formation (e.g. clustered and distributed). Comparing the observed surface density distribution to previously reported surface density threshold definitions of clusters, we find that the fraction of stars in clusters is crucially dependent on the adopted definitions, ranging from 40 to 90%. However, we find that only a low fraction (< 26%) of stars are formed in dense environments where their formation/evolution (along with their circumstellar disks and/or planets) may be affected by the close proximity of their low-mass neighbours.

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 Star Formation Rate and Gas Surface Density Relation in the Milky Way: Implications for Extragalactic Studies

We investigate the relation between star formation rate (SFR) and gas surface densities in Galactic star forming regions using a sample of YSOs and massive clumps. Our YSO sample consists of objects located in 20 molecular clouds from the Spitzer cores to disks and Gould's Belt surveys. We estimate the gas surface density (Sigma_gas) from Av maps and YSO SFR surface densities (Sigma_SFR) from the number of YSOs, assuming a mean mass and lifetime. We also divide the clouds into contour levels of Av, counting only the youngest Class I and Flat SED YSOs. For a sample of massive star forming clumps, we derive SFRs from the infrared luminosity and use HCN gas maps to estimate Sigma_gas. We find that Galactic clouds lie above the extragalactic relations (e.g., Kennicutt-Schmidt Law) by factors up to 17. Cloud regions with high Sigma_gas lie above extragalactic relations up to a factor of 54 and overlap with massive clumps. We use 12CO and 13CO gas maps of the Perseus and Ophiuchus clouds to estimate Sigma_gas and compare to Sigma_gas from Av maps. We find that 13CO, underestimates the Av-based mass by factors of 4-5. 12CO may underestimate the total gas mass at Sigma_gas > 200 Msun pc^-2 by > 30%;however, this does not explain the large discrepancy between Galactic and extragalactic relations. We find evidence for a threshold of star formation (Sigma_th) at 129+-14 Msun pc^-2. At Sigma_gas > Sigma_th, the Galactic relation is linear. A possible reason for the difference between Galactic and extragalactic relations is that all the CO-emitting gas, including Sigma_gas below Sigma_th, is measured in extragalactic studies. If the Kennicutt-Schmidt relation (Sigma_SFR Sigma_gas^1.4) and a linear relation between dense gas and star formation is assumed, the fraction of dense star forming gas (f_dense) increases as Sigma_gas^0.4. When Sigma_gas reaches ~300Sigma_th, f_dense is 1. (Abridged)

preprint2010arXiv

The VIRUS-P Exploration of Nearby Galaxies (VENGA): Survey Design and First Results

VENGA is a large-scale extragalactic IFU survey, which maps the bulges, bars and large parts of the outer disks of 32 nearby normal spiral galaxies. The targets are chosen to span a wide range in Hubble types, star formation activities, morphologies, and inclinations, at the same time of having vast available multi-wavelength coverage from the far-UV to the mid-IR, and available CO and 21cm mapping. The VENGA dataset will provide 2D maps of the SFR, stellar and gas kinematics, chemical abundances, ISM density and ionization states, dust extinction and stellar populations for these 32 galaxies. The uniqueness of the VIRUS-P large field of view permits these large-scale mappings to be performed. VENGA will allow us to correlate all these important quantities throughout the different environments present in galactic disks, allowing the conduction of a large number of studies in star formation, structure assembly, galactic feedback and ISM in galaxies.

preprint2009arXiv

Evolutionary Signatures in the Formation of Low-Mass Protostars. II. Towards Reconciling Models and Observations

A long-standing problem in low-mass star formation is the "luminosity problem," whereby protostars are underluminous compared to the accretion luminosity expected both from theoretical collapse calculations and arguments based on the minimum accretion rate necessary to form a star within the embedded phase duration. Motivated by this luminosity problem, we present a set of evolutionary models describing the collapse of low-mass, dense cores into protostars, using the Young & Evans (2005) model as our starting point. We calculate the radiative transfer of the collapsing cores throughout the full duration of the collapse in two dimensions. From the resulting spectral energy distributions, we calculate standard observational signatures to directly compare to observations. We incorporate several modifications and additions to the original Young & Evans model in an effort to better match observations with model predictions. We find that scattering, 2-D geometry, mass-loss, and outflow cavities all affect the model predictions, as expected, but none resolve the luminosity problem. A cycle of episodic mass accretion, however, can resolve this problem and bring the model predictions into better agreement with observations. Standard assumptions about the interplay between mass accretion and mass loss in our model give star formation efficiencies consistent with recent observations that compare the core mass function (CMF) and stellar initial mass function (IMF). The combination of outflow cavities and episodic mass accretion reduce the connection between observational Class and physical Stage to the point where neither of the two common observational signatures (bolometric temperature and ratio of bolometric to submillimeter luminosity) can be considered reliable indicators of physical Stage.

preprint2009arXiv

Far-infrared Observations of the Very Low-Luminosity Embedded Source L1521F-IRS in the Taurus Star-Forming Region

We investigate the environment of the very low-luminosity object L1521F IRS using data from the Taurus Spitzer Legacy Survey. The MIPS 160 micron image shows both extended emission from the Taurus cloud as well as emission from multiple cold cores over a 1 X 2 deg region. Analysis shows that the cloud dust temperature is 14.2 +- 0.4 K and the extinction ratio is A_160/A_K = 0.010 +- 0.001 up to A_V ~ 4 mag. We find kappa_160 = 0.23 +- 0.046 cm^2/g for the specific opacity of the gas-dust mixture. Therefore, for dust in the Taurus cloud we find the 160 um opacity is significantly higher than that measured for the diffuse ISM, but not too different from dense cores, even at modest extinction values. Furthermore, the 160 um image shows features that do not appear in the IRAS 100 um image. We identify six regions as cold cores, i.e. colder than 14.2 K, all of which have counterparts in extinction maps or C18O maps. We compare the effects of L1521F IRS on its natal core and find there is no evidence for dust heating at 160 or 100 um by the embedded source. From the infrared luminosity L_TIR = 0.024 Lo we find L_bol = 0.034 - 0.046L_o, thus confirming the source's low-luminosity. Comparison of L1521F IRS with theoretical simulations for the very early phases of star formation appears to rule out the first core collapse phase. The evolutionary state appears similar to or younger than the class 0 phase, and the estimated mass is likely to be substellar.

preprint2009arXiv

Fragmentation and Evolution of Molecular Clouds. II: The Effect of Dust Heating

We investigate the effect of heating by luminosity sources in a simulation of clustered star formation. Our heating method involves a simplified continuum radiative transfer method that calculates the dust temperature. The gas temperature is set by the dust temperature. We present the results of four simulations, two simulations assume an isothermal equation of state and the two other simulations include dust heating. We investigate two mass regimes, i.e., 84 Msun and 671 Msun, using these two different energetics algorithms. The mass functions for the isothermal simulations and simulations which include dust heating are drastically different. In the isothermal simulation, we do not form any objects with masses above 1 Msun. However, the simulation with dust heating, while missing some of the low-mass objects, forms high-mass objects (~20 Msun) which have a distribution similar to the Salpeter IMF. The envelope density profiles around the stars formed in our simulation match observed values around isolated, low-mass star-forming cores. We find the accretion rates to be highly variable and, on average, increasing with final stellar mass. By including radiative feedback from stars in a cluster-scale simulation, we have determined that it is a very important effect which drastically affects the mass function and yields important insights into the formation of massive stars.

preprint2009arXiv

The Bolocam Galactic Plane Survey -- II. Catalog of the Image Data

We present a catalog of 8358 sources extracted from images produced by the Bolocam Galactic Plane Survey (BGPS). The BGPS is a survey of the millimeter dust continuum emission from the northern Galactic plane. The catalog sources are extracted using a custom algorithm, Bolocat, which was designed specifically to identify and characterize objects in the large-area maps generated from the Bolocam instrument. The catalog products are designed to facilitate follow-up observations of these relatively unstudied objects. The catalog is 98% complete from 0.4 Jy to 60 Jy over all object sizes for which the survey is sensitive (<3.5'). We find that the sources extracted can best be described as molecular clumps -- large dense regions in molecular clouds linked to cluster formation. We find the flux density distribution of sources follows a power law with dN/dS ~S^(-2.4 +/- 0.1) and that the mean Galactic latitude for sources is significantly below the midplane: <b>=(-0.095 +/- 0.001) deg.

preprint2009arXiv

The Spatially Resolved Star Formation Law from Integral Field Spectroscopy: VIRUS-P Observations of NGC 5194

We investigate the relation between the star formation rate surface density (Sigma_SFR) and the mass surface density of gas (Sigma_gas) in NGC 5194. VIRUS-P integral field spectroscopy is used to measure H-alpha, H-beta [NII]6548,6584, and [SII]6717,6731 emission line fluxes for 735 regions ~170 pc in diameter, and derive extinction corrected Sigma_SFRs. HI 21cm and CO maps are used to measure the HI and H_2 gas surface density for each region. We present a new method for fitting the Star Formation Law (SFL), which includes the intrinsic scatter in the relation as a free parameter, allows the inclusion of non-detections, and is free of the systematics involved in performing linear correlations over incomplete data in logarithmic space. We use the [SII]/H-alpha ratio to separate the H-alpha flux contribution from the diffuse ionized gas (DIG). After removing the DIG, we measure a slope N=0.82+/-0.05, and an intrinsic scatter epsilon=0.43+/-0.02 dex for the molecular gas SFL. We also measure a typical depletion timescale tau~2 Gyr, in good agreement with Bigiel et al. (2008). The HI density shows no correlation with the SFR, and the total gas SFL closely follows the molecular gas SFL. We assess the validity of corrections applied in narrow-band H-alpha measurements to overcome issues related to continuum subtraction, underlying photospheric absorption, and contamination by the [NII] doublet. The disagreement with the super-linear molecular SFL measured by Kennicutt et al. (2007) is due to differences in the fitting method. Our results support a low and close to constant star formation efficiency (SFE = 1/tau) in the molecular ISM. The data shows excellent agreement with the model of Krumholz et al (2009). The large intrinsic scatter may imply the existence of other important parameters setting the SFR.

preprint2009arXiv

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

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

preprint2009arXiv

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

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

preprint2008arXiv

Properties of the Youngest Protostars in Perseus, Serpens, and Ophiuchus

We present an unbiased census of deeply embedded protostars in Perseus, Serpens, and Ophiuchus, assembled by combining large-scale 1.1 mm Bolocam continuum and Spitzer Legacy surveys. We identify protostellar candidates based on their mid-infrared properties, correlate their positions with 1.1 mm core positions, and construct well-sampled SEDs using our extensive wavelength coverage (lam=1.25-1100 micron). Source classification based on the bolometric temperature yields a total of 39 Class 0 and 89 Class I sources in the three cloud sample. We compare to protostellar evolutionary models using the bolometric temperature-luminosity diagram, finding a population of low luminosity Class I sources that are inconsistent with constant or monotonically decreasing mass accretion rates. This result argues strongly for episodic accretion during the Class I phase, with more than 50% of sources in a ``sub-Shu'' (dM/dt < 1e-6 Msun/yr) accretion state. Average spectra are compared to protostellar radiative transfer models, which match the observed spectra fairly well in Stage 0, but predict too much near-IR and too little mid-IR flux in Stage I. Finally, the relative number of Class 0 and Class I sources are used to estimate the lifetime of the Class 0 phase; the three cloud average yields a Class 0 lifetime of 1.7e5 yr, ruling out an extremely rapid early accretion phase. Correcting photometry for extinction results in a somewhat shorter lifetime (1.1e5 yr). In Ophiuchus, however, we find very few Class 0 sources (N(Class0)/N(ClassI)=0.1-0.2), similar to previous studies of that cloud. The observations suggest a consistent picture of nearly constant average accretion rate through the entire embedded phase, with accretion becoming episodic by at least the Class I stage, and possibly earlier.

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.

preprint2007arXiv

A Parameter Study of the Dust and Gas Temperature in a Field of Young Stars

We model the thermal effect of young stars on their surrounding environment in order to understand clustered star formation. We take radiative heating of dust, dust-gas collisional heating, cosmic-ray heating, and molecular cooling into account. Using Dusty, a spherical continuum radiative transfer code, we model the dust temperature distribution around young stellar objects with various luminosities and surrounding gas and dust density distributions. We have created a grid of dust temperature models, based on our modeling with Dusty, which we can use to calculate the dust temperature in a field of stars with various parameters. We then determine the gas temperature assuming energy balance. Our models can be used to make large-scale simulations of clustered star formation more realistic.

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.

preprint2006arXiv

SHARC-II Mapping of Spitzer c2d Small Clouds and Cores

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

preprint2006arXiv

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

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

preprint2006arXiv

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

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

preprint2006arXiv

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

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

preprint2006arXiv

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

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

preprint2002arXiv

The Physical Conditions for Massive Star Formation: Dust Continuum Maps and Modeling

Fifty-one dense cores associated with water masers were mapped at 350 micron. These cores are very luminous, 10^3 < Lbol/Lsun < 10^6, indicative of the formation of massive stars. Dust continuum contour maps and photometry are presented for these sources. The spectral energy distributions and normalized radial profiles of dust continuum emission were modeled for 31 sources using a one-dimensional dust radiative transfer code, assuming a power law density distribution in the envelope, n = n_f (r/r_f)^{-p}. The best fit density power law exponent, p, ranged from 0.75 to 2.5 with <p> = 1.8 +/- 0.4. The mean value of p is comparable to that found in regions forming only low mass stars. The mean p is incompatible with a logatropic sphere (p = 1), but other star formation models cannot be ruled out. Different mass estimates are compared and mean masses of gas and dust are reported within a half-power radius determined from the dust emission and within a radius where the total density exceeds 10^4 cm^3. Evolutionary indicators commonly used for low mass star formation may have some utility for regions forming massive stars. For comparison with extragalactic star formation studies, the luminosity to dust mass ratio is calculated for these sources with a method most parallel to that used in studies of distant galaxies and is found to be similar to that seen in high redshift starburst galaxies.

preprint2000arXiv

Tracing the Mass during Low-Mass Star Formation. I. Submillimeter Continuum Observations

We have obtained 850 and 450 micron continuum maps of 21 low mass cores with SED's ranging from pre-protostellar to Class I (18K < Tbol < 370K), using SCUBA at the JCMT. In this paper we present the maps, radial intensity profiles, and photometry. Pre-protostellar cores do not have power-law intensity profiles, whereas the intensity profiles of Class 0 and Class I sources can be fitted with power laws over a large range of radii. A substantial number of sources have companion sources within a few arcminutes (2 out of 5 pre-protostellar cores, 9 out of 16 Class 0/I sources). The mean separation between sources is 10800 AU. The median separation is 18000 AU including sources without companions as a lower limit. The mean value of the spectral index between 450 and 850 microns is 2.8 +- 0.4, with PPCs having slightly lower spectral indices (2.5 +- 0.4). The mean mass of the sample, based on the dust emission in a 120 arcsecond aperture, is 1.1 +- 0.9 solar masses. For the sources fitted by power-law intensity distributions (I_nu(b)/I_nu(0) = (b/b_0) ^ m), the mean value of m is 1.52 +- 0.45 for Class 0 and I sources at 850 microns and 1.44 +- 0.25 at 450 microns. Based on a simple analysis, assuming the emission is in the Rayleigh-Jeans limit and that T_d(r) ~ r ^ {-0.4}, these values of m translate into power-law density distributions (n ~ r ^ {-p}) with p ~ 2.1. However, we show that this result may be changed by more careful consideration of effects such as beam size and shape, finite outer radii, more realistic T_d(r), and failure of the Rayleigh-Jeans approximation.

preprint1996arXiv

Dense Gas and Star Formation: Characteristics of Cloud Cores Associated with Water Masers

We have observed 150 regions of massive star formation, selected originally by the presence of a water maser, in the J = 5-4, 3-2, and 2-1 transitions of CS, and 49 regions in the same transitions of C$^{34}$S. Over 90% of the 150 regions were detected in the J = 2-1 and 3-2 transitions of CS and 75% were detected in the J=5-4 transition. We have combined the data with the J = 7-6 data from our original survey (Plume et al. 1992) to determine the density by analyzing the excitation of the rotational levels. Using Large Velocity Gradient (LVG) models, we have determined densities and column densities for 71 of these regions. The gas densities are very high (the mean log of the density is 5.9), but much less than the critical density of the J=7-6 line. Small maps of 25 of the sources in the J = 5-4 line yield a mean diameter of 1.0 pc. The mean virial mass is 3800 solar masses. The mean ratio of bolometric luminosity to virial mass (L/M) is 190, about 50 times higher than estimates using CO emission, suggesting that star formation is much more efficient in the dense gas probed in this study. The gas depletion time for the dense gas is roughly 1.3 x 10^7 yr. We find no statistically significant linewidth--size or density--size relationships in our data. Instead, both linewidth and density are larger for a given size than would be predicted by the usual relationships. We find that the linewidth increases with density, the opposite of what would be predicted by the usual arguments. We estimate that the luminosity of our Galaxy (excluding the inner 400 pc) in the CS J = 5-4 transition is 15 to 23 L_sun, considerably less than the luminosity in this line within the central 100 pc of NGC 253 and M82. In addition, the ratio of far-infrared luminosity to CS luminosity is higher in M82 than in any cloud in our sample.

Neal J. Evans II

What is connected

Connect this record

See the researcher in context

Building this map preview

100 published item(s)

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

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- XI. From inflow to infall in hub-filament systems

Nobeyama Survey of Inward Motions toward Cores in Orion Identified by SCUBA-2

Slow Star Formation in the Milky Way: Theory Meets Observations

Tracers of Dense Gas in the Outer Galaxy

Atomic Shocks in the Outflow of L1551 IRS 5 Identified with SOFIA-upGREAT Observations of [OI]

Molecular Cloud Cores with High Deuterium Fractions: Nobeyama Mapping Survey

ALMA ACA and Nobeyama observations of two Orion cores in deuterated molecular lines

Large-scale molecular gas distribution in the M17 cloud complex: dense gas conditions of massive star formation?

Star Formation Occurs in Dense Gas, but What Does "Dense" Mean?

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

Probing Planet Forming Zones with Rare CO Isotopologues

The abundance and thermal history of water ice in the disk surrounding HD142527 from the DIGIT Herschel Key Program

The CDF Archive: Herschel PACS and SPIRE Spectroscopic Data Pipeline and Products for Protostars and Young Stellar Objects

The inflow signature toward different evolutionary phases of massive star formation

The Mid-Infrared Evolution of the FU Orionis Disk

"Dust, Ice, and Gas In Time" (DIGIT) Herschel Observations of GSS30-IRS1 in Ophiuchus

Detection of Infall in the Protostar B335 with ALMA

Evidence for Decay of Turbulence by MHD Shocks in Molecular Clouds via CO Emission

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

Infrared emission from tidal disruption events --- probing the pc-scale dust content around galactic nuclei

The Bolocam Galactic Plane Survey. XI. Temperatures and Substructure of Galactic Clumps Based on 350 micron Observations

The Bolocam Galactic Plane Survey. XIII. Physical Properties and Mass Functions of Dense Molecular Cloud Structures

The Gould Belt 'MISFITS' Survey: The Real Solar Neighborhood Protostars

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

Young Stellar Objects in the Gould Belt

Herschel Key Program, "Dust, Ice, and Gas In Time" (DIGIT): the origin of molecular and atomic emission in low-mass protostars in Taurus

Star Formation Relations in Nearby Molecular Clouds

The Bolocam Galactic Plane Survey. XII. Distance Catalog Expansion Using Kinematic Isolation of Dense Molecular Cloud Structures With 13CO(1-0)

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

The Gould Belt Very Large Array Survey IV: The Taurus-Auriga complex

The Star Formation Relation for Regions in the Galactic Plane: The Effect of Spatial Resolution

ALMA observations of the massive molecular outflow G331.512-0.103

An Analysis of the Environments of FU Orionis Objects with Herschel

DIGIT survey of far-infrared lines from protoplanetary discs: II. CO

DIGIT survey of far-infrared lines from protoplanetary disks I

Embedded Protostars in the Dust, Ice, and Gas In Time (DIGIT) Key Program: Continuum SEDs, and an Inventory of Characteristic Far-Infrared Lines from PACS Spectroscopy

L1448-MM observations by the Herschel Key program, "Dust, Ice, and Gas In Time" (DIGIT)

Physical characteristics of G331.5-0.1: The luminous central region of a Giant Molecular Cloud

Testing 24 micron and Infrared Luminosity as Star Formation Tracers for Galactic Star Forming Regions

The Bolocam Galactic Plane Survey IX: Data Release 2 and Outer Galaxy Extension

The Bolocam Galactic Plane Survey. VIII. A Mid-Infrared Kinematic Distance Discrimination Method

The Gould's Belt Very Large Array Survey I: The Ophiuchus complex

The Luminosities of Protostars in the Spitzer c2d and Gould Belt Legacy Clouds

The molecular ISM in luminous infrared galaxies: a 3mm line survey of Arp 157

The Star Formation Rate of Molecular Clouds

The VIRUS-P Exploration of Nearby Galaxies (VENGA): Survey Design, Data Processing, and Spectral Analysis Methods

Warm gas towards young stellar objects in Corona Australis - Herschel/PACS observations from the DIGIT key programme

A Herschel Survey of Cold Dust in Disks Around Brown Dwarfs and Low-Mass Stars

CO2 Ice toward Low-luminosity, Embedded Protostars: Evidence for Episodic Mass Accretion via Chemical History

Fragmentation and Evolution of Molecular Clouds. III: The Effect of Dust and Gas Energetics

Star Formation in the Milky Way and Nearby Galaxies

Submillimeter Follow-up of WISE-Selected Hyperluminous Galaxies

The First Hyper-Luminous Infrared Galaxy Discovered by WISE

The Herschel DIGIT Survey of Weak-line T Tauri Stars: implications for disk evolution and dissipation

The VIRUS-P Exploration of Nearby Galaxies (VENGA): The Xco Gradient in NGC 628

A Herschel Search For Cold Dust in Brown Dwarf Disks: First Results

A Mid-Infrared Census of Star Formation Activity in Bolocam Galactic Plane Survey Sources

Disentangling the Environment of the FU Orionis Candidate HBC 722 with Herschel

Herschel Observations of the T Cha Transition Disk: Constraining the Outer Disk Properties

Ices in starless and starforming cores

The Bolocam Galactic Plane Survey V: HCO+ and N2H+ Spectroscopy of 1.1 mm Dust Continuum Sources

The Bolocam Galactic Plane Survey. VII. Characterizing the Properties of Massive Star-Forming Regions

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

The Spitzer ice legacy: Ice evolution from cores to protostars

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

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

A Spitzer c2d Legacy Survey to Identify and Characterize Disks with Inner Dust Holes

c2d Spitzer IRS spectra of embedded low-mass young stars: gas-phase emission lines

Low-mass Star Formation: Observations

Observational Constraints on Submillimeter Dust Opacity

The Bolocam Galactic Plane Survey -- III. Characterizing Physical Properties of Massive Star-Forming Regions in the Gemini OB1 Molecular Cloud

The Bolocam Galactic Plane Survey IV: 1.1 and 0.35 mm Dust Continuum Emission in the Galactic Center Region

The Bolocam Galactic Plane Survey: Survey Description and Data Reduction