We have mapped the G287.84-0.82 cometary globule (with the Treasure Chest cluster embedded in it) in the South Pillars region of Carina (i) in [CII], 63micron [OI], and CO(11-10) using upGREAT on SOFIA and (ii) in J=2-1 transitions of CO, 13CO, C18O and J=3-2 transitions of H2CO using the APEX telescope in Chile. We probe the morphology, kinematics, and physical conditions of the molecular gas and the photon dominated regions (PDRs) in G287.84-0.82. The [CII] and [OI] emission suggest that the overall structure of the pillar (with red-shifted photo evaporating tails) is consistent with the effect of FUV radiation and winds from eta-Car and O stars in Trumpler 16. The gas in the head of the pillar is strongly influenced by the embedded cluster, whose brightest member is an O9.5V star, CPD-59 2661. The emission of the [CII] and [OI] lines peak at a position close to the embedded star, while all other tracers peak at another position lying to the north-east consistent with gas being compressed by the expanding PDR created by the embedded cluster. The molecular gas inside the globule is probed with the J=2-1 transitions of CO and isotopologues as well as H2CO, and analyzed using a non-LTE model (escape-probability approach), while we use PDR models to derive the physical conditions of the PDR. We identify at least two PDR gas components; the diffuse part (~10^4 cm^-3) is traced by [CII], while the dense (n~ 2-8x10^5 cm^-3) part is traced by [CII], [OI], CO(11-10). Using the F=2-1 transition of [13CII] detected at 50 positions in the region, we derive optical depths (0.9-5), excitation temperatures of [CII] (80-255 K), and N(C+) of 0.3-1x10^19 cm^-2. The total mass of the globule is ~1000 Msun, about half of which is traced by [CII]. The dense PDR gas has a thermal pressure of 10^7-10^8 K cm^-3, which is similar to the values observed in other regions.

HH 50138 is one of the brightest B[e] stars at a distance of $\sim$ 380 pc with strong infrared excess. The star was observed in [O I] 63 $\mu$m and [C II] 158 $\mu$m with high velocity resolution with upGREAT on SOFIA. The velocity resolved [O I] emission provides evidence for a large gas-disk, $\sim$ 760 au in size, around HD 50138. Whereas previous interferometric observations give strong evidence for a hot gas and dust disk in Keplerian rotation, our bservations are the first to provide unambiguous evidence for a large warm disk around the star. Herschel/PACS observations showed that the [C II] emission is extended, therefore the [C II] emission most likely originates in an ionized gas shell created by a past outflow event. We confirm the isolated nature of HD 50138. It is far from any star forming region and has low proper motion. Neither is there any sign of a remnant cloud from which it could have formed. The extended disk around the star appears carbon poor. It shows OH and [O I] emission, but no CO. The CO abundance appears to be at least an order of magnitude lower than that of OH. Furthermore $^{13}$CO is enriched by more than a factor of five, confirming that the star is not a Herbig Be star. Finally we note that our high spectral resolution [O I] and [C II] observations provide a very accurate heliocentric velocity of the star, 40.8 $\pm$ 0.2 km~s$^{-1}$.

We analyze a [C II] 158 micron map obtained with the L2 GREAT receiver on SOFIA of the emission/reflection nebula illuminated by the early B star S1 in the rho-OphA cloud core. This data set has been complemented with maps of CO(3-2), 13CO(3-2) and C18O(3-2), observed as a part of the JCMT Gould Belt Survey, with archival HCO^+(4-3) JCMT data, as well as with [O I] 63 and 145 micron imaging with Herschel/PACS. The [C II] emission is completely dominated by the strong PDR emission from the nebula surrounding S1 expanding into the dense Oph A molecular cloud west and south of S1. The [C II] emission is significantly blue shifted relative to the CO spectra and also relative to the systemic velocity, particularly in the northwestern part of the nebula. The [C II] lines are broader towards the center of the S1 nebula and narrower towards the PDR shell. The [C II] lines are strongly self-absorbed over an extended region in the S1 PDR. Based on the strength of the [13C II] F = 2-1 hyperfine component, [C II] is significantly optically thick over most of the nebula. CO and 13CO(3-2) spectra are strongly self-absorbed, while C18O(3-2) is single peaked and centered in the middle of the self-absorption. We have used a simple two-layer LTE model to characterize the background and foreground cloud contributing to the [C II] emission. From this analysis we estimate the extinction due to the foreground cloud to be ~9.9 mag, which is slightly less than the reddening estimated towards S1. Since some of the hot gas in the PDR is not traced by low J CO emission, this result appears quite plausible. Using a plane parallel PDR model with the observed [OI(145)]/[C II] brightness ratio and an estimated FUV intensity of 3100-5000 G0 suggests that the density of the [C II] emitting gas is ~3-4x10^3 cm^-3.

Determining the properties of the warm ionized medium (WIM) at the leading edge of spiral arms is important for understanding its dynamics and cloud formation. The inner edge of the Scutum arm tangency is a unique location in which to disentangle the WIM from other components. We use high spectral resolution [C II] 158 micron and [N II] 205 micron fine structure line observations taken with the upGREAT and GREAT instruments on SOFIA, along with auxiliary HI and 13CO observations. The observations were in and out of the Galactic plane along 18 lines of sight between longitude 30deg and 32deg. We detect strong [N II] emission throughout the Scutum tangency. At VLSR = 110 to 125 km/s, where there is little, if any, 13CO, we could disentangle the [N II] and [C II] emission that arises from the WIM at the inner edge. We find an average electron density, <n(e)> about 0.9 cm{-3} in the plane, and about 0.4 cm{-3} above the plane. For VLSR < 110 km/s there is [N II] emission tracing highly ionized gas throughout the arm's molecular layer. This ionized gas has a high density, n(e) ~ 30 cm{-3}, and a few percent filling factor. Thus, [N II] and [C II] at the Scutum arm tangency reveal a highly ionized gas with <n(e)> about 10 to 20 times those of the interarm WIM, which is best explained by a model in which the interarm WIM is compressed as it falls into the potential well of the arm. The widespread distribution of [N II] in the molecular layers shows that high density ionized gas is distributed throughout the Scutum arm. The n(e) derived from [N II] for these molecular cloud regions are about 30 cm{-3}, and probably arise in the ionized boundary layers of clouds. This [N II] emission from the ionized boundary layers is probably the result of the shock compression of the WIM as it impacts the arm's neutral gas, but some could arise from extended HII regions.

Herbig Ae/Be objects are pre-main sequence stars surrounded by gas- and dust-rich circumstellar discs. These objects are in the throes of star and planet formation, and their characterisation informs us of the processes and outcomes of planet formation processes around intermediate mass stars. Here we analyse the spectral energy distributions of disc host stars observed by the Herschel Open Time Key Programme `Gas in Protoplanetary Systems'. We present Herschel/PACS far-infrared imaging observations of 22 Herbig Ae/Bes and 5 debris discs, combined with ancillary photometry spanning ultraviolet to sub-millimetre wavelengths. From these measurements we determine the diagnostics of disc evolution, along with the total excess, in three regimes spanning near-, mid-, and far-infrared wavelengths. Using appropriate statistical tests, these diagnostics are examined for correlations. We find that the far-infrared flux, where the disc becomes optically thin, is correlated with the millimetre flux, which provides a measure of the total dust mass. The ratio of far-infrared to sub-millimetre flux is found to be greater for targets with discs that are brighter at millimetre wavelengths and that have steeper sub-millimetre slopes. Furthermore, discs with flared geometry have, on average, larger excesses than flat geometry discs. Finally, we estimate the extents of these discs (or provide upper limits) from the observations.

We have mapped the NGC 2023 reflection nebula in [CII] and CO(11--10) with the heterodyne receiver GREAT on SOFIA and obtained slightly smaller maps in 13CO(3--2), CO(3--2), CO(4--3), CO(6--5), and CO(7--6) with APEX in Chile. We use these data to probe the morphology, kinematics, and physical conditions of the C II region, which is ionized by FUV radiation from the B2 star HD37903. The [CII] emission traces an ellipsoidal shell-like region at a position angle of ~ -50 deg, and is surrounded by a hot molecular shell. In the southeast, where the C II region expands into a dense, clumpy molecular cloud ridge, we see narrow and strong line emission from high-J CO lines, which comes from a thin, hot molecular shell surrounding the [CII] emission. The [CII] lines are broader and show photo evaporating gas flowing into the C II region. Based on the strength of the [13CII] F=2--1 line, the [CII] line appears to be somewhat optically thick over most of the nebula with an optical depth of a few. We model the physical conditions of the surrounding molecular cloud and the PDR emission using both RADEX and simple PDR models. The temperature of the CO emitting PDR shell is ~ 90 -- 120 K, with densities of 10^5 -- 10^6 cm^-3, as deduced from RADEX modeling. Our PDR modeling indicates that the PDR layer where [CII] emission dominates has somewhat lower densities, 10^4 to a few times 10^5 cm^-3

We present dust polarization and CO molecular line images of NGC 7538 IRS1. We combined data from the SMA, CARMA and JCMT telescopes to make images with 2.5 arcsec resolution at 230 and 345 GHz. The images show a remarkable spiral pattern in both the dust polarization and molecular outflow. These data dramatically illustrate the interplay between a high infall rate onto IRS1 and a powerful outflow disrupting the dense, clumpy medium surrounding the star. The images of the dust polarization and the CO outflow presented here provide observational evidence for the exchange of energy and angular momentum between the infall and the outflow. The spiral dust pattern, which rotates through over 180 degrees from IRS1, may be a clumpy filament wound up by conservation of angular momentum in the infalling material. The redshifted CO emission ridge traces the dust spiral closely through the MM dust cores, several of which may contain protostars. We propose that the CO maps the boundary layer where the outflow is ablating gas from the dense gas in the spiral.

HD 100546 is a well-studied Herbig Be star-disk system that likely hosts a close-in companion with compelling observational evidence for an embedded protoplanet at 68 AU. We present ALMA observations of the HD 100546 disk which resolve the gas and dust structure at (sub)mm wavelengths. The CO emission (at 345.795 GHz) originates from an extensive molecular disk (390+/-20 AU in radius) whereas the continuum emission is more compact (230+/-20 AU in radius) suggesting radial drift of the mm-sized grains. The CO emission is similar in extent to scattered light images indicating well-mixed gas and um-sized grains in the disk atmosphere. Assuming azimuthal symmetry, a single-component power-law model cannot reproduce the continuum visibilities. The visibilities and images are better reproduced by a double-component model: a compact ring with a width of 21 AU centered at 26 AU and an outer ring with a width of 75+/-3 AU centered at 190+/-3 AU. The influence of a companion and protoplanet on the dust evolution is investigated. The companion at 10 AU facilitates the accumulation of mm-sized grains within a compact ring, ~ 20 - 30 AU, by ~ 10 Myr. The injection of a protoplanet at 1 Myr hastens the ring formation (~ 1.2 Myr) and also triggers the development of an outer ring (~ 100 - 200 AU). These observations provide additional evidence for the presence of a close-in companion and hint at dynamical clearing by a protoplanet in the outer disk.

The Stratospheric Observatory for Infrared Astronomy (SOFIA) has recently concluded a set of engineering flights for Observatory performance evaluation. These in-flight opportunities are viewed as a first comprehensive assessment of the Observatory's performance and are used to guide future development activities, as well as to identify additional Observatory upgrades. Pointing stability was evaluated, including the image motion due to rigid-body and flexible-body telescope modes as well as possible aero-optical image motion. We report on recent improvements in pointing stability by using an active mass damper system installed on the telescope. Measurements and characterization of the shear layer and cavity seeing, as well as image quality evaluation as a function of wavelength have also been performed. Additional tests targeted basic Observatory capabilities and requirements, including pointing accuracy, chopper evaluation and imager sensitivity. This paper reports on the data collected during these flights and presents current SOFIA Observatory performance and characterization.

We present {\lambda}1.3 mm CARMA observations of dust polarization toward 30 star-forming cores and 8 star-forming regions from the TADPOL survey. We show maps of all sources, and compare the ~2.5" resolution TADPOL maps with ~20" resolution polarization maps from single-dish submillimeter telescopes. Here we do not attempt to interpret the detailed B-field morphology of each object. Rather, we use average B-field orientations to derive conclusions in a statistical sense from the ensemble of sources, bearing in mind that these average orientations can be quite uncertain. We discuss three main findings: (1) A subset of the sources have consistent magnetic field (B-field) orientations between large (~20") and small (~2.5") scales. Those same sources also tend to have higher fractional polarizations than the sources with inconsistent large-to-small-scale fields. We interpret this to mean that in at least some cases B-fields play a role in regulating the infall of material all the way down to the ~1000 AU scales of protostellar envelopes. (2) Outflows appear to be randomly aligned with B-fields; although, in sources with low polarization fractions there is a hint that outflows are preferentially perpendicular to small-scale B-fields, which suggests that in these sources the fields have been wrapped up by envelope rotation. (3) Finally, even at ~2.5" resolution we see the so-called "polarization hole" effect, where the fractional polarization drops significantly near the total intensity peak. All data are publicly available in the electronic edition of this article.

We present spectrally resolved Herschel/HIFI observations of the young multiple system T Tau in atomic and molecular lines. While CO, H2O, [C II], and SO lines trace the envelope and the outflowing gas up to velocities of 33 km/s with respect to systemic, the CN 5-4 hyperfine structure lines at 566.7, 566.9 GHz show a narrow double-peaked profile centered at systemic velocity, consistent with an origin in the outer region of the compact disk of T Tau N. Disk modeling of the T Tau N disk with the thermo-chemical code ProDiMo produces CN line fluxes and profiles consistent with the observed ones and constrain the size of the gaseous disk (R_out = 110 (+10, -20) AU) and its inclination (i = 25 \pm 5 degree). The model indicates that the CN lines originate in a disk upper layer at 40-110 AU from the star, which is irradiated by the stellar UV field and heated up to temperatures of 50-700 K. With respect to previously observed CN 2-1 millimeter lines, the CN 5-4 lines appear to be less affected by envelope emission, due to their larger critical density and excitation temperature. Hence, high-J CN lines are a unique confusion-free tracer of embedded disks, such as the disk of T Tau N.

We report new results based on the analysis of the SMA and CARMA observations of NGC 7538\,IRS\,1 at 1.3 and 3.4 mm with sub-arcsec resolutions. With angular resolutions $\sim$ 0\farcs7, the SMA and CARMA observations show that the continuum emission at 1.3 and 3.4 mm from the hyper-compact \ion{H}{2} region IRS\,1 is dominated by a compact source with a tail-like extended structure to the southwest of IRS\,1. With a CARMA B-array image at 1.3 mm convolved to 0\farcs1, we resolve the hyper-compact \ion{H}{2} region into two components: an unresolved hyper-compact core, and a north-south extension with linear sizes of $<270$ AU and $\sim$2000 AU, respectively. The fine structure observed with CARMA is in good agreement with the previous VLA results at centimeter wavelengths, suggesting that the hyper-compact \ion{H}{2} region at the center of IRS\,1 is associated with an ionized bipolar outflow. We image the molecular lines OCS(19-18) and CH$_3$CN(12-11) as well as $^{13}$CO(2-1) surrounding IRS\,1, showing a velocity gradient along the southwest-northeast direction. The spectral line profiles in $^{13}$CO(2-1), CO(2-1), and HCN(1-0) observed toward IRS\,1 show broad redshifted absorption, providing evidence for gas infall with rates in the range of $3-10\times10^{-3}$ M$_\odot$ yr$^{-1}$ inferred from our observations.

At the distance of 99-116 pc, HD141569A is one of the nearest HerbigAe stars that is surrounded by a tenuous disk, probably in transition between a massive primordial disk and a debris disk. We observed the fine-structure lines of OI at 63 and 145 micron and the CII line at 157 micron with the PACS instrument onboard the Herschel Space Telescope as part of the open-time large programme GASPS. We complemented the atomic line observations with archival Spitzer spectroscopic and photometric continuum data, a ground-based VLT-VISIR image at 8.6 micron, and 12CO fundamental ro-vibrational and pure rotational J=3-2 observations. We simultaneously modeled the continuum emission and the line fluxes with the Monte Carlo radiative transfer code MCFOST and the thermo-chemical code ProDiMo to derive the disk gas- and dust properties assuming no dust settling. The models suggest that the oxygen lines are emitted from the inner disk around HD141569A, whereas the [CII] line emission is more extended. The CO submillimeter flux is emitted mostly by the outer disk. Simultaneous modeling of the photometric and line data using a realistic disk structure suggests a dust mass derived from grains with a radius smaller than 1 mm of 2.1E-7 MSun and from grains with a radius of up to 1 cm of 4.9E-6 MSun. We constrained the polycyclic aromatic hydrocarbons (PAH) mass to be between 2E-11 and 1..4E-10 MSun assuming circumcircumcoronene (C150H30) as the representative PAH. The associated PAH abundance relative to hydrogen is lower than those found in the interstellar medium (3E-7) by two to three orders of magnitude. The disk around HD141569A is less massive in gas (2.5 to 4.9E-4 MSun or 67 to 164 MEarth) and has a flat opening angle (<10%). [abridged]

The Herschel Space Observatory was used to observe ~ 120 pre-main-sequence stars in Taurus as part of the GASPS Open Time Key project. PACS was used to measure the continuum as well as several gas tracers such as [OI] 63 \mu m, [OI] 145 \mu m, [CII] 158 \mu m, OH, H2O and CO. The strongest line seen is [OI] at 63 \mu m. We find a clear correlation between the strength of the [OI] 63 \mu m line and the 63 \mu m continuum for disk sources. In outflow sources, the line emission can be up to 20 times stronger than in disk sources, suggesting that the line emission is dominated by the outflow. The tight correlation seen for disk sources suggests that the emission arises from the inner disk ($<$ 50 AU) and lower surface layers of the disk where the gas and dust are coupled. The [OI] 63 \mu m is fainter in transitional stars than in normal Class II disks. Simple SED models indicate that the dust responsible for the continuum emission is colder in these disks, leading to weaker line emission. [CII] 158 \mu m emission is only detected in strong outflow sources. The observed line ratios of [OI] 63 \mu m to [OI] 145 \mu m are in the regime where we are insensitive to the gas-to-dust ratio, neither can we discriminate between shock or PDR emission. We detect no Class III object in [OI] 63 \mu m and only three in continuum, at least one of which is a candidate debris disk.

Context. Circumstellar discs are the places where planets form, therefore knowledge of their evolution is crucial for our understanding of planet formation. The Herschel Space Observatory is providing valuable data for studying disc systems, thanks to its sensitivity and wavelength coverage. This paper is one of several devoted to analysing and modelling Herschel-PACS observations of various young stellar associations from the GASPS Open Time Key Programme. Aims. The aim of this paper is to elucidate the gas and dust properties of circumstellar discs in the 10 Myr TW Hya Association (TWA) using new far-infrared (IR) imaging and spectroscopy from Herschel-PACS. Methods. We obtained far-IR photometric data at 70, 100, and 160 microns of 14 TWA members; spectroscopic observations centred on the [OI] line at 63.18 microns were also obtained for 9 of the 14. The new photometry for each star was incorporated into its full spectral energy distribution (SED). Results. We detected excess IR emission that is characteristic of circumstellar discs from five TWA members, and computed upper limits for another nine. Two TWA members (TWA 01 and TWA 04B) also show [OI] emission at 63.18 microns. Discs in the TWA association display a variety of properties, with a wide range of dust masses and inner radii, based on modified blackbody modelling. Both transitional and debris discs are found in the sample. Models for sources with a detected IR excess give dust masses in the range from 0.15 Msun to 63 Msun.

GASPS is a far-infrared line and continuum survey of protoplanetary and young debris disks using PACS on the Herschel Space Observatory. The survey includes [OI] at 63 microns, as well as 70, 100 and 160um continuum, with the brightest objects also studied in [OI]145um, [CII]157um, H2O and CO. Targets included T Tauri stars and debris disks in 7 nearby young associations, and a sample of isolated Herbig AeBe stars. The aim was to study the global gas and dust content in a wide disk sample, systemically comparing the results with models. In this paper we review the main aims, target selection and observing strategy. We show initial results, including line identifications, sources detected, and a first statistical study. [OI]63um was the brightest line in most objects, by a factor of ~10. Detection rates were 49%, including 100% of HAeBe stars and 43% of T Tauri stars. Comparison with published dust masses show a dust threshold for [OI]63um detection of ~1e-5 M_solar. Normalising to 140pc distance, 32% with mass 1e-6 - 1e-5 M_solar, and a small number with lower mass were also detected. This is consistent with moderate UV excess and disk flaring. In most cases, continuum and line emission is spatially and spectrally unresolved, suggesting disk emission. ~10 objects were resolved, likely from outflows. Detection rates in [OI]145um, [CII]157um and CO J=18-17 were 20-40%, but [CII] was not correlated with disk mass, suggesting it arises instead from a compact envelope. [OI] detection rates in T Tauri associations of ages 0.3-4Myr were ~50%. ~2 stars were detectable in associations of 5-20Myr, with no detections in associations of age >20Myr. Comparing with the total number of young stars, and assuming a ISM-like gas/dust ratio, this indicates that ~18% of stars retain a gas-rich disk of total mass >1M_Jupiter for 1-4Myr, 1-7% keep such disks for 5-10Myr, and none remain beyond 10-20Myr.

We present far-IR/sub-mm imaging and spectroscopy of 49 Ceti, an unusual circumstellar disk around a nearby young A1V star. The system is famous for showing the dust properties of a debris disk, but the gas properties of a low-mass protoplanetary disk. The data were acquired with the Herschel Space Observatory PACS and SPIRE instruments, largely as part of the "Gas in Protoplanetary Systems" (GASPS) Open Time Key Programme. Disk dust emission is detected in images at 70, 160, 250, 350, and 500 \mu m; 49 Cet is significantly extended in the 70 \mu m image, spatially resolving the outer dust disk for the first time. Spectra covering small wavelength ranges centered on eight atomic and molecular emission lines were obtained, including [OI] 63 \mu m and [CII] 158 \mu m. The CII line was detected at the 5\sigma\ level - the first detection of atomic emission from the disk. No other emission lines were seen, despite the fact that the OI line is the brightest one observed in Herschel protoplanetary disk spectra (Meeus et al. 2012; Dent et al. 2013). We present an estimate of the amount of circumstellar atomic gas implied by the CII emission. The new far-IR/sub-mm data fills in a large gap in the previous spectral energy distribution (SED) of 49 Cet. A simple model of the new SED confirms the two-component structure of the disk: warm inner dust and cold outer dust that produces most of the observed excess. Finally, we discuss preliminary thermochemical modeling of the 49 Cet gas/dust disk and our attempts to match several observational results simultaneously. Although we are not yet successful in doing so, our investigations shed light on the evolutionary status of the 49 Cet gas, which might not be primordial gas but rather secondary gas coming from comets.

We present 30 and 40 micron imaging of the massive protostar G35.20-0.74 with SOFIA-FORCAST. The high surface density of the natal core around the protostar leads to high extinction, even at these relatively long wavelengths, causing the observed flux to be dominated by that emerging from the near-facing outflow cavity. However, emission from the far-facing cavity is still clearly detected. We combine these results with fluxes from the near-infrared to mm to construct a spectral energy distribution (SED). For isotropic emission the bolometric luminosity would be 3.3x10^4 Lsun. We perform radiative transfer modeling of a protostar forming by ordered, symmetric collapse from a massive core bounded by a clump with high mass surface density, Sigma_cl. To fit the SED requires protostellar masses ~20-34 Msun depending on the outflow cavity opening angle (35 - 50 degrees), and Sigma_cl ~ 0.4-1 g cm-2. After accounting for the foreground extinction and the flashlight effect, the true bolometric luminosity is ~ (0.7-2.2)x10^5 Lsun. One of these models also has excellent agreement with the observed intensity profiles along the outflow axis at 10, 18, 31 and 37 microns. Overall our results support a model of massive star formation involving the relatively ordered, symmetric collapse of a massive, dense core and the launching bipolar outflows that clear low density cavities. Thus a unified model may apply for the formation of both low and high mass stars.

Water is key in the evolution of protoplanetary disks and the formation of comets and icy/water planets. While high excitation water lines originating in the hot inner disk have been detected in several T Tauri stars (TTSs), water vapor from the outer disk, where most of water ice reservoir is stored, was only reported in the closeby TTS TW Hya. We present spectrally resolved Herschel/HIFI observations of the young TTS DG Tau in the ortho- and para- water ground-state transitions at 557, 1113 GHz. The lines show a narrow double-peaked profile, consistent with an origin in the outer disk, and are ~19-26 times brighter than in TW Hya. In contrast, CO and [C II] lines are dominated by emission from the envelope/outflow, which makes H2O lines a unique tracer of the disk of DG Tau. Disk modeling with the thermo-chemical code ProDiMo indicates that the strong UV field, due to the young age and strong accretion of DG Tau, irradiates a disk upper layer at 10-90 AU from the star, heating it up to temperatures of 600 K and producing the observed bright water lines. The models suggest a disk mass of 0.015-0.1 Msun, consistent with the estimated minimum mass of the solar nebula before planet formation, and a water reservoir of ~1e2-1e3 Earth oceans in vapour, and ~100 times larger in the form of ice. Hence, this detection supports the scenario of ocean delivery on terrestrial planets by impact of icy bodies forming in the outer disk.

Context. HD 172555 is a young A7 star belonging to the Beta Pictoris Moving Group that harbours a debris disc. The Spitzer IRS spectrum of the source showed mid-IR features such as silicates and glassy silica species, indicating the presence of a warm dust component with small grains, which places HD 172555 among the small group of debris discs with such properties. The IRS spectrum also shows a possible emission of SiO gas. Aims. We aim to study the dust distribution in the circumstellar disc of HD 172555 and to asses the presence of gas in the debris disc. Methods. As part of the GASPS Open Time Key Programme, we obtained Herschel-PACS photometric and spectroscopic observations of the source. We analysed PACS observations of HD 172555 and modelled the Spectral Energy Distribution (SED) with a modified blackbody and the gas emission with a two-level population model with no collisional de-excitation. Results. We report for the first time the detection of OI atomic gas emission at 63.18 microns in the HD 172555 circumstellar disc. We detect excesses due to circumstellar dust toward HD 172555 in the three photometric bands of PACS (70, 100, and 160 microns). We derive a large dust particle mass of 4.8e-4 Earth masses and an atomic oxygen mass of 2.5e-2*R^2 Earth masses, where R in AU is the separation between the star and the inner disc. Thus, most of the detected mass of the disc is in the gaseous phase.

Observations of the atomic and molecular line emission associated with jets and outflows emitted by young stellar objects can be used to trace the various evolutionary stages they pass through as they evolve to become main sequence stars. To understand the relevance of atomic and molecular cooling in shocks, and how accretion and ejection efficiency evolves with the source evolutionary state, we will study the far-infrared counterparts of bright optical jets associated with Class I and II sources in Taurus (T Tau, DG Tau A, DG Tau B, FS Tau A+B, and RW Aur). We have analysed Herschel/PACS observations of a number of atomic ([OI]63um, 145um, [CII]158um) and molecular (high-J CO, H2O, OH) lines, collected within the OTKP GASPS. To constrain the origin of the detected lines we have compared the FIR emission maps with the emission from optical-jets and millimetre-outflows, and the line fluxes and ratios with predictions from shock and disk models. All of the targets are associated with extended emission in the atomic lines correlated with the direction of the optical jet/mm-outflow. The atomic lines can be excited in fast dissociative J-shocks. The molecular emission, on the contrary, originates from a compact region, that is spatially and spectrally unresolved. Slow C- or J- shocks with high pre-shock densities reproduce the observed H2O and high-J CO lines; however, the disk and/or UV-heated outflow cavities may contribute to the emission. While the cooling is dominated by CO and H2O lines in Class 0 sources, [OI] becomes an important coolant as the source evolves and the environment is cleared. The cooling and mass loss rates estimated for Class II and I sources are one to four orders of magnitude lower than for Class 0 sources. This provides strong evidence to indicate that the outflow activity decreases as the source evolves.

We present Herschel PACS photometry of seventeen B- to M-type stars in the 30 Myr-old Tucana-Horologium Association. This work is part of the Herschel Open Time Key Programme "Gas in Protoplanetary Systems" (GASPS). Six of the seventeen targets were found to have infrared excesses significantly greater than the expected stellar IR fluxes, including a previously unknown disk around HD30051. These six debris disks were fitted with single-temperature blackbody models to estimate the temperatures and abundances of the dust in the systems. For the five stars that show excess emission in the Herschel PACS photometry and also have Spitzer IRS spectra, we fit the data with models of optically thin debris disks with realistic grain properties in order to better estimate the disk parameters. The model is determined by a set of six parameters: surface density index, grain size distribution index, minimum and maximum grain sizes, and the inner and outer radii of the disk. The best fitting parameters give us constraints on the geometry of the dust in these systems, as well as lower limits to the total dust masses. The HD105 disk was further constrained by fitting marginally resolved PACS 70 micron imaging.

The Stratospheric Observatory for Infrared Astronomy (SOFIA) is an airborne observatory consisting of a specially modified Boeing 747SP with a 2.7-m telescope, flying at altitudes as high as 13.7 km (45,000 ft). Designed to observe at wavelengths from 0.3 micron to 1.6 mm, SOFIA operates above 99.8 % of the water vapor that obscures much of the infrared and submillimeter. SOFIA has seven science instruments under development, including an occultation photometer, near-, mid-, and far-infrared cameras, infrared spectrometers, and heterodyne receivers. SOFIA, a joint project between NASA and the German Aerospace Center DLR, began initial science flights in 2010 December, and has conducted 30 science flights in the subsequent year. During this early science period three instruments have flown: the mid-infrared camera FORCAST, the heterodyne spectrometer GREAT, and the occultation photometer HIPO. This article provides an overview of the observatory and its early performance.

Molecular globules and pillars are spectacular features, found only in the interface region between a molecular cloud and an HII-region. Impacting Far-ultraviolet (FUV) radiation creates photon dominated regions (PDRs) on their surfaces that can be traced by typical cooling lines. With the GREAT receiver onboard SOFIA we mapped and spectrally resolved the [CII] 158 micron atomic fine-structure line and the highly excited 12CO J=11-10 molecular line from three objects in Cygnus X (a pillar, a globule, and a strong IRAS source). We focus here on the globule and compare our data with existing Spitzer data and recent Herschel Open-Time PACS data. Extended [CII] emission and more compact CO-emission was found in the globule. We ascribe this emission mainly to an internal PDR, created by a possibly embedded star-cluster with at least one early B-star. However, external PDR emission caused by the excitation by the Cyg OB2 association cannot be fully excluded. The velocity-resolved [CII] emission traces the emission of PDR surfaces, possible rotation of the globule, and high-velocity outflowing gas. The globule shows a velocity shift of ~2 km/s with respect to the expanding HII-region, which can be understood as the residual turbulence of the molecular cloud from which the globule arose. This scenario is compatible with recent numerical simulations that emphazise the effect of turbulence. It is remarkable that an isolated globule shows these strong dynamical features traced by the [CII]-line, but it demands more observational studies to verify if there is indeed an embedded cluster of B-stars.

The BD+40\degree4124 region was observed with high angular and spectral resolution with the German heterodyne instrument GREAT in CO J = 13 \rightarrow 12 and [CII] on SOFIA. These observations show that the [CII] emission is very strong in the reflection nebula surrounding the young Herbig Ae/Be star BD+40\degree4124. A strip map over the nebula shows that the [CII] emission approximately coincides with the optical nebulosity. The strongest [CII] emission is centered on the B2 star and a deep spectrum shows that it has faint wings, which suggests that the ionized gas is expanding. We also see faint CO J = 13 \rightarrow 12 at the position of BD+40\degree4124, which suggests that the star may still be surrounded by an accretion disk.We also detected [CII] emission and strong CO J = 13 \rightarrow 12 toward V1318 Cyg. Here the [CII] emission is fainter than in BD+40\degree4124 and appears to come from the outflow, since it shows red and blue wings with very little emission at the systemic velocity, where the CO emission is quite strong. It therefore appears that in the broad ISO beam the [CII] emission was dominated by the reflection nebula surrounding BD+40\degree4124, while the high J CO lines originated from the adjacent younger and more deeply embedded binary system V1318 Cyg.

We have obtained new velocity-resolved spectra of the [CII] 158 micron line towards the Ring Nebula in Lyra (NGC 6720), one of the best-studied planetary nebulae, in order to probe its controversial 3-dimensional structure and to estimate the mass of circumstellar material in this object. We used the Terahertz receiver GREAT aboard the SOFIA airborne telescope to obtain the [CII] spectra at eight locations within and outside the bright optical ring of NGC 6720. Emission was detected at all positions except for the most distant position along the nebula's minor axis, and generally covers a broad velocity range, ~50 km/s (FWZI), except at a position along the major axis located just outside the optical ring, where it is significantly narrower (~25 km/s). The one narrow spectrum appears to be probing circumstellar material lying outside the main nebular shell that has not been accelerated by past fast wind episodes from the central star, and therefore most likely comes from equatorial and/or low-latitude regions of this bipolar nebula. Along lines-of-sight passing within about 10 arcsec of the nebular center, the CII column density is a factor 46 higher than the CO column density. The total mass of gas associated with the [CII] emission inside a circular region of diameter 87.5 arcsec is at least 0.11 solar masses. A significant amount of [CII] flux arises from a photodissociation region immediately outside the bright optical ring, where we find a CII to CO ratio of 6.5, lower than that seen towards the central region. Comparing our data with lower-quality CI spectra, which indicate similarly large CI/CO ratios in NGC 6720, we conclude that the bulk of elemental carbon in NGC 6720 is divided roughly equally between CII and CI, and that the emissions from these species are far more robust tracers of circumstellar material than CO in this object and other evolved planetary nebulae.

We present high angular resolution continuum observations of the high-mass protostar NGC 7538S with BIMA and CARMA at 3 and 1.4 mm, VLA observations at 1.3, 2, 3.5 and 6 cm, and archive IRAC observations from the Spitzer Space Observatory, which detect the star at 4.5, 5.8, and 8 $\mu$m. The star looks rather unremarkable in the mid-IR. The excellent positional agreement of the IRAC source with the VLA free-free emission, the OH, CH$_3$OH, H$_2$O masers, and the dust continuum confirms that this is the most luminous object in the NGC 7538S core. The continuum emission at millimeter wavelengths is dominated by dust emission from the dense cold cloud core surrounding the protostar. Including all array configurations, the emission is dominated by an elliptical source with a size of ~ 8" x 3". If we filter out the extended emission we find three compact mm-sources inside the elliptical core. The strongest one, $S_A$, coincides with the VLA/IRAC source and resolves into a double source at 1.4 mm, where we have sub-arcsecond resolution. The measured spectral index, $\alpha$, between 3 and 1.4 mm is ~ 2.3, and steeper at longer wavelengths, suggesting a low dust emissivity or that the dust is optically thick. We argue that the dust in these accretion disks is optically thick and estimate a mass of an accretion disk or infalling envelope surrounding S$_A$ to be ~ 60 solar masses.

Line spectra of 68 Taurus T Tauri stars were obtained with the Herschel-PACS (Photodetector Array Camera & Spectrometer) instrument as part of the GASPS (Gas Evolution in Protoplanetary Systems) survey of protoplanetary discs. A careful examination of the line scans centred on the [OI] 63.18 microns fine-structure line unveiled a line at 63.32 micron in some of these spectra. We identify this line with a transition of ortho-water. It is detected confidently (i.e., >3 sigma) in eight sources, i.e., 24% of the sub-sample with gas-rich discs. Several statistical tests were used to search for correlations with other disc and stellar parameters such as line fluxes of [OI] 6300 Armstrong and 63.18 microns; X-ray luminosity and continuum levels at 63 microns and 850 microns. Correlations are found between the water line fluxes and the [OI] 63.18 microns line luminosity, the dust continuum, and possibly with the stellar X-ray luminosity. This is the first time that this line of warm water vapour has been detected in protoplanetary discs. We discuss its origins, in particular whether it comes from the inner disc and/or disc surface or from shocks in outflows and jets. Our analysis favours a disc origin, with the observed water vapour line produced within 2-3AU from the central stars, where the gas temperature is of the order of 500-600 K.

We present detailed model fits to observations of the disc around the Herbig Ae star HD 163296. This well-studied object has an age of ~ 4 Myr, with evidence of a circumstellar disc extending out to ~ 540AU. We use the radiation thermo-chemical disc code ProDiMo to model the gas and dust in the circumstellar disc of HD 163296, and attempt to determine the disc properties by fitting to observational line and continuum data. These include new Herschel/PACS observations obtained as part of the open-time key program GASPS (Gas in Protoplanetary Systems), consisting of a detection of the [OI]63mic line and upper limits for several other far infrared lines. We complement this with continuum data and ground-based observations of the 12CO 3-2, 2-1 and 13CO J=1-0 line transitions, as well as the H2 S(1) transition. We explore the effects of stellar ultraviolet variability and dust settling on the line emission, and on the derived disc properties. Our fitting efforts lead to derived gas/dust ratios in the range 9-100, depending on the assumptions made. We note that the line fluxes are sensitive in general to the degree of dust settling in the disc, with an increase in line flux for settled models. This is most pronounced in lines which are formed in the warm gas in the inner disc, but the low excitation molecular lines are also affected. This has serious implications for attempts to derive the disc gas mass from line observations. We derive fractional PAH abundances between 0.007 and 0.04 relative to ISM levels. Using a stellar and UV excess input spectrum based on a detailed analysis of observations, we find that the all observations are consistent with the previously assumed disc geometry.

Despite its importance in the thermal-balance of the gas and in the determination of primeval planetary atmospheres, the chemistry in protoplanetary discs remains poorly constrained with only a handful of detected species. We observed the emission from disc around the Herbig Be star HD 100546 with the PACS instrument in the spectroscopic mode on board the Herschel Space Telescope as part of the Gas in Protoplanetary Systems (GASPS) programme and used archival data from the DIGIT programme to search for the rotational emission of CH+. We detected in both datasets an emission line centred at 72.16 micron that most likely corresponds to the J=5-4 rotational emission of CH+. The J=3-2 and 6-5 transitions are also detected albeit with lower confidence. Other CH+ rotational lines in the PACS observations are blended with water lines. A rotational diagram analysis shows that the CH+ gas is warm at 323 (+2320/-151) K with a mass of 3e-14-5e-12 M_Sun. We modelled the CH+ chemistry with the chemo-physical code ProDiMo using a disc density structure and grain parameters that match continuum observations and near- and mid-infrared interferometric data. The model suggests that CH+ is most abundant at the location of the disc rim at 10-13 AU from the star where the gas is warm, consistent with previous observations of hot CO gas emission.

We present new continuum and line observations, along with modelling, of the faint (6-8) Myr old T Tauri star ET Cha belonging to the eta Chamaeleontis cluster. We have acquired HERSCHEL/PACS photometric fluxes at 70 mic and 160 mic, as well as a detection of the [OI] 63 mic fine-structure line in emission, and derived upper limits for some other far-IR OI, CII, CO and o-H2O lines. The HERSCHEL data is complemented by new ANDICAM B-K photometry, new HST/COS and HST/STIS UV-observations, a non-detection of CO J=3-2 with APEX, re-analysis of a UCLES high-resolution optical spectrum showing forbidden emission lines like [OI] 6300A, [SII] 6731A and 6716A, and [NII] 6583A, and a compilation of existing broad-band photometric data. We used the thermo-chemical disk code ProDiMo and the Monte-Carlo radiative transfer code MCFOST to model the protoplanetary disk around ET Cha. Based on these models we can determine the disk dust mass Mdust = (2.E-8 - 5.E-8) Msun, whereas the total disk gas mass is found to be only little constrained, Mgas = (5.E-5 - 3.E-3) Msun. In the models, the disk extends from 0.022 AU (just outside of the co-rotation radius) to only about 10 AU. Larger disks are found to be inconsistent with the CO J=3-2 non-detection. The low velocity component of the [OI] 6300A emission line is consistent with being emitted from the inner disk. The model can also reproduce the line flux of H2 v=1-0 S(1) at 2.122 mic. An additional high-velocity component of the [OI] 6300A emission line, however, points to the existence of an additional jet/outflow of low velocity (40 - 65) km/s with mass loss rate ~1.E-9 Msun/yr. In relation to our low estimations of the disk mass, such a mass loss rate suggests a disk lifetime of only ~(0.05 - 3) Myr, substantially shorter than the cluster age. The evolutionary state of this unusual protoplanetary disk is discussed.

The Stratospheric Observatory for Infrared Astronomy (SOFIA), a joint U.S./German project, is a 2.5-meter infrared airborne telescope carried by a Boeing 747-SP that flies in the stratosphere at altitudes as high as 45,000 feet (13.72 km). This facility is capable of observing from 0.3 {\mu}m to 1.6 mm with an average transmission greater than 80 percent. SOFIA will be staged out of the NASA Dryden Flight Research Center aircraft operations facility at Palmdale, CA. The SOFIA Science Mission Operations (SMO) will be located at NASA Ames Research Center, Moffett Field, CA. First science flights began in 2010 and a full operations schedule of up to one hundred 8 to 10 hour flights per year will be reached by 2014. The observatory is expected to operate until the mid 2030's. SOFIAs initial complement of seven focal plane instruments includes broadband imagers, moderate-resolution spectrographs that will resolve broad features due to dust and large molecules, and high-resolution spectrometers capable of studying the kinematics of atomic and molecular gas at sub-km/s resolution. We describe the SOFIA facility and outline the opportunities for observations by the general scientific community and for future instrumentation development. The operational characteristics of the SOFIA first-generation instruments are summarized. The status of the flight test program is discussed and we show First Light images obtained at wavelengths from 5.4 to 37 \"im with the FORCAST imaging camera. Additional information about SOFIA is available at http://www.sofia.usra.edu and http://www.sofia.usra.edu/Science/docs/SofiaScienceVision051809-1.pdf

We present high signal-to-noise, moderate spectral resolution (R ~ 2000-2500) near-infrared (0.8-5.0 micron) spectroscopy of the nearby T Tauri star TW Hya. By comparing the spectrum and the equivalent widths of several atomic and molecular features with those for stars in the IRTF near-infrared library, we revise the spectral type to M2.5V, which is later than usually adopted (K7V). This implies a substantially cooler stellar temperature than previously assumed. Comparison with various pre-main sequence models suggests that TW Hya is only ~3 Myr old; much younger than the usually adopted 8 - 10 Myr. Analysis of the relative strengths of the H lines seen in the spectrum yields estimates for the temperature and density of the emitting region of T_e > 7500 K and n_e ~ 10^{12} - 10^{13} cm^{-3}. The thickness of the emitting region is 10^2 - 10^4 km and the covering fraction is f_\ast ~ 0.04. Our derived physical parameter values agree with the predictions of the magnetospheric accretion scenario. The highest signal-to-noise H lines have profiles that indicate multiple emission components. We derive an excess spectrum (above that of the M2.5V template) that peaks in the H band. Although our derived veiling values, ~ 0.1, agree with previous estimates, the excess spectrum does not match that of current models in which this flux is generated by an inner optically thin disk. We suggest that the excess flux spectrum instead reflects the differences in atmospheric opacity, gravity, and age between TW Hya and older, higher gravity field M2.5 dwarfs.

We have acquired sub-millimeter observations of 33 fields containing 37 Herbig Ae/Be (HAEBE) stars or potential HAEBE stars, including SCUBA maps of all but two of these stars. Nine target stars show extended dust emission. The other 18 are unresolved, suggesting that the dust envelopes or disks around these stars are less than a few arcseconds in angular size. In several cases we find that the strongest sub-millimeter emission originates from younger, heavily embedded sources rather than from the HAEBE star, which means that previous models must be viewed with caution. These new data, in combination with far-infrared flux measurements available in the literature, yield SEDs from far-infrared to millimeter wavelengths for all the observed objects. Isothermal fits to these SEDs demonstrate excellent fits, in most cases, to the flux densities longward of 100 {\mu}m. We find that a smaller proportion of B-type stars than A and F-type stars are surrounded by circumstellar disks, suggesting that disks around B stars dissipate on shorter time scales than those around later spectral types. Our models also reveal that the mass of the circumstellar material and the value of beta are correlated, with low masses corresponding to low values of beta. Since low values of beta imply large grain sizes, our results suggest that a large fraction of the mass in low-beta sources is locked up in very large grains. Several of the isolated HAEBE stars have disks with very flat sub-millimeter SEDs. These disks may be on the verge of forming planetary systems.

In an effort to simultaneously study the gas and dust components of the disc surrounding the young Herbig Ae star HD 169142, we present far-IR observations obtained with the PACS instrument onboard the Herschel Space Observatory. This work is part of the Open Time Key Project GASPS, which is aimed at studying the evolution of protoplanetary discs. To constrain the gas properties in the outer disc, we observed the star at several key gas-lines, including [OI] 63.2 and 145.5 micron, [CII] 157.7 micron, CO 72.8 and 90.2 micron, and o-H2O 78.7 and 179.5 micron. We only detect the [OI] 63.2 micron line in our spectra, and derive upper limits for the other lines. We complement our data set with PACS photometry and 12/13CO data obtained with the Submillimeter Array. Furthermore, we derive accurate stellar parameters from optical spectra and UV to mm photometry. We model the dust continuum with the 3D radiative transfer code MCFOST and use this model as an input to analyse the gas lines with the thermo-chemical code ProDiMo. Our dataset is consistent with a simple model in which the gas and dust are well-mixed in a disc with a continuous structure between 20 and 200 AU, but this is not a unique solution. Our modelling effort allows us to constrain the gas-to-dust mass ratio as well as the relative abundance of the PAHs in the disc by simultaneously fitting the lines of several species that originate in different regions. Our results are inconsistent with a gas-poor disc with a large UV excess; a gas mass of 5.0 +/- 2.0 times 10^(-3) Msun is still present in this disc, in agreement with earlier CO observations.

Context - Circumstellar discs are ubiquitous around young stars, but rapidly dissipate their gas and dust on timescales of a few Myr. The Herschel space observatory allows for the study of the warm disc atmosphere, using far-infrared spectroscopy to measure gas content and excitation conditions, and far-IR photometry to constrain the dust distribution. Aims - We aim to detect and characterize the gas content of circumstellar discs in four targets as part of the Herschel science demonstration phase. Methods - We carried out sensitive medium resolution spectroscopy and high sensitivity photometry at lambda ~60-190 micron using the Photodetector Array Camera and Spectrometer instrument on the Herschel space observatory. Results - We detect [OI] 63 micron emission from the young stars HD 169142, TW Hydrae, and RECX 15, but not HD 181327. No other lines, including [CII] 158 and [OI] 145, are significantly detected. All four stars are detected in photometry at 70 and 160 micron. Extensive models are presented in associated papers.

The Herschel GASPS Key Program is a survey of the gas phase of protoplanetary discs, targeting 240 objects which cover a large range of ages, spectral types, and disc properties. To interpret this large quantity of data and initiate self-consistent analyses of the gas and dust properties of protoplanetary discs, we have combined the capabilities of the radiative transfer code MCFOST with the gas thermal balance and chemistry code ProDiMo to compute a grid of 300 000 disc models (DENT). We present a comparison of the first Herschel/GASPS line and continuum data with the predictions from the DENT grid of models. Our objective is to test some of the main trends already identified in the DENT grid, as well as to define better empirical diagnostics to estimate the total gas mass of protoplanetary discs. Photospheric UV radiation appears to be the dominant gas-heating mechanism for Herbig stars, whereas UV excess and/or X-rays emission dominates for T Tauri stars. The DENT grid reveals the complexity in the analysis of far-IR lines and the difficulty to invert these observations into physical quantities. The combination of Herschel line observations with continuum data and/or with rotational lines in the (sub-)millimetre regime, in particular CO lines, is required for a detailed characterisation of the physical and chemical properties of circumstellar discs.

Planets are formed in disks around young stars. With an age of ~10 Myr, TW Hya is one of the nearest T Tauri stars that is still surrounded by a relatively massive disk. In addition a large number of molecules has been found in the TW Hya disk, making TW Hya the perfect test case in a large survey of disks with Herschel-PACS to directly study their gaseous component. We aim to constrain the gas and dust mass of the circumstellar disk around TW Hya. We observed the fine-structure lines of [OI] and [CII] as part of the Open-time large program GASPS. We complement this with continuum data and ground-based 12CO 3-2 and 13CO 3-2 observations. We simultaneously model the continuum and the line fluxes with the 3D Monte-Carlo code MCFOST and the thermo-chemical code ProDiMo to derive the gas and dust masses. We detect the [OI] line at 63 micron. The other lines that were observed, [OI] at 145 micron and [CII] at 157 micron, are not detected. No extended emission has been found. Preliminary modeling of the photometric and line data assuming [12CO]/[13CO]=69 suggests a dust mass for grains with radius < 1 mm of ~1.9 times 10^-4 Msun (total solid mass of 3 times 10^-3 Msun) and a gas mass of (0.5--5) times 10^-3 Msun. The gas-to-dust mass may be lower than the standard interstellar value of 100.

We present deep high angular resolution observations of the high-mass protostar NGC 7538S, which is in the center of a cold dense cloud core with a radius of 0.5 pc and a mass of ~2,000 Msun. These observations show that NGC 7538S is embedded in a compact elliptical core with a mass of 85 - 115 Msun. The star is surrounded by a rotating accretion disk, which powers a very young, hot molecular outflow approximately perpendicular to the rotating accretion disk. The accretion rate is very high, ~ 1.4 - 2.8 10^-3 Msun yr^-1. Evidence for rotation of the disk surrounding the star is seen in all largely optically thin molecular tracers, H13CN J = 1-0, HN13C J = 1-0, H13CO+ J = 1-0, and DCN J = 3-2. Many molecules appear to be affected by the hot molecular outflow, including DCN and H13CO+. The emission from CH3CN, which has often been used to trace disk rotation in young high-mass stars, is dominated by the outflow, especially at higher K-levels. Our new high-angular resolution observations show that the rotationally supported part of the disk is smaller than we previously estimated. The enclosed mass of the inner, rotationally supported part of the disk (D ~ 5", i.e 14,000 AU) is ~ 14 - 24 Msun.

We investigate the young (proto)stellar population in NGC 2023 and the L 1630 molecular cloud bordering the HII region IC 434, using Spitzer IRAC and MIPS archive data, JCMT SCUBA imaging and spectroscopy as well as targeted BIMA observations of one of the Class 0 protostars, NGC 2023 MM1. We have performed photometry of all IRAC and MIPS images, and used color-color diagrams to identify and classify all young stars seen within a 22'x26' field along the boundary between IC 434 and L 1630. For some stars, which have sufficient optical, IR, and/or sub-millimeter data we have also used the online SED fitting tool for a large 2D archive of axisymmetric radiative transfer models to perform more detailed modeling of the observed SEDs. We identify 5 sub-millimeter cores in our 850 and 450 micron SCUBA images, two of which have embedded class 0 or I protostars. Observations with BIMA are used to refine the position and characteristics of the Class 0 source NGC 2023 MM 1. These observations show that it is embedded in a very cold cloud core, which is strongly enhanced in NH2D. We find that HD 37903 is the most massive member of a cluster with 20 -- 30 PMS stars. We also find smaller groups of PMS stars formed from the Horsehead nebula and another elephant trunk structure to the north of the Horsehead. We refine the spectral classification of HD 37903 to B2 Ve. Our study shows that the expansion of the IC 434 HII region has triggered star formation in some of the dense elephant trunk structures and compressed gas inside the L 1630 molecular cloud. This pre-shock region is seen as a sub-millimeter ridge in which stars have already formed. The cluster associated with NGC 2023 is very young, and has a large fraction of Class I sources.

We investigate the young stellar population in and near the cometary globule Ori,I-2. The analysis is based on deep Nordic Optical Telescope R-band and H-alpha images, JCMT SCUBA 450 and 850 micron images combined with near-infrared 2MASS photometry and mid-infrared archival Spitzer images obtained with the IRAC (3.6, 4.5, 5.8 and 8 micron), and MIPS (24 and 70 micron) instruments. We identify a total of 125 sources within the 5'x5' region imaged by IRAC. Of these sources 87 are detected in the R-band image and 51 are detected in the 2MASS survey. The detailed physical properties of the sources are explored using a combination of near/mid-infrared color-color diagrams, greybody fitting of SEDs and an online SED fitting tool that uses a library of 2D radiation transfer based accretion models of young stellar objects with disks. Ori I-2 shows clear evidence of triggered star formation with four young low luminosity pre-main sequence stars embedded in the globule. At least two, possibly as many as four, additional low-mass PMS objects, were discovered in the field which are probably part of the young sigma-Orionis cluster. Among the PMS stars which have formed in the globule, MIR-54 is a young, deeply embedded Class 0/I object, MIR-51 and 52 are young Class II sources, while MIR-89 is a more evolved, heavily extincted Class II object with its apparent colors mimicking a Class 0/I object. The Class II source, MIR-52, which is strong a H-alpha emission line star appears to drive an outflow approximately aligned with the outflow from MIR-54, and because of the proximity of the two outflows, either star could contribute. MIR-89 appears to excite a low excitation HH object, HH 992, discovered for the first time in this study.

Analysis of high spatial resolution VLA images shows that the free-free emission from NGC7538 IRS1 is dominated by a collimated ionized wind. We have re-analyzed high angular resolution VLA archive data from 6 cm to 7 mm, and measured separately the flux density from the compact bipolar core and the extended (1.5" - 3") lobes. We find that the flux density of the core is proportional to the frequency to the power of alpha, with alpha being about 0.7. The frequency dependence of the total flux density is slightly steeper with alpha = 0.8. A massive optically thick hypercompact core with a steep density gradient can explain this frequency dependence, but it cannot explain the extremely broad recombination line velocities observed in this source. Neither can it explain why the core is bipolar rather than spherical, nor the observed decrease of 4% in the flux density in less than 10 years. An ionized wind modulated by accretion is expected to vary, because the accretion flow from the surrounding cloud will vary over time. BIMA and CARMA continuum observations at 3 mm show that the free-free emission still dominates at 3 mm. HCO+ J = 1 - 0 observations combined with FCRAO single dish data show a clear inverse P Cygni profile towards IRS1. These observations confirm that IRS1 is heavily accreting with an accretion rate of about 2 times 10(-4) solar masses per year.

We aim to study dust properties of massive star forming regions in the outer Galaxy, in a direction opposite to the Galactic center. We present observations of six outer Galaxy point sources IRAS 01045+6505, 01420+6401, 05271+3059, 05345+3556, 20222+3541 and 20406+4555, taken with the Submillimeter Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT) at 450 and 850 micron. Single temperature greybody models are fitted to the Spectral Energy Distribution of the detected sub-mm cores to derive dust temperature, dust emissivity index and optical depth at 250 micron. The observed radial intensity profiles of the sub-mm cores were fitted with power laws to derive the indices describing the density distribution. At a resolution of 15" all six IRAS point sources show multiple emission peaks. Only four out of fourteen detected sub-mm cores show associated mid-infrared emission. For the sub-mm cores we derive dust temperatures of 32+-5 K and dust emissivity indices between 0.9 and 2.5. The density profiles of the sub-mm cores can be fitted by a single power law distribution with indices -1.5+-0.3, with most cores showing an index of -1.5. This is consistent with most observations of massive star forming regions and supports predictions of models of star formation which consider non-thermal support against gravitational collapse.

We present multi-wavelength, high spatial resolution imaging of the IRS 7 region in the R Corona Australis molecular cloud. Our observations include 1.1 mm continuum and HCO^+ J = $3 \to 2$ images from the SMA, ^{12}CO J = $3 \to 2$ outflow maps from the DesertStar heterodyne array receiver on the HHT, 450 $\mu$m and 850 $\mu$m continuum images from SCUBA, and archival Spitzer IRAC and MIPS 24 \micron images. The accurate astrometry of the IRAC images allow us to identify IRS 7 with the cm source VLA 10W (IRS 7A) and the X-ray source X_W. The SMA 1.1 mm image reveals two compact continuum sources which are also distinguishable at 450 $\mu$m. SMA 1 coincides with X-ray source CXOU J190156.4-365728 and VLA cm source 10E (IRS 7B) and is seen in the IRAC and MIPS images. SMA 2 has no infrared counterpart but coincides with cm source VLA 9. Spectral energy distributions constructed from SMA, SCUBA and Spitzer data yield bolometric temperatures of 83 K for SMA 1 and $\leq$70 K for SMA 2. These temperatures along with the submillimeter to total luminosity ratios indicate that SMA 2 is a Class 0 protostar, while SMA 1 is a Class 0/Class I transitional object (L=$17\pm6$ \Lsun). The ^{12}CO J = $3 \to 2$ outflow map shows one major and possibly several smaller outflows centered on the IRS 7 region, with masses and energetics consistent with previous work. We identify the Class 0 source SMA 2/VLA 9 as the main driver of this outflow. The complex and clumpy spatial and velocity distribution of the HCO^+ J = $3 \to 2$ emission is not consistent with either bulk rotation, or any known molecular outflow activity.

We present wide-field near-infrared images of the DR21/W75 high-mass star forming region, obtained with the Wide Field Camera, WFCAM, on the United Kingdom Infrared Telescope. Broad-band JHK and narrow-band H2 1-0S(1) images are compared to archival mid-IR images from the Spitzer Space Telescope, and 850 micron dust-continuum maps obtained with the Submillimeter Common User Bolometer Array (SCUBA). Together these data give a complete picture of dynamic star formation across this extensive region, which includes at least four separate star forming sites in various stages of evolution. The H2 data reveal knots and bow shocks associated with more than 50 individual flows. Most are well collimated, and at least five qualify as parsec-scale flows. Most appear to be driven by embedded, low-mass protostars. The orientations of the outflows, particularly from the few higher-mass sources in the region (DR21, DR21(OH), W75N and ERO~1), show some degree of order, being preferentially orientated roughly orthogonal to the chain of dusty cores that runs north-south through DR21. Clustering may inhibit disk accretion and therefore the production of outflows; we certainly do not see enhanced outflow activity from clusters of protostars. Finally, although the low-mass protostellar outflows are abundant and widely distributed, the current generation does not provide sufficient momentum and kinetic energy to account for the observed turbulent motions in the DR21/W75 giant molecular clouds. Rather, multiple epochs of outflow activity are required over the million-year timescale for turbulent decay.

New high resolution observations of HCO+ J=1-0, H13CN J=1-0, SO 2,2 - 1,1, and continuum with BIMA at 3.4 mm show that the NGC7538 - IRS9 cloud core is a site of active ongoing star formation. Our observations reveal at least three young bipolar molecular outflows, all ~ 10,000 -- 20,000 years old. IRS9 drives a bipolar, extreme high velocity outflow observed nearly pole on. South of IRS9 we find a cold, protostellar condensation with a size of ~ 14" x 6" with a mass > 250 Msun. This is the center of one of the outflows and shows deep, red-shifted self absorption in HCO+, suggesting that there is a protostar embedded in the core, still in a phase of active accretion. This source is not detected in the far infrared, suggesting that the luminosity < 10^4 Lsun; yet the mass of the outflow is ~ 60 Msun. The red-shifted HCO+ self-absorption profiles observed toward the southern protostar and IRS9 predict accretion rates of a few times 10^-4 to 10^-3 Msun/yr. Deep VLA continuum observations at 3.6 cm show that IRS9 coincides with a faint thermal VLA source, but no other young star in the IRS9 region has any detectable free-free emission at a level of ~ 60 microJy at 3.6 cm. The HCO+ abundance is significantly enhanced in the hot IRS9 outflow. A direct comparison of mass estimates from HCO+ and CO for the well-characterized red-shifted IRS9 outflow predicts an HCO+ enhancement of more than a factor of 30, or [HCO+/H2] >= 6 10^-8.

We present deep images of dust continuum emission at 450, 800, and 850 micron of the dark cloud LDN 1689N which harbors the low-mass young stellar objects (YSOs) IRAS 16293-2422A and B (I16293A and I16293B) and the cold prestellar object I16293E. Toward the positions of I16293A and E we also obtained spectra of CO-isotopomers and deep submillimeter observations of chemically related molecules with high critical densities. To I16293A we report the detection of the HDO 1_01 - 0_00 and H2O 1_10 - 1_01 ground-state transitions as broad self-reversed emission profiles with narrow absorption, and a tentative detection of H2D+ 1_10 - 1_11. To I16293E we detect weak emission of subthermally excited HDO 1_01 - 0_00. Based on this set of submillimeter continuum and line data we model the envelopes around I16293A and E. The density and velocity structure of I16293A is fit by an inside-out collapse model, yielding a sound speed of a=0.7 km/s, an age of t=(0.6--2.5)e4 yr, and a mass of 6.1 Msun. The density in the envelope of I16293E is fit by a radial power law with index -1.0+/-0.2, a mass of 4.4 Msun, and a constant temperature of 16K. These respective models are used to study the chemistry of the envelopes of these pre- and protostellar objects. The [HDO]/[H2O] abundance ratio in the warm inner envelope of I16293A of a few times 1e-4 is comparable to that measured in comets. This supports the idea that the [HDO]/[H2O] ratio is determined in the cold prestellar core phase and conserved throughout the formation process of low-mass stars and planets.

We report the detection of a massive rotating disk around the high-mass Class 0 candidate NGC 7538S. The disk is well-resolved with BIMA (HPBW =3.7") in 3.4 mm continuum and in H13CN J = 1 - 0. It is seen nearly edge on and has a size of ~ 30,000 AU. A young, powerful outflow perpendicular to the rotating disk is mapped in SiO J = 2 - 1 and HCO+ J = 1 - 0. The dynamical age of the outflow is < 10,000 yr. The velocity gradient seen in H13CN is consistent with Keplerian rotation. Assuming that the gas is gravitationally bound, the mass of the central object is ~ 40 Msun. The mass of the continuum ``disk'' is <= 100 Msun and has a luminosity of 10^4 Lsun. H13CN gives a mass ~ 400 Msun for the rotating disk, and ~ 1000 Msun for the extended (20") envelope. Our observations confirm that this is an extremely massive protostar in its earliest stages.

New molecular (13CO J=3-2) and dust continuum (450 and 850 micron) SCUBA maps of the NGC7129 star forming region are presented, complemented by C18O J=3-2 spectra at several positions within the mapped region. The maps include the Herbig Ae/Be star LkHalpha 234, the far-infrared source NGC 7129 FIRS2 and several other pre-stellar sources embedded within the molecular ridge. The SCUBA maps help us understand the nature of the pre-main sequence stars in this actively star forming region. A deeply embedded submillimeter source, SMM2, not clearly seen in any earlier data set, is shown to be a pre-stellar core or possibly a protostar. The highest continuum peak emission is identified with the deeply embedded source IRS6, a few arcseconds away from LkHalpha 234, and also responsible for both the optical jet and the molecular outflow. The gas and dust masses are found to be consistent, suggesting little or no CO depletion onto grains. The dust emissivity index is lower towards the dense compact sources, beta ~1 - 1.6, and higher, beta ~ 2.0, in the surrounding cloud, implying small size grains in the PDR ridge, whose mantles have been evaporated by the intense UV radiation.

This paper examines the density and velocity structure of envelopes around young stellar objects through submillimeter continuum imaging of four objects in Taurus and previously obtained molecular-line data. Observations carried out with the SCUBA on the JCMT at 850 and 450 micron of L1489 IRS, L1535 IRS, L1527 IRS, and TMC 1 reveal ~2000 AU elongated structures embedded in extended envelopes. The density distribution in these envelopes is equally well fit by a radial power-law of index p=1.0-2.0 or with a collapse model such as that of Shu (1997: ApJ, 214, 488). This inside-out collapse model predicts 13CO, C18O, HCO+, and H13CO+ line profiles which closely match observed spectra toward three of our four sources. This shows that the inside-out collapse model offers a good description of YSO envelopes, but also that reliable constraints on its parameters require independent measurements of the density and the velocity structure, e.g., through continuum and line observations. For the remaining source, L1489 IRS, we find that a model consisting of a 2000 AU radius, rotating, disk-like structure better describes the data. Possibly, this source is in transition between the embedded Class I and the optically revealed T Tauri phases. Two apparently starless cores are found at ~10,000 AU from L1489 IRS and L1535 IRS. They are cold, 10-15 K, contain 0.5-3.0 M_sol, and have flat density distributions characterized by a Gaussian of ~10,000 AU FWHM. The proximity of these cores shows that star formation in truly isolated cores is rare even in Taurus.

III Zw 2 is a spiral galaxy with an optical spectrum and faint extended radio structure typical of a Seyfert galaxy, but also with an extremely variable, blazar-like radio core. We have now discovered a new radio flare where the source has brightened more than twenty-fold within less than two years. A broad-band radio spectrum between 1.4 and 666 GHz shows a textbook-like synchrotron spectrum peaking at 43 GHz, with a self-absorbed synchrotron spectral index +2.5 at frequencies below 43 GHz and an optically thin spectral index -0.75 at frequencies above 43 GHz. The outburst spectrum can be well fitted by two homogenous, spherical components with equipartition sizes of 0.1 and 0.2 pc at 43 and 15 GHz, and with magnetic fields of 0.4 and 1 Gauss. VLBA observations at 43 GHz confirm this double structure and these sizes. Time scale arguments suggest that the emitting regions are shocks which are continuously accelerating particles. This could be explained by a frustrated jet scenario with very compact hotspots. Similar millimeter-peaked spectrum (MPS) sources could have escaped our attention because of their low flux density at typical survey frequencies and their strong variability.