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.

The density structure of the interstellar medium (ISM) determines where stars form and release energy, momentum, and heavy elements, driving galaxy evolution. Density variations are seeded and amplified by gas motion, but the exact nature of this motion is unknown across spatial scale and galactic environment. Although dense star-forming gas likely emerges from a combination of instabilities, convergent flows, and turbulence, establishing the precise origin is challenging because it requires quantifying gas motion over many orders of magnitude in spatial scale. Here we measure the motion of molecular gas in the Milky Way and in nearby galaxy NGC 4321, assembling observations that span an unprecedented spatial dynamic range ($10^{-1}{-}10^3$ pc). We detect ubiquitous velocity fluctuations across all spatial scales and galactic environments. Statistical analysis of these fluctuations indicates how star-forming gas is assembled. We discover oscillatory gas flows with wavelengths ranging from $0.3{-}400$ pc. These flows are coupled to regularly-spaced density enhancements that likely form via gravitational instabilities. We also identify stochastic and scale-free velocity and density fluctuations, consistent with the structure generated in turbulent flows. Our results demonstrate that ISM structure cannot be considered in isolation. Instead, its formation and evolution is controlled by nested, interdependent flows of matter covering many orders of magnitude in spatial scale.

We study the formation, evolution and collapse of dense cores by tracking density structures in a magnetohydrodynamic (MHD) simulation. We identify cores using the dendrogram algorithm and utilize machine learning techniques, including principal component analysis (PCA) and the k-means clustering algorithm to analyze the full density and velocity dispersion profiles of these cores. We find that there exists an evolutionary sequence consisting of three distinct phases: i) the formation of turbulent density structures (Phase I), ii) the dissipation of turbulence and the formation of coherent cores (Phase II), and iii) the transition to protostellar cores through gravitational collapse (Phase III). In dynamically evolving molecular clouds, the existence of these three phases corresponds to the coexistence of three populations of cores with distinct physical properties. The prestellar and protostellar cores frequently analyzed in previous studies of observations and simulations belong to the last phase in this evolutionary picture. We derive typical lifetimes of 1.4$\pm$1.0$\times$10$^5$ yr, 3.3$\pm$1.4$\times$10$^5$ yr and 3.3$\pm$1.4$\times$10$^5$ yr, respectively for Phase I, II and III. We find that cores can form from both converging flows and filament fragmentation and that cores may form both inside and outside the filaments. We then compare our results to previous observations of coherent cores and provide suggestions for future observations to study cores belonging to the three phases.

Star formation is a multi-scale, multi-physics problem ranging from the size scale of molecular clouds ($\sim$10s pc) down to the size scales of dense prestellar cores ($\sim$0.1 pc) that are the birth sites of stars. Several physical processes like turbulence, magnetic fields and stellar feedback, such as radiation pressure and outflows, are more or less important for different stellar masses and size scales. During the last decade a variety of technological and computing advances have transformed our understanding of star formation through the use of multi-wavelength observations, large scale observational surveys, and multi-physics multi-dimensional numerical simulations. Additionally, the use of synthetic observations of simulations have provided a useful tool to interpret observational data and evaluate the importance of various physical processes on different scales in star formation. Here, we review these recent advancements in both high- ($M \gtrsim 8 \, M_{\rm \odot}$) and low-mass star formation.

We report ALMA Band 6 continuum observations of 2000 AU resolution toward four massive molecular clouds in the Central Molecular Zone of the Galaxy. To study gas fragmentation, we use the dendrogram method to identify cores as traced by the dust continuum emission. The four clouds exhibit different fragmentation states at the observed resolution despite having similar masses at the cloud scale ($\sim$1--5 pc). Assuming a constant dust temperature of 20 K, we construct core mass functions of the clouds and find a slightly top-heavy shape as compared to the canonical initial mass function, but we note several significant uncertainties that may affect this result. The characteristic spatial separation between the cores as identified by the minimum spanning tree method, $\sim$$10^4$ AU, and the characteristic core mass, 1--7 $M_\odot$, are consistent with predictions of thermal Jeans fragmentation. The three clouds showing fragmentation may be forming OB associations (stellar mass $\sim$$10^3$ $M_\odot$). None of the four clouds under investigation seem to be currently able to form massive star clusters like the Arches and the Quintuplet ($\sim$$10^4$ $M_\odot$), but they may form such clusters by further gas accretion onto the cores.

We report the results of a pilot program for a Green Bank Telescope (GBT) MUSTANG Galactic Plane survey at 3 mm (90 GHz), MGPS90. The survey achieves a typical $1\sigma$ depth of $1-2$ mJy beam$^{-1}$ with a 9" beam. We describe the survey parameters, quality assessment process, cataloging, and comparison with other data sets. We have identified 709 sources over seven observed fields selecting some of the most prominent millimeter-bright regions between $0\deg < \ell < 50\deg$ (total area $\approx 7.5 \deg^2$). The majority of these sources have counterparts at other wavelengths. By applying flux selection criteria to these sources, we successfully recovered several known hypercompact HII (HCHII) regions, but did not confirm any new ones. We identify 126 sources that have mm-wavelength counterparts but do not have cm-wavelength counterparts and are therefore candidate HCHII regions; of these, 10 are morphologically compact and are strong candidates for new HCHII regions. Given the limited number of candidates in the extended area in this survey compared to the relatively large numbers seen in protoclusters W51 and W49, it appears that most HCHII regions exist within dense protoclusters. Comparing the counts of HCHII to ultracompact HII (UCHII) regions, we infer the HCHII region lifetime is 16-46% that of the UCHII region lifetime. We additionally separated the 3 mm emission into dust and free-free emission by comparing with archival 870 $\mu$m and 20 cm data. In the selected pilot fields, most ($\gtrsim80$%) of the 3 mm emission comes from plasma, either through free-free or synchrotron emission.

The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular clouds. We find that the results on relative alignment between cores and background magnetic field in synthetic observations slightly disagree with those measured in fully 3D simulation data, which is partly because cores identified in projected 2D maps tend to coexist within filamentary structures, while 3D cores are generally more rounded. In addition, we examine the progression of magnetic field from pc- to core-scale in the simulation, which is consistent with the anisotropic core formation model that gas preferably flow along the magnetic field toward dense cores. When comparing the observed cores identified from the GBT Ammonia Survey (GAS) and Planck polarization-inferred magnetic field orientations, we find that the relative core-field alignment has a regional dependence among different clouds. More specifically, we find that dense cores in the Taurus molecular cloud tend to align perpendicular to the background magnetic field, while those in Perseus and Ophiuchus tend to have random (Perseus) or slightly parallel (Ophiuchus) orientations with respect to the field. We argue that this feature of relative core-field orientation could be used to probe the relative significance of the magnetic field within the cloud.

The proper motions of the three stars ejected from Orion's OMC1 cloud core are combined with the requirement that their center of mass is gravitationally bound to OMC1 to show that radio source I (Src I) is likely to have a mass around 15 Solar masses consistent with recent measurements. Src I, the star with the smallest proper motion, is suspected to be either an AU-scale binary or a protostellar merger remnant produced by a dynamic interaction ~550 years ago. Near-infrared 2.2 um images spanning ~21 years confirm the ~55 km/s motion of `source x' (Src x) away from the site of stellar ejection and point of origin of the explosive OMC1 protostellar outflow. The radial velocities and masses of the Becklin-Neugebauer (BN) object and Src I constrain the radial velocity of Src x to be V_{LSR} = -28 +/-10 km/s . Several high proper-motion radio sources near BN, including Zapata 11 ([ZRK2004] 11) and a diffuse source near IRc 23, may trace a slow bipolar outflow from BN. The massive disk around Src I is likely the surviving portion of a disk that existed prior to the stellar ejection. Though highly perturbed, shocked, and re-oriented by the N-body interaction, enough time has elapsed to allow the disk to relax with its spin axis roughly orthogonal to the proper motion.

We imaged the continuum and molecular line emission from Orion Source I (SrcI) with up to 30 mas (12 AU) resolution at 43, 99, 223, and 340 GHz in an attempt to probe the structure and chemistry of the circumstellar disk and bipolar outflow associated with this high mass protostar. The continuum spectral index ranges from $\sim$2 along the midplane of the disk to $\sim$3 along the edges, consistent with dust that is optically thick in the midplane but becomes optically thin at the periphery. Salt (NaCl) emission is visible where the dust is optically thin; it provides a unique tracer of the velocity field within the disk. All other molecules that we have mapped - H$_2$O, AlO, SiO, SiS, SO, and SO$_2$ - appear to originate primarily in the bipolar outflow. The base of the outflow is corotating with the disk. SiS shows a filamentary structure that is most prominent along the edges of the outflow. The molecular distributions suggest that Si and Al released from dust grains in the disk react with oxygen derived from H$_2$O to form SiO and AlO, and with SO and SO$_2$ to form SiS.

We present evidence of 6.7 GHz methanol (CH3OH)~and 4.8 GHz formaldehyde (H2CO) absorption towards the Galactic Center (GC) point-source `N3'. Both absorption features are unresolved and spatially aligned with N3. The 6.7 GHz CH3OH contains a single velocity component (centered at ~10 km/s) while the 4.8 GHz H2CO shows two velocity components (centered at ~-3 and +8 km/s). We find that the velocity of these absorption components are similar to that of emission lines from other molecules (e.g., SiO and H3CN) detected toward this compact-source (-13 to +25 km/s; `N3 cloud'). The detection of these absorption features is a firm indication that some of the molecular gas in the N3 molecular cloud is on the near-side of the continuum source. Analysis of the CH3OH absorption kinematics shows a relatively large velocity dispersion (3.8 km/s) for the size scale of this feature (<0.1'', <0.01 pc at the GC; Ludovici et al. 2016), when compared with other similar size GC clouds in the Larson linewidth-size relationship. Further, this linewidth is closer to velocity dispersion measurements for size scales of 1.3 pc, which is roughly the width of the N3 cloud (25''; 1.0 pc). We argue that this relatively broad linewidth, over a small cross-sectional area, is due to turbulence through the depth of the cloud, where the cloud has a presumed line-of-sight thickness of ~1 pc.

We report the discovery of maser emission in the two lowest rotational transitions of CS toward the high-mass protostar W51 e2e with ALMA and the JVLA. The masers from CS J=1-0 and J=2-1 are neither spatially nor spectrally coincident (they are separated by ~150 AU and ~30 km/s), but both appear to come from the base of the blueshifted outflow from this source. These CS masers join a growing list of rarely-detected maser transitions that may trace a unique phase in the formation of high-mass protostars.

The giant molecular cloud Sagittarius B2 (hereafter SgrB2) is the most massive region with ongoing high-mass star formation in the Galaxy. In the southern region of the 40-pc large envelope of SgrB2, we encounter the SgrB2(DS) region which hosts more than 60 high-mass protostellar cores distributed in an arc shape around an extended HII region. We use the Very Large Array in its CnB and D configurations, and in the frequency bands C (4--8 GHz) and X (8--12 GHz) to observe the whole SgrB2 complex. Continuum and radio recombination line maps are obtained. We detect radio continuum emission in SgrB2(DS) in a bubble-shaped structure. From 4 to 12 GHz, we derive a spectral index between -1.2 and -0.4, indicating the presence of non-thermal emission. We decompose the contribution from thermal and non-thermal emission, and find that the thermal component is clumpy and more concentrated, while the non-thermal component is more extended and diffuse. The radio recombination lines in the region are found to be not in local thermodynamic equilibrium (LTE) but stimulated by the non-thermal emission. The thermal free-free emission is likely tracing an HII region ionized by an O7 star, while the non-thermal emission can be generated by relativistic electrons created through first-order Fermi acceleration. We have developed a simple model of the SgrB2(DS) region and found that first-order Fermi acceleration can reproduce the observed flux density and spectral index.

We present new observations of C-band continuum emission and masers to assess high-mass ($>$8 $M_\odot$) star formation at early evolutionary phases in the inner 200 pc of the Central Molecular Zone (CMZ) of the Galaxy. The continuum observation is complete to free-free emission from stars above 10-11 $M_\odot$ in 91% of the covered area. We identify 104 compact sources in the continuum emission, among which five are confirmed ultracompact H II regions, 12 are candidates of ultra-compact H II regions, and the remaining 87 sources are mostly massive stars in clusters, field stars, evolved stars, pulsars, extragalactic sources, or of unknown nature that is to be investigated. We detect class II CH$_3$OH masers at 23 positions, among which six are new detections. We confirm six known H$_2$CO masers in two high-mass star forming regions, and detect two new H$_2$CO masers toward the Sgr C cloud, making it the ninth region in the Galaxy that contains masers of this type. In spite of these detections, we find that current high-mass star formation in the inner CMZ is only taking place in seven isolated clouds. The results suggest that star formation at early evolutionary phases in the CMZ is about 10 times less efficient than expected by the dense gas star formation relation, which is in line with previous studies that focus on more evolved phases of star formation. This means that if there will be any impending, next burst of star formation in the CMZ, it has not yet begun.

We present initial results from the K-band focal plane array Examinations of Young STellar Object Natal Environments (KEYSTONE) survey, a large project on the 100-m Green Bank Telescope mapping ammonia emission across eleven giant molecular clouds at distances of $0.9-3.0$ kpc (Cygnus X North, Cygnus X South, M16, M17, MonR1, MonR2, NGC2264, NGC7538, Rosette, W3, and W48). This data release includes the NH$_3$ (1,1) and (2,2) maps for each cloud, which are modeled to produce maps of kinetic temperature, centroid velocity, velocity dispersion, and ammonia column density. Median cloud kinetic temperatures range from $11.4\pm2.2$ K in the coldest cloud (MonR1) to $23.0\pm6.5$ K in the warmest cloud (M17). Using dendrograms on the NH$_3$ (1,1) integrated intensity maps, we identify 856 dense gas clumps across the eleven clouds. Depending on the cloud observed, $40-100\%$ of the clumps are aligned spatially with filaments identified in H$_2$ column density maps derived from SED-fitting of dust continuum emission. A virial analysis reveals that 523 of the 835 clumps ($\sim63\%$) with mass estimates are bound by gravity alone. We find no significant difference between the virial parameter distributions for clumps aligned with the dust-continuum filaments and those unaligned with filaments. In some clouds, however, hubs or ridges of dense gas with unusually high mass and low virial parameters are located within a single filament or at the intersection of multiple filaments. These hubs and ridges tend to host water maser emission, multiple 70$\mu$m-detected protostars, and have masses and radii above an empirical threshold for forming massive stars.

We have investigated the formation and kinematics of sub-mm continuum cores in the Orion A molecular cloud. A comparison between sub-mm continuum and near infrared extinction shows a continuum core detection threshold of $A_V\sim$ 5-10 mag. The threshold is similar to the star formation extinction threshold of $A_V\sim$ 7 mag proposed by recent work, suggesting a universal star formation extinction threshold among clouds within 500 pc to the Sun. A comparison between the Orion A cloud and a massive infrared dark cloud G28.37+0.07 indicates that Orion A produces more dense gas within the extinction range 15 mag $\lesssim A_V \lesssim$ 60 mag. Using data from the CARMA-NRO Orion Survey, we find that dense cores in the integral-shaped filament (ISF) show sub-sonic core-to-envelope velocity dispersion that is significantly less than the local envelope line dispersion, similar to what has been found in nearby clouds. Dynamical analysis indicates that the cores are bound to the ISF. An oscillatory core-to-envelope motion is detected along the ISF. Its origin is to be further explored.

Software is a critical part of modern research, and yet there are insufficient mechanisms in the scholarly ecosystem to acknowledge, cite, and measure the impact of research software. The majority of academic fields rely on a one-dimensional credit model whereby academic articles (and their associated citations) are the dominant factor in the success of a researcher's career. In the petabyte era of astronomical science, citing software and measuring its impact enables academia to retain and reward researchers that make significant software contributions. These highly skilled researchers must be retained to maximize the scientific return from petabyte-scale datasets. Evolving beyond the one-dimensional credit model requires overcoming several key challenges, including the current scholarly ecosystem and scientific culture issues. This white paper will present these challenges and suggest practical solutions for elevating the role of software as a product of the research enterprise.

We have used an unbiased, spectral line-survey that covers the frequency range from 211 to 275 GHz and was obtained with ALMA (angular resolution of 0.4 arcsec) to study the small-scale structure of the dense gas in Sagittarius B2 (north). Eight filaments are found converging to the central hub and extending for about 0.1 pc. The spatial structure, together with the presence of the massive central region, suggest that these filaments may be associated with accretion processes. In order to derive the kinematic properties of the gas in a chemically line-rich source like Sgr B2(N), we have developed a new tool that stacks all the detected transition lines of any molecular species. This permits to increase the signal-to-noise ratio of our observations and average out line blending effects, which are a common problem in line-rich regions. We derive velocity gradients along the filaments of about 20-100 km s$^{-1}$ pc$^{-1}$, which are 10-100 times larger than those typically found on larger scales (1 pc) in other star-forming regions. The mass accretion rates of individual filaments are about 0.05 M$_\odot$ yr$^{-1}$, which result in a total accretion rate of 0.16 M$_\odot$ yr$^{-1}$. Some filaments harbor dense cores that are likely forming stars and stellar clusters. The stellar content of these dense cores is on the order of 50% of the total mass. We conclude that the cores may merge in the center when already forming stellar clusters but still containing a significant amount of gas, resulting in a "damp" merger. The high density and mass of the central region, combined with the presence of converging filaments with high mass, high accretion rates and embedded dense cores already forming stars, suggest that Sgr B2(N) may have the potential to evolve into a super stellar cluster.

We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission, especially visible as arc structures to the north and south. The Keplerian kinematics are most prominently seen in the vibrationally excited water line, H2O (Eu=3461.9K). The mass of the central source found by modelling the Keplerian rotation is consistent with 45+/-10Mo. The H30alpha (231.9GHz) radio-recombination line and the SiO (5-4) molecular line were detected at up to the 10 sigma$ level. The estimated disc mass is 0.6-2.6Mo under the optically thin assumption. Analysis of the Toomre Q parameter, in the optically thin regime, indicates that the disc stability is highly dependent on temperature. The disc currently appears stable for temperatures >150K, this does not preclude that the substructures formed earlier through disc fragmentation.

We present the observation and analysis of newly discovered coherent structures in the L1688 region of Ophiuchus and the B18 region of Taurus. Using data from the Green Bank Ammonia Survey (GAS), we identify regions of high density and near-constant, almost-thermal, velocity dispersion. Eighteen coherent structures are revealed, twelve in L1688 and six in B18, each of which shows a sharp "transition to coherence" in velocity dispersion around its periphery. The identification of these structures provides a chance to study the coherent structures in molecular clouds statistically. The identified coherent structures have a typical radius of 0.04 pc and a typical mass of 0.4 Msun, generally smaller than previously known coherent cores identified by Goodman et al. (1998), Caselli et al. (2002), and Pineda et al. (2010). We call these structures "droplets." We find that unlike previously known coherent cores, these structures are not virially bound by self-gravity and are instead predominantly confined by ambient pressure. The droplets have density profiles shallower than a critical Bonnor-Ebert sphere, and they have a velocity (VLSR) distribution consistent with the dense gas motions traced by NH3 emission. These results point to a potential formation mechanism through pressure compression and turbulent processes in the dense gas. We present a comparison with a magnetohydrodynamic simulation of a star-forming region, and we speculate on the relationship of droplets with larger, gravitationally bound coherent cores, as well as on the role that droplets and other coherent structures play in the star formation process.

We present TurbuStat (v1.0): a Python package for computing turbulence statistics in spectral-line data cubes. TurbuStat includes implementations of fourteen methods for recovering turbulent properties from observational data. Additional features of the software include: distance metrics for comparing two data sets; a segmented linear model for fitting lines with a break-point; a two-dimensional elliptical power-law model; multi-core fast-fourier-transform support; a suite for producing simulated observations of fractional Brownian Motion fields, including two-dimensional images and optically-thin HI data cubes; and functions for creating realistic world coordinate system information for synthetic observations. This paper summarizes the TurbuStat package and provides representative examples using several different methods. TurbuStat is an open-source package and we welcome community feedback and contributions.

The current generation of radio and millimeter telescopes, particularly the Atacama Large Millimeter Array (ALMA), offers enormous advances in observing capabilities. While these advances represent an unprecedented opportunity to facilitate scientific understanding, the increased complexity in the spatial and spectral structure of these ALMA data cubes lead to challenges in their interpretation. In this paper, we perform a feasibility study for applying topological data analysis and visualization techniques never before tested by the ALMA community. Through techniques based on contour trees, we seek to improve upon existing analysis and visualization workflows of ALMA data cubes, in terms of accuracy and speed in feature extraction. We review our application development process in building effective analysis and visualization capabilities for the astrophysicists. We also summarize effective design practices by identifying domain-specific needs of simplicity, integrability, and reproducibility, in order to best target and service the large astrophysics community.

It is usually assumed that the stellar initial mass function (IMF) takes a universal form and that there exists a direct mapping between this and the distribution of natal core masses (the core mass function, CMF). The IMF and CMF have been best characterized in the Solar neighborhood. Beyond 500~pc from the Sun, in diverse environments where metallicity varies and massive star feedback may dominate, the IMF has been measured only incompletely and imprecisely, while the CMF has hardly been measured at all. In order to establish if the IMF and CMF are indeed universal and related to each other, it is necessary to: 1) perform multi-wavelength large-scale imaging and spectroscopic surveys of different environments across the Galaxy; 2) require an angular resolution of < 0.1'' in the optical/near-IR for stars and < 5'' in the far-IR for cores; 3) achieve far-IR sensitivities to probe 0.1~Msol cores at 2--3 kpc.

Young massive clusters (YMCs) are the most compact, high-mass stellar systems still forming at the present day. The precursor clouds to such systems are, however, rare due to their large initial gas mass reservoirs and rapid dispersal timescales due to stellar feedback. Nonetheless, unlike their high-z counterparts, these precursors are resolvable down to the sites of individually forming stars, and hence represent the ideal environments in which to test the current theories of star and cluster formation. Using high angular resolution (1$^{\prime\prime}$ / 0.05pc) and sensitivity ALMA observations of two YMC progenitor clouds in the Galactic Centre, we have identified a suite of molecular line transitions -- e.g. c-C$_{3}$H$_{2} $($7-6$) -- that are believed to be optically thin, and reliably trace the gas structure in the highest density gas on star-forming core scales. We conduct a virial analysis of the identified core and proto-cluster regions, and show that half of the cores (5/10) and both proto-clusters are unstable to gravitational collapse. This is the first kinematic evidence of global gravitational collapse in YMC precursor clouds at such an early evolutionary stage. The implications are that if these clouds are to form YMCs, then they likely do so via the "conveyor-belt" mode, whereby stars continually form within dispersed dense gas cores as the cloud undergoes global gravitational collapse. The concurrent contraction of both the cluster-scale gas and embedded (proto)stars ultimately leads to the high (proto)stellar density in YMCs.

The closest galaxy center, our own Central Molecular Zone (CMZ; the central 500 pc of the Milky Way), is a powerful laboratory for studying the secular processes that shape galaxies across cosmic time, from large-scale gas flows and star formation to stellar feedback and interaction with a central supermassive black hole. Research over the last decade has revealed that the process of converting gas into stars in galaxy centers differs from that in galaxy disks. The CMZ is the only galaxy center in which we can identify and weigh individual forming stars, so it is the key location to establish the physical laws governing star formation and feedback under the conditions that dominate star formation across cosmic history. Large-scale surveys of molecular and atomic gas within the inner kiloparsec of the Milky Way (~10 degrees) will require efficient mapping capabilities on single-dish radio telescopes. Characterizing the detailed star formation process will require large-scale, high-resolution surveys of the protostellar populations and small-scale gas structure with dedicated surveys on the Atacama Large Millimeter/submillimeter Array, and eventually with the James Webb Space Telescope, the Next Generation Very Large Array, and the Origins Space Telescope.

We investigate star formation at very early evolutionary phases in five massive clouds in the inner 500 pc of the Galaxy, the Central Molecular Zone. Using interferometer observations of H$_2$O masers and ultra-compact H II regions, we find evidence of ongoing star formation embedded in cores of 0.2 pc scales and $\gtrsim$10$^5$ cm$^{-3}$ densities. Among the five clouds, Sgr C possesses a high (9%) fraction of gas mass in gravitationally bound and/or protostellar cores, and follows the dense ($\gtrsim$10$^4$ cm$^{-3}$) gas star formation relation that is extrapolated from nearby clouds. The other four clouds have less than 1% of their cloud masses in gravitationally bound and/or protostellar cores, and star formation rates 10 times lower than predicted by the dense gas star formation relation. At the spatial scale of these cores, the star formation efficiency is comparable to that in Galactic disk sources. We suggest that the overall inactive star formation in these Central Molecular Zone clouds could be because there is much less gas confined in gravitationally bound cores, which may be a result of the strong turbulence in this region and/or the very early evolutionary stage of the clouds when collapse has only recently started.

The fraction of stars forming in compact, gravitationally bound clusters (the `cluster formation efficiency' or CFE) is an important quantity for deriving the spatial clustering of stellar feedback and for tracing star formation using stellar clusters across the Universe. Observations of clusters in nearby galaxies have revealed a strong dependence of the CFE on the local gas density, indicating that more stars form in star clusters when the star formation rate surface density is higher. Previously, it has not been possible to test this relation at very young ages and in clusters with individual stars resolved due to the universally-low densities in the cluster-forming regions in the Local Group. This has even led to the suggestion that the CFE increases with distance from the Sun, which would suggest an observational bias. However, the Central Molecular Zone of the Milky Way hosts clouds with densities that are orders of magnitude higher than anywhere else in the Local Group. We report a measurement of the CFE in the highest-density region in the Galaxy, Sgr~B2, based on ALMA observations of high-mass young stellar objects. We find that over a third of the stars ($37\pm7\%$) in Sgr~B2 are forming in bound clusters. This value is consistent with the predictions of environmentally-dependent models for the CFE and is inconsistent with a constant CFE in the Galaxy.

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.

In this paper we provide a comprehensive description of the internal dynamics of G0.253+0.016 (a.k.a. 'the Brick'); one of the most massive and dense molecular clouds in the Galaxy to lack signatures of widespread star formation. As a potential host to a future generation of high-mass stars, understanding largely quiescent molecular clouds like G0.253+0.016 is of critical importance. In this paper, we reanalyse Atacama Large Millimeter Array cycle 0 HNCO $J=4(0,4)-3(0,3)$ data at 3 mm, using two new pieces of software which we make available to the community. First, scousepy, a Python implementation of the spectral line fitting algorithm scouse. Secondly, acorns (Agglomerative Clustering for ORganising Nested Structures), a hierarchical n-dimensional clustering algorithm designed for use with discrete spectroscopic data. Together, these tools provide an unbiased measurement of the line of sight velocity dispersion in this cloud, $\sigma_{v_{los}, {\rm 1D}}=4.4\pm2.1$ kms$^{-1}$, which is somewhat larger than predicted by velocity dispersion-size relations for the Central Molecular Zone (CMZ). The dispersion of centroid velocities in the plane of the sky are comparable, yielding $\sigma_{v_{los}, {\rm 1D}}/\sigma_{v_{pos}, {\rm 1D}}\sim1.2\pm0.3$. This isotropy may indicate that the line-of-sight extent of the cloud is approximately equivalent to that in the plane of the sky. Combining our kinematic decomposition with radiative transfer modelling we conclude that G0.253+0.016 is not a single, coherent, and centrally-condensed molecular cloud; 'the Brick' is not a \emph{brick}. Instead, G0.253+0.016 is a dynamically complex and hierarchically-structured molecular cloud whose morphology is consistent with the influence of the orbital dynamics and shear in the CMZ.

We report the discovery of the first 14NH3 (2,2) maser, seen in the Sgr B2 Main star forming region near the center of the Milky Way, using data from the Very Large Array radio telescope. The maser is seen in both lower resolution (3" or ~0.1 pc) data from 2012 and higher resolution (0''.1 or ~1000 AU) data from 2018. In the higher resolution data ammonia (2,2) maser emission is detected toward 5 independent spots. The maser spots are not spatially or kinematically coincident with any other masers in this region, or with the peaks of the radio continuum emission from the numerous ultracompact and hypercompact \hii\, regions in this area. While the (2,2) maser spots are spatially unresolved in our highest resolution observations, they have unusually broad linewidths of several kilometers per second, which suggests that each of these spots consists of multiple masers tracing unresolved velocity structure. No other ammonia lines observed in Sgr B2 Main are seen to be masers, which continues to challenge theories of ammonia, maser emission that predict simultaneous maser emission in multiple ammonia transitions.

The nuclear starburst within the central $\sim 15^{\prime\prime}$ ($\sim 250$ pc; $1^{\prime\prime} \simeq 17$ pc) of NGC253 has been extensively studied as a prototype for the starburst phase in galactic evolution. Atacama Large Millimeter/submillimeter Array (ALMA) imaging within receiver Bands 6 and 7 have been used to investigate the dense gas structure, kinetic temperature, and heating processes which drive the NGC253 starburst. Twenty-nine transitions from fifteen molecular species/isotopologues have been identified and imaged at $1.^{\prime\prime}5$ to $0.^{\prime\prime}4$ resolution, allowing for the identification of five of the previously-studied giant molecular clouds (GMCs) within the central molecular zone (CMZ) of NGC253. Ten transitions from the formaldehyde (H$_2$CO) molecule have been used to derive the kinetic temperature within the $\sim 0.^{\prime\prime}5$ to $5^{\prime\prime}$ dense-gas structures imaged. On $\sim 5^{\prime\prime}$ scales we measure $T_K \gtrsim 50$ K, while on size scales $\lesssim 1^{\prime\prime}$ we measure $T_K \gtrsim 300$ K. These kinetic temperature measurements further delineate the association between potential sources of dense gas heating. We have investigated potential heating sources by comparing our measurements to models which predict the physical conditions associated with dense molecular clouds that possess a variety of heating mechanisms. This comparison has been supplemented with tracers of recently-formed massive stars (Br$\gamma$) and shocks ([FeII]). Derived molecular column densities point to a radially-decreasing abundance of molecules with sensitivity to cosmic ray and mechanical heating within the NGC253 CMZ. These measurements are consistent with radio spectral index calculations which suggest a higher concentration of cosmic ray producing supernova remnants within the central 10 pc of NGC253.

We report the detection of NaCl, KCl, and their $^{37}$Cl and $^{41}$K isotopologues toward the disk around Orion SrcI. About 60 transitions of these molecules were identified. This is the first detection of these molecules in the interstellar medium not associated with the ejecta of evolved stars. It is also the first ever detection of the vibrationally excited states of these lines in the ISM above v = 1, with firm detections up to v = 6. The salt emission traces the region just above the continuum disk, possibly forming the base of the outflow. The emission from the vibrationally excited transitions is inconsistent with a single temperature, implying the lines are not in LTE. We examine several possible explanations of the observed high excitation lines, concluding that the vibrational states are most likely to be radiatively excited via rovibrational transitions in the 25-35 {\mu}m (NaCl) and 35-45 {\mu}m (KCl) range. We suggest that the molecules are produced by destruction of dust particles. Because these molecules are so rare, they are potentially unique tools for identifying high-mass protostellar disks and measuring the radiation environment around accreting young stars.

astroquery is a collection of tools for requesting data from databases hosted on remote servers with interfaces exposed on the internet, including those with web pages but without formal application program interfaces (APIs). These tools are built on the Python requests package, which is used to make HTTP requests, and astropy, which provides most of the data parsing functionality. astroquery modules generally attempt to replicate the web page interface provided by a given service as closely as possible, making the transition from browser-based to command-line interaction easy. astroquery has received significant contributions from throughout the astronomical community, including several significant contributions from telescope archives. astroquery enables the creation of fully reproducible workflows from data acquisition through publication. This paper describes the philosophy, basic structure, and development model of the astroquery package. The complete documentation for astroquery can be found at http://astroquery.readthedocs.io/.

We present an initial overview of the filamentary structure in the Orion A molecular cloud utilizing a high angular and velocity resolution C$^{18}$O(1-0) emission map that was recently produced as part of the CARMA-NRO Orion Survey. The main goal of this study is to build a credible method to study varying widths of filaments which has previously been linked to star formation in molecular clouds. Due to the diverse star forming activities taking place throughout its $\sim$20 pc length, together with its proximity of 388 pc, the Orion A molecular cloud provides an excellent laboratory for such an experiment to be carried out with high resolution and high sensitivity. Using the widely-known structure identification algorithm, DisPerSE, on a 3-dimensional (PPV) C$^{18}$O cube, we identified 625 relatively short (the longest being 1.74 pc) filaments over the entire cloud. We study the distribution of filament widths using FilChaP, a python package that we have developed and made publicly available. We find that the filaments identified in a 2 square degree PPV cube do not overlap spatially, except for the complex OMC-4 region that shows distinct velocity components along the line of sight. The filament widths vary between 0.02 and 0.3 pc depending on the amount of substructure that a filament possesses. The more substructure a filament has, the larger is its width. We also find that despite this variation, the filament width shows no anticorrelation with the central column density which is in agreement with previous Herschel observations.

We present combined observations of the NH3 (J,K) = (1,1) and (2,2) inversion transitions towards OMC1 in Orion A obtained by the Karl G. Jansky Very Large Array (VLA) and the 100 m Robert C. Byrd Green Bank Telescope (GBT). With an angular resolution of 6" (0.01 pc), these observations reveal with unprecedented detail the complex filamentary structure extending north of the active Orion BN/KL region in a field covering 6' x 7'. We find a 0.012 pc wide filament within OMC1, with an aspect ratio of ~37:1, that was missed in previous studies. Its orientation is directly compared to the relative orientation of the magnetic field from the James Clerk Maxwell Telescope BISTRO survey in Orion A. We find a small deviation of ~11 deg between the mean orientation of the filament and the magnetic field, suggesting that they are almost parallel to one another. The filament's column density is estimated to be 2-3 orders of magnitude larger than the filaments studied with Herschel and is possibly self-gravitating given the low values of turbulence found. We further produce maps of the gas kinematics by forward modeling the hyperfine structure of the NH3 (J,K) = (1,1) and (2,2) lines. The resulting distribution of velocity dispersions peaks at ~0.5 km/s, close to the subsonic regime of the gas. This value is about 0.2 km/s smaller than previously measured in single-dish observations of the same region, suggesting that higher angular and spectral resolution observations will identify even lower velocity dispersions that might reach the subsonic turbulence regime in dense gas filaments.

We define the molecular cloud properties of the Milky Way first quadrant using data from the JCMT CO(3-2) High Resolution Survey. We apply the Spectral Clustering for Interstellar Molecular Emission Segmentation (SCIMES) algorithm to extract objects from the full-resolution dataset, creating the first catalog of molecular clouds with a large dynamic range in spatial scale. We identify $>85\,000$ clouds with two clear sub-samples: $\sim35\,500$ well-resolved objects and $\sim540$ clouds with well-defined distance estimations. Only 35% of the cataloged clouds (as well as the total flux encompassed by them) appear enclosed within the Milky Way spiral arms. The scaling relationships between clouds with known distances are comparable to the characteristics of the clouds identified in previous surveys. However, these relations between integrated properties, especially from the full catalog, show a large intrinsic scatter ($\sim0.5$ dex), comparable to other cloud catalogs of the Milky Way and nearby galaxies. The mass distribution of molecular clouds follows a truncated-power law relationship over three orders of magnitude in mass with a form $dN/dM \propto M^{-1.7}$ with a clearly defined truncation at an upper mass of $M_0 \sim 3 \times 10^6~M_\odot$, consistent with theoretical models of cloud formation controlled by stellar feedback and shear. Similarly, the cloud population shows a power-law distribution of size with $dN/dR \propto R^{-2.8}$ with a truncation at $R_0 = 70$ pc.

We report new, $\sim$1000 AU spatial resolution observations of 225 GHz dust continuum emission towards the OB cluster-forming molecular clump G33.92+0.11. On parsec scales, this molecular clump presents a morphology with several arm-like dense gas structures surrounding the two central massive ($\gtrsim$100 $M_{\odot}$) cores. From the new, higher resolution observations, we identified 28 localized, spatially compact dust continuum emission sources, which may be candidates of young stellar objects. Only one of them is not embedded within known arm-like (or elongated) dense gas structures. The spatial separations of these compact sources can be very well explained by Jeans lengths. We found that G33.92+0.11 may be consistently described by a marginally centrifugally supported, Toomre unstable accretion flow which is approximately in a face-on projection. The arm-like overdensities are natural consequence of the Toomre instability, which can fragment to form young stellar objects in shorter time scales than the timescale of the global clump contraction. On our resolved spatial scales, there is not yet evidence that the fragmentation is halted by turbulence, magnetic field, or stellar feedback.

Outflows from AGB stars enrich the Galactic environment with metals and inject mechanical energy into the ISM. Radio spectroscopy can recover both properties through observations of molecular lines. We present results from SWAG: "Survey of Water and Ammonia in the Galactic Center". The survey covers the entire Central Molecular Zone (CMZ), the inner 3.35deg x 0.9deg (~480 x 130pc) of the Milky Way that contains ~5x10^7 Mo of molecular gas. Although our survey primarily targets the CMZ, we observe across the entire sightline through the Milky Way. AGB stars are revealed by their signature of double peaked 22 GHz water maser lines. They are distinguished by their spectral signatures and their luminosities, which reach up to 10^-7 Lo. Higher luminosities are usually associated with Young Stellar Objects located in CMZ star forming regions. We detect a population of ~600 new water masers that can likely be associated with AGB outflows.

NGC 253 hosts the nearest nuclear starburst. Previous observations show a region rich in molecular gas, with dense clouds associated with recent star formation. We used ALMA to image the 350 GHz dust continuum and molecular line emission from this region at 2 pc resolution. Our observations reveal ~14 bright, compact (~2-3 pc FWHM) knots of dust emission. Most of these sources are likely to be forming super star clusters (SSCs) based on their inferred dynamical and gas masses, association with 36 GHz radio continuum emission, and coincidence with line emission tracing dense, excited gas. One source coincides with a known SSC, but the rest remain invisible in Hubble near-infrared (IR) imaging. Our observations imply that gas still constitutes a large fraction of the overall mass in these sources. Their high brightness temperature at 350 GHz also implies a large optical depth near the peak of the IR spectral energy distribution. As a result, these sources may have large IR photospheres and the IR radiation force likely exceeds L/c. Still, their moderate observed velocity dispersions suggest that feedback from radiation, winds, and supernovae are not yet disrupting most sources. This mode of star formation appears to produce a large fraction of stars in the burst. We argue for a scenario in which this phase lasts ~1 Myr, after which the clusters shed their natal cocoons but continue to produce ionizing photons. The strong feedback that drives the observed cold gas and X-ray outflows likely occurs after the clusters emerge from this early phase.

We present MUSE integral field data of two HII region complexes in the Large Magellanic Cloud (LMC), N44 and N180. Both regions consist of a main superbubble and a number of smaller, more compact HII regions that formed on the edge of the superbubble. For a total of 11 HII regions, we systematically analyse the radiative and mechanical feedback from the massive O-type stars on the surrounding gas. We exploit the integral field property of the data and the coverage of the HeII$\lambda$5412 line to identify and classify the feedback-driving massive stars, and from the estimated spectral types and luminosity classes we determine the stellar radiative output in terms of the ionising photon flux $Q_{0}$. We characterise the HII regions in terms of their sizes, morphologies, ionisation structure, luminosity and kinematics, and derive oxygen abundances via emission line ratios. We analyse the role of different stellar feedback mechanisms for each region by measuring the direct radiation pressure, the pressure of the ionised gas, and the pressure of the shock-heated winds. We find that stellar winds and ionised gas are the main drivers of HII region expansion in our sample, while the direct radiation pressure is up to three orders of magnitude lower than the other terms. We relate the total pressure to the star formation rate per unit area, $\Sigma_{SFR}$, for each region and find that stellar feedback has a negative effect on star formation, and sets an upper limit to $\Sigma_{SFR}$ as a function of increasing pressure.

We present an analysis of gas densities in the central R=300 parsecs of the Milky Way, focusing on three clouds: GCM-0.02-0.07 (the 50 km/s cloud), GCM-0.13-0.08 (the 20 km/s cloud), and GCM0.25+0.01 (the "Brick"). Densities are determined using observations of the J=(3-2), (4-3), (5-4), (10-9), (18-17), (19-18), (21-20), and (24-23) transitions of the molecule HC3N. We find evidence of at least two excitation regimes for HC3N and constrain the low-excitation component to have a density less than 10^4 cm^-3 and the high-excitation component to have a density between 10^5 and 10^6 cm^-3. This is much less than densities of 10^7 cm^-3 that are found in Sgr B2, the most actively star-forming cloud in the Galactic center. This is consistent with the requirement of a higher density threshold for star formation in the Galactic center than is typical in the Galactic disk. We are also able to constrain the column density of each component in order to determine the mass fraction of 'dense' (n>10^5 cm^-3) gas for these clouds. We find that this is ~15% for all three clouds. Applying the results of our models to ratios of the (10-9) and (3-2) line across the entire central R=300 pc, we find that the fraction of gas with n>10^4 cm^-3 increases inward of a radius of ~140 pc, consistent with the predictions of recent models for the gas dynamics in this region. Our observations show that HC3N is an excellent molecule for probing the density structure of clouds in the Galactic center.

In this contribution we explore the new science that a Next Generation Very Large Array (ngVLA) would be able to perform on the topics of the youngest HII regions and (proto)stellar jets. Free-free continuum and radio recombination line (RRL) emission are often the only way of peering into the dense envelopes surrounding (proto)stars of all masses, and trace their initial feedback in the form of `radio jets', `hypercompact HII regions', or photoevaporating, partially-ionized flows. Properly disentangling free-free from dust emission is also mandatory in studies of protoplanetary and accretion disks. Current VLA research has reached an impasse in which a population of faint ionized radio sources, probably corresponding to the above mentioned objects, is detected, but their nature is mostly unknown. The ngVLA would allow us to resolve the density structure and kinematics of such sources, revolutionizing our knowledge of star formation across the entire stellar-mass spectrum.

Understanding the formation of massive stars is one of the unsolved problems in modern astronomy. The main difficulty is that the intense radiation from the high-luminosity stars and the thermal pressure from the resulting ionized gas (both insignificant for low-mass stars) may be able to reverse the accretion flow and prevent the star from accreting fresh material. Such feedback effects can naturally be mitigated if accretion proceeds through discs, which is the established mechanism to form sun-like stars. However, recent 3D MHD simulations have shown that accretion on 1000 au scales is through filaments rather than a large disc. This theoretical prediction has never been confirmed via observations owing to the poor linear resolution of previous studies (>1000 au). Here we present the first observational evidence that mass assembly in young high-mass stars forming in protoclusters is predominantly asymmetric and disordered. In particular, we observed the innermost regions around three deeply embedded high-mass protostars with very high spatial resolution (~100 au). We identified multiple massive (several solar masses), warm (50-150 Kelvin) filamentary streamers pointing onto the central sources, which we interpret as multi-directional accretion channels. These structures inhibit the formation of a large, steady disc. Nevertheless, the identification of fast collimated outflows in the three observed systems indicates that (non-steady) compact discs may be present (we measure upper limits on their radii of <80 for one object and <350 astronomical units for the remaining two objects). Our finding contrasts with the simplified classic paradigm of an ordered (and stable) disc/jet system and provides an experimental confirmation of a multi-directional and unsteady accretion model for massive star formation supported by recent 3D (magneto)hydrodynamic simulations.

We report ALMA long-baseline observations of Orion Source I (SrcI) with resolution 0.03-0.06" (12-24 AU) at 1.3 and 3.2 mm. We detect both continuum and spectral line emission from SrcI's disk. We also detect a central weakly resolved source that we interpret as a hot spot in the inner disk, which may indicate the presence of a binary system. The high angular resolution and sensitivity of these observations allow us to measure the outer envelope of the rotation curve of the H$_2$O $5_{5,0}-6_{4,3}$ line, which gives a mass $M_I\approx15\pm2$ Msun. We detected several other lines that more closely trace the disk, but were unable to identify their parent species. Using centroid-of-channel methods on these other lines, we infer a similar mass. These measurements solidify SrcI as a genuine high-mass protostar system and support the theory that SrcI and the Becklin Neugebauer Object were ejected from the dynamical decay of a multiple star system $\sim$500 years ago, an event that also launched the explosive molecular outflow in Orion.

We present the first results from a new, high resolution, $^{12}$CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) molecular line survey of the Orion A cloud, hereafter referred to as the CARMA-NRO Orion Survey. CARMA observations have been combined with single-dish data from the Nobeyama 45m telescope to provide extended images at about 0.01 pc resolution, with a dynamic range of approximately 1200 in spatial scale. Here we describe the practical details of the data combination in uv space, including flux scale matching, the conversion of single dish data to visibilities, and joint deconvolution of single dish and interferometric data. A $\Delta$-variance analysis indicates that no artifacts are caused by combining data from the two instruments. Initial analysis of the data cubes, including moment maps, average spectra, channel maps, position-velocity diagrams, excitation temperature, column density, and line ratio maps provides evidence of complex and interesting structures such as filaments, bipolar outflows, shells, bubbles, and photo-eroded pillars. The implications for star formation processes are profound and follow-up scientific studies by the CARMA-NRO Orion team are now underway. We plan to make all the data products described here generally accessible; some are already available at https://dataverse.harvard.edu/dataverse/CARMA-NRO-Orion

We measure H$_2$ temperatures and column densities across the Orion BN/KL explosive outflow from a set of thirteen near-IR H$_2$ rovibrational emission lines observed with the TripleSpec spectrograph on Apache Point Observatory's 3.5-meter telescope. We find that most of the region is well-characterized by a single temperature (~2000-2500 K), which may be influenced by the limited range of upper energy levels (6000-20,000 K) probed by our data set. The H$_2$ column density maps indicate that warm H$_2$ comprises 10$^{-5}$ - 10$^{-3}$ of the total H$_2$ column density near the center of the outflow. Combining column density measurements for co-spatial H$_2$ and CO at T = 2500 K, we measure a CO/H$_2$ fractional abundance of 2$\times$10$^{-3}$, and discuss possible reasons why this value is in excess of the canonical 10$^{-4}$ value, including dust attenuation, incorrect assumptions on co-spatiality of the H$_2$ and CO emission, and chemical processing in an extreme environment. We model the radiative transfer of H$_2$ in this region with UV pumping models to look for signatures of H$_2$ fluorescence from H I Ly$\alpha$ pumping. Dissociative (J-type) shocks and nebular emission from the foreground Orion H II region are considered as possible Ly$\alpha$ sources. From our radiative transfer models, we predict that signatures of Ly$\alpha$ pumping should be detectable in near-IR line ratios given a sufficiently strong source, but such a source is not present in the BN/KL outflow. The data are consistent with shocks as the H$_2$ heating source.

The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.

We report ALMA observations with resolution $\approx0.5$" at 3 mm of the extended Sgr B2 cloud in the Central Molecular Zone (CMZ). We detect 271 compact sources, most of which are smaller than 5000 AU. By ruling out alternative possibilities, we conclude that these sources consist of a mix of hypercompact HII regions and young stellar objects (YSOs). Most of the newly-detected sources are YSOs with gas envelopes which, based on their luminosities, must contain objects with stellar masses $M_*\gtrsim8$ M$_\odot$. Their spatial distribution spread over a $\sim12\times3$ pc region demonstrates that Sgr B2 is experiencing an extended star formation event, not just an isolated `starburst' within the protocluster regions. Using this new sample, we examine star formation thresholds and surface density relations in Sgr B2. While all of the YSOs reside in regions of high column density ($N(H_2)\gtrsim2\times10^{23}$ cm$^{-2}$), not all regions of high column density contain YSOs. The observed column density threshold for star formation is substantially higher than that in solar vicinity clouds, implying either that high-mass star formation requires a higher column density or that any star formation threshold in the CMZ must be higher than in nearby clouds. The relation between the surface density of gas and stars is incompatible with extrapolations from local clouds, and instead stellar densities in Sgr B2 follow a linear $\Sigma_*-\Sigma_{gas}$ relation, shallower than that observed in local clouds. Together, these points suggest that a higher volume density threshold is required to explain star formation in CMZ clouds.

We present Atacama Large Millimeter Array (ALMA) observations of $^{12}$CO2$-$1 emission from circumstellar material around the massive star $\eta$~Carinae. These observations reveal new structural details about the cool equatorial torus located $\sim$4000 au from the star. The CO torus is not a complete azimuthal loop, but rather, is missing its near side, which appears to have been cleared away. The missing material matches the direction of apastron in the eccentric binary system, making it likely that $\eta$~Car's companion played an important role in disrupting portions of the torus soon after ejection. Molecular gas seen in ALMA data aligns well with the cool dust around $\eta$~Car previously observed in mid-infrared (IR) maps, whereas hot dust resides at the inner surface of the molecular torus. The CO also coincides with the spatial and velocity structure of near-IR H$_2$ emission. Together, these suggest that the CO torus seen by ALMA is actually the pinched waist of the Homunculus polar lobes, which glows brightly because it is close to the star and warmer than the poles. The near side of the torus appears to be a blowout, associated with fragmented equatorial ejecta. We discuss implications for the origin of various features northwest of the star. CO emission from the main torus implies a total gas mass in the range of 0.2-1 $M_{\odot}$ (possibly up to 5 $M_{\odot}$ or more, although with questionable assumptions). Deeper observations are needed to constrain CO emission from the cool polar lobes.

Understanding the early stages of star formation is a research field of ongoing development, both theoretically and observationally. In this context, molecular data have been continuously providing observational constraints on the gas dynamics at different excitation conditions and depths in the sources. We have investigated the Barnard 59 core, the only active site of star formation in the Pipe Nebula, to achieve a comprehensive view of the kinematic properties of the source. These information were derived by simultaneously fitting ammonia inversion transition lines (1,1) and (2,2). Our analysis unveils the imprint of protostellar feedback, such as increasing line widths, temperature and turbulent motions in our molecular data. Combined with complementary observations of dust thermal emission, we estimate that the core is gravitationally bound following a virial analysis. If the core is not contracting, another source of internal pressure, most likely the magnetic field, is supporting it against gravitational collapse and limits its star formation efficiency.

The Galactic Center contains large amounts of molecular and ionized gas as well as a plethora of energetic objects. Water masers are an extinction-insensitive probe for star formation and thus ideal for studies of star formation stages in this highly obscured region. With the Australia Telescope Compact Array, we observed 22 GHz water masers in the entire Central Molecular Zone with sub-parsec resolution as part of the large SWAG survey: ``Survey of Water and Ammonia in the Galactic Center''. We detect of order 600 22 GHz masers with isotropic luminosities down to ~10^-7 Lo. Masers with luminosities of >~10^-6 Lo are likely associated with young stellar objects. They appear to be close to molecular gas streamers and may be due to star formation events that are triggered at pericenter passages near Sgr A*. Weaker masers are more widely distributed and frequently show double line features, a tell-tale sign for an origin in evolved star envelopes.

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.

The Survey of Water and Ammonia in the Galactic Center (SWAG) covers the Central Molecular Zone (CMZ) of the Milky Way at frequencies between 21.2 and 25.4 GHz obtained at the Australia Telescope Compact Array at $\sim 0.9$ pc spatial and $\sim 2.0$ km s$^{-1}$ spectral resolution. In this paper, we present data on the inner $\sim 250$ pc ($1.4^\circ$) between Sgr C and Sgr B2. We focus on the hyperfine structure of the metastable ammonia inversion lines (J,K) = (1,1) - (6,6) to derive column density, kinematics, opacity and kinetic gas temperature. In the CMZ molecular clouds, we find typical line widths of $8-16$ km s$^{-1}$ and extended regions of optically thick ($\tau > 1$) emission. Two components in kinetic temperature are detected at $25-50$ K and $60-100$ K, both being significantly hotter than dust temperatures throughout the CMZ. We discuss the physical state of the CMZ gas as traced by ammonia in the context of the orbital model by Kruijssen et al. (2015) that interprets the observed distribution as a stream of molecular clouds following an open eccentric orbit. This allows us to statistically investigate the time dependencies of gas temperature, column density and line width. We find heating rates between $\sim 50$ and $\sim 100$ K Myr$^{-1}$ along the stream orbit. No strong signs of time dependence are found for column density or line width. These quantities are likely dominated by cloud-to-cloud variations. Our results qualitatively match the predictions of the current model of tidal triggering of cloud collapse, orbital kinematics and the observation of an evolutionary sequence of increasing star formation activity with orbital phase.

We use Green Bank Ammonia Survey observations of NH$_3$ (1,1) and (2,2) emission with 32'' FWHM resolution from a ~ 10 pc$^{2}$ portion of the Cepheus-L1251 molecular cloud to identify hierarchical dense gas structures. Our dendrogram analysis of the NH$_3$ data results in 22 top-level structures, which reside within 13 lower-level, parent structures. The structures are compact (0.01 pc $\lesssim R_{eff} \lesssim$ 0.1 pc) and are spatially correlated with the highest H$_2$ column density portions of the cloud. We also compare the ammonia data to a catalog of dense cores identified by higher-resolution (18.2'' FWHM) Herschel Space Observatory observations of dust continuum emission from Cepheus-L1251. Maps of kinetic gas temperature, velocity dispersion, and NH$_3$ column density, derived from detailed modeling of the NH$_3$ data, are used to investigate the stability and chemistry of the ammonia-identified and Herschel-identified structures. We show that the dust and dense gas in the structures have similar temperatures, with median $T_{dust}$ and $T_K$ measurements of 11.7 $\pm$ 1.1 K and 10.3 $\pm$ 2.0 K, respectively. Based on a virial analysis, we find that the ammonia-identified structures are gravitationally dominated, yet may be in or near a state of virial equilibrium. Meanwhile, the majority of the Herschel-identified dense cores appear to be not bound by their own gravity and instead confined by external pressure. CCS $(2_0-1_0)$ and HC$_5$N $(9-8)$ emission from the region reveal broader line widths and centroid velocity offsets when compared to the NH$_3$ (1,1) emission in some cases, likely due to these carbon-based molecules tracing the turbulent outer layers of the dense cores.

STATCONT is a python-based tool designed to determine the continuum emission level in spectral data, in particular for sources with a line-rich spectrum. The tool inspects the intensity distribution of a given spectrum and automatically determines the continuum level by using different statistical approaches. The different methods included in STATCONT are tested against synthetic data. We conclude that the sigma-clipping algorithm provides the most accurate continuum level determination, together with information on the uncertainty in its determination. This uncertainty can be used to correct the final continuum emission level, resulting in the here called `corrected sigma-clipping method' or c-SCM. The c-SCM has been tested against more than 750 different synthetic spectra reproducing typical conditions found towards astronomical sources. The continuum level is determined with a discrepancy of less than 1% in 50% of the cases, and less than 5% in 90% of the cases, provided at least 10% of the channels are line free. The main products of STATCONT are the continuum emission level, together with a conservative value of its uncertainty, and datacubes containing only spectral line emission, i.e., continuum-subtracted datacubes. STATCONT also includes the option to estimate the spectral index, when different files covering different frequency ranges are provided.

We use gas temperature and velocity dispersion data from the Green Bank Ammonia Survey and core masses and sizes from the James Clerk Maxwell Telescope Gould Belt Survey to estimate the virial states of dense cores within the Orion A molecular cloud. Surprisingly, we find that almost none of the dense cores are sufficiently massive to be bound when considering only the balance between self-gravity and the thermal and non-thermal motions present in the dense gas. Including the additional pressure binding imposed by the weight of the ambient molecular cloud material and additional smaller pressure terms, however, suggests that most of the dense cores are pressure confined.

High-mass stars have generally been assumed to accrete most of their mass while already contracted onto the main sequence, but this hypothesis has not been observationally tested. We present ALMA observations of a 3 x 1.5 pc area in the W51 high-mass star-forming complex. We identify dust continuum sources and measure the gas and dust temperature through both rotational diagram modeling of CH3OH and brightness-temperature-based limits. The observed region contains three high-mass YSOs that appear to be at the earliest stages of their formation, with no signs of ionizing radiation from their central sources. The data reveal high gas and dust temperatures (T > 100 K) extending out to about 5000 AU from each of these sources. There are no clear signs of disks or rotating structures down to our 1000 AU resolution. The extended warm gas provides evidence that, during the process of forming, these high-mass stars heat a large volume and correspondingly large mass of gas in their surroundings, inhibiting fragmentation and therefore keeping a large reservoir available to feed from. By contrast, the more mature massive stars that illuminate compact H ii regions have little effect on their surrounding dense gas, suggesting that these main sequence stars have completed most or all of their accretion. The high luminosity of the massive protostars (L > 10^4 solar luminosities), combined with a lack of centimeter continuum emission from these sources, implies that they are not on the main sequence while they accrete the majority of their mass; instead, they may be bloated and cool.

We have modified the iterative procedure introduced by Lin et al. (2016), to systematically combine the submm images taken from ground based (e.g., CSO, JCMT, APEX) and space (e.g., Herschel, Planck) telescopes. We applied the updated procedure to observations of three well studied Infrared Dark Clouds (IRDCs): G11.11-0.12, G14.225-0.506 and G28.34+0.06, and then performed single-component, modified black-body fits to derive $\sim$10$"$ resolution dust temperature and column density maps. The derived column density maps show that these three IRDCs exhibit complex filamentary structures embedding with rich clumps/cores. We compared the column density probability distribution functions (N-PDFs) and two-point correlation (2PT) functions of the column density field between these IRDCs with several OB cluster-forming regions. Based on the observed correlation and measurements, and complementary hydrodynamical simulations for a 10$^{4}$ $\rm M_{\odot}$ molecular cloud, we hypothesize that cloud evolution can be better characterized by the evolution of the (column) density distribution function and the relative power of dense structures as a function of spatial scales, rather than merely based on the presence of star-forming activity. Based on the small analyzed sample, we propose four evolutionary stages, namely: {\it cloud integration, stellar assembly, cloud pre-dispersal and dispersed-cloud.} The initial {\it cloud integration} stage and the final {\it dispersed cloud} stage may be distinguished from the two intermediate stages by a steeper than $-$4 power-law index of the N-PDF. The {\it cloud integration} stage and the subsequent {\it stellar assembly} stage are further distinguished from each other by the larger luminosity-to-mass ratio ($>$40 $\rm L_{\odot}/M_{\odot}$) of the latter.

We present an overview of the first data release (DR1) and first-look science from the Green Bank Ammonia Survey (GAS). GAS is a Large Program at the Green Bank Telescope to map all Gould Belt star-forming regions with $A_V \gtrsim 7$ mag visible from the northern hemisphere in emission from NH$_3$ and other key molecular tracers. This first release includes the data for four regions in Gould Belt clouds: B18 in Taurus, NGC 1333 in Perseus, L1688 in Ophiuchus, and Orion A North in Orion. We compare the NH$_3$ emission to dust continuum emission from Herschel, and find that the two tracers correspond closely. NH$_3$ is present in over 60\% of lines-of-sight with $A_V \gtrsim 7$ mag in three of the four DR1 regions, in agreement with expectations from previous observations. The sole exception is B18, where NH$_3$ is detected toward ~ 40\% of lines-of-sight with $A_V \gtrsim 7$ mag. Moreover, we find that the NH$_3$ emission is generally extended beyond the typical 0.1 pc length scales of dense cores. We produce maps of the gas kinematics, temperature, and NH$_3$ column densities through forward modeling of the hyperfine structure of the NH$_3$ (1,1) and (2,2) lines. We show that the NH$_3$ velocity dispersion, ${\sigma}_v$, and gas kinetic temperature, $T_K$, vary systematically between the regions included in this release, with an increase in both the mean value and spread of ${\sigma}_v$ and $T_K$ with increasing star formation activity. The data presented in this paper are publicly available.

The high-mass star forming sites SgrB2(M) and SgrB2(N) have been the target of numerous studies, revealing e.g. a rich chemistry. We want to characterize their physical and chemical structure using ALMA high-angular resolution observations at mm wavelengths, reaching spatial scales of about 4000 au, and covering the whole band 6 (from 211 to 275 GHz). In order to determine the continuum emission in line-rich sources, we use a new statistical method: STATCONT. We detect 27 continuum sources in SgrB2(M) and 20 in SgrB2(N). We study the continuum emission across the ALMA band 6, and compare it with previous SMA 345 GHz and VLA 40 GHz observations, to study the nature of the sources detected. The brightest sources are dominated by (partially optically thick) dust emission, while there is an important degree of contamination from ionized gas free-free emission in weaker sources. While the total mass in SgrB2(M) is distributed in many fragments, most of the mass in SgrB2(N) arises from a single object, with filamentary-like structures converging towards the center. There seems to be a lack of low-mass dense cores in both regions. We determine H2 volume densities for the cores of about 10^5-10^7 Msun pc^-3, one to two orders of magnitude higher than the stellar densities of super star clusters. In general, SgrB2(N) is chemically richer than SgrB2(M). There seems to be a correlation between the chemical richness and the mass of the fragments, with more massive clumps being more chemically rich. Both SgrB2(N) and SgrB2(M) harbour a cluster of hot molecular cores. We compare the continuum images with predictions from a detailed 3D radiative transfer model that reproduces the structure of SgrB2 from 45 pc down to 100 au. This dataset, together with ongoing projects in the range 5 to 200 GHz, better constrain the 3D structure of SgrB2, and allow us to understand its physical and chemical structure.

We recently reported a population of protostellar candidates in the 20 km s$^{-1}$ cloud in the Central Molecular Zone of the Milky Way, traced by H$_2$O masers in gravitationally bound dense cores. In this paper, we report high-angular-resolution ($\sim$3'') molecular line studies of the environment of star formation in this cloud. Maps of various molecular line transitions as well as the continuum at 1.3 mm are obtained using the Submillimeter Array. Five NH$_3$ inversion lines and the 1.3 cm continuum are observed with the Karl G. Jansky Very Large Array. The interferometric observations are complemented with single-dish data. We find that the CH$_3$OH, SO, and HNCO lines, which are usually shock tracers, are better correlated spatially with the compact dust emission from dense cores among the detected lines. These lines also show enhancement in intensities with respect to SiO intensities toward the compact dust emission, suggesting the presence of slow shocks or hot cores in these regions. We find gas temperatures of $\gtrsim$100 K at 0.1-pc scales based on RADEX modelling of the H$_2$CO and NH$_3$ lines. Although no strong correlations between temperatures and linewidths/H$_2$O maser luminosities are found, in high-angular-resolution maps we notice several candidate shock heated regions offset from any dense cores, as well as signatures of localized heating by protostars in several dense cores. Our findings suggest that at 0.1-pc scales in this cloud star formation and strong turbulence may together affect the chemistry and temperature of the molecular gas.

The W51 cloud complex is one of the best laboratories in our Galaxy to study high-mass star formation. At a distance of about 5 kpc, it is the closest region containing a high-mass protocluster, and it has two. The cloud includes a long infrared-dark cloud, is interacting with a supernova remnant, and contains a variety of unique massive protostellar sources. This article is an observational review of the region.

Most massive stars form in dense clusters where gravitational interactions with other stars may be common. The two nearest forming massive stars, the BN object and Source I, located behind the Orion Nebula, were ejected with velocities of $\sim$29 and $\sim$13 km s$^{-1}$ about 500 years ago by such interactions. This event generated an explosion in the gas. New ALMA observations show in unprecedented detail, a roughly spherically symmetric distribution of over a hundred $^{12}$CO J=2$-$1 streamers with velocities extending from V$_{LSR}$ =$-$150 to +145 km s$^{-1}$. The streamer radial velocities increase (or decrease) linearly with projected distance from the explosion center, forming a `Hubble Flow' confined to within 50 arcseconds of the explosion center. They point toward the high proper-motion, shock-excited H$_2$ and [Fe ii ] `fingertips' and lower-velocity CO in the H$_2$ wakes comprising Orion's `fingers'. In some directions, the H$_2$ `fingers' extend more than a factor of two farther from the ejection center than the CO streamers. Such deviations from spherical symmetry may be caused by ejecta running into dense gas or the dynamics of the N-body interaction that ejected the stars and produced the explosion. This $\sim$10$^{48}$ erg event may have been powered by the release of gravitational potential energy associated with the formation of a compact binary or a protostellar merger. Orion may be the prototype for a new class of stellar explosion responsible for luminous infrared transients in nearby galaxies.

The origin and life-cycle of molecular clouds are still poorly constrained, despite their importance for understanding the evolution of the interstellar medium. We have carried out a systematic, homogeneous, spectroscopic survey of the inner Galactic plane, in order to complement the many continuum Galactic surveys available with crucial distance and gas-kinematic information. Our aim is to combine this data set with recent infrared to sub-millimetre surveys at similar angular resolutions. The SEDIGISM survey covers 78 deg^2 of the inner Galaxy (-60 deg < l < +18 deg, |b| < 0.5 deg) in the J=2-1 rotational transition of 13CO. This isotopologue of CO is less abundant than 12CO by factors up to 100. Therefore, its emission has low to moderate optical depths, and higher critical density, making it an ideal tracer of the cold, dense interstellar medium. The data have been observed with the SHFI single-pixel instrument at APEX. The observational setup covers the 13CO(2-1) and C18O(2-1) lines, plus several transitions from other molecules. The observations have been completed. Data reduction is in progress, and the final data products will be made available in the near future. Here we give a detailed description of the survey and the dedicated data reduction pipeline. Preliminary results based on a science demonstration field covering -20 deg < l < -18.5 deg are presented. Analysis of the 13CO(2-1) data in this field reveals compact clumps, diffuse clouds, and filamentary structures at a range of heliocentric distances. By combining our data with data in the (1-0) transition of CO isotopologues from the ThrUMMS survey, we are able to compute a 3D realization of the excitation temperature and optical depth in the interstellar medium. Ultimately, this survey will provide a detailed, global view of the inner Galactic interstellar medium at an unprecedented angular resolution of ~30".

We report a new analysis protocol for HCN hyperfine data, based on the PYSPECKIT package, and results of using this new protocol to analyse a sample area of seven massive molecular clumps from the Census of High- and Medium-mass Protostars (CHaMP) survey, in order to derive maps of column density for this species. There is a strong correlation between the HCN integrated intensity, IHCN, and previously reported IHCO+ in the clumps, but IN2H+ is not well correlated with either of these other two "dense gas tracers". The four fitted parameters from PYSPECKIT in this region fall in the range of VLSR = 8-10 km/s, {\sigma} V = 1.2-2.2 km/s, Tex = 4-15 K, and {\tau} = 0.2-2.5. These parameters allow us to derive a column density map of these clouds, without limiting assumptions about the excitation or opacity. A more traditional (linear) method of converting IHCN to total mass column gives much lower clump masses than our results based on the hyperfine analysis. This is primarily due to areas in the sample region of low I, low Tex, and high {\tau} . We conclude that there may be more dense gas in these massive clumps not engaged in massive star formation than previously recognized. If this result holds for other clouds in the CHaMP sample, it would have dramatic consequences for the calibration of the Kennicutt-Schmidt star formation laws, including a large increase in the gas depletion time-scale in such regions.

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.

We present an analysis of the effect of feedback from O- and B-type stars with data from the integral field spectrograph MUSE mounted on the Very Large Telescope of pillar-like structures in the Carina Nebular Complex, one of the most massive star-forming regions in the Galaxy. For the observed pillars, we compute gas electron densities and temperatures maps, produce integrated line and velocity maps of the ionised gas, study the ionisation fronts at the pillar tips, analyse the properties of the single regions, and detect two ionised jets originating from two distinct pillar tips. For each pillar tip we determine the incident ionising photon flux $Q_\mathrm{0,pil}$ originating from the nearby massive O- and B-type stars and compute the mass-loss rate $\dot{M}$ of the pillar tips due to photo-evaporation caused by the incident ionising radiation. We combine the results of the Carina data set with archival MUSE data of a pillar in NGC 3603 and with previously published MUSE data of the Pillars of Creation in M16, and with a total of 10 analysed pillars, find tight correlations between the ionising photon flux and the electron density, the electron density and the distance from the ionising sources, and the ionising photon flux and the mass-loss rate. The combined MUSE data sets of pillars in regions with different physical conditions and stellar content therefore yield an empirical quantification of the feedback effects of ionising radiation. In agreement with models, we find that $\dot{M}\propto Q_\mathrm{0,pil}^{1/2}$.

The Central Molecular Zone (CMZ) at the center of our Galaxy is the best template to study star formation processes under extreme conditions, similar to those in high-redshift galaxies. We observed on-the-fly maps of para-H$_{2}$CO transitions at 218 GHz and 291 GHz towards seven Galactic Center clouds. From the temperature-sensitive integrated intensity line ratios of H$_{2}$CO(3$_{2,1}-$2$_{2,0}$)/H$_{2}$CO(3$_{0,3}-$2$_{0,2}$) and H$_{2}$CO(4$_{2,2}-$3$_{2,1}$)/H$_{2}$CO(4$_{0,4}-$3$_{0,3}$) in combination with radiative transfer models, we produce gas temperature maps of our targets. These transitions are sensitive to gas with densities of $\sim$10$^{5}$ cm$^{-3}$ and temperatures <150 K. The measured gas temperatures in our sources are all higher (>40 K) than their dust temperatures ($\sim$25 K). Our targets have a complex velocity structure that requires a careful disentanglement of the different components. We produce temperature maps for each of the velocity components and show that the temperatures of the components differ, revealing temperature gradients in the clouds. Combining the temperature measurements with the integrated intensity line ratio of H$_{2}$CO(4$_{0,4}-$3$_{0,3}$)/H$_{2}$CO(3$_{0,3}-$2$_{0,2}$), we constrain the density of this warm gas to 10$^{4}-$10$^{6}$ cm$^{-3}$. We find a positive correlation of the line width of the main H$_{2}$CO lines with the temperature of the gas, direct evidence for gas heating via turbulence. Our data is consistent with a turbulence heating model with a density of n = 10$^5$ cm$^{-3}$.

We have developed an iterative procedure to systematically combine the millimeter and submillimeter images of OB cluster-forming molecular clouds, which were taken by ground based (CSO, JCMT, APEX, IRAM-30m) and space telescopes (Herschel, Planck). For the seven luminous ($L$$>$10$^{6}$ $L_{\odot}$) Galactic OB cluster-forming molecular clouds selected for our analyses, namely W49A, W43-Main, W43-South, W33, G10.6-0.4, G10.2-0.3, G10.3-0.1, we have performed single-component, modified black-body fits to each pixel of the combined (sub)millimeter images, and the Herschel PACS and SPIRE images at shorter wavelengths. The $\sim$10$"$ resolution dust column density and temperature maps of these sources revealed dramatically different morphologies, indicating very different modes of OB cluster-formation, or parent molecular cloud structures in different evolutionary stages. The molecular clouds W49A, W33, and G10.6-0.4 show centrally concentrated massive molecular clumps that are connected with approximately radially orientated molecular gas filaments. The W43-Main and W43-South molecular cloud complexes, which are located at the intersection of the Galactic near 3-kpc (or Scutum) arm and the Galactic bar, show a widely scattered distribution of dense molecular clumps/cores over the observed $\sim$10 pc spatial scale. The relatively evolved sources G10.2-0.3 and G10.3-0.1 appear to be affected by stellar feedback, and show a complicated cloud morphology embedded with abundant dense molecular clumps/cores. We find that with the high angular resolution we achieved, our visual classification of cloud morphology can be linked to the systematically derived statistical quantities (i.e., the enclosed mass profile, the column density probability distribution function, the two-point correlation function of column density, and the probability distribution function of clump/core separations).

We present a line survey of the ultra-luminous infrared galaxy Arp 220, taken with the newly installed SEPIA Band 5 instrument on APEX. We illustrate the capacity of SEPIA to detect the 183.3 GHz H2O 31,3-22,0 line against the atmospheric H2O absorption feature. We confirm the previous detection of the HCN(2-1) line, and detect new transitions of standard dense gas tracers such as HNC(2-1), HCO+(2-1), CS(4-3), C34S(4-3), HC3N(20-19). We also detect HCN(2-1) v2=1 and the 193.5 GHz methanol (4-3) group for the first time. The absence of time variations in the megamaser water line compared to previous observations seems to rule out an AGN nuclear origin for the line. It could, on the contrary, favor a thermal origin instead, but also possibly be a sign that the megamaser emission is associated with star-forming cores washed-out in the beam. We finally discuss how the new transitions of HCN, HNC, HCO+ refine our knowledge of the ISM physical conditions in Arp 220.

We present new JVLA observations of the high-mass cluster-forming region W51A from 2 to 16 GHz with resolution ${\theta}_{fwhm} \approx$ 0.3 - 0.5". The data reveal a wealth of observational results: (1) Currently-forming, very massive (proto-O) stars are traced by o-H2CO $2_{1,1}-2_{1,2}$ emission, suggesting that this line can be used efficiently as a massive protostar tracer. (2) There is a spatially distributed population of $\sim$mJy continuum sources, including hypercompact H ii regions and candidate colliding wind binaries, in and around the W51 proto-clusters. (3) There are two clearly detected protoclusters, W51e and W51 IRS2, that are gas-rich but may have most of their mass in stars within their inner $\sim$ 0.05 pc. The majority of the bolometric luminosity in W51 most likely comes from a third population of OB stars between these clusters. The presence of a substantial population of exposed O-stars coincident with a population of still-forming massive stars, along with a direct measurement of the low mass loss rate via ionized gas outflow from W51 IRS2, together imply that feedback is ineffective at halting star formation in massive protoclusters. Instead, feedback may shut off the large-scale accretion of diffuse gas onto the W51 protoclusters, implying that they are evolving towards a state of gas exhaustion rather than gas expulsion. Recent theoretical models predict gas exhaustion to be a necessary step in the formation of gravitationally bound stellar clusters, and our results provide an observational validation of this process.

We present near-IR (1.1-2.4 micron) position-position-velocity cubes of the 500-year-old Orion BN/KL explosive outflow with spatial resolution 1" and spectral resolution 86 km/s. We construct integrated intensity maps free of continuum sources of 15 H2 and [Fe II] lines while preserving kinematic information of individual outflow features. Included in the detected H2 lines are the 1-0 S(1) and 1-0 Q(3) transitions, allowing extinction measurements across the outflow. Additionally, we present dereddened flux ratios for over two dozen outflow features to allow for the characterization of the true excitation conditions of the BN/KL outflow. All ratios show the dominance of shock excitation of the H2 emission, although some features exhibit signs of fluorescent excitation from stellar radiation or J-type shocks. We also detect tracers of the PDR/ionization front north of the Trapezium stars in [O I] and [Fe II] and analyze other observed outflows not associated with the BN/KL outflow.

We present ALMA observations of the Orion Nebula that cover the OMC1 outflow region. Our focus in this paper is on compact emission from protoplanetary disks. We mosaicked a field containing $\sim 600$ near-IR-identified young stars, around which we can search for sub-mm emission tracing dusty disks. Approximately 100 sources are known proplyds identified with HST. We detect continuum emission at 1 mm wavelengths towards $\sim 20\%$ of the proplyd sample, and $\sim 8\%$ of the larger sample of near-IR objects. The noise in our maps allows 4$\sigma$ detection of objects brighter than $\sim 1.5$ mJy, corresponding to protoplanetary disk masses larger than 1.5 M$_{\rm J}$ (using standard assumptions about dust opacities and gas-to-dust ratios). None of these disks are detected in contemporaneous CO(2-1) or C$^{18}$O(2-1) observations, suggesting that the gas-to-dust ratios may be substantially smaller than the canonical value of 100. Furthermore, since dust grains may already be sequestered in large bodies in ONC disks, the inferred masses of disk solids may be underestimated. Our results suggest that the distribution of disk masses in this region is compatible with the detection rate of massive planets around M dwarfs, which are the dominant stellar constituent in the ONC.

Using spectral-line observations of HNCO, N2H+, and HNC, we investigate the kinematics of dense gas in the central ~250 pc of the Galaxy. We present SCOUSE (Semi-automated multi-COmponent Universal Spectral-line fitting Engine), a line fitting algorithm designed to analyse large volumes of spectral-line data efficiently and systematically. Unlike techniques which do not account for complex line profiles, SCOUSE accurately describes the {l, b, v_LSR} distribution of CMZ gas, which is asymmetric about Sgr A* in both position and velocity. Velocity dispersions range from 2.6 km/s<\sigma<53.1 km/s. A median dispersion of 9.8 km/s, translates to a Mach number, M_3D>28. The gas is distributed throughout several "streams", with projected lengths ~100-250 pc. We link the streams to individual clouds and sub-regions, including Sgr C, the 20 and 50 km/s clouds, the dust ridge, and Sgr B2. Shell-like emission features can be explained by the projection of independent molecular clouds in Sgr C and the newly identified conical profile of Sgr B2 in {l ,b, v_LSR} space. These features have previously invoked supernova-driven shells and cloud-cloud collisions as explanations. We instead caution against structure identification in velocity-integrated emission maps. Three geometries describing the 3-D structure of the CMZ are investigated: i) two spiral arms; ii) a closed elliptical orbit; iii) an open stream. While two spiral arms and an open stream qualitatively reproduce the gas distribution, the most recent parameterisation of the closed elliptical orbit does not. Finally, we discuss how proper motion measurements of masers can distinguish between these geometries, and suggest that this effort should be focused on the 20 km/s and 50 km/s clouds and Sgr C.

This paper is the third in a series of ammonia multilevel imaging studies in well-known high-mass star forming regions. Using the JVLA, we have mapped the hot and dense molecular gas in W51 Main, with about 0.2 arcsec angular resolution, in five highly-excited metastable inversion transitions of ammonia (ammonia): (J,K)=(6,6), (7,7), (9,9), (10,10), and (13,13). We have identified and characterised two main centers of high-mass star formation in W51-Main: the W51e2 complex and the W51e8 core (6 arcsec southward of W51e2). The former breaks down into three further sub-cores: W51e2-W, which surrounds the well known HC HII region, where hot ammonia is observed in absorption, and two additional dusty cores, W51e2-E (~0.8" to the East) and W51e2-NW (~1" to the North), where hot ammonia is observed in emission. The velocity maps towards the HC HII region show a clear velocity gradient that may indicate rotation, though we do not directly observe a Keplerian velocity profile. The absence of outflow and/or maser activity and the low amount of molecular gas available for accretion (~5 solar masses) with respect to the mass of the central YSO (>20 solar masses), both indicate that the central YSO has already accreted most of its final mass. On the other hand, in the nearby W51e2-E object, the relatively large amount of hot molecular gas available for accretion (~20 solar masses, within about half an arcsecond or 2500 AU), along with strong outflow and maser activity, indicates that the main accretion center in the W51e2 complex is W51e2-E rather than W51e2-W. Finally, W51e2-NW and W51e8, although less dense, are also hot cores and contain a significant amount of molecular gas (~30 and 70 solar masses, respectively). We speculate that they may host high-mass YSOs either at a previous evolutionary stage or with lower mass than W51e2-E and W51e2-W.

The formation environment of stars in massive stellar clusters is similar to the environment of stars forming in galaxies at a redshift of 1 - 3, at the peak star formation rate density of the Universe. As massive clusters are still forming at the present day at a fraction of the distance to high-redshift galaxies they offer an opportunity to understand the processes controlling star formation and feedback in conditions similar to those in which most stars in the Universe formed. Here we describe a system of massive clusters and their progenitor gas clouds in the centre of the Milky Way, and outline how detailed observations of this system may be able to: (i) help answer some of the fundamental open questions in star formation and (ii) quantify how stellar feedback couples to the surrounding interstellar medium in this high-pressure, high-redshift analogue environment.

A nebular analysis of the central Orion Nebula and its main structures is presented. We exploit MUSE integral field observations in the wavelength range 4595-9366 \r{A} to produce the first O, S and N ionic and total abundance maps of a region spanning 6' x 5' with a spatial resolution of 0.2". We use the S$_{23}$ ( = ([SII]$\lambda$6717,31+[SIII]$\lambda$9068)/H$\beta$) parameter, together with [OII]/[OIII] as an indicator of the degree of ionisation, to distinguish between the various small-scale structures. The only Orion Bullet covered by MUSE is HH 201, which shows a double component in the [FeII]$\lambda$8617 line throughout indicating an expansion, and we discuss a scenario in which this object is undergoing a disruptive event. We separate the proplyds located south of the Bright Bar into four categories depending on their S$_{23}$ values, propose the utility of the S$_{23}$ parameter as an indicator of the shock-contribution to the excitation of line-emitting atoms, and show that the MUSE data is able to identify the proplyds associated with disks and microjets. We compute the second order structure function for the H$\alpha$, [OIII]$\lambda$5007, [SII]$\lambda$6731 and [OI]$\lambda$6300 emission lines to analyse the turbulent velocity field of the region covered with MUSE. We find that the spectral and spatial resolution of MUSE is not able to faithfully reproduce the structure functions of previous works.

We sort $4683$ molecular clouds between $10^\circ< \ell <65^\circ$ from the Bolocam Galactic Plane Survey based on observational diagnostics of star formation activity: compact $70$ $\mu{\rm m}$ sources, mid-IR color-selected YSOs, ${\rm H_2O}$ and ${\rm CH_3OH}$ masers, and UCHII regions. We also present a combined ${\rm NH_3}$-derived gas kinetic temperature and ${\rm H_2O}$ maser catalog for $1788$ clumps from our own GBT 100m observations and from the literature. We identify a subsample of $2223$ ($47.5\%$) starless clump candidates, the largest and most robust sample identified from a blind survey to date. Distributions of flux density, flux concentration, solid angle, kinetic temperature, column density, radius, and mass show strong ($>1$ dex) progressions when sorted by star formation indicator. The median starless clump candidate is marginally sub-virial ($\alpha \sim 0.7$) with $>75\%$ of clumps with known distance being gravitationally bound ($\alpha < 2$). These samples show a statistically significant increase in the median clump mass of $\Delta M \sim 170-370$ M$_\odot$ from the starless candidates to clumps associated with protostars. This trend could be due to (i) mass growth of the clumps at $\dot{M}\sim200-440$ Msun Myr$^{-1}$ for an average free-fall $0.8$ Myr time-scale, (ii) a systematic factor of two increase in dust opacity from starless to protostellar phases, (iii) and/or a variation in the ratio of starless to protostellar clump lifetime that scales as $\sim M^{-0.4}$. By comparing to the observed number of ${\rm CH_3OH}$ maser containing clumps we estimate the phase-lifetime of massive ($M>10^3$ M$_\odot$) starless clumps to be $0.37 \pm 0.08 \ {\rm Myr} \ (M/10^3 \ {\rm M}_\odot)^{-1}$; the majority ($M<450$ M$_\odot$) have phase-lifetimes longer than their average free-fall time.

We have discovered a new H$_2$CO (formaldehyde) $1_{1,0}-1_{1,1}$ 4.82966 GHz maser in Galactic Center Cloud C, G0.38+0.04. At the time of acceptance, this is the eighth region containing an H$_2$CO maser detected in the Galaxy. Cloud C is one of only two sites of confirmed high-mass star formation along the Galactic Center Ridge, affirming that H$_2$CO masers are exclusively associated with high-mass star formation. This discovery led us to search for other masers, among which we found new SiO vibrationally excited masers, making this the fourth star-forming region in the Galaxy to exhibit SiO maser emission. Cloud C is also a known source of CH$_3$OH Class-II and OH maser emission. There are now two known SiO and H$_2$CO maser containing regions in the CMZ, compared to two and six respectively in the Galactic disk, while there is a relative dearth of H$_2$O and CH$_3$OH Class-II masers in the CMZ. SiO and H$_2$CO masers may be preferentially excited in the CMZ, perhaps due to higher gas-phase abundances from grain destruction and heating, or alternatively H$_2$O and CH$_3$OH maser formation may be suppressed in the CMZ. In any case, Cloud C is a new testing ground for understanding maser excitation conditions.

This paper discusses compelling science cases for a future long-baseline interferometer operating at millimeter and centimeter wavelengths, like the proposed Next Generation Vary Large Array (ngVLA). We report on the activities of the Cradle of Life science working group, which focused on the formation of low- and high-mass stars, the formation of planets and evolution of protoplanetary disks, the physical and compositional study of Solar System bodies, and the possible detection of radio signals from extraterrestrial civilizations. We propose 19 scientific projects based on the current specification of the ngVLA. Five of them are highlighted as possible Key Science Projects: (1) Resolving the density structure and dynamics of the youngest HII regions and high-mass protostellar jets, (2) Unveiling binary/multiple protostars at higher resolution, (3) Mapping planet formation regions in nearby disks on scales down to 1 AU, (4) Studying the formation of complex molecules, and (5) Deep atmospheric mapping of giant planets in the Solar System. For each of these projects, we discuss the scientific importance and feasibility. The results presented here should be considered as the beginning of a more in-depth analysis of the science enabled by such a facility, and are by no means complete or exhaustive.

The statistical description of Giant Molecular Cloud (GMC) properties relies heavily on the performance of automatic identification algorithms, which are often seriously affected by the survey design. The algorithm we designed, SCIMES (Spectral Clustering for Interstellar Molecular Emission Segmentation), is able to overcome some of these limitations by considering the cloud segmentation problem in the broad framework of the graph theory. The application of the code on the CO(3-2) High Resolution Survey (COHRS) data allowed for a robust decomposition of more than 12,000 objects in the Galactic Plane. Together with the wealth of Galactic Plane surveys of the recent years, this approach will help to open the door to a future, systematic cataloging of all discrete molecular features of our own Galaxy.

We present a generalization of the Giant Molecular Cloud (GMC) identification problem based on cluster analysis. The method we designed, SCIMES (Spectral Clustering for Interstellar Molecular Emission Segmentation) considers the dendrogram of emission in the broader framework of graph theory and utilizes spectral clustering to find discrete regions with similar emission properties. For Galactic molecular cloud structures, we show that the characteristic volume and/or integrated CO luminosity are useful criteria to define the clustering, yielding emission structures that closely reproduce "by-eye" identification results. SCIMES performs best on well-resolved, high-resolution data, making it complementary to other available algorithms. Using 12CO(1-0) data for the Orion-Monoceros complex, we demonstrate that SCIMES provides robust results against changes of the dendrogram-construction parameters, noise realizations and degraded resolution. By comparing SCIMES with other cloud decomposition approaches, we show that our method is able to identify all canonical clouds of the Orion-Monoceros region, avoiding the over-division within high resolution survey data that represents a common limitation of several decomposition algorithms. The Orion-Monoceros objects exhibit hierarchies and size-line width relationships typical to the turbulent gas in molecular clouds, although "the Scissors" region deviates from this common description. SCIMES represents a significant step forward in moving away from pixel-based cloud segmentation towards a more physical-oriented approach, where virtually all properties of the ISM can be used for the segmentation of discrete objects.

The Galactic center is the closest region in which we can study star formation under extreme physical conditions like those in high-redshift galaxies. We measure the temperature of the dense gas in the central molecular zone (CMZ) and examine what drives it. We mapped the inner 300 pc of the CMZ in the temperature-sensitive J = 3-2 para-formaldehyde (p-H$_2$CO) transitions. We used the $3_{2,1} - 2_{2,0} / 3_{0,3} - 2_{0,2}$ line ratio to determine the gas temperature in $n \sim 10^4 - 10^5 $cm$^{-3}$ gas. We have produced temperature maps and cubes with 30" and 1 km/s resolution and published all data in FITS form. Dense gas temperatures in the Galactic center range from ~60 K to > 100 K in selected regions. The highest gas temperatures T_G > 100 K are observed around the Sgr B2 cores, in the extended Sgr B2 cloud, the 20 km/s and 50 km/s clouds, and in "The Brick" (G0.253+0.016). We infer an upper limit on the cosmic ray ionization rate ${\zeta}_{CR} < 10^{-14}$ 1/s. The dense molecular gas temperature of the region around our Galactic center is similar to values found in the central regions of other galaxies, in particular starburst systems. The gas temperature is uniformly higher than the dust temperature, confirming that dust is a coolant in the dense gas. Turbulent heating can readily explain the observed temperatures given the observed line widths. Cosmic rays cannot explain the observed variation in gas temperatures, so CMZ dense gas temperatures are not dominated by cosmic ray heating. The gas temperatures previously observed to be high in the inner ~75 pc are confirmed to be high in the entire CMZ.

The ubiquity of filamentary structure at various scales through out the Galaxy has triggered a renewed interest in their formation, evolution, and role in star formation. The largest filaments can reach up to Galactic scale as part of the spiral arm structure. However, such large scale filaments are hard to identify systematically due to limitations in identifying methodology (i.e., as extinction features). We present a new approach to directly search for the largest, coldest, and densest filaments in the Galaxy, making use of sensitive Herschel Hi-GAL data complemented by spectral line cubes. We present a sample of the 9 most prominent Herschel filaments, including 6 identified from a pilot search field plus 3 from outside the field. These filaments measure 37-99 pc long and 0.6-3.0 pc wide with masses (0.5-8.3)$\times10^4 \, M_\odot$, and beam-averaged ($28"$, or 0.4-0.7 pc) peak H$_2$ column densities of (1.7-9.3)$\times 10^{22} \, \rm{cm^{-2}}$. The bulk of the filaments are relatively cold (17-21 K), while some local clumps have a dust temperature up to 25-47 K. All the filaments are located within <~60 pc from the Galactic mid-plane. Comparing the filaments to a recent spiral arm model incorporating the latest parallax measurements, we find that 7/9 of them reside within arms, but most are close to arm edges. These filaments are comparable in length to the Galactic scale height and therefore are not simply part of a grander turbulent cascade.

We present a new integral-field spectroscopic dataset of the central part of the Orion Nebula (M 42), observed with the MUSE instrument at the ESO VLT. We reduced the data with the public MUSE pipeline. The output products are two FITS cubes with a spatial size of ~5.9'x4.9' (corresponding to ~0.76 pc x 0.63 pc) and a contiguous wavelength coverage of 4595...9366 Angstrom, spatially sampled at 0.2". We provide two versions with a sampling of 1.25 Angstrom and 0.85 Angstrom in dispersion direction. Together with variance cubes these files have a size of 75 and 110 GiB on disk. They represent one of the largest integral field mosaics to date in terms of information content. We make them available for use in the community. To validate this dataset, we compare world coordinates, reconstructed magnitudes, velocities, and absolute and relative emission line fluxes to the literature and find excellent agreement. We derive a two-dimensional map of extinction and present de-reddened flux maps of several individual emission lines and of diagnostic line ratios. We estimate physical properties of the Orion Nebula, using the emission line ratios [N II] and [S III] (for the electron temperature $T_e$) and [S II] and [Cl III] (for the electron density $N_e$), and show two-dimensional images of the velocity measured from several bright emission lines.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

Integral field unit (IFU) data of the iconic Pillars of Creation in M16 are presented. The ionisation structure of the pillars was studied in great detail over almost the entire visible wavelength range, and maps of the relevant physical parameters, e.g. extinction, electron density, electron temperature, line-of-sight velocity of the ionised and neutral gas are shown. In agreement with previous authors, we find that the pillar tips are being ionised and photo-evaporated by the massive members of the nearby cluster NGC 6611. They display a stratified ionisation structure where the emission lines peak in a descending order according to their ionisation energies. The IFU data allowed us to analyse the kinematics of the photo-evaporative flow in terms of the stratified ionisation structure, and we find that, in agreement with simulations, the photo-evaporative flow is traced by a blueshift in the position-velocity profile. The gas kinematics and ionisation structure have allowed us to produce a sketch of the 3D geometry of the Pillars, positioning the pillars with respect to the ionising cluster stars. We use a novel method to detect a previously unknown bipolar outflow at the tip of the middle pillar and suggest that it has an embedded protostar as its driving source. Furthermore we identify a candidate outflow in the leftmost pillar. With the derived physical parameters and ionic abundances, we estimate a mass loss rate due to the photo-evaporative flow of 70 M$_{\odot}$ Myr$^{-1}$ which yields an expected lifetime of approximately 3 Myr.

We use the distance probability density function (DPDF) formalism of Ellsworth-Bowers et al. (2013, 2015) to derive physical properties for the collection of 1,710 Bolocam Galactic Plane Survey (BGPS) version 2 sources with well-constrained distance estimates. To account for Malmquist bias, we estimate that the present sample of BGPS sources is 90% complete above 400 $M_\odot$ and 50% complete above 70 $M_\odot$. The mass distributions for the entire sample and astrophysically motivated subsets are generally fitted well by a lognormal function, with approximately power-law distributions at high mass. Power-law behavior emerges more clearly when the sample population is narrowed in heliocentric distance (power-law index $\alpha = 2.0\pm0.1$ for sources nearer than 6.5 kpc and $\alpha = 1.9\pm0.1$ for objects between 2 kpc and 10 kpc). The high-mass power-law indices are generally $1.85 \leq \alpha \leq 2.05$ for various subsamples of sources, intermediate between that of giant molecular clouds and the stellar initial mass function. The fit to the entire sample yields a high-mass power-law $\hat{\alpha} = 1.94_{-0.10}^{+0.34}$. Physical properties of BGPS sources are consistent with large molecular cloud clumps or small molecular clouds, but the fractal nature of the dense interstellar medium makes difficult the mapping of observational categories to the dominant physical processes driving the observed structure. The face-on map of the Galactic disk's mass surface density based on BGPS dense molecular cloud structures reveals the high-mass star-forming regions W43, W49, and W51 as prominent mass concentrations in the first quadrant. Furthermore, we present a 0.25-kpc resolution map of the dense gas mass fraction across the Galactic disk that peaks around 5%.

We present 107 maps of continuum emission at 350 microns from Galactic molecular clumps. Observed sources were mainly selected from the Bolocam Galactic Plane Survey (BGPS) catalog, with 3 additional maps covering star forming regions in the outer Galaxy. The higher resolution of the SHARC-II images (8.5" beam) compared with the 1.1 mm images from BGPS (33" beam) allowed us to identify a large population of smaller substructures within the clumps. A catalog is presented for the 1386 sources extracted from the 350 micron maps. The color temperature distribution of clumps based on the two wavelengths has a median of 13.3 K and mean of 16.3 +- 0.4 K, assuming an opacity law index of 1.7. For the structures with good determination of color temperatures, the mean ratio of gas temperature, determined from NH3 observations, to dust color temperature is 0.88 and the median ratio is 0.76. About half the clumps have more than two substructures and 22 clumps have more than 10. The fraction of the mass in dense substructures seen at 350 microns compared to the mass of their parental clump is ~0.19, and the surface densities of these substructures are, on average, 2.2 times those seen in the clumps identified at 1.1 mm. For a well-characterized sample, 88 structures (31%) exceed a surface density of 0.2 g cm^(-2), and 18 (6%) exceed 1.0 g cm^(-2), thresholds for massive star formation suggested by theorists.

Aims. Adaptive optics images are used to test the hypothesis that the explosive BN/KL outflow from the Orion OMC1 cloud core was powered by the dynamical decay of a non-hierarchical system of massive stars. Methods. Narrow-band H2, [Fe II], and broad-band Ks obtained with the Gemini South multi-conjugate adaptive optics (AO) system GeMS and near-infrared imager GSAOI are presented. The images reach resolutions of 0.08 to 0.10", close to the 0.07" diffraction limit of the 8-meter telescope at 2.12 microns. Comparison with previous AO-assisted observations of sub-fields and other ground-based observations enable measurements of proper motions and the investigation of morphological changes in H2 and [Fe II] features with unprecedented precision. The images are compared with numerical simulations of compact, high-density clumps moving ~1000 times their own diameter through a lower density medium at Mach 1000. Results. Several sub-arcsecond H2 features and many [Fe II] 'fingertips' on the projected outskirts of the flow show proper motions of ~300 km/s. High-velocity, sub-arcsecond H2 knots ('bullets') are seen as far as 140" from their suspected ejection site. If these knots propagated through the dense Orion A cloud, their survival sets a lower bound on their densities of order 10^7 cm^-3, consistent with an origin within a few au of a massive star and accelerated by a final multi-body dynamic encounter that ejected the BN object and radio source I from OMC1 about 500 years ago. Conclusions. Over 120 high-velocity bow-shocks propagating in nearly all directions from the OMC1 cloud core provide evidence for an explosive origin for the BN/KL outflow triggered by the dynamic decay of a non-hierarchical system of massive stars. Such events may be linked to the origin of runaway, massive stars.

We present new 2 cm and 6 cm maps of H2CO, radio recombination lines, and the radio continuum in the W51 star forming complex acquired with Arecibo and the Green Bank Telescope at ~50" resolution. We use H2CO absorption to determine the relative line-of-sight positions of molecular and ionized gas. We measure gas densities using the H2CO densitometer, including continuous measurements of the dense gas mass fraction (DGMF) over the range $10^4$ cm$^{-3}$ < n(H$_2$) < $10^6$ cm$^{-3}$- this is the first time a dense gas mass fraction has been measured over a range of densities with a single data set. The DGMF in W51A is high,f >~70% above $n>10^4$ cm$^{-3}$, while it is low, f<20%, in W51 B. We did not detect any H2CO emission throughout the W51 GMC; all gas dense enough to emit under normal conditions is in front of bright continuum sources and therefore is seen in absorption instead. The data set has been made public at http://dx.doi.org/10.7910/DVN/26818. Conclusions. (1) The dense gas fraction in the W51 A and B clouds shows that W51 A will continue to form stars vigorously, while star formation has mostly ended in W51 B. The lack of dense, star-forming gas around W51 C indicates that collect-and-collapse is not acting or is inefficient in W51. (2) Ongoing high-mass star formation is correlated with n ~ 1x10$^5$ cm$^{-3}$ gas. Gas with n > 10$^4$ cm$^{-3}$ is weakly correlated with low and moderate mass star formation, but does not strongly correlate with high-mass star formation. (3) The nondetection of H$_2$CO emission implies that the emission detected in other galaxies, e.g. Arp 220, comes from high-density gas that is not directly affiliated with already-formed massive stars. Either the non-star-forming ISM of these galaxies is very dense, implying the star formation density threshold is higher, or H II regions have their emission suppressed.

The bandwith, sensitivity and sheer survey speed of the SKA offers unique potential for deep spectroscopic surveys of the Milky Way. Within the frequency bands available to the SKA lie many transitions that trace the ionised, radical and molecular components of the interstellar medium and which will revolutionise our understanding of many physical processes. In this chapter we describe the impact on our understanding of the Milky Way that can be achieved by spectroscopic SKA surveys, including "out of the box" early science with radio recombination lines, Phase 1 surveys of the molecular ISM using anomalous formaldehyde absorption, and full SKA surveys of ammonia inversion lines.

We present an expanded distance catalog for 1,710 molecular cloud structures identified in the Bolocam Galactic Plane Survey (BGPS) version 2, representing a nearly threefold increase over the previous BGPS distance catalog. We additionally present a new method for incorporating extant data sets into our Bayesian distance probability density function (DPDF) methodology. To augment the dense-gas tracers (e.g., HCO+(3-2), NH3(1,1)) used to derive line-of-sight velocities for kinematic distances, we utilize the Galactic Ring Survey 13CO(1-0) data to morphologically extract velocities for BGPS sources. The outline of a BGPS source is used to select a region of the GRS 13CO data, along with a reference region to subtract enveloping diffuse emission, to produce a line profile of 13CO matched to the BGPS source. For objects with a HCO+(3-2) velocity, \approx 95% of the new 13CO(1-0) velocities agree with that of the dense gas. A new prior DPDF for kinematic distance ambiguity (KDA) resolution, based on a validated formalism for associating molecular cloud structures with known objects from the literature, is presented. We demonstrate this prior using catalogs of masers with trigonometric parallaxes and HII regions with robust KDA resolutions. The distance catalog presented here contains well-constrained distance estimates for 20% of BGPS V2 sources, with typical distance uncertainties \lesssim 0.5 kpc. Approximately 75% of the well-constrained sources lie within 6 kpc of the Sun, concentrated in the Scutum-Centarus arm. Galactocentric positions of objects additionally trace out portions of the Sagittarius, Perseus, and Outer arms in the first and second Galactic quadrants, and we also find evidence for significant regions of interarm dense gas.

We present observations of near-infrared 2.12 micro-meter molecular hydrogen outflows emerging from 1.1 mm dust continuum clumps in the North America and Pelican Nebula (NAP) complex selected from the Bolocam Galactic Plane Survey (BGPS). Hundreds of individual shocks powered by over 50 outflows from young stars are identified, indicating that the dusty molecular clumps surrounding the NGC 7000 / IC 5070 / W80 HII region are among the most active sites of on-going star formation in the Solar vicinity. A spectacular X-shaped outflow, MHO 3400, emerges from a young star system embedded in a dense clump more than a parsec from the ionization front associated with the Pelican Nebula (IC 5070). Suspected to be a binary, the source drives a pair of outflows with orientations differing by 80 degrees. Each flow exhibits S-shaped symmetry and multiple shocks indicating a pulsed and precessing jet. The `Gulf of Mexico' located south of the North America Nebula (NGC 7000), contains a dense cluster of molecular hydrogen objects (MHOs), Herbig-Haro (HH) objects, and over 300 YSOs, indicating a recent burst of star formation. The largest outflow detected thus far in the North America and Pelican Nebula complex, the 1.6 parsec long MHO 3417 flow, emerges from a 500 Solar mass BGPS clump and may be powered by a forming massive star. Several prominent outflows such as MHO 3427 appear to be powered by highly embedded YSOs only visible at a wavelength > 70 micro-meters. An `activity index' formed by dividing the number of shocks by the mass of the cloud containing their source stars is used to estimate the relative evolutionary states of Bolocam clumps. Outflows can be used as indicators of the evolutionary state of clumps detected in mm and sub-mm dust continuum surveys.

ALMA HCO+ observations of the infrared dark cloud G0.253+0.016 located in the Central Molecular Zone of the Galaxy are presented. The 89 GHz emission is area-filling, optically thick, and sub-thermally excited. Two types of filaments are seen in absorption against the HCO+ emission. Broad-line absorption filaments (BLAs) have widths of less than a few arcseconds (0.07 - 0.14 pc), lengths of 30 to 50 arcseconds (1.2 - 1.8 pc), and absorption profiles extending over a velocity range larger than 20 km/sec. The BLAs are nearly parallel to the nearby G0.18 non-thermal filaments and may trace HCO+ molecules gyrating about highly ordered magnetic fields located in front of G0.253+0.016 or edge-on sheets formed behind supersonic shocks propagating orthogonal to our line-of-sight in the foreground. Narrow-line absorption filaments (NLAs) have line-widths less than 20 km/sec. Some NLAs are also seen in absorption in other species with high optical depth such as HCN and occasionally in emission where the background is faint. The NLAs, which also trace low-density, sub-thermally excited HCO+ molecules, are mostly seen on the blueshifted side of the emission from G0.253+0.016. If associated with the surface of G0.253+0.016, the kinematics of the NLAs indicate that the cloud surface is expanding. The decompression of entrained, milli-Gauss magnetic fields may be responsible for the re-expansion of the surface layers of G0.253+0.016 as it recedes from the Galactic center following a close encounter with Sgr A.

We explore the relationship between gas and dust in massive star-forming regions by comparing physical properties derived from each. We compare the temperatures and column densities in a massive star-forming Infrared Dark Cloud (IRDC, G32.02+0.05), which shows a range of evolutionary states, from quiescent to active. The gas properties were derived using radiative transfer modeling of the (1,1), (2,2), and (4,4) transitions of NH3 on the Karl G. Jansky Very Large Array (VLA), while the dust temperatures and column densities were calculated using cirrus-subtracted, modified blackbody fits to Herschel data. We compare the derived column densities to calculate an NH3 abundance, 4.6 x 10^-8. In the coldest star-forming region, we find that the measured dust temperatures are lower than the measured gas temperatures (mean and standard deviations T_dust ~ 11.6 +/- 0.2 K vs. T_gas ~ 15.2 +/- 1.5 K), which may indicate that the gas and dust are not well-coupled in the youngest regions (~0.5 Myr) or that these observations probe a regime where the dust and/or gas temperature measurements are unreliable. Finally, we calculate millimeter fluxes based on the temperatures and column densities derived from NH3 which suggest that millimeter dust continuum observations of massive star-forming regions, such as the Bolocam Galactic Plane Survey or ATLASGAL, can probe hot cores, cold cores, and the dense gas lanes from which they form, and are generally not dominated by the hottest core.

We present new NH3 (1,1), (2,2), and (4,4) observations from the Karl G. Jansky Very Large Array (VLA) compiled with work in the literature to explore the range of conditions observed in young, massive star-forming regions. To sample the effects of evolution independent from those of distance/resolution, abundance, and large-scale environment, we compare clumps in different evolutionary stages within a single Infrared Dark Cloud (IRDC), G32.02+0.06. We find that the early stages of clustered star formation are characterized by dense, parsec-scale filamentary structures interspersed with complexes of dense cores (<0.1 pc cores clustered in complexes separated by ~1 pc) with masses from about 10 to 100 Msun. The most quiescent core is the most extended while the star-forming cores are denser and more compact, showing very similar column density structure before and shortly after the onset of massive star formation, with peak surface densities >~ 1 g cm^-2. Quiescent cores and filaments show smoothly varying temperatures from 10-20 K, rising to over 40 K in star-forming cores. We calculate the virial parameters for 16 cores and find that the level of support provided by turbulence is generally insufficient to support them against gravitational collapse (alpha_vir ~ 0.6). The star-forming filaments show smooth velocity fields, punctuated by discontinuities at the sites of active star formation. We discuss the Massive Molecular Filament (MMF; M > 10^5 Msun, l > 60 pc) hosting the IRDC, hypothesizing that it may have been shaped by previous generations of massive stars.

At a distance of 8 kpc, the center of our Galaxy is the nearest galactic nucleus, and has been the subject of numerous key projects undertaken by great observatories such as Chandra, Spitzer, and Herschel. However, there are still no surveys of molecular gas properties in the Galactic center with less than 30" (1 pc) resolution. There is also no sensitive polarization survey of this region, despite numerous nonthermal magnetic features apparently unique to the central 300 parsecs. In this paper, we outline the potential the VLASS has to fill this gap. We assess multiple considerations in observing the Galactic center, and recommend a C-band survey with 10 micro-Jy continuum RMS and sensitive to molecular gas with densities greater than 10^4 cm^{-3}, covering 17 square degrees in both DnC and CnB configurations ( resolution ~5"), totaling 750 hours of observing time. Ultimately, we wish to note that the upgraded VLA is not just optimized for fast continuum surveys, but has a powerful correlator capable of simultaneously observing continuum emission and dozens of molecular and recombination lines. This is an enormous strength that should be fully exploited and highlighted by the VLASS, and which is ideally suited for surveying the center of our Galaxy.

We present optical and near-IR spectroscopic observations of the luminous blue variable SN 2009ip during its remarkable photometric evolution of 2012. The spectra sample three key points in the SN 2009ip lightcurve, corresponding to its initial brightening in August (2012-A) and its dramatic rebrightening in early October (2012-B). Based on line fluxes and velocities measured in our spectra, we find a surprisingly low I(H-alpha)/I(H-beta) ~ 1.3-1.4 in the 2012-B spectra. Such a ratio implies either a rare Case B recombination scenario where H-alpha, but not H-beta, is optically thick, or an extremely high density for the circumstellar material of n_e > 10^13 cm^(-3) The H-alpha line intensity yields a minimum radiating surface area of >~20,000 AU^2 in H-alpha at the peak of SN 2009ip's photometric evolution. Combined with the nature of this object's spectral evolution in 2012, a high circumstellar density and large radiating surface area imply the presence of a thin disk geometry around the central star (and, consequently, a possible binary companion), suggesting that the observed 2012-B rebrightening of SN 2009ip can be attributed to the illumination of the disk's inner rim by fast-moving ejecta produced by the underlying events of 2012-A.

Molecular clouds are supersonically turbulent. This turbulence governs the initial mass function and the star formation rate. In order to understand the details of star formation, it is therefore essential to understand the properties of turbulence, in particular the probability distribution of density in turbulent clouds. We present formaldehyde volume density measurements of a non-star-forming cloud along the line of sight towards W49A. We use these measurements in conjunction with total mass estimates from 13CO to infer the shape of the density probability distribution function. This method is complementary to measurements of turbulence via the column density distribution and should be applicable to any molecular cloud with detected CO. We show that turbulence in this cloud is probably compressively driven, with a compressive-to-total Mach number ratio $b = \mathcal{M}_C/\mathcal{M}>0.4$. We measure the standard deviation of the density distribution, constraining it to the range $1.5 < \sigma_s < 1.9$ assuming that the density is lognormally distributed. This measurement represents an essential input into star formation laws. The method of averaging over different excitation conditions to produce a model of emission from a turbulent cloud is generally applicable to optically thin line observations.

We present new near-IR H2, CO J=2-1, and CO J = 3-2 observations to study outflows in the massive star forming region IRAS 05358+3543. The Canada-France-Hawaii Telescope H2 images and James Clerk Maxwell Telescope CO data cubes of the IRAS 05358 region reveal several new outflows, most of which emerge from the dense cluster of sub-mm cores associated with the Sh 2-233IR NE cluster to the northeast of IRAS 05358. We used Apache Point Observatory (APO) JHK spectra to determine line of sight velocities of the outflowing material. Analysis of archival VLA cm continuum data and previously published VLBI observations reveal a massive star binary as a probable source of one or two of the outflows. We have identified probable sources for 6 outflows and candidate counterflows for 7 out of a total of 11 seen to be originating from the IRAS 05358 clusters. We classify the clumps within Sh 2-233IR NE as an early protocluster and Sh 2-233IR SW as a young cluster, and conclude that the outflow energy injection rate approximately matches the turbulent decay rate in Sh 2-233IR NE.

We present near-infrared spectroscopy and 1 mm line and continuum observations of a recently identified site of high mass star formation likely to be located in the Central Molecular Zone near Sgr C. Located on the outskirts of the massive evolved HII region associated with Sgr C, the area is characterized by an Extended Green Object measuring ~10" in size (0.4 pc), whose observational characteristics suggest the presence of an embedded massive protostar driving an outflow. Our data confirm that early-stage star formation is taking place on the periphery of the Sgr C HII region, with detections of two protostellar cores and several knots of H2 and Brackett gamma emission alongside a previously detected compact radio source. We calculate the cores' joint mass to be ~10^3 Msun, with column densities of 1-2 x 10^24 cm-2. We show the host molecular cloud to hold ~10^5 Msun of gas and dust with temperatures and column densities favourable for massive star formation to occur, however, there is no evidence of star formation outside of the EGO, indicating that the cloud is predominantly quiescent. Given its mass, density, and temperature, the cloud is comparable to other remarkable non-star-forming clouds such as G0.253 in the Eastern CMZ.

The Bolocam Galactic Plane Survey (BGPS) is a 1.1 mm continuum survey of dense clumps of dust throughout the Galaxy covering 170 square degrees. We present spectroscopic observations using the Heinrich Hertz Submillimeter Telescope of the dense gas tracers, HCO+ and N2H+ 3-2, for all 6194 sources in the Bolocam Galactic Plane Survey v1.0.1 catalog between 7.5 <= l <= 194 degrees. This is the largest targeted spectroscopic survey of dense molecular gas in the Milky Way to date. We find unique velocities for 3126 (50.5%) of the BGPS v1.0.1 sources observed. Strong N2H+ 3-2 emission (T_{mb} > 0.5 K) without HCO+ 3-2 emission does not occur in this catalog. We characterize the properties of the dense molecular gas emission toward the entire sample. HCO+ is very sub-thermally populated and the 3-2 transitions are optically thick toward most BGPS clumps. The median observed line width is 3.3 km/s consistent with supersonic turbulence within BGPS clumps. We find strong correlations between dense molecular gas integrated intensities and 1.1 mm peak flux and the gas kinetic temperature derived from previously published NH3 observations. These intensity correlations are driven by the sensitivity of the 3-2 transitions to excitation conditions rather than by variations in molecular column density or abundance. We identify a subset of 113 sources with stronger N2H+ than HCO+ integrated intensity, but we find no correlations between the N2H+ / HCO+ ratio and 1.1 mm continuum flux density, gas kinetic temperature, or line width. Self-absorbed profiles are rare (1.3%).

We present a re-reduction and expansion of the Bolocam Galactic Plane Survey, first presented by Aguirre et al. (2011) and Rosolowsky et al. (2010). The BGPS is a 1.1 mm survey of dust emission in the Northern galactic plane, covering longitudes -10 < \ell < 90 and latitudes |b| < 0.5 with a typical 1-\sigma RMS sensitivity of 30-100 mJy in a 33" beam. Version 2 of the survey includes an additional 20 square degrees of coverage in the 3rd and 4th quadrants and 2 square degrees in the 1st quadrant. The new data release has improved angular recovery, with complete recovery out to 80" and partial recovery to 300", and reduced negative bowls around bright sources resulting from the atmospheric subtraction process. We resolve the factor of 1.5 flux calibration offset between the v1.0 data release and other data sets and determine that there is no offset between v2.0 and other data sets. The v2.0 pointing accuracy is tested against other surveys and demonstrated to be accurate and an improvement over v1.0. We present simulations and tests of the pipeline and its properties, including measurements of the pipeline's angular transfer function. The Bolocat cataloging tool was used to extract a new catalog, which includes 8594 sources, with 591 in the expanded regions. We have demonstrated that the Bolocat 40" and 80" apertures are accurate even in the presence of strong extended background emission. The number of sources is lower than in v1.0, but the amount of flux and area included in identified sources is larger.

We present the first public version (v0.2) of the open-source and community-developed Python package, Astropy. This package provides core astronomy-related functionality to the community, including support for domain-specific file formats such as Flexible Image Transport System (FITS) files, Virtual Observatory (VO) tables, and common ASCII table formats, unit and physical quantity conversions, physical constants specific to astronomy, celestial coordinate and time transformations, world coordinate system (WCS) support, generalized containers for representing gridded as well as tabular data, and a framework for cosmological transformations and conversions. Significant functionality is under active development, such as a model fitting framework, VO client and server tools, and aperture and point spread function (PSF) photometry tools. The core development team is actively making additions and enhancements to the current code base, and we encourage anyone interested to participate in the development of future Astropy versions.

We present the first overview of the Herschel observations of the nearby high-mass star-forming region NGC 7538, taken as part of the Herschel imaging study of OB Young Stellar objects (HOBYS) Key Programme. These PACS and SPIRE maps cover an approximate area of one square degree at five submillimeter and far-infrared wavebands. We have identified 780 dense sources and classified 224 of those. With the intention of investigating the existence of cold massive starless or class 0-like clumps that would have the potential to form intermediate- to high-mass stars, we further isolate 13 clumps as the most likely candidates for followup studies. These 13 clumps have masses in excess of 40 M_sun and temperatures below 15 K. They range in size from 0.4 pc to 2.5 pc and have densities between 3x10^3 cm^-3 to 4x10^4 cm^-3. Spectral energy distributions are then used to characterize their energetics and evolutionary state through a luminosity-mass diagram. NGC 7538 has a highly filamentary structure, previously unseen in the dust continuum of existing submillimeter surveys. We report the most complete imaging to date of a large, evacuated ring of material in NGC 7538 which is bordered by many cool sources.

The 2012 explosion of SN2009ip raises questions about our understanding of the late stages of massive star evolution. Here we present a comprehensive study of SN2009ip during its remarkable re-brightening(s). High-cadence photometric and spectroscopic observations from the GeV to the radio band obtained from a variety of ground-based and space facilities (including the VLA, Swift, Fermi, HST and XMM) constrain SN2009ip to be a low energy (E~ 10^50 erg for an ejecta mass ~ 0.5 Msun) and likely asymmetric explosion in a complex medium shaped by multiple eruptions of the restless progenitor star. Most of the energy is radiated as a result of the shock breaking out through a dense shell of material located at 5x10^14 cm with M~0.1 Msun, ejected by the precursor outburst ~40 days before the major explosion. We interpret the NIR excess of emission as signature of dust vaporization of material located further out (R>4x 10^15 cm), the origin of which has to be connected with documented mass loss episodes in the previous years. Our modeling predicts bright neutrino emission associated with the shock break-out if the cosmic ray energy is comparable to the radiated energy. We connect this phenomenology with the explosive ejection of the outer layers of the massive progenitor star, that later interacted with material deposited in the surroundings by previous eruptions. Future observations will reveal if the luminous blue variable (LBV) progenitor star survived. Irrespective of whether the explosion was terminal, SN2009ip brought to light the existence of new channels for sustained episodic mass-loss, the physical origin of which has yet to be identified.

We present a new distance estimation method for dust-continuum-identified molecular cloud clumps. Recent (sub-)millimeter Galactic plane surveys have cataloged tens of thousands of these objects, but detailed study of their physical properties requires robust distance determinations. We derive Bayesian distance probability density functions (DPDFs) for 770 objects from the Bolocam Galactic Plane Survey in the longitude range 7.5 < l < 65. The DPDF formalism is based on kinematic distances, and uses external data sets to place prior distance probabilities to resolve the kinematic distance ambiguity (KDA) for objects in the inner Galaxy. We present priors related to the mid-infrared absorption of dust in dense molecular regions and the Galactic distribution of molecular gas. By assuming a numerical model of Galactic mid-infrared emission and simple radiative transfer, we match the morphology of millimeter thermal dust emission with mid-infrared absorption to compute a prior DPDF for distance discrimination. Selecting objects first from (sub-)millimeter source catalogs avoids a bias towards the darkest infrared dark clouds (IRDCs) and extends the range of heliocentric distance and contrast range probed by mid-infrared extinction. We derive well-constrained KDA resolutions for 618 molecular cloud clumps, with approximately 15% placed at or beyond the tangent distance. Objects sufficiently dark to be cataloged as IRDCs are generally placed at the near kinematic distance. Distance comparisons with Galactic Ring Survey KDA resolutions yield a 92% agreement, and a Galactic face-on view reveals sections of the Sagittarius and Scutum-Centaurus Arms. This KDA-resolution method for large catalogs of sources through the combination of (sub-)millimeter and mid-infrared observations of molecular cloud clumps is generally applicable to other dust-continuum Galactic plane surveys.

We have mapped the Auriga/California molecular cloud with the Herschel PACS and SPIRE cameras and the Bolocam 1.1 mm camera on the Caltech Submillimeter Observatory (CSO) with the eventual goal of quantifying the star formation and cloud structure in this Giant Molecular Cloud (GMC) that is comparable in size and mass to the Orion GMC, but which appears to be forming far fewer stars. We have tabulated 60 compact 70/160um sources that are likely pre-main-sequence objects and correlated those with Spitzer and WISE mid-IR sources. At 1.1 mm we find 18 cold, compact sources and discuss their properties. The most important result from this part of our study is that we find a modest number of additional compact young objects beyond those identified at shorter wavelengths with Spitzer. We also describe the dust column density and temperature structure derived from our photometric maps. The column density peaks at a few x 10^22 cm^-2 (N_H2) and is distributed in a clear filamentary structure along which nearly all the pre-main-sequence objects are found. We compare the YSO surface density to the gas column density and find a strong non-linear correlation between them. The dust temperature in the densest parts of the filaments drops to ~10K from values ~ 14--15K in the low density parts of the cloud. We also derive the cumulative mass fraction and probability density function of material in the cloud which we compare with similar data on other star-forming clouds.

SBW1 is a B supergiant surrounded by a ring nebula that is a nearby twin of SN 1987A. We present images and spectra of SBW1 obtained with HST, Spitzer, and Gemini South. HST images of SBW1 do not exhibit long Rayleigh-Taylor fingers, which are presumed to cause the hotspots in the SN1987A ring, but instead show a geometrically thin clumpy ring. The radial mass distribution and size scales of inhomogeneities in SBW1's ring closely resemble those in the SN1987A ring, but the more complete disk expected to reside at the base of the RT fingers is absent. This structure may explain why portions of the SN1987A ring between the hotspots have not yet brightened after more than 15 years. The model we suggest does not require a fast wind colliding with a previous red supergiant wind. More surprisingly, images of SBW1 also reveal diffuse emission filling the interior of the ring in H-alpha and thermal-IR emission; 190K dust dominates the 8-20 micron luminosity (but contains only 1e-5Msun). Cooler (85K dust resides in the equatorial ring itself (dust mass of 5e-3Msun). Diffuse emission extends inward to 1 arcsecond from the central star, where a paucity of emission suggests an inner hole excavated by the wind. We propose that diffuse emission inside the ring arises from an ionized flow of material photoevaporated from the dense ring, and it prevents the supergiant wind from advancing in the equator. This inner emission could correspond to a structure hypothesized to reside around SN1987A that was never directly detected. We suggest that photoionization can play an important dynamical role in shaping the ring nebula, and we speculate that this might help explain the origin of the polar rings around SN1987A. abridged.

We propose that bound, young massive stellar clusters form from dense clouds that have escape speeds greater than the sound speed in photo-ionized gas. In these clumps, radiative feedback in the form of gas ionization is bottled up, enabling star formation to proceed to sufficiently high efficiency so that the resulting star cluster remains bound even after gas removal. We estimate the observable properties of the massive proto-clusters (MPCs) for existing Galactic plane surveys and suggest how they may be sought in recent and upcoming extragalactic observations. These surveys will potentially provide a significant sample of MPC candidates that will allow us to better understand extreme star-formation and massive cluster formation in the Local Universe.

We search the lambda = 1.1 mm Bolocam Galactic Plane Survey for clumps containing sufficient mass to form ~10^4 M\odot star clusters. 18 candidate massive proto-clusters are identified in the first Galactic quadrant outside of the central kiloparsec. This sample is complete to clumps with mass M(clump) > 10^4 M_sun and radius r < 2.5 pc. The overall Galactic massive cluster formation rate is CFR(M_cluster > 10^4) ~ 5 Myr^-1, which is in agreement with the rates inferred from Galactic open clusters and M31 massive clusters. We find that all massive proto-clusters in the first quadrant are actively forming massive stars and place an upper limit of t_starless < 0.5 Myr on the lifetime of the starless phase of massive cluster formation. If massive clusters go through a starless phase with all of their mass in a single clump, the lifetime of this phase is very short.

We present a pilot survey of 21 lines of sight towards UCHII regions in the H2CO 1-1 (6cm) and 2-2 (2cm) transitions, using the H2CO centimeter lines as a molecular gas densitometer. Using Arecibo and Green Bank beam-matched observations, we measure the density of 51 detected H2CO line pairs and present upper limits on density for an additional 24 detected 1-1 lines. We analyze the systematic uncertainties in the H2CO densitometer, achieving H2 density measurements with accuracies ~ 0.1-0.3 dex. The densities measured are not correlated with distance, implying that it is possible to make accurate density measurements throughout the galaxy without a distance bias. We confirm that ultracompact HII regions are associated with gas at densities n(H2)~10^5cm^-3. The densities measured in line-of-sight molecular clouds suggest that they consist of low volume filling factor (f~10^-2) gas at high (n(H2)>10^4cm^-3) density, which is inconsistent with purely supersonic turbulence and requires high-density clumping greater than typically observed in gravoturbulent simulations. We observe complex line morphologies that indicate density variations with velocity around UCHII regions, and we classify a subset of the UCHII molecular envelopes as collapsing or expanding. We compare these measurements to Bolocam Galactic Plane Survey 1.1mm observations, and note that most UCHII regions have 1.1mm emission consisting of significant (5-70%) free-free emission and are therefore not necessarily dominated by optically thin dust emission. A comparison of our data with the Mangum et al. starburst sample shows that the area filling factor of dense (n(H2)~10^5cm^-3) molecular gas in typical starburst galaxies is <~0.01, but in extreme starburst galaxies like Arp 220, is ~0.1, suggesting that Arp 220 is physically similar to an oversized UCHII region.

New JCMT HARP CO 3-2 observations of the W5 star forming complex are presented, totaling an area of 12000 arcmin^2 with sensitivity better than 0.1 K per 0.4 km/s channel. We discovered 55 CO outflow candidates, of which 40 are associated with W5 and 15 are more distant than the Perseus arm. Most of the outflows are located on the periphery of the W5 HII region. However, two outflow clusters are > 5 pc from the ionization fronts, indicating that their driving protostars formed without directly being triggered by the O-stars in W5. We compare the derived outflow properties to those in Perseus and find that the total W5 outflow mass is surprisingly low given the cloud masses. The outflow mass deficiency in the more massive W5 cloud (M(H2) ~ 5 \times 10^4 Msun) can be explained if ionizing radiation dissociates molecules as they break out of their host cloud cores. Although CO J=3-2 is a good outflow tracer, it is likely to be a poor mass tracer because of sub-thermal line excitation and high opacity, which may also contribute to the outflow mass discrepancy. It is unlikely that outflows could provide the observed turbulent energy in the W5 molecular clouds even accounting for undetected outflow material. Many cometary globules have been observed with velocity gradients from head to tail, displaying strong interaction with the W5 HII region and exhibiting signs of triggered or revealed star formation in their heads. Because it is observed face-on, W5 is an excellent region to study feedback effects, both positive and negative, of massive stars on star formation.

We observed the Sh 2-233IR (S233IR) region with better sensitivity in near-infrared than previous studies for this region. By applying statistical subtraction of the back- ground stars, we identified member sources and derived the age and mass of three distinguishable sub-groups in this region: Sh 2-233IR NE, Sh 2-233IR SW, and the "distributed stars" over the whole cloud. Star formation may be occurring sequentially with a relatively small age difference (\sim 0.2-0.3 Myrs) between subclusters. We found that the slopes for initial mass function (Gamma \sim -0.5) of two subclusters are flatter than that of Salpeter, which suggests that more massive stars were preferentially formed in those clusters compared to other Galactic star-forming regions. These subclusters may not result from the overall collapse of the whole cloud, but have formed by triggering before the previous star formation activities disturbed the natal molecular cloud. Addi- tionally, high star formation efficiency (&40%) of the subclusters may also suggest that stars form very efficiently in the center of NE.

We present the results of observations of 1882 sources in the Bolocam Galactic Plane Survey (BGPS) at 1.1 mm with the 10m Heinrich Hertz Telescope simultaneously in HCO+ J=3-2 and N2H+ J=3-2. We detect 77% of these sources in HCO^+ and 51% in N2H+ at greater than 3$\sigma$. We find a strong correlation between the integrated intensity of both dense gas tracers and the 1.1 mm dust emission of BGPS sources. We determine kinematic distances for 529 sources (440 in the first quadrant breaking the distance ambiguity and 89 in the second quadrant) We derive the size, mass, and average density for this subset of clumps. The median size of BGPS clumps is 0.75 pc with a median mass of 330 M$_{\odot}$ (assuming T_{Dust}=20 K). The median HCO+ linewidth is 2.9 km s$^{-1}$ indicating that BGPS clumps are dominated by supersonic turbulence or unresolved kinematic motions. We find no evidence for a size-linewidth relationship for BGPS clumps. We analyze the effects of the assumed dust temperature on the derived clump properties with a Monte Carlo simulation and we find that changing the temperature distribution will change the median source properties (mass, volume-averaged number density, surface density) by factors of a few. The observed differential mass distribution has a power-law slope that is intermediate between that observed for diffuse CO clouds and the stellar IMF. BGPS clumps represent a wide range of objects (from dense cores to more diffuse clumps) and are typically characterized by larger sizes and lower densities than previously published surveys of high-mass star-forming regions. This collection of objects is a less-biased sample of star-forming regions in the Milky Way that likely span a wide range of evolutionary states.

The detected variety in chemistry and circumstellar shell morphology of the limited sample of Galactic post-AGB stars is so large that there is no consensus yet on how the different objects are linked by evolutionary channels. The evaluation is complicated by the fact that their distances and hence luminosities remain largely unknown. Via cross-correlation of the Spitzer SAGE catalogue with optical catalogues we selected a sample of LMC post-AGB candidates based on their [8]-[24] colour index and estimated luminosity. We determined the fundamental properties of the central stars of 105 of these objects using low-resolution, optical spectra that we obtained at Siding Spring Observatory and SAAO, and constructed a catalogue of 70 high probability and 1337 candidate post-AGB stars that is available at the CDS. The sample forms an ideal testbed for stellar evolution theory predictions of the final phase of low- and intermediate-mass stars, because the distance and hence luminosity and also the current and initial mass of these objects is well constrained. About half of the objects in our sample of post-AGB candidates show a spectral energy distribution (SED) that is indicative of a disc rather than an expanding and cooling AGB remnant. Like in the Galaxy, the disc sources are likely associated with binary evolution. Important side products of this research are catalogues of candidate young stellar objects, candidate supergiants with circumstellar dust, and discarded objects for which a spectrum was obtained. These too are available at the CDS.

Context: Despite their profound effect on the universe, the formation of massive stars and stellar clusters remains elusive. In the past decade, compelling evidence has emerged that suggests Infrared Dark Clouds (IRDCs) may be precursors to stellar clusters. However, the usual method for identifying IRDCs is biased by the requirement that they are seen in absorption against bright mid-IR emission, whereas dust continuum observations allow cold, dense pre-stellar-clusters to be identified anywhere. Aims: We aim to understand what physical properties characterize IRDCs, to explore the population of dust continuum sources that are not IRDCs, and to roughly characterize the star formation activity in dust continuum sources. Results: We present temperature and column density maps in the Hi-GAL l=30 and l=59 SDP fields, as well as a robust algorithm for cirrus subtraction and source identification using Hi-GAL. We report on the fraction of Hi-GAL sources which are mid-IR-dark, -neutral, or -bright in both fields. We find significant trends in column density and temperature between mid-IR-dark and -bright pixels; mid-IR-dark pixels are about 10 K colder and have a factor of 2 higher column density on average than mid-IR-bright pixels. We find that Hi-GAL dust continuum sources span a range of evolutionary states from pre- to star-forming, and that warmer sources are associated with more star formation tracers. There is a trend of increasing temperature with tracer type from mid-IR-dark at the coldest, to outflow/maser sources in the middle, and finally to 8 and 24 micron bright sources at the warmest. Finally, we identify five candidate IRDC-like sources on the far-side of the Galaxy. These are cold (~ 20 K), high column density (N(H2) > 10^22 cm^-2) clouds identified with Hi-GAL which, despite bright surrounding mid-IR emission, show little to no absorption at 8 micron. (abridged)

The Bolocam Galactic Plane Survey (BGPS) data for a six square degree region of the Galactic plane containing the Galactic center is analyzed and compared to infrared and radio continuum data. The BGPS 1.1 mm emission consists of clumps interconnected by a network of fainter filaments surrounding cavities, a few of which are filled with diffuse near-IR emission indicating the presence of warm dust or with radio continuum characteristic of HII regions or supernova remnants. New 350 {\mu}m images of the environments of the two brightest regions, Sgr A and B, are presented. Sgr B2 is the brightest mm-emitting clump in the Central Molecular Zone and may be forming the closest analog to a super star cluster in the Galaxy. The Central Molecular Zone (CMZ) contains the highest concentration of mm and sub-mm emitting dense clumps in the Galaxy. Most 1.1 mm features at positive longitudes are seen in silhouette against the 3.6 to 24 {\mu}m background observed by the Spitzer Space Telescope. However, only a few clumps at negative longitudes are seen in absorption, confirming the hypothesis that positive longitude clumps in the CMZ tend to be on the near-side of the Galactic center, consistent with the suspected orientation of the central bar in our Galaxy. Some 1.1 mm cloud surfaces are seen in emission at 8 {\mu}m, presumably due to polycyclic aromatic hydrocarbons (PAHs). A ~0.2\degree (~30 pc) diameter cavity and infrared bubble between l \approx 0.0\degree and 0.2\degree surrounds the Arches and Quintuplet clusters and Sgr A. The bubble contains several clumpy dust filaments that point toward Sgr A\ast; its potential role in their formation is explored. [abstract truncated]

We present the Bolocam Galactic Plane Survey (BGPS), a 1.1 mm continuum survey at 33" effective resolution of 170 square degrees of the Galactic Plane visible from the northern hemisphere. The survey is contiguous over the range -10.5 < l < 90.5, |b| < 0.5 and encompasses 133 square degrees, including some extended regions |b| < 1.5. In addition to the contiguous region, four targeted regions in the outer Galaxy were observed: IC1396, a region towards the Perseus Arm, W3/4/5, and Gem OB1. The BGPS has detected approximately 8400 clumps over the entire area to a limiting non-uniform 1-sigma noise level in the range 11 to 53 mJy/beam in the inner Galaxy. The BGPS source catalog is presented in a companion paper (Rosolowsky et al. 2010). This paper details the survey observations and data reduction methods for the images. We discuss in detail the determination of astrometric and flux density calibration uncertainties and compare our results to the literature. Data processing algorithms that separate astronomical signals from time-variable atmospheric fluctuations in the data time-stream are presented. These algorithms reproduce the structure of the astronomical sky over a limited range of angular scales and produce artifacts in the vicinity of bright sources. Based on simulations, we find that extended emission on scales larger than about 5.9' is nearly completely attenuated (> 90%) and the linear scale at which the attenuation reaches 50% is 3.8'. Comparison with other millimeter-wave data sets implies a possible systematic offset in flux calibration, for which no cause has been discovered. This presentation serves as a companion and guide to the public data release through NASA's Infrared Processing and Analysis Center (IPAC) Infrared Science Archive (IRSA). New data releases will be provided through IPAC IRSA with any future improvements in the reduction.

We examine the physical properties and evolutionary stages of a sample of 17 clumps within 8 Infrared Dark Clouds (IRDCs) by combining existing infrared, millimeter, and radio data with new Bolocam Galactic Plane Survey (BGPS) 1.1 mm data, VLA radio continuum data, and HHT dense gas (HCO+ and N2H+) spectroscopic data. We combine literature studies of star formation tracers and dust temperatures within IRDCs with our search for ultra-compact (UC) HII regions to discuss a possible evolutionary sequence for IRDC clumps. In addition, we perform an analysis of mass tracers in IRDCs and find that 8 micron extinction masses and 1.1 mm Bolocam Galactic Plane Survey (BGPS) masses are complementary mass tracers in IRDCs except for the most active clumps (notably those containing UCHII regions), for which both mass tracers suffer biases. We find that the measured virial masses in IRDC clumps are uniformly higher than the measured dust continuum masses on the scale of ~1 pc. We use 13CO, HCO+, and N2H+ to study the molecular gas properties of IRDCs and do not see any evidence of chemical differentiation between hot and cold clumps on the scale of ~1 pc. However, both HCO+ and N2H+ are brighter in active clumps, due to an increase in temperature and/or density. We report the identification of four UCHII regions embedded within IRDC clumps and find that UCHII regions are associated with bright (>1 Jy) 24 micron point sources, and that the brightest UCHII regions are associated with "diffuse red clumps" (an extended enhancement at 8 micron). The broad stages of the discussed evolutionary sequence (from a quiescent clump to an embedded HII region) are supported by literature dust temperature estimates; however, no sequential nature can be inferred between the individual star formation tracers.

We present the 1.1 millimeter Bolocam Galactic Plane Survey (BGPS) observations of the Gemini OB1 molecular cloud complex, and targeted ammonia observations of the BGPS sources. When paired with molecular spectroscopy of a dense gas tracer, millimeter observations yield physical properties such as masses, radii, mean densities, kinetic temperatures and line widths. We detect 34 distinct BGPS sources above 5-sigma=0.37 Jy/beam with corresponding 5-sigma detections in the ammonia (1,1) transition. Eight of the objects show water maser emission (20%). We find a mean millimeter source FWHM of 1.12 pc, and a mean kinetic temperature of 20 K for the sample of 34 BGPS sources. The observed ammonia line widths are dominated by non-thermal motions, typically found to be a few times the thermal sound speed expected for the derived kinetic temperature. We calculate the mass for each source from the millimeter flux assuming the sources are isothermal and find a mean isothermal mass within a 120" aperture of 230 +/- 180 solar masses. We find a total mass of 8,400 solar masses for all BGPS sources in the Gemini OB1 molecular cloud, representing 6.5% of the cloud mass. By comparing the millimeter isothermal mass to the virial mass within a radius equal to the mm source size calculated from the ammonia line widths, we find a mean virial parameter (M_vir/M_iso) of 1.0 +/- 0.9 for the sample. We find mean values for the distributions of column densities of 10^22 cm^-2 for H_2, and 3.0x10^14 cm^-2 for ammonia, giving a mean ammonia abundance of 3.0x10^-8 relative to H_2. We find volume-averaged densities on the order of 10^3-10^4 cm^-3. The sizes and densities suggest that in the Gem OB1 region the BGPS is detecting the clumps from which stellar clusters form, rather than smaller, higher density cores where single stars or small multiple systems form.

Aims: To explore the infrared and radio properties of one of the closest Galactic starburst regions. Methods: Images obtained with the Herschel Space Observatory at wavelengths of 70, 160, 250, 350, and 500 microns using the PACS and SPIRE arrays are analyzed and compared with radio continuum VLA data and 8 micron images from the Spitzer Space Telescope. The morphology of the far-infrared emission is combined with radial velocity measurements of millimeter and centimeter wavelength transitions to identify features likely to be associated with the W43 complex. Results: The W43 star-forming complex is resolved into a dense cluster of protostars, infrared dark clouds, and ridges of warm dust heated by massive stars. The 4 brightest compact sources with L > 1.5 x 10^4 Lsun embedded within the Z-shaped ridge of bright dust emission in W43 remain single at 4" (0.1 pc) resolution. These objects, likely to be massive protostars or compact clusters in early stages of evolution are embedded in clumps with masses of 10^3 to 10^4 Msun, but contribute only 2% to the 3.6 x 10^6 Lsun far-IR luminosity of W43 measured in a 16 by 16 pc box. The total mass of gas derived from the far-IR dust emission inside this region is ~10^6 Msun. Cometary dust clouds, compact 6 cm radio sources, and warm dust mark the locations of older populations of massive stars. Energy release has created a cavity blowing-out below the Galactic plane. Compression of molecular gas in the plane by the older HII region near G30.684-0.260 and the bipolar structure of the resulting younger W43 HII region may have triggered the current mini-star burst.

We present a catalog of 8358 sources extracted from images produced by the Bolocam Galactic Plane Survey (BGPS). The BGPS is a survey of the millimeter dust continuum emission from the northern Galactic plane. The catalog sources are extracted using a custom algorithm, Bolocat, which was designed specifically to identify and characterize objects in the large-area maps generated from the Bolocam instrument. The catalog products are designed to facilitate follow-up observations of these relatively unstudied objects. The catalog is 98% complete from 0.4 Jy to 60 Jy over all object sizes for which the survey is sensitive (<3.5'). We find that the sources extracted can best be described as molecular clumps -- large dense regions in molecular clouds linked to cluster formation. We find the flux density distribution of sources follows a power law with dN/dS ~S^(-2.4 +/- 0.1) and that the mean Galactic latitude for sources is significantly below the midplane: <b>=(-0.095 +/- 0.001) deg.

(abridged) We study photoelectric heating throughout the Large Magellanic Cloud. We quantify the importance of the [CII] cooling line and the photoelectric heating process of various environments in the LMC and investigate which parameters control the extent of photoelectric heating. We use the BICE [CII] map and the Spitzer/SAGE infrared maps. We examine the spatial variations in the efficiency of photoelectric heating: photoelectric heating rate over power absorbed by grains. We correlate the photoelectric heating efficiency and the emission from various dust constituents and study the variations as a function of H\alpha emission, dust temperatures, and the total infrared luminosity. From this we estimate radiation field, gas temperature, and electron density. We find systematic variations in photoelectric efficiency. The highest efficiencies are found in the diffuse medium, while the lowest coincide with bright star-forming regions (~1.4 times lower). The [CII] line emission constitutes 1.32% of the far infrared luminosity across the whole of the LMC. We find correlations between the [CII] emission and ratios of the mid infrared and far infrared bands, which comprise various dust constituents. The correlations are interpreted in light of the spatial variations of the dust abundance and by the local environmental conditions that affect the dust emission properties. As a function of the total infrared surface brightness, S_{TIR}, the [CII] surface brightness can be described as: S_{[CII]}=1.25 S_{TIR}^{0.69} [10^{-3} erg s^{-1} cm^{-2} sr^{-1}]. The [CII] emission is well-correlation with the 8 micrometer emission, suggesting that the polycyclic aromatic hydrocarbons play a dominant role in the photoelectric heating process.

context: The spectrum of IRAS 16594-4656 shows shock excited H_2 emission and collisionally excited emission lines such as[O I],[C I],and [Fe II]. aim: The goal is to determine the location of the H_2 and [Fe II] shock emission, to determine the shock velocities,and constrain the physical properties in the shock. methods: High resolution spectra of the H_2 1-0 S(1),H_2 2-1 S(1), [Fe II], and Pa$\beta$ emission lines were obtained with the near infrared spectrograph Phoenix on Gemini South. results: The position-velocity diagrams of H_2 1-0 S(1), H_2 2-1 S(1), and [Fe II] are presented. The H_2 and [Fe II] emission is spatially extended. The collisionally excited [O I] and [C I] optical emission lines have a similar double peaked profile compared to the extracted H_2 profile and appear to be produced in the same shock. They all indicate an expansion velocity of ~8 km/s and the presence of a neutral, very high density region with $n_{\rm e}$ about 3 x 10^6 to 5 x10^7 cm$^{-3}$. The [Fe II] emission however is single peaked. It has a gaussian FWHM of 30 km/s and a total width of 62 km/s at 1% of the peak. The Pa$\beta$ profile is even wider with a gaussian FWHM of 48 km/s and a total width of 75 km/s at 1% of the peak. conclusions: The H$_2$ emission is excited in a slow 5 to 20 km/s shock into dense material at the edge of the lobes, caused by the interaction of the AGB ejecta and the post-AGB wind. The 3D representation of the H_2 data shows a hollow structure with less H_2 emission in the equatorial region. The [Fe II] emission is not present in the lobes, but originates close to the central star in fast shocks in the post-AGB wind or in a disk. The Pa$\beta$ emission also appears to originate close to the star.