Arxiv:1304.7389V3 [Astro-Ph.GA] 9 May 2013 Rcdtruhawat Fmlclr Tmc N Ionic Millimeter-Wave and and Telescopes
Total Page:16
File Type:pdf, Size:1020Kb
Draft version June 9, 2018 Preprint typeset using LATEX style emulateapj v. 5/2/11 EMBEDDED PROTOSTARS IN THE DUST, ICE, AND GAS IN TIME (DIGIT) HERSCHEL KEY PROGRAM: CONTINUUM SEDS, AND AN INVENTORY OF CHARACTERISTIC FAR-INFRARED LINES FROM PACS SPECTROSCOPY Joel D. Green 1, Neal J. Evans II1, Jes K. Jørgensen2,3, Gregory J. Herczeg4,5, Lars E. Kristensen6,7, Jeong-Eun Lee8, Odysseas Dionatos3,2,9, Umut A. Yildiz6, Colette Salyk10, Gwendolyn Meeus11, Jeroen Bouwman12, Ruud Visser13, Edwin A. Bergin13, Ewine F. van Dishoeck6,5, Michelle R. Rascati1, Agata Karska5, Tim A. van Kempen6,14, Michael M. Dunham15, Johan E. Lindberg3,2, Davide Fedele5, & the DIGIT Team 1. The University of Texas at Austin, Department of Astronomy, 2515 Speedway, Stop C1400, Austin, TX 78712-1205, USA; [email protected] 2. Niels Bohr Institute, University of Copenhagen. Denmark 3. Centre for Star and Planet Formation, Natural History Museum of Denmark, University of Copenhagen, Denmark 4. Kavli Institute for Astronomy and Astrophysics, Peking University, Beijing, 100871, PR China 5. Max-Planck Institute for Extraterrestrial Physics, Postfach 1312, 85741, Garching, Germany 6. Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands 7. Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA, 02183, USA 8. Department of Astronomy & Space Science, Kyung Hee University, Gyeonggi, 446-701, Korea 9. University of Vienna, Department of Astronomy, T¨urkenschanzstrasse 17, 1180 Vienna, Austria 10. National Optical Astronomy Observatory, 950 N Cherry Ave Tucson, AZ 85719, USA 11. Universidad Autonoma de Madrid, Dpt. Fisica Teorica, Campus Cantoblanco, Spain 12. Max Planck Institute for Astronomy, Heidelberg, Germany 13. Department of Astronomy, University of Michigan, 500 Church Street, Ann Arbor, MI 48109-1042, USA 14. Joint ALMA offices, Av. Alonso de Cordova, Santiago, Chile 15. Dept. of Astronomy, Yale University, New Haven, CT, USA Draft version June 9, 2018 ABSTRACT We present 50-210 µm spectral scans of 30 Class 0/I protostellar sources, obtained with Herschel- PACS, and 0.5-1000 µm 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 µm [O I] flux ratio is 17.2 ± 9.2. The [O I] 63 µm 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 (hN i = (0.70 ± 1.12) × 1049 total particles). N CO correlates strongly with Lbol, but neither Trot nor N CO(warm)/N CO(hot) cor- relates 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 (hTroti = 194 ± 85 K) and OH (hTroti = 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. 1. INTRODUCTION Large samples of protostars in relatively nearby (d ≤ arXiv:1304.7389v3 [astro-ph.GA] 9 May 2013 Embedded protostars represent a transition from the 300 pc) clouds have been developed through recent sur- core collapse within molecular clouds to the eventual veys with the Spitzer Space Telescope (e.g, Evans et al. young star and protoplanetary disk system, revealed as 2009, Allen et al. in prep.), along with ground-based sur- the surrounding envelope is accreted and/or removed. veys (e.g., Jørgensen et al. 2009, Dunham et al. 2013). Processes occurring in the protostellar stage of evolu- We have selected a sample of well-studied protostars from tion determine the conditions in the disks, including the these studies, including both Class 0 and Class I objects. amount and nature of material available for planet forma- Class 0 and Class I sources are characterized observa- tion in later stages of evolution. In a picture of episodic tionally by rising spectral energy distributions (SEDs) accretion (e.g., Dunham et al. 2010), the disk mass at the between near-infrared and mid-infrared wavelengths. end of the protostellar phase may be determined by the In addition to the continuum emission, the far- phasing of the last burst of accretion onto the star and infrared/submillimeter bands contain numerous pure ro- the end of infall. Tracing the physical processes in these tational transitions of the CO ladder, as well as OH systems will lead to a better understanding of the con- and H2O, and several fine structure lines, all useful straints on planet formation. Physical processes can be tracers of gas content and properties. The transitions traced through a wealth of molecular, atomic, and ionic and collisional rates of these simple molecules are well- tracers available to optical, infrared, and millimeter-wave understood (see, e.g., Yang et al. 2010; Neufeld 2012, telescopes. for a recent update on CO). Thus these lines make 2 Green et al. excellent diagnostics of opacity, density, temperature, Herschel data calibration will follow. In §2 we describe and shock velocities (e.g., Kaufman & Neufeld 1996; the sample in the context of other surveys, and describe Flower & Pineau Des Forˆets 2010) of the gas surround- the observations. In §3 we present our detailed data re- ing these systems. duction method for acquiring accurate continuum SEDs ISO/LWS covered the 20-200 µm spectral region and linefluxes. In §4 we present results and statistics for and was well-suited to study the warm (T > 100 the sample. In §5 we analyze rotational diagrams for all K) region of protostellar envelopes, distinguishable detected molecular species. In §6 we present correlations from the ambient cloud typically probed in ground- amongst the detected and derived parameters. In §7 we based millimeter studies. ISO-LWS detected gas consider the origin of the line emission. We conclude in phase H2O, high-J CO rotational transitions, and §8. fine structure emission lines toward protostars and re- 2. lated sources (e.g., Lorenzetti et al. 1999; Giannini et al. OBSERVATIONS 1999; Ceccarelli et al. 1999; Lorenzetti et al. 2000; The DIGIT program uses Herschel to achieve far- Giannini et al. 2001; Nisini et al. 2002). These lines are infrared spectral coverage of protostars and more evolved indicative of the innermost regions of the protostellar en- young stellar objects in order to make a full inven- velope, exposed to heating by the central object, and the tory of their chemical tracers. The full DIGIT sam- outflow cavity region, where winds and jets may interact ple consists of 94 sources: 24 Herbig Ae/Be stars, 9 with the envelope and the surrounding cloud. T Tauri stars, 30 protostars observed with PACS and The Herschel Space Observatory (Pilbratt et al. 2010) HIFI spectroscopy, and an additional 31 weak-line T provides greater spatial and spectral resolution, as well as Tauri stars observed with PACS imaging. Cieza et al. sensitivity, allowing study of sources with weaker contin- (2013) presented the results for the weak-line T Tauri uum and line emission. Herschel is an ESA space-based stars; various papers are summarizing the results on submillimeter telescope with a 3.5-meter primary mir- disks in the Herbig and T Tauri star spectroscopic ror. The Herschel-PACS spectrograph (Photodetector survey (Sturm et al. 2013, Meeus et al., Fedele et al., Array Camera and Spectrometer, Poglitsch et al. 2010) subm.) In this paper we focus on the embedded pro- provides a generational improvement over the spectral tostar sample. We presented the first DIGIT obser- range covered by ISO-LWS, and is complementary with vations of a protostar, DK Cha, in van Kempen et al. the Infrared Spectrograph (Houck et al. 2004) on the (2010a). The 30 object sample is drawn from previ- Spitzer Space Telescope (Werner et al. 2004). Spitzer- ous studies (summarized by Lahuis et al. 2006), focus- IRS revealed a wealth of diagnostics in the mid-infrared, ing on protostars with high-quality Spitzer-IRS 5-40 µm at relatively low spectral (R ∼ 60-600) and spatial (2– and UV, optical, infrared, and submillimeter complemen- ′′ 5 ) resolution. For example, H2O and OH were detected tary data. Thus all of the sources in the DIGIT sam- in some protostars (Watson et al. 2007) but in transi- ple are well-studied at other wavelengths, and many are tions with typically greater critical densities than those part of ongoing studies with Herschel in other programs, detectable by Herschel. Combined with enhanced spatial especially in the Water in Star-Forming Regions with resolution, this difference in excitation requirements can Herschel (WISH) program (van Dishoeck et al. 2011; lead to significantly different conclusions on the origin of see also Nisini et al. (2010); Kristensen et al. (2012); the emission without the context of the Herschel lines Karska et al. (2013); Wampfler et al. (2013)). Specific (Herczeg et al. 2012). Similarly while pure rotational outflow positions that are well separated from the cen- lines of H2 were detected with Spitzer in many regions tral source positions are observed by, e.g., Lefloch et al.