The tori around active galactic nuclei (AGN) are potential formation sites for large amounts of dust, and they may help resolve the so-called dust budget crisis at high redshift. We investigate the dust composition in 53 of the 87 Palomar Green (PG) quasars showing the 9.7 micron silicate feature in emission. By simultaneously fitting the mid-infrared spectroscopic features and the underlying continuum, we estimate the mass fraction in various amorphous and crystalline dust species. We find that the dust consists predominantly of alumina and amorphous silicates, with a small fraction in crystalline form. The mean crystallinity is 8 +/- 6%, with more than half of the crystallinities greater than 5%, well above the upper limit determined for the Galaxy. Higher values of crystallinity are found for higher oxide fractions and for more luminous sources.

We present 350 micron polarization observations of four low-mass cores containing Class 0 protostars: L483, L1157, L1448-IRS2, and Serp-FIR1. This is the second paper in a larger survey aimed at testing magnetically regulated models for core-collapse. One key prediction of these models is that the mean magnetic field in a core should be aligned with the symmetry axis (minor axis) of the flattened YSO inner envelope (aka pseudodisk). Furthermore, the field should exhibit a pinched or hour-glass shaped morphology as gravity drags the field inward towards the central protostar. We combine our results for the four cores with results for three similar cores that were published in the first paper from our survey. An analysis of the 350 micron polarization data for the seven cores yields evidence of a positive correlation between mean field direction and pseudodisk symmetry axis. Our rough estimate for the probability of obtaining by pure chance a correlation as strong as the one we found is about 5%. In addition, we combine together data for multiple cores to create a source-averaged magnetic field map having improved signal-to-noise ratio, and this map shows good agreement between mean field direction and pseudodisk axis (they are within 15 degrees). We also see hints of a magnetic pinch in the source-averaged map. We conclude that core-scale magnetic fields appear to be strong enough to guide gas infall, as predicted by the magnetically regulated models. Finally, we find evidence of a positive correlation between core magnetic field direction and bipolar outflow axis.

We present a multi-wavelength observational study towards the high-mass young stellar object G8.68-0.37. A single massive gas-and-dust core is observed in the (sub)millimeter continuum and molecular line emissions. We fitted the spectral energy distribution (SED) from the dust continuum emission. The best-fit SED suggests the presence of two components with temperature of $T_{\rm d}=20$ K and 120 K, respectively. The core has a total mass of up to $1.5\times10^3$ $M_{\odot}$ and bolometric luminosity of $2.3\times10^4 L_{\odot}$. Both the mass and luminosity are dominated by the cold component ($T_{\rm d}=20$ K). The molecular lines of C$^{18}$O, C$^{34}$S, DCN, and thermally excited CH$_3$OH are detected in this core. Prominent infall signatures are observed in the $^{12}$CO $(1-0)$ and $(2-1)$. We estimated an infall velocity of 0.45 km s$^{-1}$ and mass infall rate of $7\times10^{-4} M_{\odot}$ year$^{-1}$. From the molecular lines, we have found a high DCN abundance and relative abundance ratio to HCN. The overabundant DCN may originate from a significant deuteration in the previous cold pre-protostellar phase. And the DCN should now be rapidly sublimated from the grain mantles to maintain the overabundance in the gas phase.

We report the first detection of hydrogen fluoride (HF) toward a high redshift quasar. Using the Caltech Submillimeter Observatory (CSO) we detect the HF J = 1 - 0 transition in absorption toward the Cloverleaf, a broad absorption line (BAL) quasi-stellar object (QSO) at z=2.56. The detection is statistically significant at the ~ 6 sigma level. We estimate a lower limit of 4 \times 1014 cm-2 for the HF column density and using a previous estimate of the hydrogen column density, we obtain a lower limit of 1.7 \times 10-9 for the HF abundance. This value suggests that, assuming a Galactic N(HF)/NH ratio, HF accounts for at least ~10% of the fluorine in the gas phase along the line of sight to the Cloverleaf quasar. This observation corroborates the prediction that HF should be a good probe of the molecular gas at high redshift. Measurements of the HF abundance as a function of redshift are urgently needed to better constrain the fluorine nucleosynthesis mechanism(s).

The nuclei of active comets emit molecules anisotropically from discrete vents. As the nucleus rotates, we expect to observe periodic variability in the molecular emission line profiles, which can be studied through mm/submm spectroscopy. Using this technique we investigated the HCN atmosphere of comet 103P/Hartley 2, the target of NASA's EPOXI mission, which had an exceptionally favorable apparition in late 2010. We detected short-term evolution of the spectral line profile, which was stimulated by the nucleus rotation, and which provides evidence for rapid deceleration and excitation of the rotation state. The measured rate of change in the rotation period is +1.00 \pm 0.15 min per day and the period itself is 18.32 \pm 0.03 hr, both applicable at the epoch of the EPOXI encounter. Surprisingly, the spin-down efficiency is lower by two orders of magnitude than the measurement in comet 9P/Tempel 1 and the best theoretical prediction. This secures rotational stability of the comet's nucleus during the next few returns, although we anticipate a catastrophic disruption from spin-up as its ultimate fate.

We report new observations of the fundamental $J=1-0$ transition of HCl (at 625.918GHz) toward a sample of 25 galactic star-forming regions, molecular clouds, and evolved stars, carried out using the Caltech Submillimeter Observatory. Fourteen sources in the sample are also observed in the corresponding H\tscl\ $J=1-0$ transition (at 624.978GHz). We have obtained clear detections in all but four of the targets, often in emission. Absorptions against bright background continuum sources are also seen in nine cases, usually involving a delicate balance between emission and absorption features. From RADEX modeling, we derive gas densities and HCl column densities for sources with HCl emission. HCl is found in a wide range of environments, with gas densities ranging from $10^5$ to $10^7$~cm$^{-3}$. The HCl abundance relative to H$_2$ is in the range of $(3-30)\times10^{-10}$. Comparing with the chlorine abundance in the solar neighborhood, this corresponds to a chlorine depletion factor of up to $\sim$400, assuming that HCl accounts for one third of the total chlorine in the gas phase. The [\tfcl]/[\tscl] isotopic ratio is rather varied, from unity to $\sim$5, mostly lower than the terrestrial value of 3.1. Such variation is highly localized, and could be generated by the nucleosynthesis in supernovae, which predicts a \tscl\ deficiency in most models. The lower ratios seen in W3IRS4 and W3IRS5 likely confine the progenitors of the supernovae to stars with relatively large mass ($\ga$25M$_\sun$) and high metallicity (Z$\sim$0.02).

We used the Extended Submillimeter Array (eSMA) in its most extended configuration to investigate the innermost (within a radius of 290 R* from the star) circumstellar envelope (CSE) of IRC+10216. We imaged the CSE using HCN and other molecular lines with a beam size of 0."22 x 0."46, deeply into the very inner edge (15 R*) of the envelope where the expansion velocity is only 3 km/s. The excitation mechanism of hot HCN and KCl maser lines is discussed. HCN maser components are spatially resolved for the first time on an astronomical object. We identified two discrete regions in the envelope: a region with a radius of . 15 R*, where molecular species have just formed and the gas has begun to be accelerated (region I) and a shell region (region II) with a radius of 23 R* and a thickness of 15 R*, whose expansion velocity has reached up to 13 km/s, nearly the terminal velocity of 15 km/s. The Si$^{34}$S line detected in region I shows a large expansion velocity of 16 km/s due to strong wing components, indicating that the emission may arise from a shock region in the innermost envelope. In region II, the P.A. of the most copious mass loss direction was found to be 120 +/- 10 degrees, which may correspond to the equatorial direction of the star. Region II contains a torus-like feature. These two regions may have emerged due to significant differences in the size distributions of the dust particles in the two regions.

Context. The high degree of deuteration observed in some prestellar cores depends on the ortho-to-para H2 ratio through the H3+ fractionation. Aims. We want to constrain the ortho/para H2 ratio across the L183 prestellar core. This is mandatory to correctly describe the deuter- ation amplification phenomenon in depleted cores such as L183 and to relate the total (ortho+para) H2D+ abundance to the sole ortho-H2D+ column density measurement. Methods. To constrain this ortho/para H2 ratio and derive its profile, we make use of the N2D+ /N2H+ ratio and of the ortho-H2D+ observations performed across the prestellar core. We use two simple chemical models limited to an almost totally depleted core description. New dissociative recombination and trihydrogen cation-dihydrogen reaction rates (including all isotopologues) are presented in this paper and included in our models. Results. We estimate the H2D+ ortho/para ratio in the L183 cloud, and constrain the H2 ortho/para ratio : we show that it is varying across the prestellar core by at least an order of magnitude being still very high (~0.1) in most of the cloud. Our time-dependent model indicates that the prestellar core is presumably older than 1.5-2 x 10^5 years but that it may not be much older. We also show that it has reached its present density only recently and that its contraction from a uniform density cloud can be constrained. Conclusions. A proper understanding of deuteration chemistry cannot be attained without taking into account the whole ortho/para family of molecular hydrogen and trihydrogen cation isotopologues as their relations are of utmost importance in the global scheme. Tracing the ortho/para H2 ratio should also give useful constraints on the dynamical evolution of prestellar cores.

To study the evolution of high mass cores, we have searched for evidence of collapse motions in a large sample of starless cores in the Orion molecular cloud. We used the Caltech Submillimeter Observatory telescope to obtain spectra of the optically thin (\H13CO+) and optically thick (\HCO+) high density tracer molecules in 27 cores with masses $>$ 1 \Ms. The red- and blue-asymmetries seen in the line profiles of the optically thick line with respect to the optically thin line indicate that 2/3 of these cores are not static. We detect evidence for infall (inward motions) in 9 cores and outward motions for 10 cores, suggesting a dichotomy in the kinematic state of the non-static cores in this sample. Our results provide an important observational constraint on the fraction of collapsing (inward motions) versus non-collapsing (re-expanding) cores for comparison with model simulations.

We present a new method for the simultaneous calculation of the cosmic ray ionization rate, zeta(H2), and the ionization fraction, chi(e), in dense molecular clouds. A simple network of chemical reactions dominant in the creation and destruction of HCNH+ and HCO+ is used in conjunction with observed pairs of rotational transitions of several molecular species in order to determine the electron abundance and the H3+ abundance. The cosmic ray ionization rate is then calculated by taking advantage of the fact that, in dark clouds, it governs the rate of creation of H3+. We apply this technique to the case of the star-forming region DR21(OH), where we successfully detected the (J=3-2) and (J=4-3) rotational transitions of HCNH+. We also determine the C and O isotopic ratios in this source to be 12C/13C=63+-4 and 16O/18O=318+-64, which are in good agreement with previous measurements in other clouds. The significance of our method lies in the ability to determine N(H3+) and chi(e) directly from observations, and estimate zeta(H2) accordingly. Our results, zeta(H2)=3.1x10^(-18) 1/s and chi(e)=3.2x10^(-8), are consistent with recent determinations in other objects.

Asteroid 3200 Phaethon resembles a comet in some ways, including a highly-eccentric orbit (e=0.89) and a strong associated meteor shower (the Geminids). Yet this object has never been observed to exhibit any cometary activity, i.e., gas production. We observed 3200 Phaethon with the Caltech Submillimeter Observatory on two occasions, once while it was near its closest approach to Earth as it neared perihelion, and another while it was further from Earth post-perihelion. Observations of the J=2-1 and J=3-2 rotational transitions of 12CO, typically strong lines in comets and indicative of gas production, yielded no detection. Upper limits on the 12CO production of 1.8e28 molecules/s and 7.6e28 molecules/s for Phaethon were determined on these two occasions.

Prestellar cores are unique laboratories for studies of the chemical and physical conditions preceding star formation. We observed the prestellar core L1544 in the fundamental transition of ortho-H2D+ (1_1,0-1_1,1) at different positions over 100", and found a strong correlation between its abundance and the CO depletion factor. We also present a tentative detection of the fundamental transition of para-D2H+ (1_1,0-1_0,1) at the dust emission peak. Maps in N2H+, N2D+, HCO+ and DCO+ are used, and interpreted with the aid of a spherically symmetric chemical model that predicts the column densities and abundances of these species as a function of radius. The correlation between the observed deuterium fractionation of H3+, N2H+ and HCO+ and the observed integrated CO depletion factor across the core can be reproduced by this chemical model. In addition a simpler model is used to study the H2D+ ortho-to-para ratio. We conclude that, in order to reproduce the observed ortho-H2D+ observations, the grain radius should be larger than 0.3 microns.

We use extensive 350 um polarimetry and continuum maps obtained with Hertz and SHARC II along with HCN and HCO+ spectroscopic data to trace the orientation of the magnetic field in the Orion A star-forming region. Using the polarimetry data, we find that the direction of the projection of the magnetic field in the plane of the sky relative to the orientation of the integral-shaped filament varies considerably as one moves from north to south. While in IRAS 05327-0457 and OMC-3 MMS 1-6 the projection of the field is primarily perpendicular to the filament it becomes better aligned with it at OMC-3 MMS 8-9 and well aligned with it at OMC-2 FIR 6. The OMC-2 FIR 4 cloud, located between the last two, is a peculiar object where we find almost no polarization. The projected angle of the field is more complicated in OMC-1 where it exhibits smooth variations in its orientation across the face of this massive complex. By combining the polarimetry and spectroscopic data we were able to measure a set of average values for the inclination angle of the magnetic field relative to the line of sight. We find that the field is oriented quite close to the plane of the sky in most places. More precisely, the inclination of the magnetic field is ~73 deg. around OMC-3 MMS 6, ~74 deg. at OMC-3 MMS 8-9, ~80 deg. at OMC-2 FIR 4, ~65 deg. in the northeastern part of OMC-1, and ~49 deg. in the Bar. We also present polarimetry data for the OMC-4 region located some 13 arcminutes south of OMC-1.

We demonstrate that the combination of Zeeman, polarimetry and ion-to-neutral molecular line width ratio measurements permits the determination of the magnitude and orientation of the magnetic field in the weakly ionized parts of molecular clouds. Zeeman measurements provide the strength of the magnetic field along the line of sight, polarimetry measurements give the field orientation in the plane of the sky and the ion-to-neutral molecular line width ratio determines the angle between the magnetic field and the line of sight. We apply the technique to the M17 star-forming region using a HERTZ 350 um polarimetry map and HCO+-to-HCN molecular line width ratios to provide the first three-dimensional view of the magnetic field in M17.

The magnetic field is believed to play a crucial role in the process of star formation. From the support it provides during the initial collapse of molecular clouds to the creation of strong collimated jets responsible for large mass losses, current theories predict its importance in many different stages during the formation of stars. Here we report on observational evidence which tests one aspect that can be inferred from these theories: the alignment between the local magnetic field and collimated bipolar outflows in such environments. There is good evidence of an alignment in the case of NGC 2071.

We present further observational evidence in support of our earlier proposal (Houde et al. 2000) for detecting the presence of the magnetic field in molecular clouds by comparing spectra of molecular ions with those of neutral molecules. The ion lines tend to be narrower and do not show the wings due to flows, when the magnetic field is sufficiently strong. We obtained spectra for the optically thin lines of the H13CN and H13CO+ species in a sample of ten molecular clouds and found the results to be in agreement with our previous observations of the main isotopic species, HCN and HCO+, made in OMC-1, OMC-2, OMC-3 and DR21OH, thus eliminating the possibility of optical depth effects playing a role in the ion line narrowing. HCS+ was also detected in four of these star forming regions. We also discuss previously published results by (Benson et al. 1998) of N2H+ detections in a large sample of dark clouds. We show that the similarity in line widths between ion and neutral species in their sample is consistent with the relatively small amount of turbulence and other flows observed in these clouds.

Observations of the effect of the magnetic field on its environment are usually achieved with techniques which rely on the interaction with the spin of the particles under study. Because of the relative weakness of this effect, extraction of the field characteristics proves to be a most challenging task. We take a totally different approach to the problem and show that the manifestation of the magnetic field can be directly observed by means of a comparison of the spectra of molecular ions with those of neutral molecules. This takes advantage of the strong cyclotron interaction between the ions and the field, but requires the presence of flows or turbulent motion in the gas. We compare our theory to data obtained on the OMC-1, OMC-2, OMC-3 and DR21(OH) molecular clouds.

The cumulene carbenes are important components of hydrocarbon chemistry in low mass star forming cores. Here we report the first astronomical detection of the long chain cumulene carbene H_2C_6 in the interstellar cloud TMC1, from observations of two of its rotational transitions: $J_{K,K'} = 7_{1,7} \rightarrow 6_{1,6}$ at 18.8 GHz and $8_{1,8} \rightarrow 7_{1,7}$ at 21.5 GHz, using NASA's Deep Space Network 70 m antenna at Goldstone, California. In addition we also observed the shorter cumulene carbene, H_2C_4 at the same position. The fractional abundance of H_2C_6 relative to H_2 is about $4.7 \times 10^{-11}$ and H_2C_4 is about $1.1 \times 10^{-9}$. The abundance of H_2C_6 is in fairly good agreement with gas phase chemical models for young molecular cloud cores, but the abundance of H_2C_4 is significantly larger than predicted.