We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations of the tadpole, a small globule in the Carina Nebula that hosts the HH 900 jet+outflow system. Our data include $^{12}$CO, $^{13}$CO, C$^{18}$O J=2-1, $^{13}$CO, C$^{18}$O J=3-2, and serendipitous detections of DCN J=3-2 and CS J=7-6. With angular resolution comparable to the Hubble Space Telescope (HST), our data reveal for the first time the bipolar molecular outflow in CO, seen only inside the globule, that is launched from the previously unseen jet-driving protostar (the HH 900 YSO). The biconical morphology joins smoothly with the externally irradiated outflow seen in ionized gas tracers outside the globule, tracing the overall morphology of a jet-driven molecular outflow. Continuum emission at the location of the HH 900 YSO appears to be slightly flattened perpendicular to outflow axis. Model fits to the continuum have a best-fit spectral index of $\sim 2$, suggesting cold dust and the onset of grain growth. In position-velocity space, $^{13}$CO and C$^{18}$O gas kinematics trace a C-shaped morphology, similar to infall profiles seen in other sources, although the global dynamical behaviour of the gas remains unclear. Line profiles of the CO isotopologues display features consistent with externally heated gas. We estimate a globule mass of $\sim 1.9$ M$_{\odot}$, indicating a remaining lifetime of $\sim 4$ Myr, assuming a constant photoevaporation rate. This long globule lifetime will shield the disk from external irradiation perhaps prolonging its life and enabling planet formation in regions where disks are typically rapidly destroyed.

Complex organic molecules (COMs) are thought to form on icy dust grains in the earliest phase of star formation. The evolution of these COMs from the youngest Class 0/I protostellar phases toward the more evolved Class II phase is still not fully understood. Since planet formation seems to start early, and mature disks are too cold for characteristic COM emission lines, studying the inventory of COMs on solar system scales in the Class 0/I stage is relevant. ALMA Band 3 (3 mm) and Band 6 (1 mm) observations are obtained of seven Class 0 protostars in the Perseus and Serpens star-forming regions. By modeling the inner protostellar region using 'LTE' models, the excitation temperature and column densities are determined for several O-bearing COMs. B1-c, B1-bS, and Serpens S68N show COM emission, i.e, three out of the seven sources. No clear correlation seems to exist between the occurrence of COMs and source luminosity. The abundances of several COMs with respect to CH3OH are remarkably similar for the three COM-rich sources, and to IRAS 16293-2422B and HH 212. For other COMs the abundances differ by up to an order of magnitude, indicating that local source conditions are case determining. B1-c hosts a cold ($T_{ex}\approx60$ K), more extended component of COM emission with a column density of typically a few % of the warm/hot ($T_{ex}\sim 200$ K), central component. A D/H ratio of 1-3 % is derived based on the CH2DOH/CH3OH ratio suggesting a temperature of $\sim$15~K during the formation of methanol. This ratio is consistent with other low-mass protostars. Future mid-infrared facilities such as JWST/MIRI will be essential to directly observe COM ices. Combining this with a larger sample of COM-rich sources with ALMA will allow for directly linking ice and gas-phase abundances in order to constrain the routes that produce and maintain chemical complexity during the star formation process.

CMB-HD is a proposed ultra-deep (0.5 uk-arcmin), high-resolution (15 arcseconds) millimeter-wave survey over half the sky that would answer many outstanding questions in both fundamental physics of the Universe and astrophysics. This survey would be delivered in 7.5 years of observing 20,000 square degrees, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors.

Our aim was to directly detect the thermal emission of the putative exoring system responsible for the complex deep transits observed in the light curve for the young Sco-Cen star 1SWASP J140747.93-394542.6 (V1400 Cen, hereafter J1407), confirming it as the occulter seen in May 2007, and to determine its orbital parameters with respect to the star. We used the Atacama Large Millimeter/submillimeter Array (ALMA) to observe the field centred on J1407 in the 340 GHz (Band 7) continuum in order to determine the flux and astrometric location of the ring system relative to the star. We used the VLT/NACO camera to observe the J1407 system in March 2019 and to search for the central planetary mass object at thermal infrared wavelengths. We detect no point source at the expected location of J1407, and derive an upper limit $3\sigma$ level of $57.6~\mu\rm{Jy}$. There is a point source detected at an angular separation consistent with the expected location for a free-floating ring system that occulted J1407 in May 2007, with a flux of $89~\mu\rm{Jy}$ consistent with optically thin dust surrounding a massive substellar companion. At 3.8 microns with the NACO camera, we detect the star J1407 but no other additional point sources within 1.3 arcseconds of the star, with a lower bound on the sensitivity of $6M_{Jup}$ at the location of the ALMA source, and down to $4M_{Jup}$ in the sky background limit. The ALMA upper limit at the location of J1407 implies that a hypothesised bound ring system is composed of dust smaller than $1\rm{~mm}$ in size, implying a young ring structure. The detected ALMA source has multiple interpretations, including: (i) it is an unbound substellar object surrounded by warm dust in Sco-Cen with an upper mass limit of $6M_{Jup}$, or (ii) it is a background galaxy.

Scattered light high-resolution imaging of the proto-planetary disc orbiting HD100453 shows two symmetric spiral arms, possibly launched by an external stellar companion. In this paper we present new, sensitive high-resolution ($\sim$30 mas) Band 7 ALMA observations of this source. This is the first source where we find counterparts in the sub-mm continuum to both scattered light spirals. The CO J=3-2 emission line also shows two spiral arms; in this case they can be traced over a more extended radial range, indicating that the southern spiral arm connects to the companion position. This is clear evidence that the companion is responsible for launching the spirals. The pitch angle of the sub-millimeter continuum spirals ($\sim 6 ^{\circ}$) is lower than the one in scattered light ($\sim 16 ^{\circ}$). We show that hydrodynamical simulations of binary-disc interaction can account for the difference in pitch angle only if one takes into account that the midplane is colder than the upper layers of the disc, as expected for the case of externally irradiated discs.

Forming high-mass stars have a significant effect on their natal environment. Their feedback pathways, including winds, outflows, and ionising radiation, shape the evolution of their surroundings which impacts the formation of the next generation of stars. They create or reveal dense pillars of gas and dust towards the edges of the cavities they clear. They are modelled in feedback simulations, and the sizes and shapes of the pillars produced are consistent with those observed. However, these models predict measurably different kinematics which provides testable discriminants. Here we present the first ALMA Compact Array (ACA) survey of 13 pillars in Carina, observed in $^{12}$CO, $^{13}$CO and C$^{18}$O J=2-1, and the 230 GHz continuum. The pillars in this survey were chosen to cover a wide range in properties relating to the amount and direction of incident radiation, proximity to nearby irradiating clusters and cloud rims, and whether they are detached from the cloud. With these data, we are able to discriminate between models. We generally find pillar velocity dispersions of $<$ 1 km s$^{-1}$ and that the outer few layers of molecular emission in these pillars show no significant offsets from each other, suggesting little bulk internal motions within the pillars. There are instances where the pillars are offset in velocity from their parental cloud rim, and some with no offset, hinting at a stochastic development of these motions.

We present a detailed observational and modeling study of the hot core VLA 3 in the high-mass star-forming region AFGL 2591, which is a target region of the NOrthern Extended Millimeter Array (NOEMA) large program CORE. Using NOEMA observations at 1.37 mm with an angular resolution of ~0."42 (1 400 au at 3.33 kpc), we derived the physical and chemical structure of the source. We modeled the observed molecular abundances with the chemical evolution code MUSCLE (MUlti Stage ChemicaL codE). Results. With the kinetic temperature tracers CH3CN and H2CO we observe a temperature distribution with a power-law index of q = 0.41+-0.08. Using the visibilities of the continuum emission we derive a density structure with a power-law index of p = 1.7+-0.1. The hot core spectra reveal high molecular abundances and a rich diversity in complex molecules. The majority of the molecules have an asymmetric spatial distribution around the forming protostar(s), which indicates a complex physical structure on scales < 1 400 au. Using MUSCLE, we are able to explain the observed molecular abundance of 10 out of 14 modeled species at an estimated hot core chemical age of ~21 100 years. In contrast to the observational analysis, our chemical modeling predicts a lower density power-law index of p < 1.4. Reasons for this discrepancy are discussed. Conclusions. Combining high spatial resolution observations with detailed chemical modeling allows us to derive a concise picture of the physical and chemical structure of the famous AFGL 2591 hot core. The next steps are to conduct a similar analysis for the whole CORE sample, and then use this analysis to constrain the chemical diversity in high-mass star formation to a much greater depth.

We present new MUSE/VLT observations of a small globule in the Carina H II region that hosts the HH 900 jet+outflow system. Data were obtained with the GALACSI ground-layer adaptive optics system in wide-field mode, providing spatially-resolved maps of diagnostic emission lines. These allow us to measure the variation of the physical properties in the globule and jet+outflow system. We find high temperatures ($T_e \approx 10^4$ K), modest extinction ($A_V \approx 2.5$ mag), and modest electron densities ($n_e \approx 200$ cm$^{-3}$) in the ionized gas. Higher excitation lines trace the ionized outflow; both the excitation and ionization in the outflow increase with distance from the opaque globule. In contrast, lower excitation lines that are collisionally de-excited at densities $\gtrsim 10^4$ cm$^{-3}$ trace the highly collimated protostellar jet. Assuming the globule is an isothermal sphere confined by the pressure of the ionization front, we compute a Bonnor-Ebert mass of $\sim 3.7$ M$_{\odot}$. This is two orders of magnitude higher than previous mass estimates, calling into question whether small globules like the Tadpole contribute to the bottom of the IMF. Derived globule properties are consistent with a cloud that has been and/or will be compressed by the ionization front on its surface. At the estimated globule photoevaporation rate of $\sim 5 \times 10^{-7}$ M$_{\odot}$ yr$^{-1}$, the globule will be completely ablated in $\sim 7$ Myr. Stars that form in globules like the Tadpole will emerge into the H II later and may help resolve some of the temporal tension between disk survival and enrichment.

The sub-mm sky is a unique window for probing the architecture of the Universe and structures within it. From the discovery of dusty sub-mm galaxies, to the ringed nature of protostellar disks, our understanding of the formation, destruction, and evolution of objects in the Universe requires a comprehensive view of the sub-mm sky. The current generation single-dish sub-mm facilities have shown of the potential for discovery, while interferometers have presented a high resolution view into the finer details. However, our understanding of large-scale structure and our full use of these interferometers is now hampered by the limited sensitivity of our sub-mm view of the universe at larger scales. Thus, now is the time to start planning the next generation of sub-mm single dish facilities, to build on these revolutions in our understanding of the sub-mm sky. Here we present the case for the Atacama Large Aperture Submillimeter Telescope (AtLAST), a concept for a 50m class single dish telescope. We envision AtLAST as a facility operating as an international partnership with a suite of instruments to deliver the transformative science described in many Astro2020 science white papers. A 50m telescope with a high throughput and 1$^\circ$ FoV with a full complement of advanced instrumentation, including highly multiplexed high-resolution spectrometers, continuum cameras and Integral Field Units, AtLAST will have mapping speeds thousands of times greater than any current or planned facility. It will reach confusion limits below $L_*$ in the distant universe and resolve low-mass protostellar cores at the distance of the Galactic Center, providing synergies with upcoming facilities across the spectrum. Located on the Atacama plateau, to observe frequencies un-obtainable by other observatories, AtLAST will enable a fundamentally new understanding of the sub-mm universe at unprecedented depths.

The formation process of high-mass stars (>8M$_\odot$) is poorly constrained, particularly, the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500M$_\odot$) star-forming region IRAS 23033+5951. Using the Northern Extended Millimeter Array (NOEMA) in three configurations and the IRAM 30-m single-dish telescope at 220GHz, we probe the gas and dust emission at an angular resolution of ~0.45arcsec, corresponding to 1900au. In the mm continuum emission, we identify a protostellar cluster with at least four mm-sources, where three of them show a significantly higher peak intensity well above a signal-to-noise ratio of 100. Hierarchical fragmentation from large to small spatial scales is discussed. Two fragments are embedded in rotating structures and drive molecular outflows, traced by $^{13}$CO (2-1) emission. The velocity profiles across two of the cores are similar to Keplerian but are missing the highest velocity components close to the center of rotation, which is a common phenomena from observations like these, and other rotation scenarios are not excluded entirely. Position-velocity diagrams suggest protostellar masses of ~6 and 19M$_\sun$. Rotational temperatures from fitting CH$_3$CN ($12_K-11_K$) spectra are used for estimating the gas temperature and by that the disk stability against gravitational fragmentation, utilizing Toomre's $Q$ parameter. [We] identify only one candidate disk to be unstable against gravitational instability caused by axisymmetric perturbations. The dominant sources cover different evolutionary stages within the same maternal gas clump. The appearance of rotation and outflows of the cores are similar to those found in low-mass star-forming regions.

(abridged) The Atacama Large Aperture Submillimeter Telescope (AtLAST) project aims to build a 50-m-class submm telescope with $>1^\circ$ field of view, high in the Atacama Desert, providing fast and detailed mapping of the mm/submm sky. It will thus serve as a strong complement to existing facilities such as ALMA. ALMA's small field of view ($<15^{\prime\prime}$ at 350 GHz) limits its mapping speed for large surveys. Instead, a single dish with a large field of view such as the AtLAST concept can host large multi-element instruments that can more efficiently map large portions of the sky. Small aperture survey instruments (typically much smaller than $<3\times$ the size of an interferometric array element) can mitigate this somewhat but lack the resolution for accurate recovery of source location and have small collecting areas. Furthermore, small aperture survey instruments do not provide sufficient overlap in the spatial scales they sample to provide a complete reconstruction of extended sources (i.e.\ the zero-spacing information is incomplete in $u,v$-space.) The heterodyne instrumentation for the AtLAST telescope that we consider here will take advantage of extensive developments in the past decade improving the performance and pixel count of heterodyne focal plane arrays. Such instrumentation, with higher pixel counts, has alredy begun to take advantage of integration in the focal planes to increase packaging efficiency over simply stacking modular mixer blocks in the focal plane. We extrapolate from the current state-of-the-art to present concept first-generation heterodyne designs for AtLAST.

A millimeter-wave survey over half the sky, that spans frequencies in the range of 30 to 350 GHz, and that is both an order of magnitude deeper and of higher-resolution than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. By providing such a deep, high-resolution millimeter-wave survey (about 0.5 uK-arcmin noise and 15 arcsecond resolution at 150 GHz), CMB-HD will enable major advances. It will allow 1) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10/hMpc), which probes dark matter particle properties. It will also allow 2) measurements of the thermal and kinetic Sunyaev-Zel'dovich effects on small scales to map the gas density and gas pressure profiles of halos over a wide field, which probes galaxy evolution and cluster astrophysics. In addition, CMB-HD would allow us to cross critical thresholds in fundamental physics: 3) ruling out or detecting any new, light (< 0.1eV), thermal particles, which could potentially be the dark matter, and 4) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe. Such a survey would also 5) monitor the transient sky by mapping the full observing region every few days, which opens a new window on gamma-ray bursts, novae, fast radio bursts, and variable active galactic nuclei. Moreover, CMB-HD would 6) provide a census of planets, dwarf planets, and asteroids in the outer Solar System, and 7) enable the detection of exo-Oort clouds around other solar systems, shedding light on planet formation. CMB-HD will deliver this survey in 5 years of observing half the sky, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. The telescopes will field about 2.4 million detectors (600,000 pixels) in total.

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.

This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additional single-dish maps, merged with the interferometric data to recover the extended emission. Our correlator setup covered a number of lines from well-known hot core tracers and a few outflow tracers. The angular (~0.45"$-$0.9") and spectral (0.5 km/s) resolutions were sufficient to resolve the clump in IRAS23385+6053 and investigate the existence of large-scale motions due to rotation, infall, or expansion. We find that the clump splits into six distinct cores when observed at sub-arcsecond resolution. These are identified through their 1.4 mm continuum and molecular line emission. We produce maps of the velocity, line width, and rotational temperature from the methanol and methyl cyanide lines, which allow us to investigate the cores and reveal a velocity and temperature gradient in the most massive core. We also find evidence of a bipolar outflow, possibly powered by a low-mass star. We present the tentative detection of a circumstellar self-gravitating disk lying in the most massive core and powering a large-scale outflow previously known in the literature. In our scenario, the star powering the flow is responsible for most of the luminosity of IRAS23385+6053 (~$3000~L_\odot$). The other cores, albeit with masses below the corresponding virial masses, appear to be accreting material from their molecular surroundings and are possibly collapsing or on the verge of collapse. We conclude that we are observing a sample of star-forming cores that is bound to turn into a cluster of massive stars.

The majority of the ultimate main-sequence mass of a star is assembled in the protostellar phase, where a forming star is embedded in an infalling envelope and encircled by a protoplanetary disk. Studying mass accretion in protostars is thus a key to understanding how stars gain their mass and ultimately how their disks and planets form and evolve. At this early stage, the dense envelope reprocesses most of the luminosity generated by accretion to far-infrared and submillimeter wavelengths. Time-domain photometry at these wavelengths is needed to probe the physics of accretion onto protostars, but variability studies have so far been limited, in large part because of the difficulty in accessing these wavelengths from the ground. We discuss the scientific progress that would be enabled with far-infrared and submillimeter programs to probe protostellar variability in the nearest kiloparsec.

We present the science case for mapping several thousand galaxy (proto)clusters at z=1-10 with a large aperture single dish sub-mm facility, producing a high-redshift counterpart to local large surveys of rich clusters like the well-studied Abell catalogue. Principal goals of a large survey of distant clusters are the evolution of galaxy clusters over cosmic time and the impact of environment on the evolution and formation of galaxies. To make a big leap forward in this emerging research field, the community would benefit from a large-format, wide-band, direct-detection spectrometer (e.g., based on MKID technology), covering a wide field of ~1 square degree and a frequency coverage from 70 to 700 GHz.

The Sloan Digital Sky Survey (SDSS) was revolutionary because of the extraordinary breadth and ambition of its optical imaging and spectroscopy. We argue that a 'sub-millimeter SDSS' - a sensitive large-area imaging+spectroscopic survey in the sub-mm window - will revolutionize our understanding of galaxy evolution in the early Universe. By detecting the thermal dust continuum emission and atomic and molecular line emission of galaxies out to z~10 it will be possible to measure the redshifts, star formation rates, dust and gas content of hundreds of thousands of high-z galaxies down to ~L*. Many of these galaxies will have counterparts visible in the deep optical imaging of the Large Synoptic Survey Telescope. This 3D map of galaxy evolution will span the peak epoch of galaxy formation all the way back to cosmic dawn, measuring the co-evolution of the star formation rate density and molecular gas content of galaxies, tracking the production of metals and charting the growth of large-scale structure.

Sunyaev-Zeldovich (SZ) effects were first proposed in the 1970s as tools to identify the X-ray emitting hot gas inside massive clusters of galaxies and obtain their velocities relative to the cosmic microwave background (CMB). Yet it is only within the last decade that they have begun to significantly impact astronomical research. Thanks to the rapid developments in CMB instrumentation, measurement of the dominant thermal signature of the SZ effects has become a routine tool to find and characterize large samples of galaxy clusters and to seek deeper understanding of several important astrophysical processes via high-resolution imaging studies of many targets. With the notable exception of the Planck satellite and a few combinations of ground-based observatories, much of this "SZ revolution" has happened in the photometric mode, where observations are made at one or two frequencies in the millimeter regime to maximize the cluster detection significance and minimize the foregrounds. Still, there is much more to learn from detailed and systematic analyses of the SZ spectra across multiple wavelengths, specifically in the submillimeter (>300 GHz) domain. The goal of this Science White Paper is to highlight this particular aspect of SZ research, point out what new and potentially groundbreaking insights can be obtained from these studies, and emphasize why the coming decade can be a golden era for SZ spectral measurements.

The past two decades have seen extensive surveys of the far-infrared to submillimeter continuum emission in the plane of our Galaxy. We line out prospects for the coming decade for corresponding molecular and atomic line surveys which are needed to fully understand the formation of the dense structures that give birth to clusters and stars out of the diffuse interstellar medium. We propose to work towards Galaxy wide surveys in mid-J CO lines to trace shocks from colliding clouds, Galaxy-wide surveys for atomic Carbon lines in order to get a detailed understanding of the relation of atomic and molecular gas in clouds, and to perform extensive surveys of the structure of the dense parts of molecular clouds to understand the importance of filaments/fibers over the full range of Galactic environments and to study how dense cloud cores are formed from the filaments. This work will require a large (50m) Single Dish submillimeter telescope equipped with massively multipixel spectrometer arrays, such as envisaged by the AtLAST project.

The cycling of baryons in and out of galaxies is what ultimately drives galaxy formation and evolution. The circumgalactic medium (CGM) represents the interface between the interstellar medium and the cosmic web, hence its properties are directly shaped by the baryon cycle. Although traditionally the CGM is thought to consist of warm and hot gas, recent breakthroughs are presenting a new scenario according to which an important fraction of its mass may reside in the cold atomic and molecular phase. This would represent fuel that is readily available for star formation, with crucial implications for feeding and feedback processes in galaxies. However, such cold CGM, especially in local galaxies where its projected size on sky is expected to be of several arcminutes, cannot be imaged by ALMA due to interferometric spatial scale filtering of large-scale structures. We show that the only way to probe the multiphase CGM including its coldest component is through a large (e.g. 50-m) single dish (sub-)mm telescope.

Opening up a new window of millimeter-wave observations that span frequency bands in the range of 30 to 500 GHz, survey half the sky, and are both an order of magnitude deeper (about 0.5 uK-arcmin) and of higher-resolution (about 10 arcseconds) than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. In particular, such a survey would allow for major advances in measuring the distribution of dark matter and gas on small-scales, and yield needed insight on 1.) dark matter particle properties, 2.) the evolution of gas and galaxies, 3.) new light particle species, 4.) the epoch of inflation, and 5.) the census of bodies orbiting in the outer Solar System.

The Sunyaev-Zeldovich (SZ) effect was first predicted nearly five decades ago, but has only recently become a mature tool for performing high resolution studies of the warm and hot ionized gas in and between galaxies, groups, and clusters. Galaxy groups and clusters are powerful probes of cosmology, and they also serve as hosts for roughly half of the galaxies in the Universe. In this white paper, we outline the advances in our understanding of thermodynamic and kinematic properties of the warm-hot universe that can come in the next decade through spatially and spectrally resolved measurements of the SZ effects. Many of these advances will be enabled through new/upcoming millimeter/submillimeter (mm/submm) instrumentation on existing facilities, but truly transformative advances will require construction of new facilities with larger fields of view and broad spectral coverage of the mm/submm bands.

We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. We detect SiO (5-4) emission that is marginally resolved and elongated in a direction perpendicular to the large-scale outflow seen in the 13CO (2-1) line using the main ALMA array in conjunction with the Atacama Compact Array (ACA). Morphologically, the SiO appears to represent a disc-like structure. Using parametric models we show that the position-velocity profile of the SiO is consistent with the Keplerian rotation of a disc around an object between 10-30Mo in mass, only if there is also radial expansion from a separate structure. The radial motion component can be interpreted as a disc wind from the disc surface. Models with a central stellar object mass between 20 and 30Mo are the most consistent with the stellar luminosity (100000 Lo) and indicative of an O-type star. The H30a millimetre recombination line (231.9GHz) is also detected, but spatially unresolved, and is indicative of a very compact, hot, ionised region co-spatial with the dust continuum core. Accounting for all observables, we suggest that G17.64 is consistent with a O-type young stellar object in the final stages of protostellar assembly, driving a wind, but that has not yet developed into a compact HII region. The existance and detection of the disc in G17.64 is likely related to its isolated and possibly more evolved nature, traits which may underpin discs in similar sources.

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

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

This work aims to understand which midplane conditions are probed by the DCO$^+$ emission in the disk around the Herbig Ae star HD 169142. We explore the sensitivity of the DCO$^+$ formation pathways to the gas temperature and the CO abundance. The DCO$^+$ $J$=3-2 transition was observed with ALMA at a spatial resolution of 0.3". The HD 169142 DCO$^+$ radial intensity profile reveals a warm, inner component at radii <30 AU and a broad, ring-like structure from ~50-230 AU with a peak at 100 AU just beyond the millimeter grain edge. We modeled DCO$^+$ emission in HD 169142 with a physical disk structure adapted from the literature, and employed a simple deuterium chemical network to investigate the formation of DCO$^+$ through the cold deuterium fractionation pathway via H$_2$D$^+$. Contributions from the warm deuterium fractionation pathway via CH$_2$D$^+$ are approximated using a constant abundance in the intermediate disk layers. Parameterized models show that alterations to the midplane gas temperature and CO abundance of the literature model are both needed to recover the observed DCO$^+$ radial intensity profile. The best-fit model contains a shadowed, cold midplane in the region z/r < 0.1 with an 8 K decrease in gas temperature and a factor of five CO depletion just beyond the millimeter grain edge, and a 2 K decrease in gas temperature for r > 120 AU. The warm deuterium fractionation pathway is implemented as a constant DCO$^+$ abundance of 2.0$\times$10$^{-12}$ between 30-70 K. The DCO$^+$ emission probes a reservoir of cold material in the HD 169142 outer disk that is not revealed by the millimeter continuum, the SED, nor the emission from the 12CO, 13CO, or C18O $J$=2-1 lines.

Highly collimated parsec-scale jets, generally linked to the presence of an accretion disk, are a commonly observed phenomenon from revealed low-mass young stellar objects. In the past two decades, only a very few of these objects have been directly (or indirectly) observed towards high-mass (M > 8 M$_{\odot}$) young stellar objects, adding to the growing evidence that disk-mediated accretion is a phenomenon also occurring in high-mass stars, the formation mechanism of which is still poorly understood. Of the observed jets from massive young stars, none is in the optical regime (due to these being typically highly obscured by their native material), and none are found outside of the Milky Way. Here, we report the detection of HH 1177, the first extragalactic optical ionized jet originating from a massive young stellar object located in the Large Magellanic Cloud. The jet is highly collimated over the entire measured extent of at least 10 pc, and has a bipolar geometry. The presence of a jet indicates ongoing, disk-mediated accretion, and together with the high degree of collimation, this system is therefore likely to be an up-scaled version of low-mass star formation. We conclude that the physics governing jet launching and collimation is independent of stellar mass.

High-mass stars form in much richer environments than those associated with isolated low-mass stars, and once they reach a certain mass, produce ionised (HII) regions. The formation of these pockets of ionised gas are unique to the formation of high-mass stars (M $>8$ M$_\odot$), and present an excellent opportunity to study the final stages of accretion, which could include accretion through the HII region itself. This study of the dynamics of the gas on both sides of these ionisation boundaries in very young HII regions aims to quantify the relationship between the HII regions and their immediate environments.We present high-resolution ($\sim$ 0.5$"$) ALMA observations of nine HII regions selected from the Red MSX Source (RMS) survey with compact radio emission and bolometric luminosities greater than 10$^4$ L$_\odot$. We focus on the initial presentation of the data, including initial results from the radio recombination line H29$\alpha$, some complementary molecules, and the 256 GHz continuum emission. Of the six (out of nine) regions with H29$\alpha$ detections, two appear to have cometary morphologies with velocity gradients across them, and two appear more spherical with velocity gradients suggestive of infalling ionised gas. The remaining two were either observed at low resolution or had signals that were too weak to draw robust conclusions. We also present a description of the interactions between the ionised and molecular gas (as traced by CS (J=5-4)), often (but not always) finding theHII region had cleared its immediate vicinity of molecules. Of our sample of nine, the observations of the two clusters expected to have the youngest HII regions (from previous radio observations) are suggestive of having infalling motions in the H29$\alpha$ emission, which could be indicative of late stage accretion onto the stars despite the presence of an HII region.

The protoplanetary system HD 169142 is one of the few cases where a potential candidate protoplanet has been recently detected via direct imaging. To study the interaction between the protoplanet and the disk itself observations of the gas and dust surface density structure are needed. This paper reports new ALMA observations of the dust continuum at 1.3\,mm, $^{12}$CO, $^{13}$CO and C$^{18}$O $J=2-1$ emission from the system HD 169142 at angular resolution of $\sim 0".18 - 0".28$ ($\sim 20\,$au$ - 33\,$au). The dust continuum emission reveals a double-ring structure with an inner ring between $0".17-0".28$ ($\sim 20 - 35\,$au) and an outer ring between $0".48-0".64$ ($\sim 56 - 83\,$au). The size and position of the inner ring is in good agreement with previous polarimetric observations in the near-infrared and is consistent with dust trapping by a massive planet. No dust emission is detected inside the inner dust cavity ($R \lesssim 20\,$au) or within the dust gap ($\sim 35 - 56\,$au). In contrast, the channel maps of the $J=2-1$ line of the three CO isotopologues reveal the presence of gas inside the dust cavity and dust gap. The gaseous disk is also much larger than the compact dust emission extending to $\sim 1'.5$ ($\sim 180\,$au) in radius. This difference and the sharp drop of the continuum emission at large radii point to radial drift of large dust grains ($>$ \micron-size). Using the thermo-chemical disk code \textsc{dali}, the continuum and the CO isotopologues emission are modelled to quantitatively measure the gas and dust surface densities. The resulting gas surface density is reduced by a factor of $\sim 30-40$ inward of the dust gap. The gas and dust distribution hint at the presence of multiple planets shaping the disk structure via dynamical clearing (dust cavity and gap) and dust trapping (double ring dust distribution).

The Mid-Infrared Instrument (MIRI) Medium Resolution Spectrometer (MRS) is the only mid-IR Integral Field Spectrometer on board James Webb Space Telescope. The complexity of the MRS requires a very specialized pipeline, with some specific steps not present in other pipelines of JWST instruments, such as fringe corrections and wavelength offsets, with different algorithms for point source or extended source data. The MRS pipeline has also two different variants: the baseline pipeline, optimized for most foreseen science cases, and the optimal pipeline, where extra steps will be needed for specific science cases. This paper provides a comprehensive description of the MRS Calibration Pipeline from uncalibrated slope images to final scientific products, with brief descriptions of its algorithms, input and output data, and the accessory data and calibration data products necessary to run the pipeline.

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}$.

Outflowing motions, whether a wind launched from the disk, a jet launched from the protostar, or the entrained molecular outflow, appear to be an ubiquitous feature of star formation. These outwards motions have a number of root causes, and how they manifest is intricately linked to their environment as well as the process of star formation itself. Using the ALMA Science Verification data of HL Tau, we investigate the high velocity molecular gas being removed from the system as a result of the star formation process. We aim to place these motions in context with the optically detected jet, and the disk. With these high resolution ($\sim 1"$) ALMA observations of CO (J=1-0), we quantify the outwards motions of the molecular gas. We find evidence for a bipolar outwards flow, with an opening angle, as measured in the red-shifted lobe, starting off at 90$^\circ$, and narrowing to 60$^\circ$ further from the disk, likely because of magnetic collimation. Its outwards velocity, corrected for inclination angle is of order 2.4 km s$^{-1}$.

We use a simple, self-consistent, self-gravitating semi-analytic disc model to conduct an examination of the parameter space in which self-gravitating discs may exist. We then use Monte-Carlo radiative transfer to generate synthetic ALMA images of these self-gravitating discs to determine the subset of this parameter space in which they generate non-axisymmetric structure that is potentially detectable by ALMA. Recently, several transition discs have been observed to have non-axisymmetric structure that extends out to large radii. It has been suggested that one possible origin of these asymmetries could be spiral density waves induced by disc self-gravity. We use our simple model to see if these discs exist in the region of parameter space where self-gravity could feasibly explain these spiral features. We find that for self-gravity to play a role in these systems typically requires a disc mass around an order of magnitude higher than the observed disc masses for the systems. The spiral amplitudes produced by self-gravity in the local approximation are relatively weak when compared to amplitudes produced by tidal interactions, or spirals launched at Lindblad resonances due to embedded planets in the disc. As such, we ultimately caution against diagnosing spiral features as being due to self-gravity, unless the disc exists in the very narrow region of parameter space where the spiral wave amplitudes are large enough to produce detectable features, but not so large as to cause the disc to fragment.

The gas-solid budget of carbon in protoplanetary disks is related to the composition of the cores and atmospheres of the planets forming in them. The key gas-phase carbon carriers CO, C$^{0}$ and C$^{+}$ can now be observed in disks. The gas-phase carbon abundance in disks has not yet been well characterized, we aim to obtain new constraints on the [C]/[H] ratio in a sample of disks, and to get an overview of the strength of [CI] and warm CO emission. We carried out a survey of the CO$\,6$--$5$ and [CI]$\,1$--$0$ and $2$--$1$ lines towards $37$ disks with APEX, and supplemented it with [CII] data from the literature. The data are interpreted using a grid of models produced with the DALI code. We also investigate how well the gas-phase carbon abundance can be determined in light of parameter uncertainties. The CO$\,6$--$5$ line is detected in $13$ out of $33$ sources, the [CI]$\,1$--$0$ in $6$ out of $12$, and the [CI]$\,2$--$1$ in $1$ out of $33$. With deep integrations, the first unambiguous detections of [CI]~$1$--$0$ in disks are obtained, in TW~Hya and HD~100546. Gas-phase carbon abundance reductions of a factor $5$--$10$ or more can be identified robustly based on CO and [CI] detections. The atomic carbon detection in TW~Hya confirms a factor $100$ reduction of [C]/[H]$_{\rm gas}$ in that disk, while the data are consistent with an ISM-like carbon abundance for HD~100546. In addition, BP~Tau, T~Cha, HD~139614, HD~141569, and HD~100453 are either carbon-depleted or gas-poor disks. The low [CI]~$2$--$1$ detection rates in the survey mostly reflect insufficient sensitivity to detect T~Tauri disks. The Herbig~Ae/Be disks with CO and [CII] upper limits below the models are debris disk like systems. A roughly order of magnitude increase in sensitivity compared to our survey is required to obtain useful constraints on the gas-phase [C]/[H] ratio in most of the targeted systems.

Intermediate mass protostarsprovide a bridge between theories of low- and high-mass star formation. Emerging molecular outflows can be used to determine the influence of fragmentation and multiplicity on protostellar evolution through the correlation of outflow forces of intermediate mass protostars with the luminosity. The aim of this paper is to derive outflow forces from outflows of six intermediate mass protostellar regions and validate the apparent correlation between total luminosity and outflow force seen in earlier work, as well as remove uncertainties caused by different methodology. By comparing CO 6--5 observations obtained with APEX with non-LTE radiative transfer model predictions, optical depths, temperatures, densities of the gas of the molecular outflows are derived. Outflow forces, dynamical timescales and kinetic luminosities are subsequently calculated. Outflow parameters, including the forces, were derived for all sources. Temperatures in excess of 50 K were found for all flows, in line with recent low-mass results. However, comparison with other studies could not corroborate conclusions from earlier work on intermediate mass protostars which hypothesized that fragmentation enhances outflow forces in clustered intermediate mass star formation. Any enhancement in comparison with the classical relation between outflow force and luminosity can be attributed the use of a higher excitation line and improvement in methods; They are in line with results from low-mass protostars using similar techniques. The role of fragmentation on outflows is an important ingredient to understand clustered star formation and the link between low and high-mass star formation. However, detailed information on spatial scales of a few 100 AU, covering all individual members is needed to make the necessary progress.

Context: Jets and outflows are key ingredients in the formation of stars across the mass spectrum. In clustered regions, understanding powering sources and outflow components poses a significant problem. Aims: To understand the dynamics in the outflow(s) from a cluster in the process of forming massive stars. Methods: We use new VLA observations of the molecular gas (SiO, CS, OCS and \molec) in the massive star forming region IRAS 17233-3606 which contains a number of HII regions. We compare these observations to previously published molecular data for this source in order to get a holistic view of the outflow dynamics. Results:We find that the dynamics of the various species can be explained by a single large scale ($\sim 0.15$ pc) outflow when compared to the sizes of the HII regions, with the different morphologies of the blue and red outflow components explained with respect to the morphology of the surrounding envelope. We further find that the direction of the velocity gradients seen in OCS and \molec are suggestive of a combination of rotation and outflow motions in the warm gas surrounding the HII regions near the base of the large scale outflow. Conclusions: Our results show that the massive protostars forming within this region appear to be contributing to a single outflow on large scales. This single large scale outflow is traced by a number of different species as the outflow interacts with its surroundings. On the small scales, there appear to be multiple mechanisms contributing to the dynamics which could be a combination of either a small scale outflow or rotation with the dynamics of the large scale outflow.

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.

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

The formation of high-mass stars is usually accompanied by powerful protostellar outflows. Such high-mass outflows are not simply scaled-up versions of their lower-mass counterparts, since observations suggest that the collimation degree degrades with stellar mass. Theoretically, the origins of massive outflows remain open to question because radiative feedback and fragmentation of the accretion flow around the most massive stars, with M > 15 M_Sun, may impede the driving of magnetic disk winds. We here present a three-dimensional simulation of the early stages of core fragmentation and massive star formation that includes a subgrid-scale model for protostellar outflows. We find that stars that form in a common accretion flow tend to have aligned outflow axes, so that the individual jets of multiple stars can combine to form a collective outflow. We compare our simulation to observations with synthetic H_2 and CO observations and find that the morphology and kinematics of such a collective outflow resembles some observed massive outflows, such as Cepheus A and DR 21. We finally compare physical quantities derived from simulated observations of our models to the actual values in the models to examine the reliability of standard methods for deriving physical quantities, demonstrating that those methods indeed recover the actual values to within a factor of 2-3.

The major physical processes responsible for shaping and sculpting pillars in the clouds surrounding massive stars (i.e. the `Pillars of Creation') are now being robustly incorporated into models quantifying the ionizing radiation from massive stars. The detailed gas dynamics within these pillars can now be compared with observations. Our goal is to quantify the gas dynamics in a pillar being sculpted by a nearby massive star. To do this, we use the CO, $^{13}$CO, and C$^{18}$O J=1-0 emission towards a pillar in the Vulpecula Rift. These data are a combination of CARMA and FCRAO observations providing high resolution ($\sim5''$) imaging of large scale pillar structures ($>100''$). We find that this cold ($\sim18$ K), low density material ($8\times10^3$ cm$^{-3}$) material is fragmenting on Jeans scales, has very low velocity dispersions ($\sim0.5$ km s$^{-1}$), and appears to be moving away from the ionizing source. We are able to draw direct comparisons with three models from the literature, and find that those with lower velocity dispersions best fit our data, although the dynamics of any one model do not completely agree with our observations. We do however, find that our observed pillar exhibits many of the characteristics expected from simulations.

Herbig Ae/Be stars are intermediate-mass pre-main sequence stars surrounded by circumstellar dust disks. Some are observed to produce jets, whose appearance as a sequence of shock fronts (knots) suggests a past episodic outflow variability. This "jet fossil record" can be used to reconstruct the outflow history. We present the first optical to near-infrared (NIR) VLT/X-shooter spectra of the jet from the Herbig Ae star HD 163296. We determine physical conditions in the knots, as well as their kinematic "launch epochs". Knots are formed simultaneously on either side of the disk, with a regular interval of ~16 yr. The velocity dispersion versus jet velocity and the energy input are comparable in both lobes. However, the mass loss rate, velocity, and shock conditions are asymmetric. We find Mjet/Macc ~ 0.01-0.1, consistent with magneto-centrifugal jet launching models. No evidence for dust is found in the high-velocity jet, suggesting it is launched within the sublimation radius (<0.5 au). The jet inclination measured from proper motions and radial velocities confirms it is perpendicular to the disk. A tentative relation is found between the structure of the jet and the photometric variability of the source. Episodes of NIR brightening were previously detected and attributed to a dusty disk wind. We report for the first time significant optical fadings lasting from a few days up to a year, coinciding with the NIR brightenings. These are likely caused by dust lifted high above the disk plane; this supports the disk wind scenario. The disk wind is launched at a larger radius than the high-velocity atomic jet, although their outflow variability may have a common origin. No significant relation between outflow and accretion variability could be established. Our findings confirm that this source undergoes periodic ejection events, which may be coupled with dust ejections above the disk plane.

The formation of stars is usually accompanied by the launching of protostellar outflows. Observations with the Atacama Large Millimetre/sub-millimetre Array (ALMA) will soon revolutionalise our understanding of the morphologies and kinematics of these objects. In this paper, we present synthetic ALMA observations of protostellar outflows based on numerical magnetohydrodynamic collapse simulations. We find significant velocity gradients in our outflow models and a very prominent helical structure within the outflows. We speculate that the disk wind found in the ALMA Science Verification Data of HD 163296 presents a first instance of such an observation.

The high spatial and line sensitivity of ALMA opens the possibility of resolving emission from molecules in circumstellar disks. With an understanding of physical conditions under which molecules have high abundance, they can be used as direct tracers of distinct physical regions. In particular, DCO+ is expected to have an enhanced abundance within a few Kelvin of the CO freezeout temperature of 19 K, making it a useful probe of the cold disk midplane. We compare ALMA line observations of HD 163296 to a grid of models. We vary the upper- and lower-limit temperatures of the region in which DCO+ is present as well as the abundance of DCO+ in order to fit channel maps of the DCO+ J=5-4 line. To determine the abundance enhancement compared to the general interstellar medium, we carry out similar fitting to HCO+ J=4-3 and H13CO+ J=4-3 observations. ALMA images show centrally peaked extended emission from HCO+ and H13CO+. DCO+ emission lies in a resolved ring from ~110 to 160 AU. The outer radius approximately corresponds to the size of the CO snowline as measured by previous lower resolution observations of CO lines in this disk. The ALMA DCO+ data now resolve and image the CO snowline directly. In the best fitting models, HCO+ exists in a region extending from the 19 K isotherm to the photodissociation layer with an abundance of 3x10^-10 relative to H2. DCO+ exists within the 19-21 K region of the disk with an abundance ratio [DCO+] / [HCO+] = 0.3. This represents a factor of 10^4 enhancement of the DCO+ abundance within this narrow region of the HD 163296 disk. Such a high enhancement has only previously been seen in prestellar cores. The inferred abundances provide a lower limit to the ionization fraction in the midplane of the cold outer disk (approximately greater than 4x10^-10), and suggest the utility of DCO+ as a tracer of its parent molecule H2D+. Abridged

Aims: The aim of this work is to study the structure of the protoplanetary disk surrounding the Herbig Ae star HD 163296. Methods: We have used high-resolution and high-sensitivity ALMA observations of the CO(3-2) emission line and the continuum at 850 microns, as well as the 3- dimensional radiative transfer code MCFOST to model the data presented in this work. Results: The CO(3-2) emission unveils for the first time at sub-millimeter frequencies the vertical structure details of a gaseous disk in Keplerian rotation, showing the back- and the front-side of a flared disk. Continuum emission at 850 microns reveals a compact dust disk with a 240 AU outer radius and a surface brightness profile that shows a very steep decline at radius larger than 125 AU. The gaseous disk is more than two times larger than the dust disk, with a similar critical radius but with a shallower radial profile. Radiative transfer models of the continuum data confirms the need for a sharp outer edge to the dust disk. The models for the CO(3-2) channel map require the disk to be slightly more geometrically thick than previous models suggested, and that the temperature at which CO gas becomes depleted (frozen-out) from the outer regions of the disk midplane is T < 20 K, in agreement with previous studies.

Whether high mass stars continue to accrete material beyond the formation of an HII region is still an open question. Ionized infall and outflow have been seen in some sources, but their ties to the surrounding molecular gas are not well constrained. We aim to quantify the ionized and molecular gas dynamics in a high mass star forming region (K3-50A) and their interaction. We present CARMA observations of the 3mm continuum, H41\alpha, and HCO+ emission, and VLA continuum observations at 23 GHz and 14.7 GHz to quantify the gas and its dynamics in K3-50A. We find large scale dynamics consistent with previous observations. On small scales, we find evidence for interaction between the ionized and molecular gas which suggests the ionized outflow is entraining the molecular one. This is the first time such an outflow entrained by photo ionized gas has been observed. Accretion may be ongoing in K3-50A because an ionized bipolar outflow is still being powered, which is in turn entraining part of the surrounding molecular gas. This outflow scenario is similar to that predicted by ionization feedback models.

Disk winds have been postulated as a mechanism for angular momentum release in protostellar systems for decades. HD 163296 is a Herbig Ae star surrounded by a disk and has been shown to host a series of HH knots (HH 409) with bow shocks associated with the farthest knots. Here we present ALMA Science Verification data of CO J=2-1 and J=3-2 emission which are spatially coincident with the blue shifted jet of HH knots, and offset from the disk by -18.6 km/s. The emission has a double corkscrew morphology and extends more than 10" from the disk with embedded emission clumps coincident with jet knots. We interpret this double corkscrew as emission from material in a molecular disk wind, and that the compact emission near the jet knots is being heated by the jet which is moving at much higher velocities. We show that the J=3-2 emission is likely heavily filtered by the interferometer, but the J=2-1 emission suffers less due to the larger beam and measurable angular scales. Excitation analysis suggests temperatures exceeding 900 K in these compact features, with the wind mass, momentum and energy being of order 10^{-5} M_sun, 10^{-4} M_sun km/s and 10^{40} erg respectively. The high mass loss rate suggests that this star is dispersing the disk faster than it is funneling mass onto the star.

Observations of higher-excited transitions of abundant molecules such as CO are important for determining where energy in the form of shocks is fed back into the parental envelope of forming stars. The nearby prototypical and protobinary low-mass hot core, IRAS16293-2422 (I16293) is ideal for such a study. The source was targeted with ALMA for science verification purposes in band 9, which includes CO J=6-5 (E_up/k_B ~ 116 K), at an unprecedented spatial resolution (~0.2", 25 AU). I16293 itself is composed of two sources, A and B, with a projected distance of 5". CO J=6-5 emission is detected throughout the region, particularly in small, arcsecond-sized hotspots, where the outflow interacts with the envelope. The observations only recover a fraction of the emission in the line wings when compared to data from single-dish telescopes, with a higher fraction of emission recovered at higher velocities. The very high angular resolution of these new data reveal that a bow shock from source A coincides, in the plane of the sky, with the position of source B. Source B, on the other hand, does not show current outflow activity. In this region, outflow entrainment takes place over large spatial scales, >~ 100 AU, and in small discrete knots. This unique dataset shows that the combination of a high-temperature tracer (e.g., CO J=6-5) and very high angular resolution observations is crucial for interpreting the structure of the warm inner environment of low-mass protostars.

The formation of massive stars exceeding 10 solar masses usually results in large-scale molecular outflows. Numerical simulations, including ionization, of the formation of such stars show evidence for ionization-driven molecular outflows. We here examine whether the outflows seen in these models reproduce the observations. We compute synthetic ALMA and CARMA maps of CO emission lines of the outflows, and compare their signatures to existing single-dish and interferometric data. We find that the ionization-driven models can only reproduce weak outflows around high-mass star-forming regions. We argue that expanding H II regions probably do not represent the dominant mechanism for driving observed outflows. We suggest instead that observed outflows are driven by the collective action of the outflows from the many lower-mass stars that inevitably form around young massive stars in a cluster.

(Context) Many physical parameters change with time in star forming regions. Here we attempt to correlate changes in infall and outflow motions in high mass star forming regions with evolutionary stage using JCMT observations. (Aims) From a sample of 45 high mass star forming regions in three phases of evolution, we investigate the presence of established infall and outflow tracers to determine whether there are any trends attributable to the age of the source. (Methods) We obtained JCMT observations of HCO+/H13CO+ J=4-3 to trace large scale infall, and SiO J=8-7 to trace recent outflow activity. We compare the infall and outflow detections to the evolutionary stage of the host source (high mass protostellar objects, hypercompact HII regions and ultracompact HII regions). We also note that the integrated intensity of SiO varies with the full width at half maximum of the H13CO+. (Results) We find a surprising lack of SiO detections in the middle stage (Hypercompact HII regions), which may be due to an observational bias. When SiO is detected, we find that the integrated intensity of the line increases with evolutionary stage. We also note that all of the sources with infall signatures onto Ultracompact HII regions have corresponding outflow signatures as well.

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

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

We present maps of seven sources selected from Klaassen & Wilson (2007a) in SiO (J=8-7) and HCO$^+$ and H$^{13}$CO$^+$ (J=4-3) which were obtained using HARP-B on the James Clerk Maxwell Telescope. We find that four out of our seven sources have infall signatures based on their HCO$^+$ emission profiles. From our maps, we have determined the extent of both the outflowing and infalling regions towards these sources, and constrained the amount of infalling and outflowing mass as well as the mass infall rate for each massive star forming region. From our SiO observations, we estimate the source luminosity required to shock the surroundings of these massive star forming regions and find luminosities similar to those of the HII regions themselves. We find that the ratio between our infall and outflow masses is less than one, suggesting high mass entrainment rates in the molecular outflows. We also find that the large scale molecular infall rate towards G10.6-0.4 is comparable to the small scale molecular infall rate found in previous studies.

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

In order to distinguish between the various components of massive star forming regions (i.e. infalling, outflowing and rotating gas structures) within our own Galaxy, we require high angular resolution observations which are sensitive to structures on all size scales. To this end, we present observations of the molecular and ionized gas towards massive star forming regions at 230 GHz from the SMA (with zero spacing from the JCMT) and at 22 and 23 GHz from the VLA at arcsecond or better resolution. These observations (of sources such as NGC7538, W51e2 and K3-50A) form an integral part of a multi-resolution study of the molecular and ionized gas dynamics of massive star forming regions (i.e. Klaassen & Wilson 2007). Through comparison of these observations with 3D radiative transfer models, we hope to be able to distinguish between various modes of massive star formation, such as ionized or halted accretion (i.e Keto 2003 or Klaassen et al. 2006 respectively).

We present single pointing observations of SiO, HCO$^+$ and H$^{13}$CO$^+$ from the James Clerk Maxwell Telescope towards 23 massive star forming regions previously known to contain molecular outflows and ultracompact HII regions. We detected SiO towards 14 sources and suggest that the non-detections in the other nine sources could be due to those outflows being older and without ongoing shocks to replenish the SiO. We serendipitously detected SO$_2$ towards 17 sources in the same tuning as HCO$^+$. We detected HCO$^+$ towards all sources, and suggest that it is tracing infall in nine cases. For seven infall candidates, we estimate mass infall rates between 1$\times10^{-2}$ and 2$\times10^{-5}$ M$_{\odot}$ yr$^{-1}$. Seven sources show both SiO detections (young outflows) and HCO$^+$ infall signatures. We also find that the abundance of H$^{13}$CO$^+$ tends to increase along with the abundance of SiO in sources for which we could determine abundances. We discuss these results with respect to current theories of massive star formation via accretion. From this survey, we suggest that perhaps both models of ionized accretion and halted accretion may be important in describing the evolution of a massive protostar (or protostars) beyond the formation of an HII region.

We have obtained maps of the large scale outflow associated with the UCHII regionG5.89-0.39 in CO and 13CO (J=3-2), SiO (J=8-7,J=5-4), SO2 (13(2,12)-13(1,13)) and H13CO+(J=4-3). From these maps we have been able to determine the mass (3.3 M_sun), momentum (96 M_sun km/s), energy (3.5x10^{46} erg), mechanical luminosity (141 L_sun), and mass loss rate (~1x10^{-3}M_sun/yr) in the large scale outflow. The observationally derived parameters were used toguide 3D magnetohydrodynamic models of the jet entrained outflow. Through the combination of observations and simulations, we suggest that the large scale outflowmay be inclined by approximately 45 deg to the line of sight, and that the jet entraining the observed molecular outflow may have been active for as little as 1000 years, half the kinematic age of the outflow.

We have observed 12CO J = 2-1 and J = 1-0, and 13CO J = 1-0 emission in two regions of HI Self-Absorption (HISA) in Perseus: a small, isolated HISA feature called the globule and a more extended HISA cloud called the complex. Using both Large Velocity Gradient and Monte Carlo radiative transfer codes we found that, in the globule, N(12CO) < 6.0x10^15 cm-2 which, using PDR models, implies that N(H_2) < 9.9x10^20 cm-2. In the complex we found that the H_2 column densities ranged from 1.2 - 2.2 x 10^21 cm-2. By comparing the HISA and CO observations we are able to constrain the physical conditions and atomic gas fraction (f). In the globule, 8 K < T_spin < 22 K and 0.02 < f < 0.2 depending on whether the (unknown) gas density is 10^2, 10^3, or 10^4 cm-3. In the complex, 12 K < T_spin < 24 K, 0.02 < f < 0.05, and we were also able to constrain the gas density (100 < n < 1200 cm-3). These results imply that the gas in the HISA clouds is colder and denser than that usually associated with the atomic ISM and, indeed, is similar to that seen in molecular clouds. The small atomic gas fractions also imply that there is a significant molecular component in these HISA clouds, even when little or no 12CO is detected. The level of 12CO detected and the visual extinction due to dust is consistent with the idea that these HISA clouds are undergoing a transition from the atomic to molecular phase.