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.

Understanding planet formation requires to discern how dust grows in protoplanetary disks. An important parameter to measure in disks is the maximum dust grain size present. This is usually estimated through measurements of the dust opacity at different millimeter wavelengths assuming optically thin emission and dust opacity dominated by absorption. However, ALMA observations have shown that these assumptions might not be correct in the case of protoplanetary disks, leading to overestimation of particle sizes and to underestimation of the disk's mass. Here, we present an analysis of high quality ALMA and VLA images of the HL Tau protoplanetary disk, covering a wide range of wavelengths, from 0.8 mm to 1 cm, and with a physical resolution of $\sim$7.35 au. We describe a procedure to analyze a set of millimeter images without any assumption about the optical depth of the emission, and including the effects of absorption and scattering in the dust opacity. This procedure allows us to obtain the dust temperature, the dust surface density and the maximum particle size at each radius. In the HL Tau disk, we found that particles have already grown up to a few millimeters in size. We detect differences in the dust properties between dark and bright rings, with dark rings containing low dust density and small dust particles. Different features in the HL Tau disk seem to have different origins. Planet-disk interactions can explain substructure at the external half of the disk, but the internal rings seem to be associated to the presence of snow lines of several molecules.

We report new Atacama Large Millimeter/submillimeter Array Band 3 (86-100 GHz; $\sim$80 mas angular resolution) and Band 4 (146-160 GHz; $\sim$50 mas angular resolution) observations of the dust continuum emission towards the archetypal and ongoing accretion burst young stellar object FU Ori, which simultaneously covered its companion, FU Ori S. In addition, we present near-infrared (2-2.45 $\mu$m) observations of FU Ori taken with the General Relativity Analysis via VLT InTerferometrY (GRAVITY; $\sim$1 mas angular resolution) instrument on the Very Large Telescope Interferometer (VLTI). We find that the emission in both FU Ori and FU Ori S at (sub)millimeter and near infrared bands is dominated by structures inward of $\sim$10 au radii. We detected closure phases close to zero from FU Ori with VLTI/GRAVITY, which indicate the source is approximately centrally symmetric and therefore is likely viewed nearly face-on. Our simple model to fit the GRAVITY data shows that the inner 0.4 au radii of the FU Ori disk has a triangular spectral shape at 2-2.45 $\mu$m, which is consistent with the H$_{2}$O and CO absorption features in a $\dot{M}\sim$10$^{-4}$ $M_{\odot}\,yr^{-1}$, viscously heated accretion disk. At larger ($\sim$0.4-10 au) radii, our analysis shows that viscous heating may also explain the observed (sub)millimeter and centimeter spectral energy distribution when we assume a constant, $\sim$10$^{-4}$ $M_{\odot}\,yr^{-1}$ mass inflow rate in this region. This explains how the inner 0.4 au disk is replenished with mass at a modest rate, such that it neither depletes nor accumulates significant masses over its short dynamic timescale. Finally, we tentatively detect evidence of vertical dust settling in the inner 10 au of the FU Ori disk, but confirmation requires more complete spectral sampling in the centimeter bands.

In previous contributions, we have presented an analytical model describing the evolution of molecular clouds (MCs) undergoing hierarchical gravitational contraction. The cloud's evolution is characterized by an initial increase in its mass, density, and star formation rate (SFR) and efficiency (SFE) as it contracts, followed by a decrease of these quantities as newly formed massive stars begin to disrupt the cloud. The main parameter of the model is the maximum mass reached by the cloud during its evolution. Thus, specifying the instantaneous mass and some other variable completely determines the cloud's evolutionary stage. We apply the model to interpret the observed scatter in SFEs of the cloud sample compiled by Lada et al.\ as an evolutionary effect so that, although clouds such as California and Orion A have similar masses, they are in very different evolutionary stages, causing their very different observed SFRs and SFEs. The model predicts that the California cloud will eventually reach a significantly larger total mass than the Orion A cloud. Next, we apply the model to derive estimated ages of the clouds since the time when approximately 25\% of their mass had become molecular. We find ages from $\sim 1.5$ to 27 Myr, with the most inactive clouds being the youngest. Further predictions of the model are that clouds with very low SFEs should have massive atomic envelopes constituting the majority of their gravitational mass, and that low-mass clouds ($M \sim 10^3$-$10^4 \, M_\odot$) end their lives with a mini-burst of star formation, reaching SFRs $\sim 300$-$500\, M_\odot$ Myr$^{-1}$. By this time, they have contracted to become compact ($\sim 1$ pc) massive star-forming clumps, in general embedded within larger GMCs.

In recent years, the disk populations in a number of young star-forming regions have been surveyed with ALMA. Understanding the disk properties and their correlation with those of the central star is critical to understand planet formation. In particular, a decrease of the average measured disk dust mass with the age of the region has been observed. We conducted high-sensitivity continuum ALMA observations of 43 Class II young stellar objects in CrA at 1.3 mm (230 GHz). The typical spatial resolution is 0.3". The continuum fluxes are used to estimate the dust masses of the disks, and a survival analysis is performed to estimate the average dust mass. We also obtained new VLT/X-Shooter spectra for 12 of the objects in our sample. 24 disks are detected, and stringent limits have been put on the average dust mass of the non-detections. Accounting for the upper limits, the average disk mass in CrA is $6\pm3\,\rm M_\oplus$, significantly lower than that of disks in other young (1-3 Myr) star forming regions (e.g. Lupus) and appears consistent with the 5-10 Myr old Upper Sco. The position of the stars in our sample on the HR diagram, however, seems to confirm that that CrA has age similar to Lupus. Neither external photoevaporation nor a lower than usual stellar mass distribution can explain the low disk masses. On the other hand, a low-mass disk population could be explained if the disks are small, which could happen if the parent cloud has a low temperature or intrinsic angular momentum, or if the the angular momentum of the cloud is removed by some physical mechanism such as magnetic braking. In order to fully explain and understand the dust mass distribution of protoplanetary disks and their evolution, it may also be necessary to take into consideration the initial conditions of star and disk formation process, which may vary from region to region, and affect planet formation.

We used ALMA to observe the star-forming region GGD27 at 1.14 mm with an unprecedented angular resolution, 40 mas (56 au) and sensitivity (0.002 Msun). We detected a cluster of 25 continuum sources, most of which are likely tracing disks around Class 0/I protostars. Excluding the two most massive objects, disks masses are in the range 0.003-0.05 Msun. The analysis of the cluster properties indicates that GGD27 displays moderate subclustering. This result combined with the dynamical timescale of the radio jet (10000 years) suggests the youthfulness of the cluster. The lack of disk mass segregation signatures may support this too. We found a clear paucity of disks with Rdisk >100 au. The median value of the radius is 34 au, smaller than the median of 92 au for Taurus but comparable to the value found in Ophiuchus and in the Orion Nebula Cluster. In GGD27 there is no evidence of a distance-dependent disk mass distribution (i. e., disk mass depletion due to external photoevaporation), most likely due to the cluster youth. There is a clear deficit of disks for distances <0.02 pc. Only for distances >0.04 pc stars can form larger and more massive disks, suggesting that dynamical interactions far from the cluster center are weaker, although the small disks found could be the result of disk truncation. This work demonstrates the potential to characterize disks from low-mass YSOs in distant and massive (still deeply embedded) clustered environments.

Context. Colliding winds in massive binaries are able to accelerate particles up to relativistic speeds as the result of the interaction between the winds of the different stellar components. HD 167971 exhibits this phenomenology which makes it a strong radio source. Aims. We aim at characterizing the morphology of the radio emission and its dependence on the orbital motion, traced independently by NIR-interferometry, of the spectroscopic binary and the tertiary component that conforms HD 167971. Methods. We analyze 2006 and 2016 very long baseline interferometric data at C and X bands. We complement our analysis with a geometrical model of the wind-wind collision region, and with an astrometric description of the system. Results. We confirm that the detected non-thermal radio emission is associated with the wind-wind collision region of the spectroscopic binary and the tertiary component in HD 167971. The wind-wind collision region changes orientation in agreement with the orbital motion of the tertiary around the spectroscopic binary. The total intensity also changes between the two observing epochs in a way inversely proportional to the separation between the SB and T, with a negative-steep spectral index typical of an optically thin synchrotron emission possibly steepened by an inverse Compton cooling effect. The wind-wind collision bow-shock shape and its position with respect to the stars indicates that the wind momentum from the spectroscopic binary is stronger than that of the tertiary. Finally, the astrometric solution derived for the stellar system and the wind-wind collision region is consistent with independent Gaia data.

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.

The filamentary infrared dark cloud L1287 is actively forming a dense cluster of low-mass YSOs at its inner $\sim$0.1 pc region. To help understand the origin of this low-mass YSO cluster, the present work aims at resolving the gas structures and kinematics. We have performed $\sim$1$"$ angular resolution ($\sim$930 AU) SMA observations at $\sim$1.3 mm. From a $\sim$2$"$ resolution 1.3 mm continuum image we identified six dense cores, namely SMA1-6 with masses in the range of $\sim0.4-4$ M$_\odot$. From a $\sim$1$"$ resolution 1.3 mm continuum image, we find a high fragmentation level, with 14 compact millimeter sources within 0.1 pc (two of them associated with the known accretion outburst YSOs RNO 1C and RNO 1B). The dense gas tracer DCN (3--2) traces well the dust continuum emission and shows a clear velocity gradient along the NW-SE direction centered at SMA3. There is another velocity gradient with opposite direction around the most luminous YSO IRAS 00338+6312. The fragmentation within 0.1 pc in L1287 is very high compared to other regions at the same spatial scales. The incoherent motions of dense gas flows are sometimes interpreted by being influenced by (proto)stellar feedback (e.g., outflows), which is not yet ruled out in this particular target source. The directions of the velocity gradients traced by DCN are approximately perpendicular to those of the dominant CO outflow(s). Therefore, we alternatively hypothesize that the velocity gradients revealed by DCN trace the convergence from the $\gtrsim$0.1 pc scales infalling motion towards the rotational motions around the more compact ($\sim0.02$ pc) sources. This global molecular gas converging flow may feed the formation of the dense low-mass YSO cluster. IRAS 00338+6312 is the most likely powering source of the dominant CO outflow. A compact blue-shifted outflow from RNO 1C is also identified.

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.

In this paper we explore the effects of self-obscuration in protostellar disks with a radially decreasing temperature gradient and a colder midplane. We are motivated by recent reports of resolved dark lanes (`hamburgers') and (sub)mm spectral indices systematically below the ISM value for optically thin dust $\alpha_{\rm ISM} =3.7$. We explore several model grids, scaling disk mass and varying inclination angle $i$ and observing frequency $\nu$ from the VLA Ka band ($\sim 37$ GHz) to ALMA Band 8 ($\sim 405$ GHz). We also consider the effects of decreasing the index of the (sub-)mm dust opacity power law $\beta$ from 1.7 to 1. We find that a distribution of disk masses in the range $M_{\rm disk} = 0.01-2~M_\odot$ is needed to reproduce the observed distribution of spectral indices, and that assuming a fixed $\beta =1.7$ gives better results than $\beta=1$. A wide distribution of disk masses is also needed to produce some cases with $\alpha <2$, as reported for some sources in the literature. Such extremely low spectral indices arise naturally when the selected observing frequencies sample the appropriate change in the temperature structure of the optically thick model disk. Our results show that protostellar disk masses could often be underestimated by $> \times10$, and are consistent with recent hydrodynamical simulations. Although we do not rule out the possibility of some grain growth occurring within the short protostellar timescales, we conclude that self-obscuration needs to be taken into account.

A number of ultracompact H II regions in Galactic star forming environments have been observed to vary significantly in radio flux density on timescales of 10-20 years. Theory predicted that such variations should occur when the accretion flow that feeds a young massive star becomes unstable and clumpy. We have targeted the massive star-forming region W49A with the Karl G. Jansky Very Large Array (VLA) for observations at 3.6 cm with the B-configuration at 0.8'' resolution, to compare to nearly identical observations taken almost 21 years earlier (February 2015 and August 1994). Most of the sources in the crowded field of ultracompact and hypercompact H II regions exhibit no significant changes over this time period. However, one source, W49A/G2, decreased by 20% in peak intensity (from 71+/-4 mJy/beam to 57+/-3 mJy/beam), and 40% in integrated flux (from 0.109+/-0.011 Jy to 0.067+/-0.007 Jy), where we cite 5 sigma errors in peak intensity, and 10% errors in integrated flux. We present the radio images of the W49A region at the two epochs, the difference image that indicates the location of the flux density decrease, and discuss explanations for the flux density decrease near the position of W49A/G2.

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.

As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41+0.31 with ALMA at 1.4 mm with an angular resolution of~0.22" (~1700 au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is associated with two free-free continuum sources. The Main core looks featureless and homogeneous in dust continuum emission and does not present any hint of fragmentation. Each transition of CH3CN and CH3OCHO, both ground and vibrationally excited, as well as those of CH3CN isotopologues, shows a clear velocity gradient along the NE-SW direction, with velocity linearly increasing with distance from the center, consistent with solid-body rotation. However, when comparing the velocity field of transitions with different upper level energies, the rotation velocity increases with increasing energy of the transition, which suggests that the rotation speeds up towards the center. Spectral lines towards the dust continuum peak show an inverse P-Cygni profile that supports the existence of infall in the core. The infall velocity increases with the energy of the transition suggesting that the infall is accelerating towards the center of the core, consistent with gravitational collapse. Despite the monolithic appearance of the Main core, the presence of red-shifted absorption, the existence of two embedded free-free sources at the center, and the rotational spin-up are consistent with an unstable core undergoing fragmentation with infall and differential rotation due to conservation of angular momentum. Therefore, the most likely explanation for the monolithic morphology is that the large opacity of the dust emission prevents the detection of any inhomogeneity in the core.

The magnetic field plays an important role in every stage of the star-formation process from the collapse of the initial protostellar core to the star's arrival on the main sequence. Consequently, the goal of this science case is to explore a wide range of magnetic phenomena that can be investigated using the polarization capabilities of the Next Generation Very Large Array (ngVLA). These include (1) magnetic fields in protostellar cores via polarized emission from aligned dust grains, including in regions optically thick at wavelengths observable by the Atacama Large Millimeter/submillimeter Array (ALMA); (2) magnetic fields in both protostellar cores and molecular outflows via spectral-line polarization from the Zeeman and Goldreich-Kylafis effects; (3) magnetic fields in protostellar jets via polarized synchrotron emission; and (4) gyrosynchrotron emission from magnetospheres around low-mass stars.

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

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.

Young stellar objects (YSOs) may undergo periods of active accretion (outbursts), during which the protostellar accretion rate is temporarily enhanced by a few orders of magnitude. Whether or not these accretion outburst YSOs possess similar dust/gas reservoirs to each other, and whether or not their dust/gas reservoirs are similar as quiescent YSOs, are issues not yet clarified. The aim of this work is to characterize the millimeter thermal dust emission properties of a statistically significant sample of long and short duration accretion outburst YSOs (i.e., FUors and EXors) and the spectroscopically identified candidates of accretion outbursting YSOs (i.e., FUor-like objects). We have carried out extensive Submillimeter Array (SMA) observations mostly at $\sim$225 GHz (1.33 mm) and $\sim$272 GHz (1.10 mm), from 2008 to 2017. We covered accretion outburst YSOs located at $<$1 kpc distances from the solar system. We analyze all the existing SMA data of such objects, both published and unpublished, in a coherent way to present a millimeter interferometric database of 29 objects. We obtained 21 detections at $>$3-$\sigma$ significance. Detected sources except for the two cases of V883 Ori and NGC 2071 MM3 were observed with $\sim$1$"$ angular resolution. Overall our observed targets show a systematically higher millimeter luminosity distribution than those of the $M_{*}>$0.3 $M_{\odot}$ Class II YSOs in the nearby ($\lesssim$400 pc) low-mass star-forming molecular clouds (e.g., Taurus, Lupus, Upp Scorpio, and Chameleon I). In addition, at 1 mm our observed confirmed binaries or triple-system sources are systematically fainter than the rest of the sources even though their 1 mm fluxes are broadly distributed. We may have detected $\sim$30-60\% millimeter flux variability from V2494 Cyg and V2495 Cyg, from the observations separated by $\sim$1 year.

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

We resolved FU Ori at 29-37 GHz using the JVLA with $\sim$0$''$.07 resolution, and performed the complementary JVLA 8-10 GHz observations, the SMA 224 GHz and 272 GHz observations, and compared with archival ALMA 346 GHz observations to obtain the SEDs. Our 8-10 GHz observations do not find evidence for the presence of thermal radio jets, and constrain the radio jet/wind flux to at least 90 times lower than the expected value from the previously reported bolometric luminosity-radio luminosity correlation. The emission at $>$29 GHz may be dominated by the two spatially unresolved sources, which are located immediately around FU Ori and its companion FU Ori S, respectively. Their deconvolved radii at 33 GHz are only a few au. The 8-346 GHz SEDs of FU Ori and FU Ori S cannot be fit by constant spectral indices (over frequency). The more sophisticated models for SEDs suggest that the $>$29 GHz emission is contributed by a combination of free-free emission from ionized gas, and thermal emission from optically thick and optically thin dust components. We hypothesize that dust in the innermost parts of the disks ($\lesssim$0.1 au) has been sublimated, and thus the disks are no more well shielded against the ionizing photons. The estimated overall gas and dust mass based on SED modeling, can be as high as a fraction of a solar mass, which is adequate for developing disk gravitational instability. Our present explanation for the observational data is that the massive inflow of gas and dust due to disk gravitational instability or interaction with a companion/intruder, was piled up at the few au scale due to the development of a deadzone with negligible ionization. The piled up material subsequently triggered the thermal and the MRI instabilities when the ionization fraction in the inner sub-au scale region exceeded a threshold value, leading to the high protostellar accretion rate.

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

Supermassive black holes in galaxy centres can grow by the accretion of gas, liberating energy that might regulate star formation on galaxy-wide scales. The nature of the gaseous fuel reservoirs that power black hole growth is nevertheless largely unconstrained by observations, and is instead routinely simplified as a smooth, spherical inflow of very hot gas. Recent theory and simulations instead predict that accretion can be dominated by a stochastic, clumpy distribution of very cold molecular clouds - a departure from the "hot mode" accretion model - although unambiguous observational support for this prediction remains elusive. Here we report observations that reveal a cold, clumpy accretion flow towards a supermassive black hole fuel reservoir in the nucleus of the Abell 2597 Brightest Cluster Galaxy (BCG), a nearby (redshift z=0.0821) giant elliptical galaxy surrounded by a dense halo of hot plasma. Under the right conditions, thermal instabilities can precipitate from this hot gas, producing a rain of cold clouds that fall toward the galaxy's centre, sustaining star formation amid a kiloparsec-scale molecular nebula that inhabits its core. The observations show that these cold clouds also fuel black hole accretion, revealing "shadows" cast by the molecular clouds as they move inward at about 300 kilometres per second towards the active supermassive black hole in the galaxy centre, which serves as a bright backlight. Corroborating evidence from prior observations of warmer atomic gas at extremely high spatial resolution, along with simple arguments based on geometry and probability, indicate that these clouds are within the innermost hundred parsecs of the black hole, and falling closer towards it.

Young massive clusters (YMCs) have central stellar mass surface densities exceeding $10^{4} M_{\odot} pc^{-2}$. It is currently unknown whether the stars formed at such high (proto)stellar densities. We compile a sample of gas clouds in the Galaxy which have sufficient gas mass within a radius of a few parsecs to form a YMC, and compare their radial gas mass distributions to the stellar mass distribution of Galactic YMCs. We find that the gas in the progenitor clouds is distributed differently than the stars in YMCs. The mass surface density profiles of the gas clouds are generally shallower than the stellar mass surface density profiles of the YMCs, which are characterised by prominent dense core regions with radii ~ 0.1 pc, followed by a power-law tail. On the scale of YMC core radii, we find that there are no known clouds with significantly more mass in their central regions when compared to Galactic YMCs. Additionally, we find that models in which stars form from very dense initial conditions require surface densities that are generally higher than those seen in the known candidate YMC progenitor clouds. Our results show that the quiescent, less evolved clouds contain less mass in their central regions than in the highly star-forming clouds. This suggests an evolutionary trend in which clouds continue to accumulate mass towards their centres after the onset of star formation. We conclude that a conveyor-belt scenario for YMC formation is consistent with the current sample of Galactic YMCs and their progenitor clouds.

We report Submillimeter Array (SMA) 1.3 mm high angular resolution observations towards the four EXor type outbursting young stellar objects (YSOs) VY Tau, V1118 Ori, V1143 Ori, and NY Ori. The data mostly show low dust masses $M_{dust}$ in the associated circumstellar disks. Among the sources, NY Ori possesses a relatively massive disk with $M_{dust} \sim 9 \times 10^{-4}$ $M_{\odot}$. V1118 Ori has a marginal detection equivalent to $M_{dust} \sim 6 \times 10^{-5}$ $M_{\odot}$. V1143 Ori has a non-detection also equivalent to $M_{dust} < 6 \times 10^{-5}$ $M_{\odot}$. For the nearest source VY Tau, we get a surprising non-detection which provides a stringent upper limit $M_{dust} < 6 \times 10^{-6}$ $M_{\odot}$. We interpret our findings as suggesting that the gas and dust reservoirs that feed the short duration, repetitive optical outbursts seen in some EXors may be limited to the small scale, innermost region of their circumstellar disks. This hot dust may have escaped our detection limits. Follow-up, more sensitive millimeter observations are needed to improve our understanding of the triggering mechanisms of EXor type outbursts.

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

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.

There is increasing evidence that episodic accretion is a common phenomenon in Young Stellar Objects (YSOs). Recently, the source HOPS 383 in Orion was reported to have a $\times 35$ mid-infrared -- and bolometric -- luminosity increase between 2004 and 2008, constituting the first clear example of a class 0 YSO (a protostar) with a large accretion burst. The usual assumption that in YSOs accretion and ejection follow each other in time needs to be tested. Radio jets at centimeter wavelengths are often the only way of tracing the jets from embedded protostars. We searched the Very Large Array archive for the available observations of the radio counterpart of HOPS 383. The data show that the radio flux of HOPS 383 varies only mildly from January 1998 to December 2014, staying at the level of $\sim 200$ to 300 $\mu$Jy in the X band ($\sim 9$ GHz), with a typical uncertainty of 10 to 20 $\mu$Jy in each measurement. We interpret the absence of a radio burst as suggesting that accretion and ejection enhancements do not follow each other in time, at least not within timescales shorter than a few years. Time monitoring of more objects and specific predictions from simulations are needed to clarify the details of the connection betwen accretion and jets/winds in YSOs.

How rapidly collapsing parsec-scale massive molecular clumps feed high-mass stars, and how they fragment to form OB clusters, have been outstanding questions in the field of star-formation. In this work, we report the resolved structures and kinematics of the approximately face-on, rotating massive molecular clump, G33.92+0.11. Our high resolution Atacama Large Millimeter/submillimeter Array (ALMA) images show that the spiral arm-like gas overdensities form in the eccentric gas accretion streams. First, we resolved that the dominant part of the $\sim$0.6 pc scale massive molecular clump (3.0$^{+2.8}_{-1.4}$$\cdot$10$^{3}$ $M_{\odot}$) G33.92+0.11 A is tangled with several 0.5-1 pc size molecular arms spiraling around it, which may be connected further to exterior gas accretion streams. Within G33.92+0.11 A, we resolved the $\sim$0.1 pc width gas mini-arms connecting with the two central massive (100-300 $M_{\odot}$) molecular cores. The kinematics of arms and cores elucidate a coherent accretion flow continuing from large to small scales. We demonstrate that the large molecular arms are indeed the cradles of dense cores, which are likely current or future sites of high-mass star formation. Since these deeply embedded massive molecular clumps preferentially form the highest mass stars in the clusters, we argue that dense cores fed by or formed within molecular arms play a key role in making the upper end of the stellar and core mass functions.

Photoevaporation due to high-energy stellar photons is thought to be one of the main drivers of protoplanetary disk dispersal. The fully or partially ionized disk surface is expected to produce free-free continuum emission at centimeter (cm) wavelengths that can be routinely detected with interferometers such as the upgraded Very Large Array (VLA). We use deep (rms noise down to 8 $\mu$Jy beam$^{-1}$ in the field of view center) 3.5 cm maps of the nearby (130 pc) Corona Australis (CrA) star formation (SF) region to constrain disk photoevaporation models. We find that the radio emission from disk sources in CrA is surprisingly faint. Only 3 out of 10 sources within the field of view are detected, with flux densities of order $10^2$ $\mu$Jy. However, a significant fraction of their emission is non-thermal. Typical upper limits for non-detections are $3\sigma\sim 60~\mu$Jy beam$^{-1}$. Assuming analytic expressions for the free-free emission from extreme-UV (EUV) irradiation, we derive stringent upper limits to the ionizing photon luminosity impinging on the disk surface $\Phi_\mathrm{EUV}<1-4\times10^{41}$ s$^{-1}$. These limits constrain $\Phi_\mathrm{EUV}$ to the low end of the values needed by EUV-driven photoevaporation models to clear protoplanetary disks in the observed few Myr timescale. Therefore, at least in CrA, EUV-driven photoevaporation is not likely to be the main agent of disk dispersal. We also compare the observed X-ray luminosities $L_X$ of disk sources with models in which photoevaporation is driven by such photons. Although predictions are less specific than for the EUV case, most of the observed fluxes (upper limits) are roughly consistent with the (scaled) predictions. Deeper observations, as well as predictions spanning a wider parameter space, are needed to properly test X-ray driven photoevaporation.

The object W33 is a giant molecular cloud that contains star forming regions at various evolutionary stages from quiescent clumps to developed H II regions. Since its star forming regions are located at the same distance and the primary material of the birth clouds is probably similar, we conducted a comparative chemical study to trace the chemical footprint of the different phases of evolution. We observed six clumps in W33 with the Atacama Pathfinder Experiment (APEX) telescope at 280 GHz and the Submillimeter Array (SMA) at 230 GHz. We detected 27 transitions of 10 different molecules in the APEX data and 52 transitions of 16 different molecules in the SMA data. The chemistry on scales larger than $\sim$0.2 pc, which are traced by the APEX data, becomes more complex and diverse the more evolved the star forming region is. On smaller scales traced by the SMA data, the chemical complexity and diversity increase up to the hot core stage. In the H II region phase, the SMA spectra resemble the spectra of the protostellar phase. Either these more complex molecules are destroyed or their emission is not compact enough to be detected with the SMA. Synthetic spectra modelling of the H$_{2}$CO transitions, as detected with the APEX telescope, shows that both a warm and a cold component are needed to obtain a good fit to the emission for all sources except for W33 Main1. The temperatures and column densities of the two components increase during the evolution of the star forming regions. The integrated intensity ratios N$_{2}$H$^{+}$(3$-$2)/CS(6$-$5) and N$_{2}$H$^{+}$(3$-$2)/H$_{2}$CO(4$_{2,2}$$-$3$_{2,1}$) show clear trends as a function of evolutionary stage, luminosity, luminosity-to-mass ratio, and H$_{2}$ peak column density of the clumps and might be usable as chemical clocks.

We summarize our comprehensive gas surveys of some of the most luminous, deeply embedded (optically obscured) star formation regions in the Milky Way, which are the local cases of massive star clusters and/or associations in the making. Our approach emphasizes multi-scale, multi-resolution imaging in dust and free-free continuum, as well as in molecular- and hydrogen recombination lines, to trace the multiple gas components from 0.1 pc (core scale) all the way up to the scales of the entire giant molecular cloud (GMC), or $\sim 100$ pc. We highlight our results in W49A, the most luminous Galactic star formation region ($L\sim 10^7~L_\odot$), which appears to be forming a young massive cluster (or a binary star cluster) with $M_\star > 5\times10^4~M_\odot$ that may remain bound after gas dispersal. The surveyed sources share elements in common, in particular, the 10-100 pc scale GMCs are filamentary but have one or two central condensations (clumps) far denser than the surrounding filaments and that host the (forming) massive stars.

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

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

We report JVLA 8-10 GHz ($\lambda$=3.0-3.7 cm) monitoring observations toward the YSO cluster R Coronae Australis (R\,CrA), taken in 2012, from March 15 to September 12. These observations were planned to measure the radio flux variabilities in timescales from 0.5 hours to several days, to tens of days, and up to $\sim$200 days. We found that among the YSOs detectable in individual epochs, in general, the most reddened objects in the \textit{Spitzer} observations show the highest mean 3.5 cm Stokes \textit{I} emission, and the lowest fractional variabilities on $<$200-day timescales. The brightest radio flux emitters in our observations are the two reddest sources IRS7W and IRS7E. In addition, by comparing with observations taken in 1996-1998 and 2005, we found that the radio fluxes of these two sources have increased by a factor $\sim$1.5. The mean 3.5-cm fluxes of the three Class I/II sources IRSI, IRS2, and IRS6 appear to be correlated with their accretion rates derived by a previous near infrared line survey. The weakly accreting Class I/II YSOs, or those in later evolutionary stages, present radio flux variability on $<$0.5-hour timescales. Some YSOs were detected only during occasional flaring events. The source R\,CrA went below our detection limit during a few fading events.

We present sensitive high angular resolution ($\sim$ 1$"$) millimeter continuum and line observations from the massive star forming region DR21(OH) located in the Cygnus X molecular cloud. Within the well-known dusty MM1-2 molecular cores, we report the detection of a new cluster of about ten compact continuum millimeter sources with masses between 5 and 24 M$_\odot$, and sizes of a few thousands of astronomical units. These objects are likely to be large dusty envelopes surrounding massive protostars, some of them most probably driving several of the outflows that emanate from this region. Additionally, we report the detection of strong millimeter emission of formaldehyde (H$_2$CO) and methanol (CH$_3$OH) near 218 GHz as well as compact emission from the typical outflow tracers carbon monoxide and silicon monoxide (CO and SiO) toward this massive star-forming region. The H$_2$CO and CH$_3$OH emission is luminous ($\sim$ 10$^{-4}$ L$_{\odot}$), well resolved, and found along the collimated methanol maser outflow first identified at centimeter wavelengths and in the sources SMA6 and SMA7. Our observations suggest that this maser outflow might be energized by a millimeter source called SMA4 located in the MM2 dusty core. The CO and SiO emission traces some other collimated outflows that emanate from MM1-2 cores, and are not related with the low velocity maser outflow.

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

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.

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

We present high angular resolution (FWHM_beam < 0.2") observations of the 23.1-GHz methanol (CH_3OH) transition toward the massive-star forming region NGC 7538 IRS 1. The two velocity components previously reported by Wilson et al. are resolved into distinct spatial features with brightness temperatures (T_B) greater than 10^4 K, proving their maser nature. Thus, NGC 7538 IRS 1 is the third region confirmed to show methanol maser emission at this frequency. The brighter 23.1-GHz spot coincides in position with a rare formaldehyde (H_2CO) maser, and marginally with a 22.2-GHz water (H_2O) maser, for which we report archival observations. The weaker CH_3OH spot coincides with an H_2O maser. The ratio of T_B for the 23.1-GHz masers to that of the well-known 12.2-GHz CH_3OH masers in this region roughly agrees with model predictions. However, the 23.1-GHz spots are offset in position from the CH_3OH masers at other frequencies. This is difficult to interpret in terms of models that assume that all the masers arise from the same clumps, but it may result from turbulent conditions within the gas or rapid variations in the background radiation field.

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

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

We present measurements of the mean dense core lifetimes in numerical simulations of magnetically supercritical, turbulent, isothermal molecular clouds, in order to compare with observational determinations. "Prestellar" lifetimes (given as a function of the mean density within the cores, which in turn is determined by the density threshold n_thr used to define them) are consistent with observationally reported values, ranging from a few to several free-fall times. We also present estimates of the fraction of cores in the "prestellar", "stellar'', and "failed" (those cores that redisperse back into the environment) stages as a function of n_thr. The number ratios are measured indirectly in the simulations due to their resolution limitations. Our approach contains one free parameter, the lifetime of a protostellar object t_yso (Class 0 + Class I stages), which is outside the realm of the simulations. Assuming a value t_yso = 0.46 Myr, we obtain number ratios of starless to stellar cores ranging from 4-5 at n_thr = 1.5 x 10^4 cm^-3 to 1 at n_thr = 1.2 x 10^5 cm^-3, again in good agreement with observational determinations. We also find that the mass in the failed cores is comparable to that in stellar cores at n_thr = 1.5 x 10^4 cm^-3, but becomes negligible at n_thr = 1.2 x 10^5 cm^-3, in agreement with recent observational suggestions that at the latter densities the cores are in general gravitationally dominated. We conclude by noting that the timescale for core contraction and collapse is virtually the same in the subcritical, ambipolar diffusion-mediated model of star formation, in the model of star formation in turbulent supercritical clouds, and in a model intermediate between the previous two, for currently accepted values of the clouds' magnetic criticality.