We examine the role of the interstellar magnetic field to modulate the orientation of turbulent flows within the Taurus molecular cloud using spatial gradients of thin velocity slices of 12CO and 13CO antenna temperatures. Our analysis accounts for the random errors of the gradients that arise from the thermal noise of the spectra. The orientations of the vectors normal to the antenna temperature gradient vectors are compared to the magnetic field orientations that are calculated from Planck 353~GHz polarization data. These relative orientations are parameterized with the projected Rayleigh statistic and mean resultant vector. For 12CO, 28 percent and 39 percent of the cloud area exhibit strongly parallel or strongly perpendicular relative orientations respectively. For the lower opacity 13CO emission, strongly parallel and strongly perpendicular orientations are found in 7 percent and 43 percent of the cloud area respectively. For both isotopologues, strongly parallel or perpendicular alignments are restricted to localized regions with low levels of turbulence. If the relative orientations serve as an observational proxy to the Alfvenic Mach number then our results imply local variations of the Alfvenic Mach number throughout the cloud.

We explore the relation between the stellar mass surface density and the mass surface density of molecular hydrogen gas in twelve nearby molecular clouds that are located at $<$1.5 kpc distance. The sample clouds span an order of magnitude range in mass, size, and star formation rates. We use thermal dust emission from $Herschel$ maps to probe the gas surface density and the young stellar objects from the most recent $Spitzer$ Extended Solar Neighborhood Archive (SESNA) catalog to probe the stellar surface density. Using a star-sampled nearest neighbor technique to probe the star-gas surface density correlations at the scale of a few parsecs, we find that the stellar mass surface density varies as a power-law of the gas mass surface density, with a power-law index of $\sim$2 in all the clouds. The consistent power-law index implies that star formation efficiency is directly correlated with gas column density, and no gas column density threshold for star formation is observed. We compare the observed correlations with the predictions from an analytical model of thermal fragmentation, and with the synthetic observations of a recent hydrodynamic simulation of a turbulent star-forming molecular cloud. We find that the observed correlations are consistent for some clouds with the thermal fragmentation model and can be reproduced using the hydrodynamic simulations.

We report results of a project to map HCN and HCO+ J = 1-0 emission toward a sample of molecular clouds in the inner Galaxy, all containing dense clumps that are actively engaged in star formation. We compare these two molecular line tracers with millimeter continuum emission and extinction, as inferred from 13CO, as tracers of dense gas in molecular clouds. The fraction of the line luminosity from each tracer that comes from the dense gas, as measured by AV > 8 mag, varies substantially from cloud to cloud. In all cases, a substantial fraction (in most cases, the majority) of the total luminosity arises in gas below the AV > 8 mag threshold and outside the region of strong mm continuum emission. Measurements of the luminosity of HCN toward other galaxies will likely be dominated by such gas at lower surface density. Substantial, even dominant, contributions to the total line luminosity can arise in gas with densities typical of the cloud as a whole (densities about 100 per cubic cm). Defining the dense clump from the HCN or HCO+ emission itself, similarly to previous studies, leads to a wide range of clump properties, with some being considerably larger and less dense than in previous studies. HCN and HCO+ have similar ability to trace dense gas for the clouds in this sample. For the two clouds with low virial parameters, the 13CO is definitely a worse tracer of the dense gas, but for the other four, it is equally good (or bad) at tracing dense gas.

With the $Karl~G.~Jansky$ Very Large Array (VLA) in C-configuration, we observed a large portion of the first Galactic quadrant ($l=14.0-67.4^\circ $ and $\lvert b \rvert \leq 1.25^\circ $) achieving an angular resolution of $\leq 40^{\prime\prime}$. At $L$ Band, the WIDAR correlator at the VLA was set to cover the 21~cm HI line, four OH transitions, a series of H$n\alpha$ radio recombination lines (RRLs; $n=151$ to 186), and eight 128~MHz wide continuum spectral windows (SPWs) simultaneously. The HI emission shows clear filamentary substructures at negative velocities with low velocity crowding. The emission at positive velocities is more smeared-out likely due to higher spatial and velocity crowding of structures at the positive velocities. Comparing to the spiral arm model of the Milky Way, the atomic gas follows the Sagittarius and Perseus Arm well but with significant material in the inter-arm regions. With the C-configuration-only HI+continuum data, we produced a HI optical depth map of the THOR areal coverage from 228 absorption spectra with the nearest-neighbor method. With this $\tau$ map, we corrected the HI emission for optical depth and the derived column density is 38% higher than the column density with optically thin assumption. The total HI mass with optical depth correction in the survey region is 4.7$\times10^8~M_\odot$, 31% more than the mass derived assuming the emission is optically thin. If we apply this 31% correction to the whole Milky Way, the total atomic gas mass would be 9.4-10.5$\times 10^9~M_\odot$. Comparing the HI with existing CO data, we find a significant increase in the atomic-to-molecular gas ratio from the spiral arms to the inter-arm regions.

In this work we investigate the interplay between gravity and turbulence at different spatial scales and in different density regimes. We analyze a sample of 70 $\mu$m quiet clumps that are divided into three surface density bins and we compare the dynamics of each group with the dynamics of their respective filaments. The densest clumps form within the densest filaments on average, and they have the highest value of the velocity dispersion. The kinetic energy is transferred from the filaments down to the clumps most likely through a turbulent cascade, but we identify a critical value of the surface density, $\Sigma\simeq0.1$ g cm$^{2}$, above which the dynamics changes from being mostly turbulent-driven to mostly gravity-driven. The scenario we obtain from our data is a continuous interplay between turbulence and gravity, where the former creates structures at all scales and the latter takes the lead when the critical surface density threshold is reached. In the densest filaments this transition can occur at the parsec, or even larger scales, leading to a global collapse of the whole region and most likely to the formation of the massive objects.

Aims:We take advantage of the second data release of the Gaia space mission and the state-of-the-art astrometry delivered from very long baseline interferometry observations to revisit the structure and kinematics of the nearby Taurus star-forming region. Methods: We apply a hierarchical clustering algorithm for partitioning the stars in our sample into groups (i.e., clusters) that are associated with the various molecular clouds of the complex, and derive the distance and spatial velocity of individual stars and their corresponding molecular clouds. Results: We show that the molecular clouds are located at different distances and confirm the existence of important depth effects in this region reported in previous studies. For example, we find that the L 1495 molecular cloud is located at $d=129.9^{+0.4}_{-0.3}$ pc, while the filamentary structure connected to it (in the plane of the sky) is at $d=160.0^{+1.2}_{-1.2}$ pc. We report B 215 and L 1558 as the closest ($d=128.5^{+1.6}_{-1.6}$ pc) and most remote ($d=198.1^{+2.5}_{-2.5}$ pc) substructures of the complex, respectively. The median inter-cloud distance is 25 pc and the relative motion of the subgroups is on the order of a few km/s. We find no clear evidence for expansion (or contraction) of the Taurus complex, but signs of the potential effects of a global rotation. Finally, we compare the radial velocity of the stars with the velocity of the underlying $^{13}$CO molecular gas and report a mean difference of $0.04\pm0.12$ km/s (with r.m.s. of 0.63 km/s) confirming that the stars and the gas are tightly coupled.

We present images of radio recombination lines (RRLs) at wavelengths around 17 cm from the star-forming region W49A to determine the kinematics of ionized gas in the THOR survey (The HI/OH/Recombination line survey of the inner Milky Way) at an angular resolution of 16.8"x13.8". The distribution of ionized gas appears to be affected by feedback processes from the star clusters in W49A. The velocity structure of the RRLs shows a complex behavior with respect to the molecular gas. We find a shell-like distribution of ionized gas as traced by RRL emission surrounding the central cluster of OB stars in W49A. We describe the evolution of the shell with the recent feedback model code WARPFIELD that includes the important physical processes and has previously been applied to the 30 Doradus region in the Large Magellanic Cloud. The cloud structure and dynamics of W49A are in agreement with a feedback-driven shell that is re-collapsing. The shell may have triggered star formation in other parts of W49A. We suggest that W49A is a potential candidate for star formation regulated by feedback-driven and re-collapsing shells.

We introduce the histogram of oriented gradients (HOG), a tool developed for machine vision that we propose as a new metric for the systematic characterization of observations of atomic and molecular gas and the study of molecular cloud formation models. In essence, the HOG technique takes as input extended spectral-line observations from two tracers and provides an estimate of their spatial correlation across velocity channels. We characterize HOG using synthetic observations of HI and $^{13}$CO(J=1-0) emission from numerical simulations of MHD turbulence leading to the formation of molecular gas after the collision of two atomic clouds. We find a significant spatial correlation between the two tracers in velocity channels where $v_{HI}\approx v_{^{13}CO}$, independent of the orientation of the collision with respect to the line of sight. We use HOG to investigate the spatial correlation of the HI, from the THOR survey, and the $^{13}$CO(J=1-0) emission, from the GRS, toward the portion of the Galactic plane 33.75$\lt l\lt$35.25$^{o}$ and $|b|\lt$1.25$^{o}$. We find a significant spatial correlation between the tracers in extended portions of the studied region. Although some of the regions with high spatial correlation are associated with HI self-absorption features, suggesting that it is produced by the cold atomic gas, the correlation is not exclusive to this kind of region. The HOG results also indicate significant differences between individual regions: some show spatial correlation in channels around $v_{HI}\approx v_{^{13}CO}$ while others present this correlation in velocity channels separated by a few km/s. We associate these velocity offsets to the effect of feedback and to the presence of physical conditions that are not included in the atomic-cloud-collision simulations, such as more general magnetic field configurations, shear, and global gas infall.

We present a 1.1~mm census of dense cores in the Mon~R2 Giant Molecular Cloud with the AzTEC instrument on the Large Millimeter Telescope (LMT). We detect 295 cores (209 starless, and 86 with protostars) in a two square degree shallow survey. We also carry out a deep follow-up survey of 9 regions with low to intermediate ($3<A_V<7$) gas column densities and detect 60 new cores in the deeper survey which allows us to derive a completeness limit. After performing corrections for low signal-to-noise cores, we find a median core mass of $\sim 2.1 \text{M}_{\odot}$ and a median size of $ 0.08$~pc. $46\%$ of the cores (141) have masses exceeding the local Bonor-Ebert mass for cores with T=12K, suggesting that in the absence of supporting non-thermal pressure, these regions are unstable to gravitational collapse. We present the core mass function (CMF) for various subdivisions of the core sample. We find that cores with masses $>$10~$M_{\odot}$ are exclusively found in regions with high core number densities and that the CMF of the starless cores has an excess of low-mass cores ($<$5~$M_{\odot}$) compared to the CMF of protostellar cores. We report a power law correlation of index $1.99 \pm 0.03$ between local core mass density and gas column density (as traced by Herschel) over a wide range of size scales (0.3-5~pc). This power law is consistent with that predicted for thermal fragmentation of a self-gravitating sheet.

We present new trigonometric parallaxes and proper motions of young stellar objects in the Taurus molecular cloud complex from observations collected with the Very Long Baseline Array as part of the Gould's Belt Distances Survey (GOBELINS). We detected 26 young stellar objects and derived trigonometric parallaxes for 18 stars with an accuracy of 0.3$\%$ to a few percent. We modeled the orbits of six binaries and determined the dynamical masses of the individual components in four of these systems (V1023 Tau, T Tau S, V807 Tau and V1000 Tau). Our results are consistent with the first trigonometric parallaxes delivered by the Gaia satellite and reveal the existence of significant depth effects. We find that the central portion of the dark cloud Lynds 1495 is located at $d=129.5\pm 0.3$ pc while the B 216 clump in the filamentary structure connected to it is at $d=158.1\pm 1.2$ pc. The closest and remotest stars in our sample are located at $d=126.6\pm 1.7$ pc and $d=162.7\pm 0.8$ pc yielding a distance difference of about 36 pc. We also provide a new distance estimate for HL Tau that was recently imaged. Finally, we compute the spatial velocity of the stars with published radial velocity and investigate the kinematic properties of the various clouds and gas structures in this region.

We investigate the relation between gas and star formation in sub-galactic regions, ~360 pc to ~1.5 kpc in size, within the nearby starburst dwarf NGC4449, in order to separate the underlying relation from the effects of sampling at varying spatial scales. Dust and gas mass surface densities are derived by combining new observations at 1.1 mm, obtained with the AzTEC instrument on the Large Millimeter Telescope, with archival infrared images in the range 8-500 micron from the Spitzer Space Telescope and the Herschel Space Observatory. We extend the dynamic range of our mm (and dust) maps at the faint end, using a correlation between the far-infrared/millimeter colors F(70)/F(1100) [and F(160)/F(1100)] and the mid-infrared color F(8)/F(24) that we establish for the first time for this and other galaxies. Supplementing our data with maps of the extinction-corrected star formation rate (SFR) surface density, we measure both the SFR-molecular gas and the SFR-total gas relations in NGC4449. We find that the SFR-molecular gas relation is described by a power law with exponent that decreases from ~1.5 to ~1.2 for increasing region size, while the exponent of the SFR-total gas relation remains constant with value ~1.5 independent of region size. We attribute the molecular law behavior to the increasingly better sampling of the molecular cloud mass function at larger region sizes; conversely, the total gas law behavior likely results from the balance between the atomic and molecular gas phases achieved in regions of active star formation. Our results indicate a non-linear relation between SFR and gas surface density in NGC4449, similar to what is observed for galaxy samples.

Sensitive, imaging observations of the 1.1 mm dust continuum emission from a 1 deg^2 area collected with the AzTEC bolometer camera on the Large Millimeter Telescope are presented. A catalog of 1545 compact sources is constructed based on a Wiener-optimization filter. These sources are linked to larger clump structures identified in the Bolocam Galactic Plane Survey. Hydrogen column densities are calculated for all sources and mass and mean volume densities are derived for the subset of sources for which kinematic distances can be assigned. The AzTEC sources are localized, high density peaks within the massive clumps of molecular clouds and comprise 5-15% of the clump mass. We examine the role of the gravitational instability in generating these fragments by comparing the mass of embedded AzTEC sources to the Jeans' mass of the parent BGPS object. For sources with distances less than 6 kpc the fragment masses are comparable to the clump Jeans' mass, despite having isothermal Mach numbers between 1.6 and 7.2. AzTEC sources linked to ultra-compact HII regions have mass surface densities greater than the critical value implied by the mass-size relationship of infrared dark clouds with high mass star formation while AzTEC sources associated with Class II methanol masers have mass surface densities greater than 0.7 g cm^{-2} that approaches the proposed threshold required to form massive stars.

There is a considerable deficiency in the number of known supernova remnants (SNRs) in the Galaxy compared to that expected. Searches for extended low-surface brightness radio sources may find new Galactic SNRs, but confusion with the much larger population of HII regions makes identifying such features challenging. SNRs can, however, be separated from HII regions using their significantly lower mid-infrared (MIR) to radio continuum intensity ratios. We use the combination of high-resolution 1-2 GHz continuum data from The HI, OH, Recombination line survey of the Milky Way (THOR) and lower-resolution VLA 1.4 GHz Galactic Plane Survey (VGPS) continuum data, together with MIR data from the Spitzer GLIMPSE, Spitzer MIPSGAL, and WISE surveys to identify SNR candidates. To ensure that the candidates are not being confused with HII regions, we exclude radio continuum sources from the WISE Catalog of Galactic HII Regions, which contains all known and candidate H II regions in the Galaxy. We locate 76 new Galactic SNR candidates in the THOR and VGPS combined survey area of 67.4deg>l>17.5deg, |b|<1.25deg and measure the radio flux density for 52 previously-known SNRs. The candidate SNRs have a similar spatial distribution to the known SNRs, although we note a large number of new candidates near l=30deg, the tangent point of the Scutum spiral arm. The candidates are on average smaller in angle compared to the known regions, 6.4'+/-4.7' versus 11.0'+/-7.8', and have lower integrated flux densities. If the 76 candidates are confirmed as true SNRs, for example using radio polarization measurements or by deriving radio spectral indices, this would more than double the number of known Galactic SNRs in the survey area. This large increase would still, however, leave a discrepancy between the known and expected SNR populations of about a factor of two.

Context: The past decade has witnessed a large number of Galactic plane surveys at angular resolutions below 20". However, no comparable high-resolution survey exists at long radio wavelengths around 21cm in line and continuum emission. Methods: Employing the Very Large Array (VLA) in the C-array configuration and a large program, we observe the HI 21cm line, four OH lines, nineteen Halpha radio recombination lines as well as the continuum emission from 1 to 2GHz in full polarization over a large part of the first Galactic quadrant. Results: Covering Galactic longitudes from 14.5 to 67.4deg and latitudes between +-1.25deg, we image all of these lines and the continuum at ~20" resolution. These data allow us to study the various components of the interstellar medium (ISM): from the atomic phase, traced by the HI line, to the molecular phase, observed by the OH transitions, to the ionized medium, revealed by the cm continuum and the Halpha radio recombination lines. Furthermore, the polarized continuum emission enables magnetic field studies. In this overview paper, we discuss the survey outline and present the first data release as well as early results from the different datasets. We now release the first half of the survey; the second half will follow later after the ongoing data processing has been completed. The data in fits format (continuum images and line data cubes) can be accessed through the project web-page http://www.mpia.de/thor. Conclusions: The HI/OH/Recombination line survey of the Milky Way (THOR) opens a new window to the different parts of the ISM. It enables detailed studies of molecular cloud formation, conversion of atomic to molecular gas, and feedback from HII regions as well as the magnetic field in the Milky Way. It is highly complementary to other surveys of our Galaxy, and comparing different datasets allows us to address many open questions.

The relations between star formation and properties of molecular clouds are studied based on a sample of star forming regions in the Galactic Plane. Sources were selected by having radio recombination lines to provide identification of associated molecular clouds and dense clumps. Radio continuum and mid-infrared emission were used to determine star formation rates, while 13CO and submillimeter dust continuum emission were used to obtain masses of molecular and dense gas, respectively. We test whether total molecular gas or dense gas provides the best predictor of star formation rate. We also test two specific theoretical models, one relying on the molecular mass divided by the free-fall time, the other using the free-fall time divided by the crossing time. Neither is supported by the data. The data are also compared to those from nearby star forming regions and extragalactic data. The star formation "efficiency," defined as star formation rate divided by mass, spreads over a large range when the mass refers to molecular gas; the standard deviation of the log of the efficiency decreases by a factor of three when the mass of relatively dense molecular gas is used rather than the mass of all the molecular gas.

The origin of striations aligned along the local magnetic field direction in the translucent envelope of the Taurus molecular cloud is examined with new observations of 12CO and 13CO J=2-1 emission obtained with the 10~m submillimeter telescope of the Arizona Radio Observatory. These data identify a periodic pattern of excess blue and redshifted emission that is responsible for the striations. For both 12CO and 13CO, spatial variations of the J=2-1 to J=1-0 line ratio are small and are not spatially correlated with the striation locations. A medium comprised of unresolved CO emitting substructures (cells) with a beam area filling factor less than unity at any velocity is required to explain the average line ratios and brightness temperatures. We propose that the striations result from the modulation of velocities and the beam filling factor of the cells as a result of either the Kelvin-Helmholtz instability or magnetosonic waves propagating through the envelope of the Taurus molecular cloud. Both processes are likely common features in molecular clouds that are sub-Alfvenic and may explain low column density, cirrus-like features similarly aligned with the magnetic field observed throughout the interstellar medium in far-infrared surveys of dust emission.

We analyze radial and azimuthal variations of the phase balance between the molecular and atomic ISM in the Milky Way. In particular, the azimuthal variations -- between spiral arm and interarm regions -- are analyzed without any explicit definition of spiral arm locations. We show that the molecular gas mass fraction, i.e., fmol=H2/ (HI+H2) in mass, varies predominantly in the radial direction: starting from ~100% at the center, remaining ~>50% (~>60%) to R~6kpc, and decreasing to ~10-20% (~50%) at R=8.5 kpc when averaged over the whole disk thickness (in the mid plane). Azimuthal, arm-interarm variations are secondary: only ~20%, in the globally molecule-dominated inner MW, but becoming larger, ~40-50%, in the atom-dominated outskirts. This suggests that in the inner MW, the gas stays highly molecular (fmol>50%) as it goes from an interarm region, into a spiral arm, and back into the next interarm region. Stellar feedback does not dissociate molecules much, and the coagulation and fragmentation of molecular clouds dominate the evolution of the ISM at these radii. The trend differs in the outskirts, where the gas phase is globally atomic (fmol<50%). The HI and H2 phases cycle through spiral arm passage there. These different regimes of ISM evolution are also seen in external galaxies (e.g., LMC, M33, and M51). We explain the radial gradient of fmol by a simple flow continuity model. The effects of spiral arms on this analysis are illustrated in Appendix.

We have carried out a search for filamentary structures in the Taurus molecular cloud using $\rm^{13}CO$ line emission data from the FCRAO survey of $\rm \sim100 \, deg^2$. We have used the topological analysis tool, DisPerSe, and post-processed its results to include a more strict definition of filaments that requires an aspect ratio of at least 3:1 and cross section intensity profiles peaked on the spine of the filament. In the velocity-integrated intensity map only 10 of the hundreds of filamentary structures identified by DisPerSe comply with our criteria. Unlike Herschel analyses, which find a characteristic width for filaments of $\rm \sim0.1 \, pc$, we find a much broader distribution of profile widths in our structures, with a peak at 0.4 pc. Furthermore, even if the identified filaments are cylindrical objects, their complicated velocity structure and velocity dispersions imply that they are probably gravitationally unbound. Analysis of velocity channel maps reveals the existence of hundreds of `velocity-coherent' filaments. The distribution of their widths is peaked at lower values (0.2 pc) while the fluctuation of their peak intensities is indicative of stochastic origin. These filaments are suppressed in the integrated intensity map due to the blending of diffuse emission from different velocities. Conversely, integration over velocities can cause filamentary structures to appear. Such apparent filaments can also be traced, using the same methodology, in simple simulated maps consisting of randomly placed cores. They have profile shapes similar to observed filaments and contain most of the simulated cores.

Newborn stars form within the localized, high density regions of molecular clouds. The sequence and rate at which stars form in dense clumps and the dependence on local and global environments are key factors in developing descriptions of stellar production in galaxies. We seek to observationally constrain the rate and latency of star formation in dense massive clumps that are distributed throughout the Galaxy and to compare these results to proposed prescriptions for stellar production. A sample of 24 micron-based Class~I protostars are linked to dust clumps that are embedded within molecular clouds selected from the APEX Telescope Large Area Survey of the Galaxy. We determine the fraction of star-forming clumps, f*, that imposes a constraint on the latency of star formation in units of a clump's lifetime. Protostellar masses are estimated from models of circumstellar environments of young stellar objects from which star formation rates are derived. Physical properties of the clumps are calculated from 870 micron dust continuum emission and NH_3 line emission. Linear correlations are identified between the star formation rate surface density, Sigma_{SFR}, and the quantities Sigma_{H2}/tau_{ff} and Sigma_{H2}/tau_{cross}, suggesting that star formation is regulated at the local scales of molecular clouds. The measured fraction of star forming clumps is 23%. Accounting for star formation within clumps that are excluded from our sample due to 24 micron saturation, this fraction can be as high as 31%. Dense, massive clumps form primarily low mass (< 1-2 msun) stars with emergent 24 micron fluxes below our sensitivity limit or are incapable of forming any stars for the initial 70% of their lifetimes. The low fraction of star forming clumps in the Galactic center relative to those located in the disk of the Milky Way is verified.

We carried out a large program with the Karl G. Jansky Very Large Array (VLA): "THOR: The HI, OH, Recombination line survey of the Milky Way". We observed a significant portion of the Galactic plane in the first quadrant of the Milky Way in the 21cm HI line, 4 OH transitions, 19 radio recombination lines, and continuum from 1 to 2 GHz. In this paper we present a catalog of the continuum sources in the first half of the survey (l=14.0-37.9deg and l=47.1-51.2deg, |b|<1.1deg) at a spatial resolution of 10-25", with a spatially varying noise level of ~0.3-1 mJy/beam. The catalog contains ~4400 sources. Around 1200 of these are spatially resolved, and ~1000 are possible artifacts, given their low signal-to-noise ratios. Since the spatial distribution of the unresolved objects is evenly distributed and not confined to the Galactic plane, most of them are extragalactic. Thanks to the broad bandwidth of the observations from 1 to 2 GHz, we are able to determine a reliable spectral index for ~1800 sources. The spectral index distribution reveals a double-peaked profile with maxima at spectral indices of alpha = -1 and alpha = 0 , corresponding to steep declining and flat spectra, respectively. This allows us to distinguish between thermal and non-thermal emission, which can be used to determine the nature of each source. We examine the spectral index of ~300 known HII regions, for which we find thermal emission with spectral indices around alpha = 0. In contrast, supernova remnants (SNR) show non-thermal emission with alpha = -0.5 and extragalactic objects generally have a steeper spectral index of alpha = -1. Using the spectral index information of the THOR survey, we investigate potential SNR candidates. We classify the radiation of four SNR candidates as non-thermal, and for the first time, we provide strong evidence for the SNR origin of these candidates.

Emission from carbon monoxide (CO) is ubiquitously used as a tracer of dense star forming molecular clouds. There is, however, growing evidence that a significant fraction of CO emission originates from diffuse molecular gas. Quantifying the contribution of diffuse CO-emitting gas is vital for understanding the relation between molecular gas and star formation. We examine the Galactic distribution of two CO-emitting gas components, a high column density component detected in 13CO and 12CO, and a low column density component detected in 12CO, but not in 13CO. The "diffuse" and "dense" components are identified using a combination of smoothing, masking, and erosion/dilation procedures, making use of three large-scale 12CO and 13CO surveys of the Inner and Outer Milky Way. The diffuse component, which globally represents 25 (1.5x1e8 Mo) of the total molecular gas mass (6.5x1e8 Mo), is more extended perpendicular to the Galactic plane. The fraction of diffuse gas increases from 15% at a galactocentric radius of 3 kpc to 50% at 15 kpc, and increases with decreasing surface density. In the Inner Galaxy, a yet denser component traced by CS emission represents 14% of the total molecular gas mass traced by 12CO emission. Only 14% of the molecular gas mass traced by 12CO emission is identified as part of molecular clouds in 12CO surveys by cloud identification algorithms. This study indicates that CO emission not only traces star forming clouds, but also a significant diffuse molecular ISM component.

We present a detailed analysis of narrow of NaI and KI absorption resonance lines toward nearly 40 T Tauri stars in Taurus with the goal of clarifying their origin. The NaI 5889.95 angstrom line is detected toward all but one source, while the weaker KI 7698.96 angstrom line in about two thirds of the sample. The similarity in their peak centroids and the significant positive correlation between their equivalent widths demonstrate that these transitions trace the same atomic gas. The absorption lines are present towards both disk and diskless young stellar objects, which excludes cold gas within the circumstellar disk as the absorbing material. A comparison of NaI and CO detections and peak centroids demonstrates that the atomic and molecular gas are not co-located, the atomic gas is more extended than the molecular gas. The width of the atomic lines corroborates this finding and points to atomic gas about an order of magnitude warmer than the molecular gas. The distribution of NaI radial velocities shows a clear spatial gradient along the length of the Taurus molecular cloud filaments. This suggests that absorption is associated with the Taurus molecular cloud. Assuming the gradient is due to cloud rotation, the rotation of the atomic gas is consistent with differential galactic rotation while the rotation of the molecular gas, although with the same rotation axis, is retrograde. Our analysis shows that narrow NaI and KI absorption resonance lines are useful tracers of the atomic envelope of molecular clouds. In line with recent findings from giant molecular clouds, our results demonstrate that the velocity fields of the atomic and molecular gas are misaligned. The angular momentum of a molecular cloud is not simply inherited from the rotating Galactic disk from which it formed but may be redistributed by cloud-cloud interactions.

Measuring redshifted CO line emission is an unambiguous method for obtaining an accurate redshift and total cold gas content of optically faint, dusty starburst systems. Here, we report the first successful spectroscopic redshift determination of AzTEC J095942.9+022938 ("COSMOS AzTEC-1"), the brightest 1.1mm continuum source found in the AzTEC/JCMT survey (Scott et al. 2008), through a clear detection of the redshifted CO (4-3) and CO (5-4) lines using the Redshift Search Receiver on the Large Millimeter Telescope. The CO redshift of $z=4.3420\pm0.0004$ is confirmed by the detection of the redshifted 158 micron [C II] line using the Submillimeter Array. The new redshift and Herschel photometry yield $L_{FIR}=(1.1\pm0.1)\times 10^{13} L_\odot$ and $SFR = 1300\, M_\odot$ yr$^{-1}$. Its molecular gas mass derived using the ULIRG conversion factor is $1.4\pm0.2 \times 10^{11} M_\odot$ while the total ISM mass derived from the 1.1mm dust continuum is $3.7\pm0.7 \times 10^{11} M_\odot$ assuming dust temperature of 35 K. Our dynamical mass analysis suggests that the compact gas disk ($r\approx 1.1$ kpc, inferred from dust continuum and SED analysis) has to be nearly face-on, providing a natural explanation for the uncommonly bright, compact stellar light seen by the HST. The [C II] line luminosity $L_{[C~II]} = 7.8\pm1.1 \times 10^9 L_\odot$ is remarkably high, but it is only 0.04 per cent of the total IR luminosity. AzTEC COSMOS-1 and other high redshift sources with a spatially resolved size extend the tight trend seen between [C II]/FIR ratio and $\Sigma_{FIR}$ among IR-bright galaxies reported by Diaz-Santos et al. (2013) by more than an order of magnitude, supporting the explanation that the higher intensity of the IR radiation field is responsible for the "[C II] deficiency" seen among luminous starburst galaxies.

To study the atomic, molecular and ionized emission of Giant Molecular Clouds (GMCs), we have initiated a Large Program with the VLA: 'THOR - The HI, OH, Recombination Line survey of the Milky Way'. We map the 21cm HI line, 4 OH lines, 19 H_alpha recombination lines and the continuum from 1 to 2 GHz of a significant fraction of the Milky Way (l=15-67deg, |b|<1deg) at ~20" resolution. In this paper, we focus on the HI emission from the W43 star-formation complex. Classically, the HI 21cm line is treated as optically thin with column densities calculated under this assumption. This might give reasonable results for regions of low-mass star-formation, however, it is not sufficient to describe GMCs. We analyzed strong continuum sources to measure the optical depth, and thus correct the HI 21cm emission for optical depth effects and weak diffuse continuum emission. Hence, we are able to measure the HI mass of W43 more accurately and our analysis reveals a lower limit of M=6.6x10^6 M_sun, which is a factor of 2.4 larger than the mass estimated with the assumption of optically thin emission. The HI column densities are as high as N(HI)~150 M_sun/pc^2 ~ 1.9x10^22 cm^-2, which is an order of magnitude higher than for low mass star formation regions. This result challenges theoretical models that predict a threshold for the HI column density of ~10 M_sun/pc^2, at which the formation of molecular hydrogen should set in. By assuming an elliptical layered structure for W43, we estimate the particle density profiles. While at the cloud edge atomic and molecular hydrogen are well mixed, the center of the cloud is strongly dominated by molecular hydrogen. We do not identify a sharp transition between hydrogen in atomic and molecular form. Our results are an important characterization of the atomic to molecular hydrogen transition in an extreme environment and challenges current theoretical models.

In this contribution, we describe the applied methods to construct a 24 micron based point source catalog derived from the image data of the MIPSGAL 24 micron Galactic Plane Survey and the corresponding data products. The high quality catalog product contains 933,818 sources, with a total of 1,353,228 in the full archive catalog. The source tables include positional and photometric informa- tion derived from the 24 micron images, source quality and confusion flags and coun- terpart photometry from matched 2MASS, GLIMPSE, and WISE point sources. Completeness decay data cubes are constructed at 1 arcminute angular resolution that describe the varying background levels over the MIPSGAL field and the ability to extract sources of a given magnitude from this background. The completeness decay cubes are included in the set of data products. We present the results of our efforts to verify the astrometric and photometric calibration of the catalog, and present several analyses of minor anomalies in these measurements to justify adopted mitigation strategies.

Giant molecular clouds (GMCs) are the primary reservoirs of cold, star-forming molecular gas in the Milky Way and similar galaxies, and thus any understanding of star formation must encompass a model for GMC formation, evolution, and destruction. These models are necessarily constrained by measurements of interstellar molecular and atomic gas, and the emergent, newborn stars. Both observations and theory have undergone great advances in recent years, the latter driven largely by improved numerical simulations, and the former by the advent of large-scale surveys with new telescopes and instruments. This chapter offers a thorough review of the current state of the field.

The mass fraction of dense gas within giant molecular clouds (GMCs) of the Milky Way is investigated using 13CO data from the FCRAO Galactic Plane Surveys and the Bolocam Galactic Plane Survey (BGPS) of 1.1 mm dust continuum emission. A sample of 860 compact dust sources are selected from the BGPS catalog and kinematically linked to 344 clouds of extended (>3') 13CO J=1-0 emission. Gas masses are tabulated for the full dust source and subregions within the dust sources with mass surface densities greater than 200 M_solar/pc^2, which are assumed to be regions of enhanced volume density. Masses of the parent GMCs are calculated assuming optically thin 13CO J=1-0 emission and LTE conditions. The mean fractional mass of dust sources to host GMC mass is 0.11 (-0.06, +0.12). The high column density subregions comprise 0.07 (-0.05, +0.13) of the mass of the cloud. Owing to our assumptions, these values are upper limits to the true mass fractions. The fractional mass of dense gas is independent of GMC mass and gas surface density. The low dense gas mass fraction suggests that the formation of dense structures within GMCs is the primary bottleneck for star formation. The distribution of velocity differences between the dense gas and the low density material along the line of sight is also examined. We find a strong, centrally peaked distribution centered on zero velocity displacement. This distribution of velocity differences is modeled with radially converging flows towards the dense gas position that are randomly oriented with respect to the observed line of sight. These models constrain the infall velocities to be 2-4 km/s for various flow configurations.

Principal component analysis is a powerful statistical system to investigate the structure and dynamics of the molecular interstellar medium, with particular emphasis on the study of turbulence, as revealed by spectroscopic imaging of molecular line emission. To-date, the method to retrieve the power law index of the velocity structure function or power spectrum has relied on an empirical calibration and testing with model turbulent velocity fields, while lacking a firm theoretical basis. In this paper, we present an analytic formulation that reveals the detailed mechanics of the method and confirms previous empirical calibrations of its recovery of the scale dependence of turbulent velocity fluctuations.

Infrared stellar photometry from 2MASS and spectral line imaging observations of 12CO and 13CO J = 1-0 line emission from the FCRAO 14m telescope are analysed to assess the variation of the CO abundance with physical conditions throughout the Orion A and Orion B molecular clouds. Three distinct Av regimes are identified in which the ratio between the 13CO column density and visual extinction changes corresponding to the photon dominated envelope, the strongly self-shielded interior, and the cold, dense volumes of the clouds. Within the strongly self-shielded interior of the Orion A cloud, the 13CO abundance varies by 100% with a peak value located near regions of enhanced star formation activity. The effect of CO depletion onto the ice mantles of dust grains is limited to regions with AV > 10 mag and gas temperatures less than 20 K as predicted by chemical models that consider thermal-evaporation to desorb molecules from grain surfaces. Values of the molecular mass of each cloud are independently derived from the distributions of Av and 13CO column densities with a constant 13CO-to-H2 abundance over various extinction ranges. Within the strongly self-shielded interior of the cloud (Av > 3 mag), 13CO provides a reliable tracer of H2 mass with the exception of the cold, dense volumes where depletion is important. However, owing to its reduced abundance, 13CO does not trace the H2 mass that resides in the extended cloud envelope, which comprises 40-50% of the molecular mass of each cloud. The implied CO luminosity to mass ratios, M/L_{CO}, are 3.2 and 2.9 for Orion A and Orion B respectively, which are comparable to the value (2.9), derived from gamma-ray observations of the Orion region. Our results emphasize the need to consider local conditions when applying CO observations to derive H2 column densities.

The Galactic supernova remnant (SNR) IC443 is one of the most studied core-collapse SNRs for its interaction with molecular clouds. However, the ambient molecular clouds with which IC443 is interacting have not been thoroughly studied and remain poorly understood. Using Five College Radio Astronomy Observatory 14m telescope, we obtained fully sampled maps of ~ 1{\deg} \times 1{\deg} region toward IC443 in the 12CO J=1-0 and HCO+ J=1-0 lines. In addition to the previously known molecular clouds in the velocity range v_lsr = -6 to -1 km/s (-3 km/s clouds), our observations reveal two new ambient molecular cloud components: small (~ 1') bright clouds in v_lsr = -8 to -3 km/s (SCs), and diffuse clouds in v_lsr = +3 to +10 km/s (+5 km/s clouds). Our data also reveal the detailed kinematics of the shocked molecular gas in IC443, however the focus of this paper is the physical relationship between the shocked clumps and the ambient cloud components. We find strong evidence that the SCs are associated with the shocked clumps. This is supported by the positional coincidence of the SCs with shocked clumps and other tracers of shocks. Furthermore, the kinematic features of some shocked clumps suggest that these are the ablated material from the SCs upon the impact of the SNR shock. The SCs are interpreted as dense cores of parental molecular clouds that survived the destruction by the pre-supernova evolution of the progenitor star or its nearby stars. We propose that the expanding SNR shock is now impacting some of the remaining cores and the gas is being ablated and accelerated producing the shocked molecular gas. The morphology of the +5 km/s clouds suggests an association with IC443. On the other hand, the -3 km/s clouds show no evidence for interaction.

Magnetically aligned velocity anisotropy over varying physical conditions and environments within the Taurus Molecular Cloud is evaluated from analysis of wide field spectroscopic imaging of 12CO and 13CO J=1-0 emission. Such anisotropy is a result of MHD turbulence in the strong magnetic field regime and provides an indirect measure of the role of magnetic fields upon the gas. Velocity anisotropy aligned with the local, projected mean magnetic field direction is limited to fields with low surface brightness 12CO emission corresponding to regions of low visual extinction and presumably, low gas volume density. The more optically thin 13CO J=1-0 emission shows little evidence for velocity anisotropy. We compare our results with computational simulations with varying degrees of magnetic field strength and Alfvenic Mach numbers. In the diffuse, molecular envelope of the cloud, a strong magnetic field and sub-Alfvenic turbulent motions are inferred. Super-Alfvenic motions are present within the high column density filaments of the Taurus cloud. From this trans-Alfvenic flow, we constrain the scaling exponent, kappa, of the magnetic field density relation (B ~ n^kappa) to be near zero as expected for ambipolar diffusion or material loading of magnetic flux tubes.

Turbulence plays a major role in the formation and evolution of molecular clouds. The problem is that turbulent velocities are convolved with the density of an observed region. To correct for this convolution, we investigate the relation between the turbulence spectrum of model clouds, and the statistics of their synthetic observations obtained from Principal Component Analysis (PCA). We apply PCA to spectral maps generated from simulated density and velocity fields, obtained from hydrodynamic simulations of supersonic turbulence, and from fractional Brownian motion fields with varying velocity, density spectra, and density dispersion. We examine the dependence of the slope of the PCA structure function, alpha_PCA, on intermittency, on the turbulence velocity (beta_v) and density (beta_n) spectral indexes, and on density dispersion. We find that PCA is insensitive to beta_n and to the log-density dispersion sigma_s, provided sigma_s < 2. For sigma_s>2, alpha_PCA increases with sigma_s due to the intermittent sampling of the velocity field by the density field. The PCA calibration also depends on intermittency. We derive a PCA calibration based on fBms with sigma_s<2 and apply it to 367 CO spectral maps of molecular clouds in the Galactic Ring Survey. The average slope of the PCA structure function, <alpha_PCA>=0.62\pm0.2, is consistent with the hydrodynamic simulations and leads to a turbulence velocity exponent <beta_v>=2.06\pm0.6 for a non-intermittent, low density dispersion flow. Accounting for intermittency and density dispersion, the coincidence between the PCA slope of the GRS clouds and the hydrodynamic simulations suggests beta_v~1.9, consistent with both Burgers and compressible intermittent turbulence.

The feedback from galactic outflows is thought to play an important role in shaping the gas content, star formation history, and ultimately the stellar mass function of galaxies. Here we present evidence for massive molecular outflows associated with ultra-luminous infrared galaxies (ULIRGs) in the coadded Redshift Search Receiver 12CO(1-0) spectrum. Our stacked spectrum of 27 ULIRGs at z = 0.043-0.11 (freq_rest = 110-120 GHz) shows broad wings around the CO line with delta_V(FWZI)~2000 km/s. Its integrated line flux accounts for up to 25+/-5% of the total CO line luminosity. When interpreted as a massive molecular outflow wind, the associated mechanical energy can be explained by a concentrated starburst with SFR \geq 100 M_sun/yr, which agrees well with their SFR derived from the FIR luminosity. Using the high signal-to-noise stacked composite spectrum, we also probe 13CO and 12CN emission in the sample and discuss how the chemical abundance of molecular gas may vary depending on the physical conditions of the nuclear region.

The 3 mm wavelength spectra of 10 galaxies have been obtained at the Five College Radio Astronomy Observatory using a new, very broadband receiver and spectrometer, called the Redshift Search Receiver (RSR). The RSR has an instantaneous bandwidth of 37 GHz covering frequencies from 74 to 111 GHz, and has a spectral resolution of 31 MHz (~100 km/s). During tests of the RSR on the FCRAO 14 m telescope the complete 3 mm spectra of the central regions of NGC 253, Maffei 2, NGC1068, IC 342, M82, NGC 3079, NGC 3690, NGC 4258, Arp 220 and NGC 6240 were obtained. Within the wavelength band covered by the RSR, 20 spectral lines from 14 different atomic and molecular species were detected. Based on simultaneous fits to the spectrum of each galaxy, a number of key molecular line ratios are derived. A simple model which assumes the emission arises from an ensemble of Milky Way-like Giant Nolecular Cloud cores can adequately fit the observed line ratios using molecular abundances based on Galactic molecular cloud cores. Variations seen in some line ratios, such as 13CO/HCN and HCO+/HCN, can be explained if the mean density of the molecular gas varies from galaxy to galaxy. However, NGC 3690, NGC 4258 and NGC 6240 show very large HCO+/HCN ratios and require significant abundance enhancement of HCO+ over HCN, possible due to the proximity to active galactic nucleus activity. Finally, the mass of dense molecular gas is estimated and we infer that 25-85 % of the total molecular gas in the central regions of these galaxies must have densities greater than 10^4 cm^-3.

We derive the physical properties of 580 molecular clouds based on their 12CO and 13CO line emission detected in the University of Massachusetts-Stony Brook (UMSB) and Galactic Ring surveys. We provide a range of values of the physical properties of molecular clouds, and find a power-law correlation between their radii and masses, suggesting that the fractal dimension of the ISM is around 2.36. This relation, M = (228 +/- 18) R^{2.36+/-0.04}, allows us to derive masses for an additional 170 GRS molecular clouds not covered by the UMSB survey. We derive the Galactic surface mass density of molecular gas and examine its spatial variations throughout the Galaxy. We find that the azimuthally averaged Galactic surface density of molecular gas peaks between Galactocentric radii of 4 and 5 kpc. Although the Perseus arm is not detected in molecular gas, the Galactic surface density of molecular gas is enhanced along the positions of the Scutum-Crux and Sagittarius arms. This may indicate that molecular clouds form in spiral arms and are disrupted in the inter-arm space. Last, we find that the CO excitation temperature of molecular clouds decreases away from the Galactic center, suggesting a possible decline in the star formation rate with Galactocentric radius. There is a marginally significant enhancement in the CO excitation temperature of molecular clouds at a Galactocentric radius of about 6 kpc, which in the longitude range of the GRS corresponds to the Sagittarius arm. This temperature increase could be associated with massive star formation in the Sagittarius spiral arm.

We investigate the relation between star formation rate (SFR) and gas surface densities in Galactic star forming regions using a sample of YSOs and massive clumps. Our YSO sample consists of objects located in 20 molecular clouds from the Spitzer cores to disks and Gould's Belt surveys. We estimate the gas surface density (Sigma_gas) from Av maps and YSO SFR surface densities (Sigma_SFR) from the number of YSOs, assuming a mean mass and lifetime. We also divide the clouds into contour levels of Av, counting only the youngest Class I and Flat SED YSOs. For a sample of massive star forming clumps, we derive SFRs from the infrared luminosity and use HCN gas maps to estimate Sigma_gas. We find that Galactic clouds lie above the extragalactic relations (e.g., Kennicutt-Schmidt Law) by factors up to 17. Cloud regions with high Sigma_gas lie above extragalactic relations up to a factor of 54 and overlap with massive clumps. We use 12CO and 13CO gas maps of the Perseus and Ophiuchus clouds to estimate Sigma_gas and compare to Sigma_gas from Av maps. We find that 13CO, underestimates the Av-based mass by factors of 4-5. 12CO may underestimate the total gas mass at Sigma_gas > 200 Msun pc^-2 by > 30%;however, this does not explain the large discrepancy between Galactic and extragalactic relations. We find evidence for a threshold of star formation (Sigma_th) at 129+-14 Msun pc^-2. At Sigma_gas > Sigma_th, the Galactic relation is linear. A possible reason for the difference between Galactic and extragalactic relations is that all the CO-emitting gas, including Sigma_gas below Sigma_th, is measured in extragalactic studies. If the Kennicutt-Schmidt relation (Sigma_SFR Sigma_gas^1.4) and a linear relation between dense gas and star formation is assumed, the fraction of dense star forming gas (f_dense) increases as Sigma_gas^0.4. When Sigma_gas reaches ~300Sigma_th, f_dense is 1. (Abridged)

We have performed an On-The-Fly (OTF) mapping survey of ${\rm ^{12}{CO(J=1-0)}}$ emission in 28 Virgo cluster spiral galaxies using the Five College Radio Astronomy Observatory (FCRAO) 14-m telescope. This survey aims to characterize the CO distribution, kinematics, and luminosity of a large sample of galaxies covering the full extents of stellar disks, rather than sampling only the inner disks or the major axis as was done by many previous single dish and interferometric CO surveys. CO emission is detected in 20 galaxies among the 28 Virgo spirals observed. An atlas consisting of global measures, radial measures, and maps, is presented for each detected galaxy. A note summarizing the CO data is also presented along with relevant information from the literature. The CO properties derived from our OTF observations are presented and compared with the results from the FCRAO Extragalactic CO Survey by Young et al. (1995) which utilized position-switching observations along the major axis and a model fitting method. We find that our OTF derived CO properties agree well with the Young et al. results in many cases, but the Young et al. measurements are larger by a factor of 1.4 - 2.4 for seven (out of 18) cases. We will explore further the possible causes for the discrepancy in the analysis paper currently under preparation.

Supersonic turbulence in molecular clouds is a dominant agent that strongly affects the clouds' evolution and star formation activity. Turbulence may be initiated and maintained by a number of processes, acting at a wide range of physical scales. By examining the dynamical state of molecular clouds, it is possible to assess the primary candidates for how the turbulent energy is injected. The aim of this paper is to constrain the scales at which turbulence is driven in the molecular interstellar medium, by comparing simulated molecular spectral line observations of numerical magnetohydrodynamic (MHD) models and molecular spectral line observations of real molecular clouds. We use principal component analysis, applied to both models and observational data, to extract a quantitative measure of the driving scale of turbulence. We find that only models driven at large scales (comparable to, or exceeding, the size of the cloud) are consistent with observations. This result applies also to clouds with little or no internal star formation activity. Astrophysical processes acting on large scales, including supernova-driven turbulence, magnetorotational instability, or spiral shock forcing, are viable candidates for the generation and maintenance of molecular cloud turbulence. Small scale driving by sources internal to molecular clouds, such as outflows, can be important on small scales, but cannot replicate the observed large-scale velocity fluctuations in the molecular interstellar medium.

We present 12CO J=1-->0 observations of ultraluminous infrared galaxies (ULIRGs) obtained using the Redshift Search Receiver (RSR) on the 14-m telescope of the Five College Radio Astronomy Observatory. The RSR is a novel, dual-beam, dual-polarization receiver equipped with an ultra-wideband spectrometer backend that is being built as a facility receiver for the Large Millimeter Telescope. Our sample consists of 29 ULIRGs in the redshift range of 0.04-0.11, including 10 objects with no prior 12CO measurements. We have detected 27 systems (a detection rate of 93%), including 9 ULIRGs that are detected in CO for the first time. Our study has increased the number of local ULIRGs with CO measurements by ~15%. The CO line luminosity L'_CO, correlates well with far-infrared luminosity L_FIR, following the general trend of other local ULIRGs. However, compared to previous surveys we probe deeper into the low CO luminosity end of ULIRG population as a single study by including a number of CO faint objects in the sample. As a result, we find 1) a smoother transition between the ULIRG population and local QSOs in L_FIR-L'_CO ("star formation efficiency") space, and 2) a broader range of L_FIR/L'_CO flux ratio (~60--1000 L_sun/[K km/s/pc^2]) than previously reported. In our new survey, we also have found a small number of ULIRGs with extreme L_FIR/L'_CO which had been known to be rare. The mid-IR color and radio-excess of 56 local ULIRGs as a function of FIR-to-CO flux ratio is examined and compared with those of spirals/starburst galaxies and low-z QSOs. In this paper, using a large sample of local ULIRGs we explore the origin of their current power source and potential evolution to QSOs.

Supernova explosions inject a considerable amount of energy into the interstellar medium (ISM) in regions with high to moderate star formation rates. In order to assess whether the driving of turbulence by supernovae is also important in the outer Galactic disk, where the star formation rates are lower, we study the spatial distribution of molecular cloud (MC) inclinations with respect to the Galactic plane. The latter contains important information on the nature of the mechanism of energy injection into the ISM. We analyze the spatial correlations between the position angles (PAs) of a selected sample of MCs (the largest clouds in the catalogue of the outer Galaxy published by Heyer et al. 2001). Our results show that when the PAs of the clouds are all mapped to values into the [0,90]degrees interval, there is a significant degree of spatial correlation between the $PA$s on spatial scales in the range of 100-800 pc. These scales are of the order of the sizes of individual SN shells in low density environments such as those prevailing in the outer Galaxy and where the metallicity of the ambient gas is of the order of the solar value or smaller. These findings suggest that individual SN explosions, occurring in the outer regions of the Galaxy and in likewise spiral galaxies, albeit at lower rates, continue to play an important role in shaping the structure and dynamics of the ISM in those regions. The SN explosions we postulate here are likely associated with the existence of young stellar clusters in the far outer regions of the Galaxy and the UV emission and low levels of star formation observed with the GALEX satellite in the outer regions of local galaxies.

Kinematic distances to 750 molecular clouds identified in the 13CO J=1-0 Boston University-Five College Radio Astronomy Observatory Galactic Ring Survey (BU-FCRAO GRS) are derived assuming the Clemens rotation curve of the Galaxy. The kinematic distance ambiguity is resolved by examining the presence of HI self-absorption toward the 13CO emission peak of each cloud using the Very Large Array Galactic Plane Survey (VGPS). We also identify 21 cm continuum sources embedded in the GRS clouds in order to use absorption features in the HI 21 cm continuum to distinguish between near and far kinematic distances. The Galactic distribution of GRS clouds is consistent with a four-arm model of the Milky Way. The locations of the Scutum-Crux and Perseus arms traced by GRS clouds match star count data from the Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE) star-count data. We conclude that molecular clouds must form in spiral arms and be short-lived (lifetimes < 10 Myr) in order to explain the absence of massive, 13CO bright molecular clouds in the inter-arm space.

The properties of Galactic molecular clouds tabulated by Solomon etal (1987) (SRBY) are re-examined using the Boston University-FCRAO Galactic Ring Survey of 13CO J=1-0 emission. These new data provide a lower opacity tracer of molecular clouds and improved angular and spectral resolution than previous surveys of molecular line emission along the Galactic Plane. We calculate GMC masses within the SRBY cloud boundaries assuming LTE conditions throughout the cloud and a constant H2 to 13CO abundance, while accounting for the variation of the 12C/13C with Galacto-centric radius. The LTE derived masses are typically five times smaller than the SRBY virial masses. The corresponding median mass surface density of molecular hydrogen for this sample is 42 Msun/pc^2, which is significantly lower than the value derived by SRBY (median 206 Msun/pc^2) that has been widely adopted by most models of cloud evolution and star formation. This discrepancy arises from both the extrapolation by SRBY of velocity dispersion, size, and CO luminosity to the 1K antenna temperature isophote that likely overestimates the GMC masses and our assumption of constant 13CO abundance over the projected area of each cloud. Owing to the uncertainty of molecular abundances in the envelopes of clouds, the mass surface density of giant molecular clouds could be larger than the values derived from our 13CO measurements. From velocity dispersions derived from the 13CO data, we find that the coefficient of the cloud structure functions, vo=sigma_v/R^{1/2}, is not constant, as required to satisfy Larson's scaling relationships, but rather systematically varies with the surface density of the cloud as Sigma^{0.5} as expected for clouds in self-gravitational equlibrium.

We present a new large-scale survey of J=3-2 12CO emission covering 4.8 square degrees around the Rosette Nebula. Approximately 2000 compact clumps are identified, with a spatially-invariant power law mass distribution index of -1.8. Most of the inner clumps show velocity gradients of 1-3 km/s/pc, directed away from the exciting nebula. The gradients decrease with distance from the central O stars, and are consistent with a photoionised gas acceleration model, assuming clump lifetimes of a few 10^5 yrs. However, in one clear case, the observed near-constant velocity gradient is difficult to explain with simple models. Most blue-shifted but very few of the red-shifted clumps are associated with dark absorbing optical globules, confirming that the dominant molecular gas motion is expansion away from the central nebula and O stars. Many clumps also lie in a molecular ring, having an expansion velocity of 30 km/s, radius 11pc, and dynamical lifetime of ~1Myr. The J=3-2/1-0 12CO line ratios of the clumps decrease with distance from the O stars, implying a gradient in their surface temperatures; the results are consistent with a simple model of clump surface heating due to the central stars. Seven high-velocity molecular flows are found in the region, with a close correspondence between these flows and embedded young clusters. These outflows are sufficiently energetic to drive local gas turbulence within each cluster. We find 14 clear examples of association between embedded young stars seen at 24um and CO clumps; these are thought to be photoevaporating molecular envelopes. The CO clumps without evidence of embedded stars tend to have lower velocity gradients, and it is suggested that the presence of the young star may extend the lifespan of the photoevaporating envelope.

Three-millimeter-wavelength spectra of a number of nearby galaxies have been obtained at the Five College Radio Astronomy Observatory (FCRAO) using a new, very broadband receiver. This instrument, which we call the Redshift Search Receiver, has an instantaneous bandwidth of 36 GHz and operates from 74 to 110.5 GHz. The receiver has been built at UMass/FCRAO to be part of the initial instrumentation for the Large Millimeter Telescope (LMT) and is intended primarily for determination of the redshift of distant, dust-obscured galaxies. It is being tested on the FCRAO 14m by measuring the 3mm spectra of a number of nearby galaxies. There are interesting differences in the chemistry of these galaxies.

We derive the molecular properties for a sample of 301 Galactic HII regions including 123 ultra compact (UC), 105 compact, and 73 diffuse nebulae. We analyze all sources within the BU-FCRAO Galactic Ring Survey (GRS) of 13CO emission known to be HII regions based upon the presence of radio continuum and cm-wavelength radio recombination line emission. Unlike all previous large area coverage 13CO surveys, the GRS is fully sampled in angle and yet covers ~75 square degrees of the Inner Galaxy. The angular resolution of the GRS 46" allows us to associate molecular gas with HII regions without ambiguity and to investigate the physical properties of this molecular gas. We find clear CO/HII morphological associations in position and velocity for ~80% of the nebular sample. Compact HII region molecular gas clouds are on average larger than UC clouds: 2.2' compared to 1.7'. Compact and UC HII regions have very similar molecular properties, with ~5K line intensities and ~4 km/s line widths. The diffuse HII region molecular gas has lower line intensities, ~3K, and smaller line widths, ~3.5 km/s. These latter characteristics are similar to those found for quiescent molecular clouds in the GRS. Our sample nebulae thus show evidence for an evolutionary sequence wherein small, dense molecular gas clumps associated with UC HII regions grow into older compact nebulae and finally fragment and dissipate into large, diffuse nebulae.

The FCRAO Survey of the Taurus Molecular Cloud observed the 12CO and 13CO J=1-0 emission from 98 square degrees of this important, nearby star forming region. This set of data with 45" resolution comprises the highest spatial dynamic range image of an individual molecular cloud constructed to date, and provides valuable insights to the molecular gas distribution, kinematics, and the star formation process. In this contribution, we describe the observations, calibration, data processing, and characteristics of the noise and line emission of the survey. The angular distribution of 12CO and 13CO emission over 1 km/s velocity intervals and the full velocity extent of the cloud are presented. These reveal a complex, dynamic medium of cold, molecular gas.

We report the results of a 100 square degree survey of the Taurus Molecular Cloud region in the J = 1-0 transition of 12CO and 13CO. The image of the cloud in each velocity channel includes ~ 3 million Nyquist sampled pixels on a 20" grid. The high sensitivity and large linear dynamic range of the maps in both isotopologues reveal a very complex, highly structured cloud morphology. There are large scale correlated structures evident in 13CO emission having very fine dimensions, including filaments, cavities, and rings. The 12CO emission shows a quite different structure, with particularly complex interfaces between regions of greater and smaller column density defining the boundaries of the largest-scale cloud structures. The axes of the striations seen in the 12CO emission from relatively diffuse gas are aligned with the direction of the magnetic field. Using a column density-dependent model for the CO fractional abundance, we derive the mass of the region mapped to be 24,000 solar masses, a factor of three greater than would be obtained with canonical CO abundance restricted to the high column density regions. We determine that half the mass of the cloud is in regions having column density below 2.1x10^{21} per square cm. The distribution of young stars in the region covered is highly nonuniform, with the probability of finding a star in a pixel with a specified column density rising sharply for N(H2) = 6x10^{21} cm^{-2}. We determine a relatively low star formation efficiency (mass of young stars/mass of molecular gas), between 0.3 and 1.2 %, and an average star formation rate during the past 3 Myr of 8x10^{-5} stars yr^{-1}.

Velocity anisotropy induced by MHD turbulence is investigated using computational simulations and molecular line observations of the Taurus molecular cloud. A new analysis method is presented to evaluate the degree and angle of velocity anisotropy using spectroscopic imaging data of interstellar clouds. The efficacy of this method is demonstrated on model observations derived from three dimensional velocity and density fields from the set of numerical MHD simulations that span a range of magnetic field strengths. The analysis is applied to 12CO J=1-0 imaging of a sub-field within the Taurus molecular cloud. Velocity anisotropy is identified that is aligned within 10 degrees of the mean local magnetic field direction derived from optical polarization measurements. Estimated values of the field strength based on velocity anisotropy are consistent with results from other methods. When combined with new column density measurements for Taurus, our magnetic field strength estimate indicates that the envelope of the cloud is magnetically subcritical. These observations favor strong MHD turbulence within the low density, sub-critical, molecular gas substrate of the Taurus cloud.

The observational record of turbulence within the molecular gas phase of the interstellar medium is summarized. We briefly review the analysis methods used to recover the velocity structure function from spectroscopic imaging and the application of these tools on sets of cloud data. These studies identify a near-invariant velocity structure function that is independent of local the environment and star formation activity. Such universality accounts for the cloud-to-cloud scaling law between the global line-width and size of molecular clouds found by Larson (1981) and constrains the degree to which supersonic turbulence can regulate star formation. In addition, the evidence for large scale driving sources necessary to sustain supersonic flows is summarized.

We present an overview of data available for the Ophiuchus and Perseus molecular clouds from ``Phase I'' of the COMPLETE Survey of Star-Forming Regions. This survey provides a range of data complementary to the Spitzer Legacy Program ``From Molecular Cores to Planet Forming Disks.'' Phase I includes: Extinction maps derived from 2MASS near-infrared data using the NICER algorithm; extinction and temperature maps derived from IRAS 60 and 100um emission; HI maps of atomic gas; 12CO and 13CO maps of molecular gas; and submillimetre continuum images of emission from dust in dense cores. Not unexpectedly, the morphology of the regions appears quite different depending on the column-density tracer which is used, with IRAS tracing mainly warmer dust and CO being biased by chemical, excitation and optical depth effects. Histograms of column-density distribution are presented, showing that extinction as derived from 2MASS/NICER gives the closest match to a log-normal distribution as is predicted by numerical simulations. All the data presented in this paper, and links to more detailed publications on their implications are publically available at the COMPLETE website.

The Boston University-Five College Radio Astronomy Observatory Galactic Ring Survey is a new survey of Galactic 13CO (1-0) emission. The survey used the SEQUOIA multi pixel array on the Five College Radio Astronomy Observatory 14 m telescope to cover a longitude range of l = 18 deg-55.7 deg and a latitude range of |b| < 1 deg, a total of 75.4 square degrees. Using both position-switching and On-The-Fly mapping modes, we achieved an angular sampling of 22 arcsec, better than half of the telescope's 46 arcsec angular resolution. The survey's velocity coverage is -5 to 135 km/s for Galactic longitudes l <= 40 deg and -5 to 85 km/s for Galactic longitudes l > 40 deg. At the velocity resolution of 0.21 km/s, the typical rms sensitivity is sigma(TA*)~0.13 K. The survey comprises a total of 1,993,522 spectra. We show integrated intensity images (zeroth moment maps), channel maps, position-velocity diagrams, and an average spectrum of the completed survey dataset. We also discuss the telescope and instrumental parameters, the observing modes, the data reduction processes, and the emission and noise characteristics of the dataset. The Galactic Ring Survey data are available to the community at www.bu.edu/galacticring or in DVD form by request.

The role of turbulent fragmentation in regulating the efficiency of star formation in interstellar clouds is examined from new wide field imaging of 12CO and 13CO J=1-0 emission from the Rosette and G216-2.5 molecular clouds. The Rosette molecular cloud is a typical star forming giant molecular cloud and G215-2.5 is a massive molecular cloud with no OB stars and very little low mass star formation. The properties of the turbulent gas flow are derived from the set of eigenvectors and eigenimages generated by Principal Component Analysis of the spectroscopic data cubes. While the two clouds represent quite divergent states of star formation activity, the velocity structure functions for both clouds are similar. The sonic scale, lambda_S, defined as the spatial scale at which turbulent velocity fluctuations are equivalent to the local sound speed, and the turbulent Mach number evaluated at 1 pc, M_{1pc}, are derived for an ensemble of clouds including the Rosette and, G216-2.5 regions that span a large range in star formation activity. We find no evidence for the positive correlations between these quantities and the star formation efficiency, that are predicted by turbulent fragmentation models. A correlation does exist between the star formation efficiency and the sonic scale for a subset of clouds with L_{FIR}/M(H_2) > 1 that are generating young stellar clusters. Turbulent fragmentation must play a limited and non-exclusive role in determining the yield of stellar masses within interstellar clouds.

The universality of interstellar turbulence is examined from observed structure functions of 27 giant molecular clouds and Monte Carlo modeling. We show that the structure functions, dv=v0 l^gamma, derived from wide field imaging of CO J=1-0 emission from individual clouds are described by a narrow range in the scaling exponent, gamma, and the scaling coefficient, v0. The similarity of turbulent structure functions emphasizes the universality of turbulence in the molecular interstellar medium and accounts for the cloud-to-cloud size-line width relationship initially identified by Larson (1981). The degree of turbulence universality is quantified by Monte Carlo simulations that reproduce the mean squared velocity residuals of the observed cloud-to-cloud relationship. Upper limits to the variation of the scaling amplitudes and exponents for molecular clouds are ~10-20%. The measured invariance of turbulence for molecular clouds with vastly different sizes, environments, and star formation activity suggests a common formation mechanism such as converging turbulent flows within the diffuse ISM and a limited contribution of energy from sources within the cloud with respect to large scale driving mechanisms.

The relationship between the star formation rate and surface density of neutral gas within the disk of M33 is examined with new imaging observations of CO J=1-0 emission gathered with the FCRAO 14m telescope and IRAS HiRes images of the 60 micron and 100 micron emission. The Schmidt law, Sigma_SFR ~ Sigma_gas^n, is constructed using radial profiles of the HI 21cm, CO, and far infrared emission. A strong correlation is identified between the star formation rate and molecular gas surface density. This suggests that the condensation of giant molecular clouds is the limiting step to star formation within the M33 disk. The corresponding molecular Schmidt index, n_{mol}, is 1.36 +/- 0.08. The star formation rate has a steep dependence on total mass gas surface density, (Sigma_{HI}+Sigma_{H_2}), owing to the shallow radial profile of the atomic gas which dominates the total gas surface density for most radii. The disk pressure of the gas is shown to play a prominent role in regulating the molecular gas fraction in M33.

Increasingly sophisticated observational tools and techniques are now being developed for probing the nature of interstellar turbulence. At the same time, theoretical advances in understanding the nature of turbulence and its effects on the structure of the ISM and on star formation are occurring at a rapid pace, aided in part by numerical simulations. These increased capabilities on both fronts open new opportunities for strengthening the links between observation and theory,and for meaningful comparisons between the two.

A large scale, high resolution map of CO(1-0) emission toward the OMC-2 and OMC-3 star forming regions is presented. The map is a mosaic of 46 fields using the Berkeley-Illinois-Maryland Array (BIMA) and covers $\sim 10'\times 15'$ at $\sim 10" $ resolution. These data are combined with singledish FCRAO observations and analyzed to identify and determine the properties of nine protostellar outflows. The BIMA data alone almost completely resolve out the cloud emission at central velocities and only recover 1/20 of the flux in the high velocity gas showing that outflows are generally broadly dispersed over $\sim 1'$ angular scales. All nine identified outflows emanate from known Class 0 or borderline Class 0/I sources, are associated with knots of shocked \h2\ emission, and have short dynamical times. It is suggested that only the youngest, most spatially compact, and energetic outflows have been found and that more distributed high velocity gas undetected by BIMA is due to older outflows continuing through the Class I phase of protostellar evolution. The mechanical energy injection rate into the cloud is estimated to be $\sim 1.5 L_\odot$ which is comparable to the turbulent energy dissipation rate. Outf`lows appear capable, therefore, of sustaining cloud turbulence but a high starformation rate is required implying a short cloud lifetime $\leq 5$ `Myr.

The multivariate tool of Principal Component Analysis (PCA) is applied to 23 fields in the FCRAO CO Survey of the Outer Galaxy. PCA enables the identification of line profile differences which are assumed to be generated from fluctuations within a turbulent velocity field. The variation of these velocity differences with spatial scale within a molecular region is described by a singular power law, delta v= c L^alpha which can be used as a powerful diagnostic to turbulent motions. For the ensemble of 23 fields, we find a mean value alpha = 0.62 +- 0.11. From a recent calibration of this method using fractal Brownian motion simulations (Brunt & Heyer 2001), the measured velocity difference-size relationship corresponds to an energy spectrum, E(k), which varies as k^-beta, where beta = 2.17 +- 0.31. We compare our results to both decaying and forced hydrodynamic simulations of turbulence. We conclude that energy must be continually injected into the regions to replenish that lost by dissipative processes such as shocks. The absence of large, widely distributed shocks within the targeted fields suggests that the energy is injected at spatial scales less than several pc.

We present spectra of the extended 12CO and 13CO J=1-0 emission along the major axes of 17 nearby galaxies. Spatial variations in the ratio of CO and 13CO integrated intensities, R, are found in nearly every galaxy observed. There is an overall variation in R of 20-40% from the inner 2 kpc to the disk. Roughly a third of the survey galaxies have such gradients in R detected above the 2-sigma confidence level. Though some galaxies show a lower central value of R, on average R inside 2 kpc is 10-30% higher than R outside of 2 kpc. An increase in the CO/13CO intensity ratio from disk to nucleus may imply that the conversion factor between CO intensity and H2 column density, X, is lower in galactic nuclei. Also variations in physical conditions, most notably the gas kinetic temperature, affect both R and X. Abundance variations probably do not cause the gradient in R, though we do not rule out a decrease in effective cloud column densities in galactic nuclei possibly caused by destructive starburst superwinds. A modest rise in temperature (less than a factor of 2 or 3) from outside a 2 kpc radius towards the nucleus can easily account for the observed gradient. These results support previous work implying that X is lower in the center of the Milky Way and probably most galactic nuclei. Therefore calculating H2 masses using the standard Galactic X-factor, especially within the central few kpc of galaxies, overestimates the true mass by factors of a few. The standard X-factor still appears to be appropriate for galactic disks.

A summary of global properties and an evaluation of the equilibrium state of molecular regions in the outer Galaxy are presented from the decomposition of the FCRAO Outer Galaxy Survey and targeted 12CO and 13CO observations of four giant molecular cloud complexes. The ensemble of identified objects includes both small, isolated clouds and clumps within larger cloud complexes. 12CO velocity dispersions show little variation with cloud sizes for radii less than 10 pc. It is demonstrated that the internal motions of regions with molecular masses greater than 10**4 msuns are bound by self gravity, yet, the constituent clumps of cloud complexes and isolated molecular clouds with M < 10**3 msuns are not in self gravitational equilibrium. The required external pressures to maintain the equilibrium of this population are (1-2)x10**4 cm-3-K.

We demonstrate the capability of Principal Component Analysis (PCA) as applied by Heyer & Schloerb (1997) to extract the statistics of turbulent interstellar velocity fields as measured by the energy spectrum, E(k)= k^-beta. Turbulent velocity and density fields with known statistics are generated from fBm simulations. These fields are translated to the observational domain, T(x,y,v), considering the excitation of molecular rotational energy levels and radiative transfer. Using PCA and the characterization of velocity and spatial scales from the eigenvectors and eigenimages respectively, a relationship is identified which describes the magnitude of line profile differences, deltav, and the scale, L, over which these differences occur, deltav=L^alpha. From a series of models with varying values of beta, we find, alpha = 0.33beta - 0.05 for 1 < beta < 3. This provides the basic calibration between the intrinsic velocity field statistics to observational measures and a diagnostic for turbulent flows in the interstellar medium. We also investigate the effects of noise, line opacity, and finite resolution on these results.

We analyze the embedded stellar content in the vicinity of the W3/W4/W5 HII regions using the FCRAO Outer Galaxy 12CO(J=1-0) Survey, the IRAS Point Source Catalog, published radio continuum surveys, and new near-infrared and molecular line observations. Thirty-four IRAS Point Sources are identified that have far-infrared colors characteristic of embedded star forming regions, and we have obtained K' mosaics and 13CO(J=1-0) maps for 32 of them. Ten of the IRAS sources are associated with an OB star and 19 with a stellar cluster, although three OB stars are not identified with a cluster. Half of the embedded stellar population identified in the K' images is found in just the 5 richest clusters, and 61% is contained in IRAS sources associated with an embedded OB star. Thus rich clusters around OB stars contribute substantially to the stellar population currently forming in the W3/W4/W5 region. Approximately 39% of the cluster population is embedded in small clouds with an average mass of ~130 Mo that are located as far as 100 pc from the W3/W4/W5 cloud complex. We speculate that these small clouds are fragments of a cloud complex dispersed by previous episodes of massive star formation. Finally, we find that 4 of the 5 known embedded massive star forming sites in the W3 molecular cloud are found along the interface with the W4 HII region despite the fact that most of the molecular mass is contained in the interior regions of the cloud. These observations are consistent with the classical notion that the W4 HII region has triggered massive star formation along the eastern edge of the W3 molecular cloud.