We have used the new 90GHz MUSTANG camera on the Robert C. Byrd Green Bank Telescope (GBT) to map the bright Huygens region of the star-forming region M42 with a resolution of 9" and a sensitivity of 2.8mJy/beam. 90GHz is an interesting transition frequency, as MUSTANG detects both the free-free emission characteristic of the HII region created by the Trapezium stars, normally seen at lower frequencies, and thermal dust emission from the background OMC1 molecular cloud, normally mapped at higher frequencies. We also present similar data from the 150GHz GISMO camera taken on the IRAM telescope. This map has 15" resolution. By combining the MUSTANG data with 1.4, 8, and 21GHz radio data from the VLA and GBT, we derive a new estimate of the emission measure (EM) averaged electron temperature of Te = 11376K by an original method relating free-free emission intensities at optically thin and optically thick frequencies. Combining ISO-LWS data with our data, we derive a new estimate of the dust temperature and spectral emissivity index within the 80" ISO-LWS beam toward OrionKL/BN, Td = 42K and Beta=1.3. We show that both Td and Beta decrease when going from the HII region and excited OMC1 interface to the denser UV shielded part of OMC1 (OrionKL/BN, Orion S). With a model consisting of only free-free and thermal dust emission we are able to fit data taken at frequencies from 1.5GHz to 854GHz.

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 present new (April 2000) MERLIN observations of the H2O masers located near the protostar Cepheus A HW2. The MERLIN observations detect many of the structures found in earlier (1996) VLBA observations of Torrelles and collaborators, and the changed positions of these structures are compatible with the VLBA proper motions and astrometric uncertainties. The radius of curvature of the R4 structure of maser arcs appears to have grown by a factor of two, and the displacement of the arcs between 1996 and 2000 are compatible with expansion about a common center. In addition, the MERLIN observations detect red-shifted masers not previously found; taken with the newly discovered masers, the R4 structure now resembles patchy emission from an elliptical ring. We demonstrate that a simple bow-shock model cannot simultaneously account for the shape and the velocity gradient of the R4 structure. A model involving a slow, hydromagnetic shock propagating into a rotating, circumstellar disk better describes the maser spot kinematics and luminosities. In this model, the central mass is 3 solar masses, and we demonstrate that the mass of the disk is negligible in comparison. The expansion velocity of the post-shock gas, roughly 5 km / s, is slow compared to the average shock velocity (roughly 13 km / s) suggesting that the post-shock gas is magnetically supported with a characteristic field strength of roughly 30 mG. We speculate that the expanding maser rings R4 and R5 may be generated by periodic, instability-driven winds from young stars that periodically send spherical shocks into the surrounding circumstellar material.

Standard stellar evolution theory is inconsistent with the observed isotopic carbon ratio, 12C/13C, in evolved stars. This theory is also inconsistent with the 3He/H abundance ratios observed in Galactic HII regions, when combined with chemical evolution theory. These discrepancies have been attributed to an extra, non-standard mixing which further processes material during the RGB and should lower both the 12C/13C and 3He/H abundance ratios for stars with masses < 2 solar masses. Measurements of isotopic ratios in planetary nebulae probe material which escapes the star to be further processed by future generations of stars. We have measured the carbon isotopic abundance ratio, 12C/13C, in 11 planetary nebulae (PNe) by observing the J=2-->1 and J=3-->2 millimeter transitions of 12CO and 13CO in molecular clouds associated with the PNe. A large velocity gradient (LVG) model has been used to determine the physical conditions for each PNe where both transitions have been detected. We detect both 12CO and 13CO in 9 PNe. If 12CO/13CO = 12C/13C, the range of 12C/13C is 2.2--31. Our results support theories which include some form of extra mixing.

We report the first extragalactic detection of the neutral carbon [CI] 3P2-3P1 fine structure line at 809 GHz. The line was observed towards M82 simultaneously with the 3P1-3P0 line at 492 GHz, providing a precise measurement of the J=2-1/J=1-0 integrated line ratio of 0.96 (on a [K km s^-1] -scale). This ratio constrains the [CI] emitting gas to have a temperature of at least 50 K and a density of at least 10^4 cm^-3. Already at this minimum temperature and density, the beam averaged CI-column density is large, 2.1 10^18 cm^-2, confirming the high CI/CO abundance ratio of approximately 0.5 estimated earlier from the 492 GHz line alone. We argue that the [CI] emission from M82 most likely arises in clouds of linear size around a few pc with a density of about 10^4 cm^-3 or slightly higher and temperatures of 50 K up to about 100 K.