Martin Houde

Dispersion of Magnetic Fields in Molecular Clouds

We describe a method for determining the dispersion of magnetic field vectors about large-scale fields in turbulent molecular clouds. The method is designed to avoid inaccurate estimates of magnetohydrodynamic or turbulent dispersion - and help avoiding inaccurate estimates of field strengths - due to large-scale, non-turbulent field structure when using the well-known method of Chandrasekhar and Fermi. Our method also provides accurate, independent estimates of the turbulent to large-scale magnetic field strength ratio. We discuss applications to the molecular clouds OMC-1, M17, and DR21(Main).

The James Clerk Maxwell telescope legacy survey of nearby star-forming regions in the gould belt

This paper describes a James Clerk Maxwell Telescope (JCMT) legacy survey that has been awarded roughly 500 hr of observing time to be carried out from 2007 to 2009. In this survey, we will map with SCUBA-2 (Submillimetre Common-User Bolometer Array 2) almost all of the well-known low-mass and intermediate-mass star-forming regions within 0.5 kpc that are accessible from the JCMT. Most of these locations are associated with the Gould Belt. From these observations, we will produce a flux-limited snapshot of star formation near the Sun, providing a legacy of images, as well as point-source and extended-source catalogs, over almost 700 deg(2) of sky. The resulting images will yield the first catalog of prestellar and protostellar sources selected by submillimeter continuum emission, and should increase the number of known sources by more than an order of magnitude. We will also obtain with the array receiver HARP (Heterodyne Array Receiver Program) CO maps, in three CO isotopologues, of a large typical sample of prestellar and protostellar sources. We will then map the brightest hundred sources with the SCUBA-2 polarimeter (POL-2), producing the first statistically significant set of polarization maps in the submillimeter. The images and source catalogs will be a powerful reference set for astronomers, providing a detailed legacy archive for future telescopes, including ALMA, Herschel, and JWST.

ALIGNMENT BETWEEN FLATTENED PROTOSTELLAR INFALL ENVELOPES AND AMBIENT MAGNETIC FIELDS

We present 350 μm 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 young stellar object inner envelope (aka pseudodisk). Furthermore, the field should exhibit a pinched or hourglass-shaped morphology as gravity drags the field inward toward 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 μm 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°). 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.

PROBING THE TURBULENCE DISSIPATION RANGE AND MAGNETIC FIELD STRENGTHS IN MOLECULAR CLOUDS

We study the turbulent velocity dispersion spectra of the coexistent HCN and HCO+ molecular species as a function of length scale in the M17 star-forming molecular cloud. We show that the observed downward shift of the ions spectrum relative to that of the neutral is readily explained by the existence of an ambipolar diffusion range within which ion and neutral turbulent energies dissipate differently. We use these observations to evaluate this decoupling scale and show how to estimate the strength of the plane-of-the-sky component of the embedded magnetic field in a completely novel way.

PROBING THE MAGNETIC FIELD WITH MOLECULAR ION SPECTRA. II.

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.

On the Measurement of the Magnitude and Orientation of the Magnetic Field in Molecular Clouds

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 μm polarimetry map and HCO+-to-HCN molecular line width ratios to provide the first three-dimensional view of the magnetic field in M17.

DISPERSION OF MAGNETIC FIELDS IN MOLECULAR CLOUDS. III.

We apply our technique on the dispersion of magnetic fields in molecular clouds to high spatial resolution Submillimeter Array polarization data obtained for Orion KL in OMC-1, IRAS 16293, and NGC 1333 IRAS 4A. We show how one can take advantage of such high resolution data to characterize the magnetized turbulence power spectrum in the inertial and dissipation ranges. For Orion KL we determine that in the inertial range the spectrum can be approximately fitted with a power law k^-(2.9\pm0.9) and we report a value of 9.9 mpc for {\lambda}_AD, the high spatial frequency cutoff presumably due to turbulent ambipolar diffusion. For the same parameters we have \sim k^-(1.4\pm0.4) and a tentative value of {\lambda}_AD \simeq 2.2 mpc for NGC 1333 IRAS 4A, and \sim k^-(1.8\pm0.3) with an upper limit of {\lambda}_AD < 1.8 mpc for IRAS 16293. We also discuss the application of the technique to interferometry measurements and the effects of the inherent spatial filtering process on the interpretation of the results.

DAMPING OF MAGNETOHYDRODYNAMIC TURBULENCE IN PARTIALLY IONIZED GAS AND THE OBSERVED DIFFERENCE OF VELOCITIES OF NEUTRALS AND IONS

Theoretical and observational studies on the turbulence of the interstellar medium developed fast in the past decades. The theory of supersonic magnetized turbulence, as well as the understanding of the projection effects of observed quantities, is still in progress. In this work, we explore the characterization of the turbulent cascade and its damping from observational spectral line profiles. We address the difference of ion and neutral velocities by clarifying the nature of the turbulence damping in the partially ionized. We provide theoretical arguments in favor of the explanation of the larger Doppler broadening of lines arising from neutral species compared to ions as arising from the turbulence damping of ions at larger scales. Also, we compute a number of MHD numerical simulations for different turbulent regimes and explicit turbulent damping, and compare both the three-dimensional distributions of velocity and the synthetic line profile distributions. From the numerical simulations, we place constraints on the precision with which one can measure the three-dimensional dispersion depending on the turbulence sonic Mach number. We show that no universal correspondence between the three-dimensional velocity dispersions measured in the turbulent volume and minima of the two-dimensional velocity dispersions available through observations exist. For instance, for subsonic turbulence the correspondence is poor at scales much smaller than the turbulence injection scale, while for supersonic turbulence the correspondence is poor for the scales comparable with the injection scale. We provide a physical explanation of the existence of such a two-dimensional to three-dimensional correspondence and discuss the uncertainties in evaluating the damping scale of ions that can be obtained from observations. However, we show that the statistics of velocity dispersion from observed line profiles can provide the spectral index and the energy transfer rate of turbulence. Also, by comparing two similar simulations with different viscous coefficients, it was possible to constrain the turbulent cut-off scale. This may especially prove useful since it is believed that ambipolar diffusion may be one of the dominant dissipative mechanisms in star-forming regions. In this case, the determination of the ambipolar diffusion scale may be used as a complementary method for the determination of magnetic field intensity in collapsing cores. We discuss the implications of our findings in terms of a new approach to magnetic field measurement proposed by Li & Houde.

The Alignment of the Magnetic Field and Collimated Outflows in Star-forming Regions: The Case of NGC 2071

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.

SIMULTANEOUS DETERMINATION OF THE COSMIC RAY IONIZATION RATE AND FRACTIONAL IONIZATION IN DR 21(OH)

We present a new method for the simultaneous calculation of the cosmic ray ionization rate, ζH2, and the ionization fraction, χ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 DR 21(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 Xe directly from observations, and estimate ζH2 accordingly. Our results, ζH2 = 3.1 x 10^-18 s^-1 and χe = 3.2 x 10^-8, are consistent with recent determinations in other objects.