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 pr
obe 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
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] 15
8 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.
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 unresolve
d, 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.
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 co
mpact 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.
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.
