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
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 22x26 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.
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.
