We present an overview of the CMZoom survey and its first data release. CMZoom is the first blind, high-resolution survey of the Central Molecular Zone (CMZ; the inner 500 pc of the Milky Way) at wavelengths sensitive to the pre-cursors of high-mass stars. CMZoom is a 500-hour Large Program on the Submillimeter Array (SMA) that mapped at 1.3 mm all of the gas and dust in the CMZ above a molecular hydrogen column density of 10^23 cm^-2 at a resolution of ~3" (0.1 pc). In this paper, we focus on the 1.3 mm dust continuum and its data release, but also describe CMZoom spectral line data which will be released in a forthcoming publication. While CMZoom detected many regions with rich and complex substructure, its key result is an overall deficit in compact substructures on 0.1 - 2 pc scales (the compact dense gas fraction: CDGF). In comparison with clouds in the Galactic disk, the CDGF in the CMZ is substantially lower, despite having much higher average column densities. CMZ clouds with high CDGFs are well-known sites of active star formation. The inability of most gas in the CMZ to form compact substructures is likely responsible for the dearth of star formation in the CMZ, surprising considering its high density. The factors responsible for the low CDGF are not yet understood but are plausibly due to the extreme environment of the CMZ, having far-reaching ramifications for our understanding of the star formation process across the cosmos.

Bondi (1952) and Parker (1958} derived a steady-state solution for Bernouilli's equation in spherical symmetry around a point mass for two cases, respectively, an inward accretion flow and an outward wind. Left unanswered were the stability of the steady-state solution, the solution itself of time-dependent flows, whether the time-dependent flows would evolve to the steady-state, and under what conditions a transonic flow would develop. In a Hamiltonian description, we find that the steady state solution is equivalent to the Lagrangian implying that time-dependent flows evolve to the steady state. We find that the second variation is definite in sign for isothermal and adiabatic flows, implying at least linear stability. We solve the partial differential equation for the time-dependent flow as an initial-value problem and find that a transonic flow develops under a wide range of realistic initial conditions. We present some examples of time-dependent solutions.

The kinetic energy of supersonic turbulence within interstellar clouds is subject to cooling by dissipation in shocks and subsequent line radiation. The clouds are therefore susceptible to a condensation process controlled by the specific entropy. In a form analogous to the thermodynamic entropy, the entropy for supersonic turbulence is proportional to the log of the product of the mean turbulent velocity and the size scale. We derive a dispersion relation for the growth of entropic instabilities in a spherical self-gravitating cloud and find that there is a critical maximum dissipation time scale, about equal to the crossing time, that allows for fragmentation and subsequent star formation. However, the time scale for the loss of turbulent energy may be shorter or longer, for example with rapid thermal cooling or the injection of mechanical energy. Differences in the time scale for energy loss in different star-forming regions may result in differences in the outcome, for example, in the initial mass function.

The formation of massive stars is a long standing problem. Although a number of theories of massive star formation exist, ideas appear to converge to a disk-mediated accretion scenario. Here we present radiative hydrodynamic simulations of a star accreting mass via a disk embedded in a torus. We use a Monte Carlo based radiation hydrodynamics code to investigate the impact that ionizing radiation has on the torus. Ionized regions in the torus midplane are found to be either gravitationally trapped or in pressure driven expansion depending on whether or not the size of the ionized region exceeds a critical radius. Trapped Hii regions in the torus plane allow accretion to progress, while expanding Hii regions disrupt the accretion torus preventing the central star from aggregating more mass, thereby setting the star's final mass. We obtain constraints for the luminosities and torus densities that lead to both scenarios.

This paper investigates the gravitational trapping of HII regions predicted by steady-state analysis using radiation hydrodynamical simulations. We present idealised spherically symmetric radiation hydrodynamical simulations of the early evolution of HII regions including the gravity of the central source. As with analytic steady state solutions of spherically symmetric ionised Bondi accretion flows, we find gravitationally trapped HII regions with accretion through the ionisation front onto the source. We found that, for a constant ionising luminosity, fluctuations in the ionisation front are unstable. This instability only occurs in this spherically symmetric accretion geometry. In the context of massive star formation, the ionising luminosity increases with time as the source accretes mass. The maximum radius of the recurring HII region increases on the accretion timescale until it reaches the sonic radius, where the infall velocity equals the sound speed of the ionised gas, after which it enters a pressure-driven expansion phase. This expansion prevents accretion of gas through the ionisation front, the accretion rate onto the star decreases to zero, and it stops growing from accretion. Because of the time required for any significant change in stellar mass and luminosity through accretion our simulations keep both mass and luminosity constant and follow the evolution from trapped to expanding in a piecewise manner. Implications of this evolution of HII regions include a continuation of accretion of material onto forming stars for a period after the star starts to emit ionising radiation, and an extension of the lifetime of ultracompact HII regions.

We investigate the general stability of 1D spherically symmetric ionized Bondi accretion onto a massive object in the specific context of accretion onto a young stellar object. We first derive a new analytic expression for a steady state two temperature solution that predicts the existence of compact and hypercompact HII regions. We then show that this solution is only marginally stable if ionization is treated self-consistently. This leads to a recurring collapse of the HII region over time. We derive a semi-analytic model to explain this instability, and test it using spatially converged 1D radiation hydrodynamical simulations. We discuss the implications of the 1D instability on 3D radiation hydrodynamics simulations of supersonic accreting flows.

Stars like our Sun form in self-gravitating dense and cold structures within interstellar clouds, called pre-stellar cores. Although much is known about the physical structure of dense clouds just before and soon after the switch-on of a protostar, the central few thousand astronomical units (au) of pre-stellar cores are unexplored. It is within these central regions that stellar systems assemble and fragmentation may take place, with the consequent formation of binaries and multiple systems. We present ALMA Band 6 observations (ACA and 12m array) of the dust continuum emission of the 8 Msun pre-stellar core L1544, with angular resolution of 2'' x 1.6'' (linear resolution 270 au x 216 au). Within the primary beam, a compact region of 0.1 Msun, which we call a "kernel", has been unveiled. The kernel is elongated, with a central flat zone with radius Rker ~ 10'' (~ 1400 au). The average number density within Rker is ~1 x 10^6 cm^{-3}, with possible local density enhancements. The region within Rker appears to have fragmented, but detailed analysis shows that similar substructure can be reproduced by synthetic interferometric observations of a smooth centrally concentrated dense core with a similar central flat zone. The presence of a smooth kernel within a dense core is in agreement with non-ideal magneto-hydro-dynamical simulations of a contracting cloud core with a peak number density of 1 x 10^7 cm^{-3}. Dense cores with lower central densities are completely filtered out when simulated 12m-array observations are carried out. These observations demonstrate that the kernel of dynamically evolved dense cores can be investigated at high angular resolution with ALMA.

The evolution of molecular clouds in galactic centres is thought to differ from that in galactic discs due to a significant influence of the external gravitational potential. We present a set of numerical simulations of molecular clouds orbiting on the 100-pc stream of the Central Molecular Zone (the central $\sim500$ pc of the Galaxy) and characterise their morphological and kinematic evolution in response to the background potential and eccentric orbital motion. We find that the clouds are shaped by strong shear and torques, by tidal and geometric deformation, and by their passage through the orbital pericentre. Within our simulations, these mechanisms control cloud sizes, aspect ratios, position angles, filamentary structure, column densities, velocity dispersions, line-of-sight velocity gradients, spin angular momenta, and kinematic complexity. By comparing these predictions to observations of clouds on the Galactic Centre 'dust ridge', we find that our simulations naturally reproduce a broad range of key observed morphological and kinematic features, which can be explained in terms of well-understood physical mechanisms. We argue that the accretion of gas clouds onto the central regions of galaxies, where the rotation curve turns over and the tidal field is fully compressive, is accompanied by transformative dynamical changes to the clouds, leading to collapse and star formation. This can generate an evolutionary progression of cloud collapse with a common starting point, which either marks the time of accretion onto the tidally-compressive region or of the most recent pericentre passage. Together, these processes may naturally produce the synchronised starbursts observed in numerous (extra)galactic nuclei.

The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.37 mm wavelength regime in both line and dust continuum emission. We present the spectral line setup of the CORE survey and a case study for W3(H2O). At ~0.35" (700 AU at 2 kpc) resolution, the W3(H2O) clump fragments into two cores (West and East), separated by ~2300 AU. Velocity shifts of a few km/s are observed in the dense-gas tracer, CH3CN, across both cores, consistent with rotation and perpendicular to the directions of two bipolar outflows, one emanating from each core. The kinematics of the rotating structure about W3(H2O) W shows signs of differential rotation of material, possibly in a disk-like object. The observed rotational signature around W3(H2O) E may be due to a disk-like object, an unresolved binary (or multiple) system, or a combination of both. We fit the emission of CH3CN (12-11) K = 4-6 and derive a gas temperature map with a median temperature of ~165 K across W3(H2O). We create a Toomre Q map to study the stability of the rotating structures against gravitational instability. The rotating structures appear to be Toomre unstable close to their outer boundaries, with a possibility of further fragmentation in the differentially-rotating core W3(H2O) W. Rapid cooling in the Toomre-unstable regions supports the fragmentation scenario. Combining millimeter dust continuum and spectral line data toward the famous high-mass star-forming region W3(H2O), we identify core fragmentation on large scales, and indications for possible disk fragmentation on smaller spatial scales.

Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high angular resolution (~0.4'') for a sample of 20 well-known high-mass star-forming regions with distances below 5.5kpc and luminosities larger than 10^4Lsun. Results: We present the overall survey scope, the selected sample, the observational setup and the main goals of CORE. Scientifically, we concentrate on the mm continuum emission on scales on the order of 1000AU. We detect strong mm continuum emission from all regions, mostly due to the emission from cold dust. The fragmentation properties of the sample are diverse. We see extremes where some regions are dominated by a single high-mass core whereas others fragment into as many as 20 cores. A minimum-spanning-tree analysis finds fragmentation at scales on the order of the thermal Jeans length or smaller suggesting that turbulent fragmentation is less important than thermal gravitational fragmentation. The diversity of highly fragmented versus singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths. Conclusions: The smallest observed separations between cores are found around the angular resolution limit which indicates that further fragmentation likely takes place on even smaller spatial scales. The CORE project with its numerous spectral line detections will address a diverse set of important physical and chemical questions in the field of high-mass star formation.

We present a composite model and radiative transfer simulations of the massive star forming core W33A MM1. The model was tailored to reproduce the complex features observed with ALMA at $\approx 0.2$ arcsec resolution in CH$_3$CN and dust emission. The MM1 core is fragmented into six compact sources coexisting within $\sim 1000$ au. In our models, three of these compact sources are better represented as disc-envelope systems around a central (proto)star, two as envelopes with a central object, and one as a pure envelope. The model of the most prominent object (Main) contains the most massive (proto)star ($M_\star\approx7~M_\odot$) and disc+envelope ($M_\mathrm{gas}\approx0.4~M_\odot$), and is the most luminous ($L_\mathrm{Main} \sim 10^4~L_\odot$). The model discs are small (a few hundred au) for all sources. The composite model shows that the elongated spiral-like feature converging to the MM1 core can be convincingly interpreted as a filamentary accretion flow that feeds the rising stellar system. The kinematics of this filament is reproduced by a parabolic trajectory with focus at the center of mass of the region. Radial collapse and fragmentation within this filament, as well as smaller filamentary flows between pairs of sources are proposed to exist. Our modelling supports an interpretation where what was once considered as a single massive star with a $\sim 10^3$ au disc and envelope, is instead a forming stellar association which appears to be virialized and to form several low-mass stars per high-mass object.

Although, from a theoretical point of view, magnetic fields are believed to have a significant role during the early stages of star formation, especially during the main accretion phase, the magnetic field strength and topology is poorly constrained in the youngest accreting Class 0 protostars that lead to the formation of solar-type stars. We carried out observations of the polarized dust continuum emission with the SMA interferometer at 0.87mm, in order to probe the structure of the magnetic field in a sample of 12 low-mass Class 0 envelopes, including both single protostars and multiple systems, in nearby clouds. Our SMA observations probe the envelope emission at scales of 600-5000 au with a spatial resolution ranging from 600 to 1500 au depending on the source distance. We report the detection of linearly polarized dust continuum emission in all of our targets, with average polarization fractions ranging from 2% to 10% in these protostellar envelopes. The polarization fraction decreases with the continuum flux density, which translates into a decrease with the H2 column density within an individual envelope. Our analysis show that the envelope-scale magnetic field is preferentially observed either aligned or perpendicular to the outflow direction. Interestingly, our results suggest for the first time a relation between the orientation of the magnetic field and the rotational energy of envelopes, with a larger occurrence of misalignment in sources where strong rotational motions are detected at hundreds to thousands of au scales. We also show that the best agreement between the magnetic field and outflow orientation is found in sources showing no small-scale multiplicity and no large disks at about 100 au scales.

The role of the magnetic field during protostellar collapse is poorly constrained from an observational point of view, although it could be significant if we believe state-of-the-art models of protostellar formation. We present polarimetric observations of the 233 GHz thermal dust continuum emission obtained with ALMA in the B335 Class 0 protostar. Linearly polarized dust emission arising from the circumstellar material in the envelope of B335 is detected at all scales probed by our observations, from radii of 50 to 1000 au. The magnetic field structure producing the dust polarization has a very ordered topology in the inner envelope, with a transition from a large-scale poloidal magnetic field, in the outflow direction, to strongly pinched in the equatorial direction. This is probably due to magnetic field lines being dragged along the dominating infall direction since B335 does not exhibit prominent rotation. Our data and their qualitative comparison to a family of magnetized protostellar collapse models show that, during the magnetized collapse in B335, the magnetic field is maintaining a high level of organization from scales 1000 au to 50 au: this suggests the field is dynamically relevant and capable of influencing the typical outcome of protostellar collapse, such as regulating the disk size in B335.

We present intensity-corrected Herschel maps at 100 um, 160 um, 250 um, 350 um, and 500 um for 56 isolated low-mass clouds. We determine the zero-point corrections for Herschel PACS and SPIRE maps from the Herschel Science Archive (HSA) using Planck data. Since these HSA maps are small, we cannot correct them using typical methods. Here, we introduce a technique to measure the zero-point corrections for small Herschel maps. We use radial profiles to identify offsets between the observed HSA intensities and the expected intensities from Planck. Most clouds have reliable offset measurements with this technique. In addition, we find that roughly half of the clouds have underestimated HSA-SPIRE intensities in their outer envelopes relative to Planck, even though the HSA-SPIRE maps were previously zero-point corrected. Using our technique, we produce corrected Herschel intensity maps for all 56 clouds and determine their line-of-sight average dust temperatures and optical depths from modified black body fits. The clouds have typical temperatures of ~ 14-20 K and optical depths of ~ 1e-5 - 1e-3. Across the whole sample, we find an anti-correlation between temperature and optical depth. We also find lower temperatures than what was measured in previous Herschel studies, which subtracted out a background level from their intensity maps to circumvent the zero-point correction. Accurate Herschel observations of clouds are key to obtain accurate density and temperature profiles. To make such future analyses possible, intensity-corrected maps for all 56 clouds are publicly available in the electronic version.

The star formation rate in the Central Molecular Zone (CMZ) is an order of magnitude lower than predicted according to star formation relations that have been calibrated in the disc of our own and nearby galaxies. Understanding how and why star formation appears to be different in this region is crucial if we are to understand the environmental dependence of the star formation process. Here, we present the detection of a sample of high-mass cores in the CMZ's "dust ridge" that have been discovered with the Submillimeter Array as part of the CMZoom survey. These cores range in mass from ~ 50 - 2150 Msun within radii of 0.1 - 0.25 pc. All appear to be young (pre-UCHII), meaning that they are prime candidates for representing the initial conditions of high-mass stars and sub-clusters. We report that at least two of these cores ('c1' and 'e1') contain young, high-mass protostars. We compare all of the detected cores with high-mass cores in the Galactic disc and find that they are broadly similar in terms of their masses and sizes, despite being subjected to external pressures that are several orders of magnitude greater - ~ 10^8 K/cm^3, as opposed to ~ 10^5 K/cm^3. The fact that > 80% of these cores do not show any signs of star-forming activity in such a high-pressure environment leads us to conclude that this is further evidence for an increased critical density threshold for star formation in the CMZ due to turbulence.

Herschel has revolutionized our ability to measure column densities (N$_{\rm H}$) and temperatures (T) of molecular clouds thanks to its far infrared multiwavelength coverage. However, the lack of a well defined background intensity level in the Herschel data limits the accuracy of the N$_{\rm H}$ and T maps. We provide a method that corrects the missing Herschel background intensity levels using the Planck model for foreground Galactic thermal dust emission. We present a Fourier method that combines the publicly available Planck model on large angular scales with the Herschel images on smaller angular scales. We apply our method to two regions spanning a range of Galactic environments: Perseus and the Galactic plane region around $l = 11\deg$ (HiGal--11). We post-process the combined dust continuum emission images to generate column density and temperature maps. We compare these to previously adopted constant--offset corrections. We find significant differences ($\gtrsim$20\%) over significant ($\sim$15\%) areas of the maps, at low column densities ($N_{\rm H}\lesssim10^{22}$\,cm$^{-2}$) and relatively high temperatures ($T\gtrsim20$\,K). We also apply our method to synthetic observations of a simulated molecular cloud to validate our method. Our method successfully corrects the Herschel images, including both the constant--offset intensity level and the scale-dependent background variations measured by Planck. Our method improves the previous constant--offset corrections, which did not account for variations in the background emission levels.

Pre-stellar cores represent the initial conditions in the process of star and planet formation, therefore it is important to study their physical and chemical structure. Because of their volatility, nitrogen-bearing molecules are key to study the dense and cold gas present in pre-stellar cores. The NH_3 rotational transition detected with Herschel-HIFI provides a unique combination of sensitivity and spectral resolution to further investigate physical and chemical processes in pre-stellar cores. Here we present the velocity-resolved Herschel-HIFI observations of the ortho-NH_3(1_0-0_0) line at 572 GHz and study the abundance profile of ammonia across the pre-stellar core L1544 to test current theories of its physical and chemical structure. Recently calculated collisional coefficients have been included in our non-LTE radiative transfer code to reproduce Herschel observations. A gas-grain chemical model, including spin-state chemistry and applied to the (static) physical structure of L1544 is also used to infer the abundance profile of ortho-NH_3 . The hyperfine structure of ortho-NH_3(1_0-0_0) is resolved for the first time in space. All the hyperfine components are strongly self-absorbed. The profile can be reproduced if the core is contracting in quasi-equilibrium, consistent with previous work, and if the NH_3 abundance is slightly rising toward the core centre, as deduced from previous interferometric observations of para-NH_3(1,1). The chemical model overestimates the NH_3 abundance at radii between ~ 4000 and 15000 AU by about two orders of magnitude and underestimates the abundance toward the core centre by more than one order of magnitude. Our observations show that chemical models applied to static clouds have problems in reproducing NH_3 observations.

We targeted the massive star forming region W33A using the Atacama Large Sub/Millimeter Array (ALMA) in band 6 (230 GHz) and 7 (345 GHz) to search for a sub-1000au disc around the central O-type massive young stellar object (MYSO) W33A MM1-Main. Our data achieves a resolution of ~0.2" (~500au) and resolves the central core, MM1, into multiple components and reveals complex and filamentary structures. There is strong molecular line emission covering the entire MM1 region. The kinematic signatures are inconsistent with only Keplerian rotation although we propose that the shift in the emission line centroids within ~1000au of MM1-Main could hint at an underlying compact disc with Keplerian rotation. We cannot however rule out the possibility of an unresolved binary or multiple system. A putative smaller disc could be fed by the large scale spiral `feeding filament' we detect in both gas and dust emission. We also discuss the nature of the now-resolved continuum sources.

The inner few hundred parsecs of the Milky Way, the Central Molecular Zone (CMZ), is our closest laboratory for understanding star formation in the extreme environments (hot, dense, turbulent gas) that once dominated the universe. We present an update on the first large-area survey to expose the sites of star formation across the CMZ at high-resolution in submillimeter wavelengths: the CMZoom survey with the Submillimeter Array (SMA). We identify the locations of dense cores and search for signatures of embedded star formation. CMZoom is a three-year survey in its final year and is mapping out the highest column density regions of the CMZ in dust continuum and a variety of spectral lines around 1.3 mm. CMZoom combines SMA compact and subcompact configurations with single-dish data from BGPS and the APEX telescope, achieving an angular resolution of about 4" (0.2 pc) and good image fidelity up to large spatial scales.

Molecules with hyperfine splitting of their rotational line spectra are useful probes of optical depth, via the relative line strengths of their hyperfine components.The hyperfine splitting is particularly advantageous in interpreting the physical conditions of the emitting gas because with a second rotational transition, both gas density and temperature can be derived. For HCN however, the relative strengths of the hyperfine lines are anomalous. They appear in ratios which can vary significantly from source to source, and are inconsistent with local thermodynamic equilibrium. This is the HCN hyperfine anomaly, and it prevents the use of simple LTE models of HCN emission to derive reliable optical depths. In this paper we demonstrate how to model HCN hyperfine line emission, and derive accurate line ratios, spectral line shapes and optical depths. We show that by carrying out radiative transfer calculations over each hyperfine level individually, as opposed to summing them over each rotational level, the anomalous hyperfine emission emerges naturally. To do this requires not only accurate radiative rates between hyperfine states, but also accurate collisional rates. We investigate the effects of different sets of hyperfine collisional rates, derived via the 'proportional method' and through direct recoupling calculations. Through an extensive parameter sweep over typical low mass star forming conditions, we show the HCN line ratios to be highly variable to optical depth. We also reproduce an observed effect whereby the red-blue asymmetry of the hyperfine lines (an infall signature) switches sense within a single rotational transition.

We present new spectral line observations of the CH3CN molecule in the accretion disk around the massive protostar IRAS 20126+4104 with the Submillimeter Array that for the first time measure the disk density, temperature, and rotational velocity with sufficient resolution (0.37", equivalent to ~600 AU) to assess the gravitational stability of the disk through the Toomre-Q parameter. Our observations resolve the central 2000 AU region that shows steeper velocity gradients with increasing upper state energy, indicating an increase in the rotational velocity of the hotter gas nearer the star. Such spin-up motions are characteristics of an accretion flow in a rotationally supported disk. We compare the observed data with synthetic image cubes produced by three-dimensional radiative transfer models describing a thin flared disk in Keplerian motion enveloped within the centrifugal radius of an angular-momentum-conserving accretion flow. Given a luminosity of 1.3x10^4 Lsun, the optimized model gives a disk mass of 1.5 Msun and a radius of 858 AU rotating about a 12.0 Msun protostar with a disk mass accretion rate of 3.9x10^{-5} Msun/yr. Our study finds that, in contrast to some theoretical expectations, the disk is hot and stable to fragmentation with Q > 2.8 at all radii which permits a smooth accretion flow. These results put forward the first constraints on gravitational instabilities in massive protostellar disks, which are closely connected to the formation of companion stars and planetary systems by fragmentation.

We have recently published observations of significant flux density variations at 1.3 cm in HII regions in the star forming regions Sgr B2 Main and North (De Pree et al. 2014). To further study these variations, we have made new 7 mm continuum and recombination line observations of Sgr B2 at the highest possible angular resolution of the Karl G. Jansky Very Large Array (VLA). We have observed Sgr B2 Main and North at 42.9 GHz and at 45.4 GHz in the BnA configuration (Main) and the A configuration (North). We compare these new data to archival VLA 7 mm continuum data of Sgr B2 Main observed in 2003 and Sgr B2 North observed in 2001. We find that one of the 41 known ultracompact and hypercompact HII regions in Sgr B2 (K2-North) has decreased $\sim$27% in flux density from 142$\pm$14 mJy to 103$\pm$10 mJy (2.3$\sigma$) between 2001 and 2012. A second source, F3c-Main has increased $\sim$30% in flux density from 82$\pm$8 mJy to 107 $\pm$11 mJy (1.8$\sigma$) between 2003 and 2012. F3c-Main was previously observed to increase in flux density at 1.3 cm over a longer time period between 1989 and 2012 (De Pree et al. 2014). An observation of decreasing flux density, such as that observed in K2-North, is particularly significant since such a change is not predicted by the classical hypothesis of steady expansion of HII regions during massive star accretion. Our new observations at 7 mm, along with others in the literature, suggest that the formation of massive stars occurs through time-variable and violent accretion.

Ground-based observations at terahertz (THz) frequencies are a newly explorable area of astronomy for the next ten years. We discuss science cases for a first-generation 10-m class THz telescope, focusing on the Greenland Telescope as an example of such a facility. We propose science cases and provide quantitative estimates for each case. The largest advantage of ground-based THz telescopes is their higher angular resolution (~ 4 arcsec for a 10-m dish), as compared to space or airborne THz telescopes. Thus, high-resolution mapping is an important scientific argument. In particular, we can isolate zones of interest for Galactic and extragalactic star-forming regions. The THz windows are suitable for observations of high-excitation CO lines and [N II] 205 um lines, which are scientifically relevant tracers of star formation and stellar feedback. Those lines are the brightest lines in the THz windows, so that they are suitable for the initiation of ground-based THz observations. THz polarization of star-forming regions can also be explored since it traces the dust population contributing to the THz spectral peak. For survey-type observations, we focus on ``sub-THz'' extragalactic surveys, whose uniqueness is to detect galaxies at redshifts z ~ 1--2, where the dust emission per comoving volume is the largest in the history of the Universe. Finally we explore possibilities of flexible time scheduling, which enables us to monitor active galactic nuclei, and to target gamma-ray burst afterglows. For these objects, THz and submillimeter wavelength ranges have not yet been explored.

We present mid-IR (19 - 37 microns) imaging observations of S106 from SOFIA/FORCAST, complemented with IR observations from Spitzer/IRAC (3.6 - 8.0 microns), IRTF/MIRLIN (11.3 and 12.5 microns), and Herschel/PACS (70 and 160 microns). We use these observations, observations in the literature, and radiation transfer modeling to study the heating and composition of the warm (~ 100 K) dust in the region. The dust is heated radiatively by the source S106 IR, with little contributions from grain-electron collisions and Ly-alpha radiation. The dust luminosity is >~ (9.02 +/- 1.01) x 10^4 L_sun, consistent with heating by a mid- to late-type O star. We find a temperature gradient (~ 75 - 107 K) in the lobes, which is consistent with a dusty equatorial geometry around S106 IR. Furthermore, the SOFIA observations resolve several cool (~ 65 - 70 K) lanes and pockets of warmer (~ 75 - 90 K) dust in the ionization shadow, indicating that the environment is fragmented. We model the dust mass as a composition of amorphous silicates, amorphous carbon, big grains, very small grains, and PAHs. We present the relative abundances of each grain component for several locations in S106.

We present Atacama Large Millimeter/submillimeter Array (ALMA) line and continuum observations at 1.2mm with ~0.3" resolution that uncover a Keplerian-like disk around the forming O-type star AFGL 4176. The continuum emission from the disk at 1.21 mm (source mm1) has a deconvolved size of 870+/-110 AU x 330+/-300 AU and arises from a structure ~8 M_sun in mass, calculated assuming a dust temperature of 190 K. The first-moment maps, pixel-to-pixel line modeling, assuming local thermodynamic equilibrium (LTE), and position-velocity diagrams of the CH3CN J=13-12 K-line emission all show a velocity gradient along the major axis of the source, coupled with an increase in velocity at small radii, consistent with Keplerian-like rotation. The LTE line modeling shows that where CH3CN J=13-12 is excited, the temperatures in the disk range from ~70 to at least 300 K and that the H2 column density peaks at 2.8x10^24 cm^-2. In addition, we present Atacama Pathfinder Experiment (APEX) 12CO observations which show a large-scale outflow from AFGL 4176 perpendicular to the major axis of mm1, supporting the disk interpretation. Finally, we present a radiative transfer model of a Keplerian disk surrounding an O7 star, with a disk mass and radius of 12 M_sun and 2000 AU, that reproduces the line and continuum data, further supporting our conclusion that our observations have uncovered a Keplerian disk around an O-type star.

We present a method for modelling star-forming clouds that combines two different models of the thermal evolution of the interstellar medium (ISM). In the combined model, where the densities are low enough that at least some part of the spectrum is optically thin, a model of the thermodynamics of the diffuse ISM is more significant in setting the temperatures. Where the densities are high enough to be optically thick across the spectrum, a model of flux limited diffusion is more appropriate. Previous methods either model the low-density interstellar medium and ignore the thermal behaviour at high densities (e.g. inside collapsing molecular cloud cores), or model the thermal behaviour near protostars but assume a fixed background temperature (e.g. approximately 10 K) on large-scales. Our new method treats both regimes. It also captures the different thermal evolution of the gas, dust, and radiation separately. We compare our results with those from the literature, and investigate the dependence of the thermal behaviour of the gas on the various model parameters. This new method should allow us to model the ISM across a wide range of densities and, thus, develop a more complete and consistent understanding of the role of thermodynamics in the star formation process.

Low-mass stars are generally understood to form by the gravitational collapse of the dense molecular clouds known as starless cores. Continuum observations have not been able to distinguish among the several different hypotheses that describe the collapse because the predicted density distributions are the almost the same, as they are for all spherical self-gravitating clouds. However, the predicted contraction velocities are different enough that the models can be discriminated by comparing the velocities at large and small radii. This can be done by observing at least two different molecular line transitions that are excited at different densities. For example, the spectral lines of the H2O (110 - 101) and C18O (1-0) have critical densities for collisional de-excitation that differ by 5 orders of magnitude. We compare observations of these lines from the contracting starless core L1544 against the spectra predicted for several different hypothetical models of contraction including the Larson-Penston flow, the inside-out collapse of the singular isothermal sphere, the quasi-equilibrium contraction of an unstable Bonnor-Ebert sphere, and the non-equilibrium collapse of an over-dense Bonnor-Ebert sphere. Only the model of the unstable quasi-equilibrium Bonnor-Ebert sphere is able to produce the observed shapes of both spectral lines. This model allows us to interpret other observations of molecular lines in L1544 to find that the inward velocities seen in observations of CS(2-1) and N2H+ are located within the starless core itself, in particular in the region where the density profile follows an inverse square law. If these conclusions were to hold in the analysis of other starless cores, this would imply that the formation of hydrostatic clouds within the turbulent interstellar medium is not only possible but not exceptional and may be an evolutionary phase in low-mass star formation.

Massive stars ($M > 8$ \msun) typically form in parsec-scale molecular clumps that collapse and fragment, leading to the birth of a cluster of stellar objects. We investigate the role of magnetic fields in this process through dust polarization at 870 $\mu$m obtained with the Submillimeter Array (SMA). The SMA observations reveal polarization at scales of $\lsim$ 0.1 pc. The polarization pattern in these objects ranges from ordered hour-glass configurations to more chaotic distributions. By comparing the SMA data with the single dish data at parsec scales, we found that magnetic fields at dense core scales are either aligned within $40^\circ$ of or perpendicular to the parsec-scale magnetic fields. This finding indicates that magnetic fields play an important role during the collapse and fragmentation of massive molecular clumps and the formation of dense cores. We further compare magnetic fields in dense cores with the major axis of molecular outflows. Despite a limited number of outflows, we found that the outflow axis appears to be randomly oriented with respect to the magnetic field in the core. This result suggests that at the scale of accretion disks ($\lsim 10^3$ AU), angular momentum and dynamic interactions possibly due to close binary or multiple systems dominate over magnetic fields. With this unprecedentedly large sample massive clumps, we argue on a statistical basis that magnetic fields play an important role during the formation of dense cores at spatial scale of 0.01 - 0.1 pc in the context of massive star and cluster star formation.

A 12-m diameter radio telescope will be deployed to the Summit Station in Greenland to provide direct confirmation of a Super Massive Black Hole (SMBH) by observing its shadow image in the active galaxy M87. The telescope (Greenland Telescope: GLT) is to become one of the Very Long Baseline Interferometry (VLBI) stations at sub-millimeter (submm) regime, providing the longest baseline > 9,000 km to achieve an exceptional angular resolution of 20 micro arc sec at 350 GHz, which will enable us to resolve the shadow size of ~40 micro arc sec. The triangle with the longest baselines formed by the GLT, the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile, and the Submillimeter Array (SMA) in Hawaii will play a key role for the M87 observations. We have been working on the image simulations based on realistic conditions for a better understanding of the possible observed images. In parallel, retrofitting of the telescope and the site developments are in progress. Based on three years of opacity monitoring at 225 GHz, our measurements indicate that the site is excellent for submm observations, comparable to the ALMA site. The GLT is also expected to make single-dish observations up to 1.5 THz.

We examine the role of the gravitational instability in an isothermal, self-gravitating layer threaded by magnetic fields on the formation of filaments and dense cores. Using numerical simulation we follow the non-linear evolution of a perturbed equilibrium layer. The linear evolution of such a layer is described in the analytic work of Nagai et al (1998). We find that filaments and dense cores form simultaneously. Depending on the initial magnetic field, the resulting filaments form either a spiderweb-like network (for weak magnetic fields) or a network of parallel filaments aligned perpendicular to the magnetic field lines (for strong magnetic fields). Although the filaments are radially collapsing, the density profile of their central region (up to the thermal scale height) can be approximated by a hydrodynamical equilibrium density structure. Thus, the magnetic field does not play a significant role in setting the density distribution of the filaments. The density distribution outside of the central region deviates from the equilibrium. The radial column density distribution is then flatter than the expected power law of $r^{-4}$ and similar to filament profiles observed with Herschel. Our results does not explain the near constant filament width of $\sim 0.1$pc. However, our model does not include turbulent motions. It is expected that accretion-driven amplification of these turbulent motions provides additional support within the filaments against gravitational collapse. Finally, we interpret the filamentary network of the massive star forming complex G14.225-0.506 in terms of the gravitational instability model and find that the properties of the complex are consistent with being formed out of an unstable layer threaded by a strong, parallel magnetic field.

The Herschel Space Observatory's recent detections of water vapor in the cold, dense cloud L1544 allow a direct comparison between observations and chemical models for oxygen species in conditions just before star formation. We explain a chemical model for gas phase water, simplified for the limited number of reactions or processes that are active in extreme cold ($<$ 15 K). In this model, water is removed from the gas phase by freezing onto grains and by photodissociation. Water is formed as ice on the surface of dust grains from O and OH and released into the gas phase by photodesorption. The reactions are fast enough with respect to the slow dynamical evolution of L1544 that the gas phase water is in equilibrium for the local conditions thoughout the cloud. We explain the paradoxical radiative transfer of the H$_2$O ($1_{10}-1_{01}$) line. Despite discouragingly high optical depth caused by the large Einstein A coefficient, the subcritical excitation in the cold, rarefied H$_2$ causes the line brightness to scale linearly with column density. Thus the water line can provide information on the chemical and dynamical processes in the darkest region in the center of a cold, dense cloud. The inverse P-Cygni profile of the observed water line generally indicates a contracting cloud. This profile is reproduced with a dynamical model of slow contraction from unstable quasi-static hydrodynamic equilibrium (an unstable Bonnor-Ebert sphere).

Long wavelength sonic oscillations are observed or inferred in many of the small, dark molecular clouds, the starless cores, that are the precursors to protostars. The oscillations provide significant internal energy and the time scale for their dissipation may control the rate of star formation in the starless cores. Despite their importance, their origin is unknown. We explore one hypothesis that the oscillations develop as a starless core forms from a filament. We model this process with a numerical hydrodynamic simulation and generate synthetic molecular line observations with a radiative transfer simulation. Comparison with actual observations shows general agreement suggesting the proposed mechanism is viable.

The Multi-Scale Continuum and Line Exploration of W49 (MUSCLE W49) is a comprehensive gas and dust survey of the giant molecular cloud (GMC) of W49A, the most luminous star-formation region in the Milky Way. It covers the entire GMC at different scales and angular resolutions. In this paper we present: 1) an all-configuration SMA mosaic in the 230-GHz band covering the central 3 arcmin (10 pc, known as W49N), with most of the embedded massive stars; and 2) PMO 14m telescope observations in the 90-GHz band, covering the entire GMC with maps up to 35 arcmin in size, or 113 pc. We also make use of archival data from the VLA, JCMT-SCUBA, IRAM 30m, and the CSO BOLOCAM GPS. Our main findings are: 1) The W49 GMC is one of the most massive in the Galaxy, with a total mass ~1.1x10^6 Msun within a radius of 60 pc. Within a radius of 6 pc, the total gas mass is ~2x10^5 Msun. At these scales only 1% of the material is photoionized. The mass reservoir is sufficient to form several young massive clusters (YMCs) as massive as a globular cluster. 2) The mass of the GMC is distributed in a hierarchical network of filaments. At scales <10 pc, a triple, centrally condensed structure peaks toward the ring of HC HII regions in W49N. This structure extends to scales from ~10 to 100 pc. The W49A starburst most likely formed from global gravitational contraction with localized collapse in a "hub-filament" geometry. 3) Currently, feedback from the central YMCs (with a present mass Mcl > 5x10^4 Msun) is still not enough to entirely disrupt the GMC, but further stellar mass growth could be enough to allow radiation pressure to clear the cloud and halt star formation. 4) The resulting stellar content will probably remain as a gravitationally bound massive star cluster, or a small system of bound clusters. (ABRIDGED)

Accretion flows onto massive stars must transfer mass so quickly that they are themselves gravitationally unstable, forming dense clumps and filaments. These density perturbations interact with young massive stars, emitting ionizing radiation, alternately exposing and confining their HII regions. As a result, the HII regions are predicted to flicker in flux density over periods of decades to centuries rather than increasing monotonically in size as predicted by simple Spitzer solutions. We have recently observed the Sgr B2 region at 1.3 cm with the VLA in its three hybrid configurations (DnC, CnB and BnA) at a resolution of 0.25''. These observations were made to compare in detail with matched continuum observations from 1989. At 0.25'' resolution, Sgr B2 contains 41 UC HII regions, 6 of which are hypercompact. The new observations of Sgr B2 allow comparison of relative peak flux densites for the HII regions in Sgr B2 over a 23 year time baseline (1989-2012) in one of the most source-rich massive star forming regions in the Milky Way. The new 1.3 cm continuum images indicate that four of the 41 UC HII regions exhibit significant changes in their peak flux density, with one source (K3) dropping in peak flux density, and the other 3 sources (F10.303, F1 and F3) increasing in peak flux density. The results are consistent with statistical predictions from simulations of high mass star formation, suggesting that they offer a solution to the lifetime problem for ultracompact HII regions.

Many ultracompact HII regions exhibit a cometary morphology in radio continuum emission. In such regions, a young massive star is probably ablating, through its ultraviolet radiation, the molecular cloud clump that spawned it. On one side of the star, the radiation drives an ionization front that stalls in dense molecular gas. On the other side, ionized gas streams outward into the more rarefied environment. This wind is underpressured with respect to the neutral gas. The difference in pressure draws in more cloud material, feeding the wind until the densest molecular gas is dissipated. Recent, time-dependent simulations of massive stars turning on within molecular gas show the system evolving in a direction similar to that just described. Here, we explore a semi-analytic model in which the wind is axisymmetric and has already achieved a steady state. Adoption of this simplified picture allows us to study the dependence of both the wind and its bounding ionization front on the stellar luminosity, the peak molecular density, and the displacement of the star from the center of the clump. For typical parameter values, the wind accelerates transonically to a speed of about 15 km / s, and transports mass outward at a rate of 10^-4 solar masses per year. Stellar radiation pressure acts to steepen the density gradient of the wind.

Array configurations built on a hierarchy of simple elements have excellent properties for cross-correlation imaging interferometers including a smooth distribution of measured Fourier components, high angular resolution, low side lobes, and compact array size. Compared to arrays with a Gaussian distribution of antenna separations, hierarchical arrays (H-arrays) produce beams with higher angular resolution and a tighter concentration of the total power (encircled energy) within a smaller area around the main beam. An attractive feature of H-arrays is their simplicity. The relationships between the Fourier coverage and the array configuration are easy enough to understand that they can be adjusted to achieve different design goals without the need for numerical optimization. H-arrays will be useful for future multi-element interferometers.

Water is a crucial molecule in molecular astrophysics as it controls much of the gas/grain chemistry, including the formation and evolution of more complex organic molecules in ices. Pre-stellar cores provide the original reservoir of material from which future planetary systems are built, but few observational constraints exist on the formation of water and its partitioning between gas and ice in the densest cores. Thanks to the high sensitivity of the Herschel Space Observatory, we report on the first detection of water vapor at high spectral resolution toward a dense cloud on the verge of star formation, the pre-stellar core L1544. The line shows an inverse P-Cygni profile, characteristic of gravitational contraction. To reproduce the observations, water vapor has to be present in the cold and dense central few thousand AU of L1544, where species heavier than Helium are expected to freeze-out onto dust grains, and the ortho:para H2 ratio has to be around 1:1 or larger. The observed amount of water vapor within the core (about 1.5x10^{-6} Msun) can be maintained by Far-UV photons locally produced by the impact of galactic cosmic rays with H2 molecules. Such FUV photons irradiate the icy mantles, liberating water wapor in the core center. Our Herschel data, combined with radiative transfer and chemical/dynamical models, shed light on the interplay between gas and solids in dense interstellar clouds and provide the first measurement of the water vapor abundance profile across the parent cloud of a future solar-type star and its potential planetary system.

Super star clusters --- extremely massive clusters found predominately in starburst environments --- are essential building blocks in the formation of galaxies and thought to dominate star formation in the high-redshift universe. However, the transformation from molecular gas into these ultra-compact star clusters is not well understood. To study this process, we used the Submillimeter Array and the Plateau de Bure Interferometer to obtain high angular resolution (~1.5" or 160 pc) images of the Antennae overlap region in CO(2--1) to search for the molecular progenitors of the super star clusters. We resolve the molecular gas distribution into a large number of clouds, extending the differential cloud mass function down to a 5\sigma completeness limit of 3.8x10^5 M_sun. We identify a distinct break in the mass function around log M_mol/M_sun ~ 6.5, which separates the molecular clouds into two distinct populations. The smaller, less massive clouds reside in more quiescent areas in the region, while the larger, more massive clouds cluster around regions of intense star formation. A broken power-law fit to the mass function yields slopes of \alpha = -1.39+/-0.10 and \alpha = -1.44+/-0.14 for the low- and high-mass cloud population, well-matched to the mass function found for super star clusters in the Antennae galaxies. We find large velocity gradients and velocity dispersions at the locations of intense star formation, suggestive of compressive shocks. It is likely that these environmental factors contribute to the formation of the observed massive molecular clouds and super star clusters in the Antennae galaxies.

Although hydrogen cyanide has become quite a common molecular tracing species for a variety of astrophysical sources, it, however, exhibits dramatic non-LTE behaviour in its hyperfine line structure. Individual hyperfine components can be strongly boosted or suppressed. If these so-called hyperfine line anomalies are present in the HCN rotational spectra towards low or high mass cores, this will affect the interpretation of various physical properties such as the line opacity and excitation temperature in the case of low mass objects and infall velocities in the case of their higher mass counterparts. This is as a consequence of the direct effects that anomalies have on the underlying line shape, be it with the line structural width or through the inferred line strength. This work involves the first observational investigation of these anomalies in two HCN rotational transitions, J=1!0 and J=3!2, towards both low mass starless cores and high mass protostellar objects. The degree of anomaly in these two rotational transitions is considered by computing the ratios of neighboring hyperfine lines in individual spectra. Results indicate some degree of anomaly is present in all cores considered in our survey, the most likely cause being line overlap effects among hyperfine components in higher rotational transitions.

Ultracompact and hypercompact HII regions appear when a star with a mass larger than about 15 solar masses starts to ionize its own environment. Recent observations of time variability in these objects are one of the pieces of evidence that suggest that at least some of them harbor stars that are still accreting from an infalling neutral accretion flow that becomes ionized in its innermost part. We present an analysis of the properties of the HII regions formed in the 3D radiation-hydrodynamic simulations presented by Peters et al. as a function of time. Flickering of the HII regions is a natural outcome of this model. The radio-continuum fluxes of the simulated HII regions, as well as their flux and size variations are in agreement with the available observations. From the simulations, we estimate that a small but non-negligible fraction (~ 10 %) of observed HII regions should have detectable flux variations (larger than 10 %) on timescales of ~ 10 years, with positive variations being more likely to happen than negative variations. A novel result of these simulations is that negative flux changes do happen, in contrast to the simple expectation of ever growing HII regions. We also explore the temporal correlations between properties that are directly observed (flux and size) and other quantities like density and ionization rates.

The recent paper by Heyer et al (2009) indicates that observations of size, linewidth and column density of interstellar clouds do not agree with simple virial equilibrium (VE) as a balance between gravitational and kinetic energies in the sense that the clouds either have too much kinetic energy or too little mass to be bound. This may be explained by violation of VE as suggested by Dobbs et al 2011, by observational underestimation of the masses as suggested by Heyer et al 2009, or by an external pressure acting as an additional confining force as suggested earlier by Heyer et al 2004. The data of Heyer et al. 2009 cannot be explained with a single value for the external pressure, but if different clouds in the sample have different external pressures in the range of Pe/k = E4 to E7 cm-3 K, then most of the clouds could be in pressure virial equilibrium (PVE). In this paper we discuss two consequences of the external pressure. First, we show that the observational data are consistent with the hypothesis (Chie\'ze 1987) that most clouds are at a critical mass for dynamical stability determined solely by the pressure. Above this mass a cloud is unstable to gravitational collapse or fragmentation. Second, we show that the external pressure modifies the well-known size-linewidth relationship first proposed by Larson (1981) so that the proportionality is no longer constant but depends on the external pressure.

In order to investigate whether massive stars form similarly to their low-mass counterparts, we have used the standard envelope plus disc geometry successfully applied to low-mass protostars to model the near-IR to sub-millimetre SED and several mid-IR images of the embedded massive star IRAS20126+4104. We have used a Monte Carlo radiative transfer dust code to model the continuum absorption, emission and scattering through two azimuthally symmetric dust geometries, the first consisting of a rotationally flattened envelope with outflow cavities, and the second which also includes a flared accretion disc. Our results show that the envelope plus disc model reproduces the observed SED and images more accurately than the model without a disc, although the latter model more closely reproduces the morphology of the mid-IR emission within a radius of 1.1" or ~1800au. We have put forward several possible causes of this discontinuity, including inner truncation of the disc due to stellar irradiation, or precession of the outflow cavity. Our best fitting envelope plus disc model has a disc radius of 9200 au. We find that it is unlikely that the outer regions of such a disc would be in hydrostatic or centrifugal equilibrium, however we calculate that the temperatures within the disc would keep it stable to fragmentation.

Context. To further understand the processes involved in the formation of massive stars, we have undertaken a study of the gas dynamics surrounding three massive star forming regions. By observing the large scale structures at high resolution, we are able to determine properties such as driving source, and spatially resolve the bulk dynamical properties of the gas such as infall and outflow. Aims. With high resolution observations, we are able to determine which of the cores in a cluster forming massive stars is responsible for the large scale structures. Methods. We present CO observations of three massive star forming regions with known HII regions and show how the CO traces both infall and outflow. By combining data taken in two SMA configurations with JCMT observations, we are able to see large scale structures at high resolution. Results. We find large (0.26-0.40 pc), massive (2-3 M_sun) and energetic (13-17 \times 10^44 erg) outflows emanating from the edges of two HII regions suggesting they are being powered by the protostar(s) within. We find infall signatures in two of our sources with mass infall rates of order 10-4 M_sun/yr. Conclusions. We suggest that star formation is ongoing in these sources despite the presence of HII regions. We further conclude that the source(s) within a single HII region are responsible for the observed large scale structures; that these large structures are not the net effect of multiple outflows from multiple HII regions and hot cores.

We report a spectral-line survey of the extreme carbon star IRC+10216 carried out between 293.9 and 354.8 GHz with the Submillimeter Array. A total of 442 lines were detected, more than 200 for the first time; 149 are unassigned. Maps at an angular resolution of ~3" were obtained for each line. A substantial new population of narrow lines with an expansion velocity of ~4 km/s (i.e. ~ 30% of the terminal velocity) was detected. Most of these are attributed to rotational transitions within vibrationally excited states, emitted from energy levels above the v=0, J=0 ground state with excitation energy of 1000-3000 K. Emission from these lines appears to be centered on the star with an angular extent of <1". We use multiple transitions detected in several molecules to derive physical conditions in this inner envelope of IRC+10216.

If star formation proceeds by thermal fragmentation and the subsequent gravitational collapse of the individual fragments, how is it possible to form fragments massive enough for O and B stars in a typical star-forming molecular cloud where the Jeans mass is about 1Msun at the typical densities (10^4 cm^-3) and temperatures (10K)? We test the hypothesis that a first generation of low-mass stars may heat the gas enough that subsequent thermal fragmentation results in fragments >=10Msun, sufficient to form B stars. We combine ATCA and SMA observations of the massive star-forming region G8.68-0.37 with radiative transfer modeling to derive the present-day conditions in the region and use this to infer the conditions in the past, at the time of core formation. Assuming the current mass/separation of the observed cores equals the fragmentation Jeans mass/length and the region's average density has not changed, requires the gas temperature to have been 100K at the time of fragmentation. The postulated first-generation of low-mass stars would still be around today, but the number required to heat the cloud exceeds the limits imposed by the observations. Several lines of evidence suggest the observed cores in the region should eventually form O stars yet none have sufficient raw material. Even if feedback may have suppressed fragmentation, it was not sufficient to halt it to this extent. To develop into O stars, the cores must obtain additional mass from outside their observationally defined boundaries. The observations suggest they are currently fed via infall from the very massive reservoir (~1500Msun) of gas in the larger pc scale cloud around the star-forming cores. This suggests that massive stars do not form in the collapse of individual massive fragments, but rather in smaller fragments that themselves continue to gain mass by accretion from larger scales.

Interferometric observations of the W33A massive star-formation region, performed with the Submillimeter Array (SMA) and the Very Large Array (VLA) at resolutions from 5 arcsec (0.1 pc) to 0.5 arcsec (0.01 pc) are presented. Our three main findings are: (1) parsec-scale, filamentary structures of cold molecular gas are detected. Two filaments at different velocities intersect in the zone where the star formation is occurring. This is consistent with triggering of the star-formation activity by the convergence of such filaments, as predicted by numerical simulations of star formation initiated by converging flows. (2) The two dusty cores (MM1 and MM2) at the intersection of the filaments are found to be at different evolutionary stages, and each of them is resolved into multiple condensations. MM1 and MM2 have markedly different temperatures, continuum spectral indices, molecular-line spectra, and masses of both stars and gas. (3) The dynamics of the "hot-core" MM1 indicates the presence of a rotating disk in its center (MM1-Main) around a faint free-free source. The stellar mass is estimated to be approximately 10 Msun. A massive molecular outflow is observed along the rotation axis of the disk.

This is the first paper of our series of high resolution (1") studies of the massive star forming region G10.6--0.4. We present the emission line observations of the hot core type tracers (O$^{13}$CS, OCS, SO$_{2}$) with $\sim$0$"$.5 resolution. By comparing the results to the high--resolution NH$_{3}$ absorption line observation, we confirm for the first time the rotationally flattened hot toroid in the central $<$0.1 pc region, which has a rotational axis perpendicular to its geometrical major axis. In addition, we present the observations of NH$_{3}$, $^{13}$CS, and CH$_{3}$CN with $\sim$1$"$ resolution, and follow the dynamics of the molecular accretion flow from the 0.3 pc radius to the inner 0.03 pc radius. With reference to the rotational axis of the hot toroid, we measure the rotational velocity from the molecular emission in the region. The results are consistent with an envelope with a rapid decrease of the specific angular momentum from the outer to the inner region. These new results improve the current understanding of the molecular accretion flow in an ultracompact (UC) H\textsc{ii} region created by the embedded O-type cluster.

We present 12mm Mopra observations of dense molecular gas towards the W28 supernova remnant (SNR) field. The focus is on the dense molecular gas towards the TeV gamma-ray sources detected by the H.E.S.S. telescopes, which likely trace the cosmic-rays from W28 and possibly other sources in the region. Using the NH3 inversion transitions we reveal several dense cores inside the molecular clouds, the majority of which coincide with high-mass star formation and HII regions, including the energetic ultra-compact HII region G5.89-0.39. A key exception to this is the cloud north east of W28, which is well-known to be disrupted as evidenced by clusters of 1720MHz OH masers and broad CO line emission. Here we detect broad NH3, up to the (9,9) transition, with linewidths up to 16 km/s. This broad NH3 emission spatially matches well with the TeV source HESS J1801-233 and CO emission, and its velocity dispersion distribution suggests external disruption from the W28 SNR direction. Other lines are detected, such as HC3N and HC5N, H2O masers, and many radio recombination lines, all of which are primarily found towards the southern high-mass star formation regions. These observations provide a new view onto the internal structures and dynamics of the dense molecular gas towards the W28 SNR field, and in tandem with future higher resolution TeV gamma-ray observations will offer the chance to probe the transport of cosmic-rays into molecular clouds.

Massive stars influence the surrounding universe far out of proportion to their numbers through ionizing radiation, supernova explosions, and heavy element production. Their formation requires the collapse of massive interstellar gas clouds with very high accretion rates. We discuss results from the first three-dimensional simulations of the gravitational collapse of a massive, rotating molecular cloud core that include heating by both non-ionizing and ionizing radiation. Local gravitational instabilities in the accretion flow lead to the build-up of a small cluster of stars. These lower-mass companions subsequently compete with the high-mass star for the same common gas reservoir and limit its overall mass growth. This process is called fragmentation-induced starvation, and explains why massive stars are usually found as members of high-order stellar systems. These simulations also show that the HII regions forming around massive stars are initially trapped by the infalling gas, but soon begin to fluctuate rapidly. Over time, the same ultracompact HII region can expand anisotropically, contract again, and take on any of the observed morphological classes. The total lifetime of HII regions is given by the global accretion timescale, rather than their short internal sound-crossing time. This solves the so-called lifetime problem of ultracompact HII region. We conclude that the the most significant differences between the formation of low-mass and high-mass stars are all explained as the result of rapid accretion within a dense, gravitationally unstable flow.

We present new images and photometry of the massive star forming complex Cygnus X obtained with IRAC and MIPS on the Spitzer Space Telescope. A combination of IRAC, MIPS, UKIDSS, and 2MASS data are used to identify and classify young stellar objects. Of the 8,231 sources detected exhibiting infrared excess in Cygnus X North, 670 are classified as Class I and 7,249 are classified as Class II. Using spectra from the FAST spectrograph at the Fred L. Whipple Observatory and Hectospec on the MMT, we spectrally typed 536 sources in the Cygnus X complex to identify the massive stars. We find that YSOs tend to be grouped in the neighborhoods of massive B stars (spectral types B0 to B9). We present a minimal spanning tree analysis of clusters in two regions in Cygnus X North. The fraction of infrared excess sources that belong to clusters with >10 members is found to be 50--70%. Most Class II objects lie in dense clusters within blown out HII regions, while Class I sources tend to reside in more filamentary structures along the bright-rimmed clouds, indicating possible triggered star formation.

SWAS and Odin provided stringent upper limits on the gas phase water abundance of dark clouds (x(H2O) < 7x10^-9). We investigate the chemistry of water vapor in starless cores beyond the previous upper limits using the highly improved angular resolution and sensitivity of Herschel and measure the abundance of water vapor during evolutionary stages just preceding star formation. High spectral resolution observations of the fundamental ortho water (o-H2O) transition (557 GHz) were carried out with Herschel HIFI toward two starless cores: B68, a Bok globule, and L1544, a prestellar core embedded in the Taurus molecular cloud complex. The rms in the brightness temperature measured for the B68 and L1544 spectra is 2.0 and 2.2 mK, respectively, in a velocity bin of 0.59 km s^-1. The continuum level is 3.5+/-0.2 mK in B68 and 11.4+/-0.4 mK in L1544. No significant feature is detected in B68 and the 3 sigma upper limit is consistent with a column density of o-H2O N(o-H2O) < 2.5x10^13 cm^-2, or a fractional abundance x(o-H2O) < 1.3x10^-9, more than an order of magnitude lower than the SWAS upper limit on this source. The L1544 spectrum shows an absorption feature at a 5 sigma level from which we obtain the first value of the o-H2O column density ever measured in dark clouds: N(o-H2O) = (8+/-4)x10^12 cm^-2. The corresponding fractional abundance is x(o-H2O) ~ 5x10^-9 at radii > 7000 AU and ~2x10^-10 toward the center. The radiative transfer analysis shows that this is consistent with a x(o-H2O) profile peaking at ~10^-8, 0.1 pc away from the core center, where both freeze-out and photodissociation are negligible. Herschel has provided the first measurement of water vapor in dark regions. Prestellar cores such as L1544 (with their high central densities, strong continuum, and large envelopes) are very promising tools to finally shed light on the solid/vapor balance of water in molecular clouds.

For a number of starless cores, self-absorbed molecular line and column density observations have implied the presence of large-amplitude oscillations. We examine the consequences of these oscillations on the evolution of the cores and the interpretation of their observations. We find that the pulsation energy helps support the cores and that the dissipation of this energy can lead toward instability and star formation. In this picture, the core lifetimes are limited by the pulsation decay timescales, dominated by non-linear mode-mode coupling, and on the order of ~few x 10^5--10^6 yr. Notably, this is similar to what is required to explain the relatively low rate of conversion of cores into stars. For cores with large-amplitude oscillations, dust continuum observations may appear asymmetric or irregular. As a consequence, some of the cores that would be classified as supercritical may be dynamically stable when oscillations are taken into account. Thus, our investigation motivates a simple hydrodynamic picture, capable of reproducing many of the features of the progenitors of stars without the inclusion of additional physical processes, such as large-scale magnetic fields.

In this paper we discuss two approximate methods previously suggested for modeling hyperfine spectral line emission for molecules whose collisional transitions rates between hyperfine levels are unknown. Hyperfine structure is seen in the rotational spectra of many commonly observed molecules such as HCN, HNC, NH3, N2H+, and C17O. The intensities of these spectral lines can be modeled by numerical techniques such as Lambda-iteration that alternately solve the equations of statistical equilibrium and the equation of radiative transfer. However, these calculations require knowledge of both the radiative and collisional rates for all transitions. For most commonly observed radio frequency spectral lines, only the net collisional rates between rotational levels are known. For such cases, two approximate methods have been suggested. The first method, hyperfine statistical equilibrium (HSE), distributes the hyperfine level populations according to their statistical weight, but allows the population of the rotational states to depart from local thermodynamic equilibrium (LTE). The second method, the proportional method approximates the collision rates between the hyperfine levels as fractions of the net rotational rate apportioned according to the statistical degeneracy of the final hyperfine levels. The second method is able to model non-LTE hyperfine emission. We compare simulations of N2H+ hyperfine lines made with approximate and more exact rates and find that satisfactory results are obtained.

We address the question of whether the formation of high-mass stars is similar to or differs from that of solar-mass stars through new molecular line observations and modeling of the accretion flow around the massive protostar IRAS20126+4104. We combine new observations of NH3(1,1) and (2,2) made at the Very Large Array, new observations of CHCN(13-12) made at the Submillimeter Array, previous VLA observations of NH(3,3), NH(4,4), and previous Plateau de Bure observations of C34S(2-1), C34S(5-4), and CHCN(12-11) to obtain a data set of molecular lines covering 15 to 419 K in excitation energy. We compare these observations against simulated molecular line spectra predicted from a model for high-mass star formation based on a scaled-up version of the standard disk-envelope paradigm developed for accretion flows around low-mass stars. We find that in accord with the standard paradigm, the observations require both a warm, dense, rapidly-rotating disk and a cold, diffuse infalling envelope. This study suggests that accretion processes around 10 M stars are similar to those of solar mass stars.

We present observations of newly discovered 24 micron circumstellar structures detected with the Multiband Imaging Photometer for Spitzer (MIPS) around three evolved stars in the Cygnus-X star forming region. One of the objects, BD+43 3710, has a bipolar nebula, possibly due to an outflow or a torus of material. A second, HBHA 4202-22, a Wolf-Rayet candidate, shows a circular shell of 24 micron emission suggestive of either a limb-brightened shell or disk seen face-on. No diffuse emission was detected around either of these two objects in the Spitzer 3.6-8 micron Infrared Array Camera (IRAC) bands. The third object is the luminous blue variable candidate G79.29+0.46. We resolved the previously known inner ring in all four IRAC bands. The 24 micron emission from the inner ring extends ~1.2 arcmin beyond the shorter wavelength emission, well beyond what can be attributed to the difference in resolutions between MIPS and IRAC. Additionally, we have discovered an outer ring of 24 micron emission, possibly due to an earlier episode of mass loss. For the two shell stars, we present the results of radiative transfer models, constraining the stellar and dust shell parameters. The shells are composed of amorphous carbon grains, plus polycyclic aromatic hydrocarbons in the case of G79.29+0.46. Both G79.29+0.46 and HBHA 4202-22 lie behind the main Cygnus-X cloud. Although G79.29+0.46 may simply be on the far side of the cloud, HBHA 4202-22 is unrelated to the Cygnus-X star formation region.

We further develop the model of molecular cloud fragmentation introduced in Field, Blackman and Keto (2007; FBK). We show that external pressure acting on fragments establishes a scale-dependent critical mass. Fragments with masses less than the critical value are confined largely by pressure, while those with masses greater than or equal to the critical value collapse under self gravitation. Both types of fragments are commonly observed. Without specifying the source of the external pressure, and without assuming any other scaling relations, we predict the power - law index in the relation between the rms velocity of supersonic motions and the size of fragments . We then investigate the possibility that the external pressure is due to the kinetic energy of H atoms released by photodissociation of hydrogen molecules in the fragment. This can account approximately for the observed values of external pressure and two additional observations: the value of the scaling coefficient in the power law mentioned above, and the observation of outflowing atomic hydrogen around molecular clouds. A further prediction is HI at fragment edges with column densities of order 1E20 per sq. cm and velocities of a few km/s that should be detectable with high resolution 21 cm observations. Finally, we predict the magnitude of the coefficient of dissipation in the observed supersonic flows.

We describe the first three-dimensional simulation of the gravitational collapse of a massive, rotating molecular cloud that includes heating by both non-ionizing and ionizing radiation. We find that as the first protostars gain sufficient mass to ionize the accretion flow, their H II regions are initially gravitationally trapped, but soon begin to rapidly fluctuate between trapped and extended states, in agreement with observations. Over time, the same ultracompact H II region can expand anisotropically, contract again, and take on any of the observed morphological classes. In their extended phases, expanding H II regions drive bipolar neutral outflows characteristic of high-mass star formation. The total lifetime of H II regions is given by the global accretion timescale, rather than their short internal sound-crossing time. The pressure of the hot, ionized gas does not terminate accretion. Instead the final stellar mass is set by fragmentation-induced starvation. Local gravitational instabilities in the accretion flow lead to the build-up of a small cluster of stars, all with relatively high masses due to heating from accretion radiation. These companions subsequently compete with the initial high-mass star for the same common gas reservoir and limit its mass growth. Our findings show that the most significant differences between the formation of low-mass and high-mass stars are all explained as the result of rapid accretion within a dense, gravitationally unstable, ionized flow.

In previous studies we identified two classes of starless cores, thermally subcritical and supercritical, distinguished by different dynamical behavior and internal structure. Here we study the evolution of the dynamically-unstable, thermally-supercritical cores by means of a numerical hydrodynamic simulation that includes radiative equilibrium and simple molecular chemistry. We use our non-LTE radiative transfer code MOLLIE to predict observable CO and N2H+ line spectra, including the non-LTE hyperfine ratios of N2H+, during the contraction. These are compared against observations of the starless core L1544.

Spectral line and continuum observations of the ionized and molecular gas in G20.08-0.14 N explore the dynamics of accretion over a range of spatial scales in this massive star-forming region. Very Large Array observations of NH_3 at 4'' angular resolution show a large-scale (0.5 pc) molecular accretion flow around and into a star cluster with three small, bright HII regions. Higher resolution (0.4'') observations with the Submillimeter Array in hot core molecules (CH_3CN, OCS, and SO_2) and the VLA in NH_3, show that the two brightest and smallest HII regions are themselves surrounded by smaller scale (0.05 pc) accretion flows. The axes of rotation of the large- and small-scale flows are aligned, and the timescale for the contraction of the cloud is short enough, 0.1 Myr, for the large-scale accretion flow to deliver significant mass to the smaller scales within the star formation timescale. The flow structure appears to be continuous and hierarchical from larger to smaller scales. Millimeter radio recombination line (RRL) observations at 0.4" angular resolution indicate rotation and outflow of the ionized gas within the brightest HII region (A). The broad recombination lines and a continuum spectral energy distribution (SED) that rises continuously from cm to mm wavelengths, are both characteristic of the class of HII regions known as "broad recombination line objects". The SED indicates a density gradient inside this HII region, and the RRLs suggest supersonic flows. These observations are consistent with photoevaporation of the inner part of the rotationally flattened molecular accretion flow. We also report the serendipitous detection of a new NH_3 (3,3) maser.

We present molecular line and 1.4 mm continuum observations towards five massive star forming regions at arcsecond resolution using the Submillimeter Array (SMA). We find that the warm molecular gas surrounding each HII region (as traced by SO_2 and OCS) appears to be undergoing bulk rotation. From the molecular line emission and thermal component of the continuum emission, we independently derived gas masses for each region which are consistent with each other. From the free-free component of the continuum emission we estimate the minimum stellar mass required to power the HII region and find that this mass, when added to the derived gas mass, is a significant fraction of the dynamical mass for that region.

An example of a cold massive core, JCMT 18354-0649S, a possible high mass analogue to a low mass star forming core is studied. Line and continuum observations from JCMT, Mopra Telescope and Spitzer are presented and modelled in detail using a 3D molecular line radiative transfer code. In almost every way JCMT 18354-0649S is a scaled-up version of a typical low mass core with similar temperatures, chemical abundances and densities. The difference is that both the infall velocity and the turbulent width of the line profiles are an order of magnitude larger. While the higher infall velocity is expected due to the large mass of JCMT 18354-0649S, we suggest that the dissipation of this highly supersonic turbulence may lead to the creation of dense clumps of gas that surround the high mass core.

We examine an Infrared Dark Cloud (IRDC) at high spatial resolution as a means to study rotation, outflow, and infall at the onset of massive star formation. Submillimeter Array observations combined with IRAM 30 meter data in 12CO(2--1) reveal the outflow orientation in the IRDC 18223-3 region, and PdBI 3 mm observations confirm this orientation in other molecular species. The implication of the outflow's presence is that an accretion disk is feeding it, so using high density tracers such as C18O, N2H+, and CH3OH, we looked for indications of a velocity gradient perpendicular to the outflow direction. Surprisingly, this gradient turns out to be most apparent in CH3OH. The large size (28,000 AU) of the flattened rotating object detected indicates that this velocity gradient cannot be due solely to a disk, but rather from inward spiraling gas within which a Keplerian disk likely exists. From the outflow parameters, we derive properties of the source such as an outflow dynamical age of ~37,000 years, outflow mass of ~13 M_sun, and outflow energy of ~1.7 x 10^46 erg. While the outflow mass and energy are clearly consistent with a high-mass star forming region, the outflow dynamical age indicates a slightly more evolved evolutionary stage than previous spectral energy distribution (SED) modeling indicates. The calculated outflow properties reveal that this is truly a massive star in the making. We also present a model of the observed methanol velocity gradient. The rotational signatures can be modeled via rotationally infalling gas. These data present evidence for one of the youngest known outflow/infall/disk systems in massive star formation. A tentative evolutionary picture for massive disks is discussed.

We present Australia Telescope Compact Array observations towards 6 massive star formation regions which, from their strong 24 GHz continuum emission but no compact 8 GHz continuum emission, appeared good candidates for hyper-compact HII regions. However, the properties of the ionised gas derived from the 19 to 93 GHz continuum emission and H70 alpha + H57 alpha radio recombination line data show the majority of these sources are, in fact, regions of spatially-extended, optically-thin free-free emission. These extended sources were missed in the previous 8 GHz observations due to a combination of spatial-filtering, poor surface brightness sensitivity and primary beam attenuation. We consider the implications that a significant number of these extended HII regions may have been missed by previous surveys of massive star formation regions. If the original sample of 21 sources is representative of the population as a whole, the fact that 6 contain previously undetected extended free-free emission suggests a large number of regions have been mis-classified. Rather than being very young objects prior to UCHII region formation, they are, in fact, associated with extended HII regions and thus significantly older. In addition, inadvertently ignoring a potentially substantial flux contribution (up to ~0.5Jy) from free-free emission has implications for dust masses derived from sub-mm flux densities. The large spatial scales probed by single-dish telescopes, which do not suffer from spatial filtering, are particularly susceptible and dust masses may be overestimated by up to a factor of ~2.

We describe a model for the thermal and dynamical equilibrium of starless cores that includes the radiative transfer of the gas and dust and simple CO chemistry. The model shows that the structure and behavior of the cores is significantly different depending on whether the central density is either above or below about 10^5 cm-3. This density is significant as the critical density for gas cooling by gas-dust collisions and also as the critical density for dynamical stability, given the typical properties of the starless cores. The starless cores thus divide into two classes that we refer to as thermally super-critical and thermally sub-critical.This two-class distinction allows an improved interpretation of the different observational data of starless cores within a single model.

In a previous paper we demonstrated that non-radial hydrodynamic oscillations of a thermally-supported (Bonnor-Ebert) sphere embedded in a low-density, high-temperature medium persist for many periods. The predicted column density variations and molecular spectral line profiles are similar to those observed in the Bok globule B68 suggesting that the motions in some starless cores may be oscillating perturbations on a thermally supported equilibrium structure. Such oscillations can produce molecular line maps which mimic rotation, collapse or expansion, and thus could make determining the dynamical state from such observations alone difficult. However, while B68 is embedded in a very hot, low-density medium, many starless cores are not, having interior/exterior density contrasts closer to unity. In this paper we investigate the oscillation damping rate as a function of the exterior density. For concreteness we use the same interior model employed in Broderick et al. (2007), with varying models for the exterior gas. We also develop a simple analytical formalism, based upon the linear perturbation analysis of the oscillations, which predicts the contribution to the damping rates due to the excitation of sound waves in the external medium. We find that the damping rate of oscillations on globules in dense molecular environments is always many periods, corresponding to hundreds of thousands of years, and persisting over the inferred lifetimes of the globules.

Previous observations show that the hypercompact HII region W51e2 is surrounded by a massive molecular accretion flow centered on the HII region. New observations of the H53alpha radio recombination line made with the VLA at 0.45 arc second angular resolution show a velocity gradient in the ionized gas within the HII region of > 500 kms-1 pc-1 comparable to the velocity gradient seen in the molecular accretion flow. New CO line observations made with the SMA at arc second angular resolution detect a molecular bipolar outflow immediately around the W51e2 HII region and extending along the axis of rotation of the molecular flow. These observations are consistent with an evolutionary phase for high mass star formation in which a newly formed massive star first begins to ionize its surroundings including its own accretion flow.

Over a timescale of a few years, an observed change in the optically thick radio continuum flux can indicate whether an unresolved H II region around a newly formed massive star is changing in size. In this Letter we report on a study of archival VLA observations of the hypercompact H II region G24.78+0.08 A1 that shows a decrease of ~ 45 % in the 6-cm flux over a five year period. Such a decrease indicates a contraction of ~ 25 % in the ionized radius and could be caused by an increase in the ionized gas density if the size of the H II region is determined by a balance between photoionization and recombination. This finding is not compatible with continuous expansion of the H II region after the end of accretion onto the ionizing star, but is consistent with the hypothesis of gravitational trapping and ionized accretion flows if the mass-accretion rate is not steady.

CO isotopes are able to probe the different components in protostellar clouds. These components, core, envelope and outflow have distinct physical conditions and sometimes more than one component contributes to the observed line profile. In this study we determine how CO isotope abundances are altered by the physical conditions in the different components. We use a 3D molecular line transport code to simulate the emission of four CO isotopomers, 12CO J=2-1, 13CO J=2-1, C18O J=2-1 and C17O J=2-1 from the Class 0/1 object L483, which contains a cold quiescent core, an infalling envelope and a clear outflow. Our models replicate JCMT (James Clerk Maxwell Telescope) line observations with the inclusion of freeze-out, a density profile and infall. Our model profiles of 12CO and 13CO have a large linewidth due to a high velocity jet. These profiles replicate the process of more abundant material being susceptible to a jet. C18O and C17O do not display such a large linewidth as they trace denser quiescent material deep in the cloud.

Velocity shifts and differential broadening of radio recombination lines are used to estimate the densities and velocities of the ionized gas in several hypercompact and ultracompact HII regions. These small HII regions are thought to be at their earliest evolutionary phase and associated with the youngest massive stars. The observations suggest that these HII regions are characterized by high densities, supersonic flows and steep density gradients, consistent with accretion and outflows that would be associated with the formation of massive stars.

We present a model in which the supersonic motions observed in molecular clouds are driven by gravitational energy released as large structures fragment into smaller ones. The fragmentation process begins in large molecular clouds, and continues down to fragments of a critical mass, at which gravitational confinement may be replaced by pressure confinement. The power laws that describe the scaling of density and mass, and number spectra and mass spectra of the fragments are given in terms of that of the observed velocity dispersion of the fragments. The results agree reasonably well with observations over the range from several to about a third of a million solar masses.

In a previous paper, we modeled the oscillations of a thermally-supported (Bonnor-Ebert) sphere as non-radial, linear perturbations following a standard analysis developed for stellar pulsations. The predicted column density variations and molecular spectral line profiles are similar to those observed in the Bok globule B68 suggesting that the motions in some starless cores may be oscillating perturbations on a thermally supported equilibrium structure. However, the linear analysis is unable to address several questions, among them the stability, and lifetime of the perturbations. In this paper we simulate the oscillations using a three-dimensional numerical hydrodynamic code. We find that the oscillations are damped predominantly by non-linear mode-coupling, and the damping time scale is typically many oscillation periods, corresponding to a few million years, and persisting over the inferred lifetime of gobules.

The hypothesis that massive stars form by accretion can be investigated by simple analytical calculations that describe the effect that the formation of a massive star has on its own accretion flow. Within a simple accretion model that includes angular momentum, that of gas flow on ballistic trajectories around a star, the increasing ionization of a massive star growing by accretion produces a three-stage evolutionary sequence. The ionization first forms a small quasi-spherical HII region gravitationally trapped within the accretion flow. At this stage the flow of ionized gas is entirely inward. As the ionization increases, the HII region transitions to a bipolar morphology in which the inflow is replaced by outflow within a narrow range of angle with about the bipolar axis. At higher rates of ionization, the opening angle of the outflow region progressively increases. Eventually, in the third stage, the accretion is confined to a thin region about an equatorial disk. Throughout this early evolution, the HII region is of hypercompact to ultracompact size depending on the mass of the enclosed star or stars. These small HII regions whose dynamics are dominated by stellar gravitation and accretion are different than compact and larger HII regions whose dynamics are dominated by the thermal pressure of the ionized gas.

New observations of CO (J=1->0) line emission from M33, using the 25 element BEARS focal plane array at the Nobeyama Radio Observatory 45-m telescope, in conjunction with existing maps from the BIMA interferometer and the FCRAO 14-m telescope, give the highest resolution (13'') and most sensitive (RMS ~ 60 mK) maps to date of the distribution of molecular gas in the central 5.5 kpc of the galaxy. A new catalog of giant molecular clouds (GMCs) has a completeness limit of 1.3 X 10^5 M_sun. The fraction of molecular gas found in GMCs is a strong function of radius in the galaxy, declining from 60% in the center to 20% at galactocentric radius R_gal ~ 4 kpc. Beyond that radius, GMCs are nearly absent, although molecular gas exists. Most (90%) of the emission from low mass clouds is found within 100 pc projected separation of a GMC. In an annulus 2.1< R_gal <4.1 kpc, GMC masses follow a power law distribution with index -2.1. Inside that radius, the mass distribution is truncated, and clouds more massive than 8 X 10^5 M_sun are absent. The cloud mass distribution shows no significant difference in the grand design spiral arms versus the interarm region. The CO surface brightness ratio for the arm to interarm regions is 1.5, typical of other flocculent galaxies.

The fine-structure lines of singly ([Ne II] 12.8 micron) and doubly ([Ne III] 15.6 micron) ionized neon are among the most prominent features in the mid-infrared spectra of star-forming regions, and have the potential to be a powerful new indicator of the star formation rate in galaxies. Using a sample of star-forming galaxies with measurements of the fine-structure lines available from the literature, we show that the sum of the [Ne II] and [Ne III] luminosities obeys a tight, linear correlation with the total infrared luminosity, over 5 orders of magnitude in luminosity. We discuss the formation of the lines and their relation with the Lyman continuum luminosity. A simple calibration between star formation rate and the [Ne II]+[Ne III] luminosity is presented.

New molecular line observations of the Bok globule Barnard 68 in HCO+ irrefutably confirm the complex pattern of red and blue asymmetric line profiles seen across the face of the cloud in previous observations of CS. The new observations thus strengthen the previous interpretation that Barnard 68 is undergoing peculiar oscillations. Furthermore, the physical chemistry of B68 indicates that the object is much older than the sound crossing time and is therefore long-lived. A model is presented for the globule in which a modest external pressure perturbation is shown to lead to oscillations about a stable equilibrium configuration. Such oscillations may be present in other stable starless cores as manifested by a similar signature of inward and outward motions.

We review the phenomenon of ultra-compact H II regions (UCHIIs) as a key phase in the early lives of massive stars. This most visible manifestation of massive star formation begins when the Lyman continuum output from the massive young stellar object becomes sufficient to ionize the surroundings from which it was born. Knowledge of this environment is gained through an understanding of the morphologies of UCHII regions and we examine the latest developments in deep radio and mid-IR imaging. SPITZER data from the GLIMPSE survey are an important new resource in which PAH emission and the ionizing stars can be seen. We review the role played by strong stellar winds from the central stars in sweeping out central cavities and causing the limb-brightened appearance. A range of evidence from velocity structure, proper motions, the molecular environment and recent hydrodynamical modeling indicates that cometary UCHII regions require a combination of champagne flow and bow shock motion. Finally, we discuss the class of hyper-compact H II regions or broad recombination line objects. They are likely to mark the transition soon after the breakout of the Lyman continuum radiation from the young star. Models for these objects are presented, including photo-evaporating disks and ionized accretion flows that are gravitationally trapped. Evolutionary scenarios tracing young massive stars passage through these ionized phases are discussed.

If the split, asymmetric molecular spectral line profiles that are seen in many starless cores are interpreted as indicative of global collapse or expansion of the core then one possible implication is that most starless cores have short lifetimes on the order of the collapse or sound crossing time scale. An alternative interpretation of the line profiles as indicative of perturbations on an underlying equilibrium structure leads to the opposite implication, that many cores have long lifetimes. While evidence suggests that some cores are collapsing on a free-fall time scale, we show that observations of some other starless cores can be reproduced by a model of non-radial oscillations about the equilibrium configuration of a pressure-bounded, thermally-supported sphere (Bonnor-Ebert sphere). We model the oscillations as linear perturbations following a standard analysis developed for stellar pulsations and compare the column densities and molecular spectral line profiles predicted from a particular model to observations of the Bok globule B68.

We report on the identification of 54 embedded clusters around 217 massive protostellar candidates of which 34 clusters are new detections. The embedded clusters are identified as stellar surface density enhancements in the 2 $\mu$m All Sky Survey (2MASS) data. Because the clusters are all associated with massive stars in their earliest evolutionary stage, the clusters should also be in an early stage of evolution. Thus the properties of these clusters should reflect properties associated with their formation rather than their evolution. For each cluster, we estimate the mass, the morphological type, the photometry and extinction. The clusters in our study, by their association with massive protostars and massive outflows, reinstate the notion that massive stars begin to form after the first generation of low mass stars have completed their accretion phase. Further, the observed high gas densities and accretion rates at the centers of these clusters is consistent with the hypothesis that high mass stars form by continuing accretion onto low mass stars.

Observations of the H66a recombination line from the ionized gas in the cluster of newly formed massive stars, G10.6-0.4, show that most of the continuum emission derives from the dense gas in an ionized accretion flow that forms an ionized disk or torus around a group of stars in the center of the cluster. The inward motion observed in the accretion flow suggests that despite the equivalent luminosity and ionizing radiation of several O stars, neither radiation pressure nor thermal pressure has reversed the accretion flow. The observations indicate why the radiation pressure of the stars and the thermal pressure of the HII region are not effective in reversing the accretion flow. The observed rate of the accretion flow, 0.001 solar masses/yr, is sufficient to form massive stars within the time scale imposed by their short main sequence lifetimes. A simple model of disk accretion relates quenched HII regions, trapped hypercompact HII regions, and photo-evaporating disks in an evolutionary sequence.

Submillimeter Array observations of Orion-KL at 1'' resolution in the 440mu/690GHz band reveal new insights about the continuum and line emission of the region. The 440mu continuum flux density measurement from source I allows us to differentiate among the various proposed physical models: Source I can be well modeled by a ``normal'' protostellar SED consisting of a proton-electron free-free emission component at low frequencies and a strong dust component in the submillimeter bands. Furthermore, we find that the protostellar object SMA1 is clearly distinct from the hot core. The non-detection of SMA1 at cm and infrared wavelengths suggests that it may be one of the youngest sources in the entire Orion-KL region. The molecular line maps show emission mainly from the sources I, SMA1 and the hot core peak position. An analysis of the CH$3CN(37_K-36_K) K-ladder (K=0...3) indicates a warm gas component of the order 600+-200K. In addition, we detect a large fraction (~58%) of unidentified lines and discuss the difficulties of line identifications at these frequencies.

The central 2x1 kpc of the starburst galaxy NGC 253 has been imaged using the Submillimeter Array at a 60 pc resolution in the J=2-1 transitions of 12CO, 13CO, and C18O as well as in the 1.3 mm continuum. Molecular gas and dust are mainly in the circumnuclear disk of ~500 pc radius, with warm (~40 K) and high area-filling factor gas in its central part. Two gas shells or cavities have been discovered in the circumnuclear disk. They have ~100 pc diameters and have large velocity widths of 80-100 km/s, suggestive of expansion at ~50 km/s. Modeled as an expanding bubble, each shell has an age of ~0.5 Myr and needed kinetic energy of ~1E46 J as well as mean mechanical luminosity of ~1E33 W for its formation. The large energy allows each to be called a superbubble. A ~10^6 Msun super star cluster can provide the luminosity, and could be a building block of the nuclear starburst in NGC 253. Alternatively, a hypernova can also be the main source of energy for each superbubble, not only because it can provide the mechanical energy and luminosity but also because the estimated rate of superbubble formation and that of hypernova explosions are comparable. Our observations indicate that the circumnuclear molecular disk harboring the starburst is highly disturbed on 100 pc or smaller scales, presumably by individual young clusters and stellar explosions, in addition to globally disturbed in the form of the well-known superwind.

We observed the nearby starburst galaxy M82 in CO in the higher frequency (2--1) transition to achieve an angular resolution below 1 arc second or 17 pc at the target. We resolved the molecular gas into a large number of compact clouds, with masses ranging from about 2x10^3 to 2x10^6 solar masses. The mass spectrum scales as N(M) ~ M^-1.5, similar to the mass spectra of young massive star clusters suggesting that individual molecular clouds are transformed in the starburst into individual star clusters. The larger clouds are surrounded by supernovae and HII regions suggesting that star formation proceeds from the outside of the clouds and progresses inward consistent with triggering by a sudden increase in external pressure. The clouds with internal star formation have velocity gradients and inverse P-Cygni spectral line profiles indicating inward motions of 35 kms consistent with shock driven compression. Diffuse free-free radio emission and X-ray emission around the clouds provides evidence for superheated ionized gas sufficient to drive the compression. Clouds with spectral lines indicating expansion show little internal star formation suggesting that the dynamics precedes and is responsible for the star formation rather than the inverse. M82 is known to be in interaction with neighboring M81. The overall picture is consistent with the formation of massive star clusters from individual giant molecular clouds crushed by a sudden galactic scale increase in external pressure generated by the changing dynamics that result from a near-collision with a neighboring galaxy. Present day globular clusters may have formed in a similar fashion in primordial galaxies.

We test the hypothesis that the starless cores may be gravitationally bound clouds supported largely by thermal pressure by comparing observed molecular line spectra to theoretical spectra produced by a simulation that includes hydrodynamics, radiative cooling, variable molecular abundance, and radiative transfer in a simple one-dimensional model. The results suggest that the starless cores can be divided into two categories: stable starless cores that are in approximate equilibrium and will not evolve to form protostars, and unstable pre-stellar cores that are proceeding toward gravitational collapse and the formation of protostars. The starless cores might be formed from the interstellar medium as objects at the lower end of the inertial cascade of interstellar turbulence. Additionally, we identify a thermal instability in the starless cores. Under par ticular conditions of density and mass, a core may be unstable to expansion if the density is just above the critical density for the collisional coupling of the gas and dust so that as the core expands the gas-dust coupling that cools the gas is reduced and the gas warms, further driving the expansion.

We present the first submm (865mum) imaging spectral line survey at one arcsecond resolution conducted with the Submillimeter Array toward Orion-KL. Within the 2x2 GHz bandpasses (lower and upper sidebands, 337.2-339.2GHz and 347.2-349.2GHz), we find about 145 spectral lines from 13 species, 6 isotopologues, and 5 vibrational excited states. Most nitrogen-bearing molecules are strong toward the hot core, whereas the oxygen-bearing molecules peak toward the south-west in the so-called compact ridge. Imaging of spectral lines is shown to be an additional tool to improve the identifications of molecular lines. Arcsecond spatial resolution allows us to distinguish the molecular line emission of the sources I and n from that of the hot core. The only molecular species detected strongly toward source I is SiO, delineating mainly the collimated north-east south-west low-velocity outflow. The two positions close to source I, which have previously been reported to show maser emission in the v=0 28SiO(1-0) and (2-1) lines, show no detectable maser emission in the v=0 28SiO(8-7) line at our spatial resolution. SiO is weak toward source n, and thus source n may not currently be driving a molecular outflow. CH$_3$OH is the molecule with the highest number of identified lines (46) in this spectral window. This ``line forest'' allows us to estimate temperatures in the region, and we find temperatures between 50 and 350K, with the peak temperatures occurring toward the hot core. The detection of strong vibrational excited line emission from the submm continuum peak SMA1 supports the interpretation that the source SMA1 is likely of protostellar nature.

We present new observations of the ultra-compact HII region G28.20-0.05 in 23 GHz continuum, and the NH3(1,1), NH3(2.2), and NH3(3,3) lines. To explain the complicated kinematics of the molecular gas, we propose a model consisting of two components. One component is an infalling, equatorial torus of molecular gas, whose dense central region has been ionized to form the ultra-compact HII region. The second component is a larger expanding molecular shell driven by some type of wide-angle outflow or wind. We estimate that the infall component includes more than 18 Msun of molecular gas. We calculate the central mass to be 79 Msun, probably comprising more than one star. The arrangement of the molecular material suggests a connection to the other disk-like structures seen around massive young stars. The central star in this case is more massive, and the whole region may be more evolved than other similar objects such as IRAS 20126+4104, IRAS 18089-1732, G192.16-3.84, and AFGL5142.

We present high resolution (0.''12 x 0.''079) observations of the ultracompact HII region G10.6-0.4 in 23 GHz radio continuum and the NH3(3,3) line. Our data show that the infall in the molecular material is largely spherical, and does not flatten into a molecular disk at radii as small as 0.03 pc. The spherical infall in the molecular gas matches in location and velocity the infall seen in the ionized gas. We use a non-detection to place a stringent upper limit on the mass of an expanding molecular shell associated with pressure driven expansion of the HII region. These data support a scenario in which the molecular accretion flow passes through an ionization front and becomes an ionized accretion flow onto one or more main sequence stars, not the classical pressure-driven expansion scenario. In the continuum emission we see evidence for externally ionized clumps of molecular gas, and cavities evacuated by an outflow from the central source.

The formation of massive stars may take place at relatively low accretion rates over a long period of time if the accretion can continue past the onset of core hydrogen ignition. The accretion may continue despite the formation of an ionized HII region around the star if the HII region is small enough that the gravitational attraction of the star dominates the thermal pressure of the HII region. The accretion may continue despite radiation pressure acting against dust grains in the molecular gas if the momentum of the accretion flow is sufficient to push the dust grains through a narrow zone of high dust opacity at the ionization boundary and into the HII region where the dust is sublimated. This model of massive star formation by continuing accretion predicts a new class of gravitationally-trapped, long-lived hypercompact HII regions. The observational characteristics of the trapped hypercompact HII regions can be predicted for comparison with observations.

We develop a method of analyzing radio frequency spectral line observations to derive data on the temperature, density, velocity, and molecular abundance of the emitting gas. The method incorporates a radiative transfer code with a new technique for handling overlapping hyperfine emission lines within the accelerated lambda iteration algorithm and a heuristic search algorithm based on simulated annnealing. We apply this method to new observations of N_2H^+ in three Lynds clouds thought to be starless cores in the first stages of star formation and determine their density structure. A comparison of the gas densities derived from the molecular line emission and the millimeter dust emission suggests that the required dust mass opacity is about kappa_{1.3mm}=0.04 cm^2/g, consistent with models of dust grains that have opacities enhanced by ice mantles and fluffy aggregrates.

Using the partially completed Submillimeter Array with five antennas, we have observed the CO J=2-1 and 3-2 emission from the envelope surrounding the carbon star V Hya. The high-angular resolution (2"-4") maps show that V Hya is powering a bipolar molecular jet having an extreme velocity of 70-185 km/s. The axis of this high velocity jet is perpendicular to the major axis of the flattened disk-like envelope, which is expanding with a velocity of ~16 km/s.There is a third kinematic component, a medium-velocity wind having a deprojected velocity of 40-120 km/s moving along the disk plane. Both the high velocity jet and the medium velocity wind have a dynamical time scale of a few hundred years. The flattened structure and the collimated jet observed in V Hya suggests that the formation of asymmetrical structure proceeds while the central star is still in the AGB phase

We present the first 865 mu continuum image with sub-arcsecond resolution obtained with the Submillimeter Array. These data resolve the Orion-KL region into the hot core, the nearby radio source I, the sub-mm counterpart to the infrared source n (radio source L), and new sub-mm continuum sources. The radio to submillimeter emission from source I may be modeled as either the result of proton-electron free-free emission that is optically thick to ~100 GHz plus dust emission that accounts for the majority of the submillimeter flux, or H- free-free emission that gives rise to a power-law spectrum with power-law index of ~1.6. The latter model would indicate similar physical conditions as found in the inner circumstellar environment of Mira variable stars. Future sub-arcsecond observations at shorter sub-mm wavelengths should easily discriminate between these two possibilities. The sub-mm continuum emission toward source n can be interpreted in the framework of emission from an accretion disk.

The search for the onset of star formation in pre-stellar cores has focussed on the identification of an infall signature in the molecular line profiles of tracer species. The classic infall signature is a double peaked line profile with an asymmetry in the strength of the peaks such that the blue peak is stronger. L1689B is a pre-stellar core and infall candidate but new JCMT HCO+ line profile data, presented here, confirms that both blue and red asymmetric line profiles are present in this source. Moreover, a dividing line can be drawn between the locations where each type of profile is found. It is argued that it is unlikely that the line profiles can be interpreted with simple models of infall or outflow and that rotation of the inner regions is the most likely explanation. A rotational model is developed in detail with a new 3D molecular line transport code and it is found that the best type of model is one in which the rotational velocity profile is in between solid body and Keplerian. It is firstly shown that red and blue asymmetric line profiles can be generated with a rotation model entirely in the absence of any infall motion. The model is then quantitively compared with the JCMT data and an iteration over a range of parameters is performed to minmize the difference between the data and model. The results indicate that rotation can dominate the line profile shape even before the onset of infall.

A plausible model is proposed for the enhancement of the abundance of molecular species in bipolar outflow sources. In this model, levels of HCO+ enhancement are considered based on previous chemical calculations, that are assumed to result from shock-induced desorption and photoprocessing of dust grain ice mantles in the boundary layer between the outflow jet and the surrounding envelope. A radiative transfer simulation that incorporates chemical variations within the flow shows that the proposed abundance enhancements in the boundary layer are capable of reproducing the observed characteristics of the outflow seen in HCO+ emission in the star forming core L1527. The radiative transfer simulation also shows that the emission lines from the enhanced molecular species that trace the boundary layer of the outflow exhibit complex line profiles indicating that detailed spatial maps of the line profiles are essential in any attempt to identify the kinematics of potential infall/outflow sources. This study is one of the first applications of a full three dimensional radiative transfer code which incorporates chemical variations within the source.

We present an analysis of the temperature, density, and velocity of the molecular gas in the star-forming core around W51 e2. A previous paper (Ho and Young 1996) describes the kinematic evidence which implies that the core around e2 is contracting onto a young massive star. The current paper presents a technique for modelling the three-dimensional structure of the core by simulating spectral line images of the source and comparing those images to observed data. The primary conclusions of this work are that the molecular gas in e2 is radially contracting at about 5 km/s and that the temperature and density of the gas decrease outward over 0.15 pc scales. The simple model of the collapse of the singular isothermal sphere for low-mass star formation (Shu 1977) is an inadequate description of this high-mass molecular core; better models have temperature decreasing outward as r^-0.6, density as r^-2, and velocity increasing as r^+0.1. The core appears to be spherical rather than disk-like at the scale of these observations, 0.3 pc. In this paper we show how a series of models of gradually increasing complexity can be used to investigate the sensitivity of the model to its parameters. Major sources of uncertainty for this method and this dataset are the interdependence of temperature and density and the assumed NH3 abundance.

We present arcsecond resolution mid-infrared (8--13 $\mu$m) images and photometry of four young stellar objects (YSOs)~--- L1551-IRS5, HL~Tau, AS~205, and AS~209 (V1121~Oph)~--- taken with the Berkeley Mid-Infrared Camera. For AS~205, a known T Tauri binary, we also present near-infrared JHK images and HKL$^{\prime}$ speckle imaging data. All three single stars are unresolved in our mid-IR images. Our data is the first to resolve in the mid-IR both components of the close binary AS~205 (projected separation $\sim$1.3$^{\prime \prime}$ (210~A.U.)). AS~205~North is the IR brighter star in our data while published observations find it to be the optically fainter star. Pre-main sequence evolutionary models suggest the AS~205 system is non-coeval; we discuss possible explanations for this result and comment on the circumstellar material and evolutionary status of this young binary. Nearly all of our objects exhibit changes in their mid-IR flux in measurements separated by intervals of days up to many years; the variations range from 30--300\%. We speculate that the cause of the variability lies in the accretion disks of these objects; the data suggest disk accretion rate fluctuations of nearly an order of magnitude. The existence of large mid-IR variability argues that simultaneous multiwavelength observations are needed for a proper analysis of YSO spectral energy distributions.