After almost 20 years of hunting, only about a dozen hot corinos, hot regions enriched in interstellar complex organic molecules (iCOMs), are known. Of them, many are binary systems with the two components showing drastically different molecular spectra. Two obvious questions arise. Why are hot corinos so difficult to find and why do their binary components seem chemically different? The answer to both questions could be a high dust opacity that would hide the molecular lines. To test this hypothesis, we observed methanol lines at centimeter wavelengths, where dust opacity is negligible, using the Very Large Array interferometer. We targeted the NGC 1333 IRAS 4A binary system, for which one of the two components, 4A1, has a spectrum deprived of iCOMs lines when observed at millimeter wavelengths, while the other component, 4A2, is very rich in iCOMs. We found that centimeter methanol lines are similarly bright toward 4A1 and 4A2. Their non-LTE analysis indicates gas density and temperature ($\geq2\times10^6$ cm$^{-3}$ and 100--190 K), methanol column density ($\sim10^{19}$ cm$^{-2}$) and extent ($\sim$35 au in radius) similar in 4A1 and 4A2, proving that both are hot corinos. Furthermore, the comparison with previous methanol line millimeter observations allows us to estimate the optical depth of the dust in front of 4A1 and 4A2, respectively. The obtained values explain the absence of iCOMs line emission toward 4A1 at millimeter wavelengths and indicate that the abundances toward 4A2 are underestimated by $\sim$30\%. Therefore, centimeter observations are crucial for the correct study of hot corinos, their census, and their molecular abundances.

We analyze multi-frequency, high-resolution continuum data obtained by ALMA and JVLA to study detailed structure of the dust distribution in the infant disk of a Class~0/I source, L1527 IRS. We find three clumps aligning in the north-south direction in the $7 {\rm \,mm}$ radio continuum image. The three clumps remain even after subtracting free-free contamination, which is estimated from the $1.3{\rm \,cm}$ continuum observations. The northern and southern clumps are located at a distance of $\sim 15{\rm \,au}$ from the central clump and are likely optically thick at $7{\rm \,mm}$ wavelength. The clumps have similar integrated intensities. The symmetric physical properties could be realized when a dust ring or spiral arms around the central protostar is projected to the plane of the sky. We demonstrates for the first time that such substructure may form even in the disk-forming stage, where the surrounding materials actively accrete toward a disk-protostar system.

Centaurs are minor planets thought to have originated in the outer Solar System region known as the Kuiper Belt. Active Centaurs enigmatically display comet-like features (e.g., tails, comae) even though they orbit in the gas giant region where it is too cold for water to readily sublimate. Only 18 active Centaurs have been identified since 1927 and, consequently, the underlying activity mechanism(s) have remained largely unknown up to this point. Here we report the discovery of activity emanating from Centaur 2014 OG392, based on archival images we uncovered plus our own new observational evidence acquired with the Dark Energy Camera (Cerro Tololo Inter-American Observatory Blanco 4 m telescope), the Inamori-Magellan Areal Camera & Spectrograph (Las Campanas Observatory 6.5 m Walter Baade Telescope) and the Large Monolithic Imager (Lowell Observatory 4.3 m Discovery Channel Telescope). We detect a coma as far as 400,000 km from 2014 OG392, and our novel analysis of sublimation processes and dynamical lifetime suggest carbon dioxide and/or ammonia are the most likely candidates for causing activity on this and other active Centaurs. We find 2014 OG392 is optically red, but CO2 and NH3 are spectrally neutral in this wavelength regime so the reddening agent is as yet unidentified.

We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and Flat Spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and Flat Spectrum protostars are 44.9$^{+5.8}_{-3.4}$, 37.0$^{+4.9}_{-3.0}$, and 28.5$^{+3.7}_{-2.3}$ au, respectively, and the mean protostellar dust disk masses are 25.9$^{+7.7}_{-4.0}$, 14.9$^{+3.8}_{-2.2}$, 11.6$^{+3.5}_{-1.9}$ Earth masses, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% out of 379 detected 0.87 mm continuum sources plus 42 non-detections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than that of Class II disks by a factor of $>$4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase.

The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($\approx$2.5"), sensitivity (a 1$\sigma$ goal of 70 $\mu$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{\circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations.

We present ALMA (0.87~mm) and VLA (9~mm) observations toward OMC2-FIR4 and OMC2-FIR3 within the Orion integral-shaped filament that are thought to be the nearest regions of intermediate mass star formation. We characterize the continuum sources within these regions on $\sim$40~AU (0\farcs1) scales and associated molecular line emission at a factor of $\sim$30 better resolution than previous observations at similar wavelengths. We identify six compact continuum sources within OMC2-FIR4, four in OMC2-FIR3, and one additional source just outside OMC2-FIR4. This continuum emission is tracing the inner envelope and/or disk emission on less than 100~AU scales. HOPS-108 is the only protostar in OMC2-FIR4 that exhibits emission from high-excitation transitions of complex organic molecules (e.g., methanol and other lines) coincident with the continuum emission. HOPS-370 in OMC2-FIR3 with L~$\sim$~360~\lsun, also exhibits emission from high-excitation methanol and other lines. The methanol emission toward these two protostars is indicative of temperatures high enough to thermally evaporate methanol from icy dust grains; overall these protostars have characteristics similar to hot corinos. We do not identify a clear outflow from HOPS-108 in \twco, but find evidence of interaction between the outflow/jet from HOPS-370 and the OMC2-FIR4 region. The multitude of observational constraints indicate that HOPS-108 is likely a low to intermediate-mass protostar in its main mass accretion phase and it is the most luminous protostar in OMC2-FIR4. The high resolution data presented here are essential for disentangling the embedded protostars from their surrounding dusty environments and characterizing them.

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 present archival observations demonstrating that main belt asteroid (6478) Gault has an extensive history of comet-like activity. Outbursts have taken place during multiple epochs since 2013 and at distances extending as far as 2.68 au, nearly aphelion. (6478) Gault is a member of the predominately S-type (i.e., volatile-poor) Phocaea family; no other main belt object of this type has ever shown more than a single activity outburst. Furthermore, our data suggest that this is the longest duration of activity caused by a body spinning near the rotational breakup barrier. If activity is indeed unrelated to volatiles, as appears to be the case, (6478) Gault represents a new class of object, perpetually active due to rotational spin-up.

We present multi-epoch continuum observations of the Class 0 protostellar system IRAS 16293-2422 taken with the Very Large Array (VLA) at multiple wavelengths between 7 mm and 15 cm (41 GHz down to 2 GHz), as well as single-epoch Atacama Large Millimeter/submillimeter Array (ALMA) continuum observations covering the range from 0.4 to 1.3 mm (700 GHz down to 230 GHz). The new VLA observations confirm that source A2 is a protostar driving episodic mass ejections, and reveal the complex relative motion between A2 and A1. The spectrum of component B can be described by a single power law ($S_\nu \propto \nu^{2.28}$) over the entire range from 3 to 700 GHz (10 cm down to 0.4 mm), suggesting that the emission is entirely dominated by dust even at $\lambda$ = 10 cm. Finally, the size of source B appears to increase with frequency up to 41 GHz, remaining roughly constant (at $0''.39$ $\equiv$ 55 AU) at higher frequencies. We interpret this as evidence that source B is a dusty structure of finite size that becomes increasingly optically thick at higher frequencies until, in the millimeter regime, the source becomes entirely optically thick. The lack of excess free-free emission at long wavelengths, combined with the absence of high-velocity molecular emission indicates that source B does not drive a powerful outflow, and might indicate that source B is at a particularly early stage of its evolution.

Disk winds are thought to play an important role in the evolution and dispersal of planet-forming disks. While high-resolution optical and infrared spectroscopy has identified several disk wind diagnostics, wind mass loss rates remain largely unconstrained mostly due to the lack of spatial resolution to measure the extent of the wind emitting region. Here, we show that the ngVLA will have the sensitivity and resolution to detect and spatially resolve the free-free emission from the fully or partially ionized component of disk winds. Hence, ngVLA observations will be critical to estimate mass loss rates and clarify the role of disk winds in the evolution and dispersal of disk mass.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of multiple protostar systems in the Perseus molecular cloud previously detected by the Karl G. Jansky Very Large Array (VLA). We observed 17 close ($<$600~AU separation) multiple systems at 1.3~mm in continuum and five molecular lines (i.e., \twco, \cateo, \thco, H$_2$CO, SO) to characterize the circum-multiple environments in which these systems are forming. We detect at least one component in the continuum for the 17 multiple systems. In three systems, one companion is not detected, and for two systems the companions are unresolved at our observed resolution. We also detect circum-multiple dust emission toward 8 out of 9 Class 0 multiples. Circum-multiple dust emission is not detected toward any of the 8 Class I multiples. Twelve systems are detected in the dense gas tracers toward their disks/inner envelopes. For these 12 systems, we use the dense gas observations to characterize their formation mechanism. The velocity gradients in the circum-multiple gas are clearly orthogonal to the outflow directions in 8 out of the 12 systems, consistent with disk fragmentation. Moreover, only two systems with separations $<$200~AU are \textit{inconsistent} with disk fragmentation, in addition to the two widest systems ($>$500~AU). Our results suggest that disk fragmentation via gravitational instability is an important formation mechanism for close multiple systems, but further statistics are needed to better determine the relative fraction formed via this method.

We present the full disk-fit results VANDAM survey of all Class 0 and I protostars in the Perseus molecular cloud. We have 18 new protostellar disk candidates around Class 0 and I sources, which are well described by a simple, parametrized disk model fit to the 8 mm VLA dust-continuum observations. 33% of Class 0 protostars and just 11% of Class I protostars have candidate disks, while 78% of Class 0 and I protostars do not have signs of disks within our 12 AU disk diameter resolution limit, indicating that at 8 mm most disks in the Class 0 and I phases are <10 AU in radius. These small radii may be a result of surface brightness sensitivity limits. Modeled 8 mm radii are similar to the radii of known Class 0 disks with detected Keplerian rotation. Since our 8 mm data trace a population of larger dust grains which radially drift towards the protostar and are lower limits on true disk sizes, large disks at early times do not seem to be particularly rare. We find statistical evidence that Class 0 and I disks are likely drawn from the same distribution, meaning disk properties may be defined early in the Class 0 phase and do not undergo large changes through the Class I phase. By combining our candidate disk properties with previous polarization observations, we find a qualitative indication that misalignment between inferred envelope-scale magnetic fields and outflows may indicate disks on smaller scales in Class 0 sources.

The Martian obliquity cycle is predominately influenced by Solar and Jovian tidal forces. The present-day axial tilt of Mars (25 degrees) is predicted to cycle between 0 degrees and 60 degrees (Mischna & Richardson 2005) with excursions up to 80 degrees(Laskar et al. 2004). We focused on Arabia Terra, a region especially well-suited for studying glacial migration due to its location in the northern mid-latitude (20-40 degrees North, 0-30 degrees East) zone, and its heavily cratered terrain. We created glacial and fluvial morphology distribution maps with the ESRI Aeronautical Reconnaissance Coverage Geographic Information System (ArcGIS). We determined the ages of our craters in order to establish a time-domain link between the Martian obliquity cycle and glacial migration; we found the average crater age to be roughly 2 billion years old.

Active asteroids behave dynamically like asteroids but display comet-like comae. These objects are poorly understood, with only about 30 identified to date. We have conducted one of the deepest systematic searches for asteroid activity by making use of deep images from the Dark Energy Camera (DECam) ideally suited to the task. We looked for activity indicators amongst 11,703 unique asteroids extracted from 35,640 images. We detected three previously-identified active asteroids ((62412), (1) Ceres and (779) Nina), though only (62412) showed signs of activity. Our activity occurrence rate of 1 in 11,703 is consistent with the prevailing 1 in 10,000 activity occurrence rate estimate. Our proof of concept demonstrates 1) our novel informatics approach can locate active asteroids and 2) DECam data are well-suited to the search for active asteroids.

Emission from protostars at centimeter radio wavelengths has been shown to trace the free-free emission arising from ionizing shocks as a result of jets and outflows driven by protostars. Therefore, measuring properties of protostars at radio frequencies can provide valuable insights into the nature of their outflows and jets. We present a C-band (4.1 cm and 6.4 cm) survey of all known protostars (Class 0 and Class I) in Perseus as part of the VLA Nascent Disk and Multiplicity (VANDAM) Survey. We examine the known correlations between radio flux density and protostellar parameters such as bolometric luminosity and outflow force, for our sample. We also investigate the relationship between radio flux density and far-infrared line luminosities from Herschel. We show that free-free emission originates most likely from J-type shocks; however, the large scatter indicates that those two types of emission probe different time and spatial scales. Using C-band fluxes, we removed an estimation of free-free contamination from the corresponding Ka-band (9 mm) flux densities that primarily probe dust emission from embedded disks. We find that the compact ($<$~1") dust emission is lower for Class I sources (median dust mass 96 M$_{\oplus}$) relative to Class 0 (248 M$_{\oplus}$), but several times higher than in Class II (5-15 M$_{\oplus}$). If this compact dust emission is tracing primarily the embedded disk, as is likely for many sources, this result provides evidence for decreasing disk masses with protostellar evolution, with sufficient mass for forming giant planet cores primarily at early times.

We present high spatial resolution observations of the continuum emission from the young multiple star system UZ Tau at frequencies from 6 to 340 GHz. To quantify the spatial variation of dust emission in the UZ Tau E circumbinary disk, the observed interferometric visibilities are modeled with a simple parametric prescription for the radial surface brightnesses at each frequency. We find evidence that the spectrum steepens with radius in the disk, manifested as a positive correlation between the observing frequency and the radius that encircles a fixed fraction of the emission ($R_{eff} \propto \nu^{0.34 \pm 0.08}$). The origins of this size--frequency relation are explored in the context of a theoretical framework for the growth and migration of disk solids. While that framework can reproduce a similar size--frequency relation, it predicts a steeper spectrum than is observed. Moreover, it comes closest to matching the data only on timescales much shorter ($\le 1$ Myr) than the putative UZ Tau age (~2-3 Myr). These discrepancies are the direct consequences of the rapid radial drift rates predicted by models of dust evolution in a smooth gas disk. One way to mitigate that efficiency problem is to invoke small-scale gas pressure modulations that locally concentrate drifting solids. If such particle traps reach high continuum optical depths at 30-340 GHz with a ~30-60% filling fraction in the inner disk ($r \lesssim20$ au), they can also explain the observed spatial gradient in the UZ Tau E disk spectrum.

We present Atacama Large Millimeter/sub-millimeter Array (ALMA) Cycle 2 observations of the 1.3 mm dust continuum emission of the protoplanetary disc surrounding the T Tauri star Elias 24 with an angular resolution of $\sim 0.2"$ ($\sim 28$ au). The dust continuum emission map reveals a dark ring at a radial distance of $0.47"$ ($\sim 65$ au) from the central star, surrounded by a bright ring at $0.58"$ ($\sim 81$ au). In the outer disc, the radial intensity profile shows two inflection points at $0.71"$ and $0.87"$ ($\sim 99$ and $121$ au respectively). We perform global three-dimensional smoothed particle hydrodynamic gas/dust simulations of discs hosting a migrating and accreting planet. Combining the dust density maps of small and large grains with three dimensional radiative transfer calculations, we produce synthetic ALMA observations of a variety of disc models in order to reproduce the gap- and ring-like features observed in Elias 24. We find that the dust emission across the disc is consistent with the presence of an embedded planet with a mass of $\sim 0.7\, \mathrm{M_{\mathrm{J}}}$ at an orbital radius of $\sim$ 60 au. Our model suggests that the two inflection points in the radial intensity profile are due to the inward radial motion of large dust grains from the outer disc. The surface brightness map of our disc model provides a reasonable match to the gap- and ring-like structures observed in Elias 24, with an average discrepancy of $\sim$ 5% of the observed fluxes around the gap region.

We report the first sub-arcsecond VLA imaging of 6 GHz continuum, methanol maser, and excited-state hydroxyl maser emission toward the massive protostellar cluster NGC6334I following the recent 2015 outburst in (sub)millimeter continuum toward MM1, the strongest (sub)millimeter source in the protocluster. In addition to detections toward the previously known 6.7 GHz Class II methanol maser sites in the hot core MM2 and the UCHII region MM3 (NGC6334F), we find new maser features toward several components of MM1, along with weaker features $\sim1''$ north, west, and southwest of MM1, and toward the non-thermal radio continuum source CM2. None of these areas have heretofore exhibited Class II methanol maser emission in three decades of observations. The strongest MM1 masers trace a dust cavity, while no masers are seen toward the strongest dust sources MM1A, 1B and 1D. The locations of the masers are consistent with a combination of increased radiative pumping due to elevated dust grain temperature following the outburst, the presence of infrared photon propagation cavities, and the presence of high methanol column densities as indicated by ALMA images of thermal transitions. The non-thermal radio emission source CM2 ($2''$ north of MM1) also exhibits new maser emission from the excited 6.035 and 6.030 GHz OH lines. Using the Zeeman effect, we measure a line-of-sight magnetic field of +0.5 to +3.7 mG toward CM2. In agreement with previous studies, we also detect numerous methanol and excited OH maser spots toward the UCHII region MM3, with predominantly negative line-of-sight magnetic field strengths of -2 to -5 mG and an intriguing south-north field reversal.

Centimeter continuum emission from protostars offers insight into the innermost part of the outflows, as shock-ionized gas produces free-free emission. We observed a complete population of Class 0 and I protostars in the Perseus molecular cloud at 4.1 cm and 6.4 cm with resolution and sensitivity superior to previous surveys. From a total of 71 detections, 8 sources exhibit resolved emission at 4.1 cm and/or 6.4 cm. In this paper we focus on this sub-sample, analyzing their spectral indices along the jet, and their alignment with respect to the large-scale molecular outflow. Spectral indices for fluxes integrated toward the position of the protostar are consistent with free-free thermal emission. The value of the spectral index along a radio jet decreases with distance from the protostar. For six sources, emission is well aligned with the outflow central axis, showing that we observe the ionized base of the jet. This is not the case for two sources, where we note misalignment of the emission with respect to the large-scale outflow. This might indicate that the emission does not originate in the radio jet, but rather in an ionized outflow cavity wall or disk surface. For five of the sources, the spectral indices along the jet decrease well below the thermal free-free limit of -0.1 with $>2\sigma$ significance. This is indicative of synchrotron emission, meaning that high energy electrons are being produced in the outflows close to the disk. This result can have far-reaching implications for the chemical composition of the embedded disks.

We present a high angular resolution ($\sim 0.2^{\prime\prime}$), high sensitivity ($\sigma \sim 0.2$ mJy) survey of the 870 $\mu$m continuum emission from the circumstellar material around 49 pre-main sequence stars in the $\rho$ Ophiuchus molecular cloud. Because most millimeter instruments have resided in the northern hemisphere, this represents the largest high-resolution, millimeter-wave survey of the circumstellar disk content of this cloud. Our survey of 49 systems comprises 63 stars; we detect disks associated with 29 single sources, 11 binaries, 3 triple systems and 4 transition disks. We present flux and radius distributions for these systems; in particular, this is the first presentation of a reasonably complete probability distribution of disk radii at millimeter-wavelengths. We also compare the flux distribution of these protoplanetary disks with that of the disk population of the Taurus-Auriga molecular cloud. We find that disks in binaries are both significantly smaller and have much less flux than their counterparts around isolated stars. We compute truncation calculations on our binary sources and find that these disks are too small to have been affected by tidal truncation and posit some explanations for this. Lastly, our survey found 3 candidate gapped disks, one of which is a newly identified transition disk with no signature of a dip in infrared excess in extant observations.

The recently discovered minor body 1I/2017 U1 (`Oumuamua) is the first known object in our Solar System that is not bound by the Sun's gravity. Its hyperbolic orbit (eccentricity greater than unity) strongly suggests that it originated outside our Solar System; its red color is consistent with substantial space weathering experienced over a long interstellar journey. We carry out an simple calculation of the probability of detecting such an object. We find that the observed detection rate of 1I-like objects can be satisfied if the average mass of ejected material from nearby stars during the process of planetary formation is ~20 Earth masses, similar to the expected value for our Solar System. The current detection rate of such interstellar interlopers is estimated to be 0.2/year, and the expected number of detections over the past few years is almost exactly one. When the Large Synoptic Survey Telescope begins its wide, fast, deep all-sky survey the detection rate will increase to 1/year. Those expected detections will provide further constraints on nearby planetary system formation through a better estimate of the number and properties of interstellar objects.

We conducted a detailed radiative transfer modeling of the dust emission from the circumstellar disk around HL Tau. The goal of our study is to derive the surface density profile of the inner disk and its structure. In addition to the Atacama Large Millimeter/submillimeter Array images at Band 3 (2.9mm), Band 6 (1.3mm), and Band 7 (0.87mm), the most recent Karl G. Jansky Very Large Array (VLA) observations at 7mm were included in the analysis. A simulated annealing algorithm was invoked to search for the optimum model. The radiative transfer analysis demonstrates that most radial components (i.e., >6AU) of the disk become optically thin at a wavelength of 7mm, which allows us to constrain, for the first time, the dust density distribution in the inner region of the disk. We found that a homogeneous grain size distribution is not sufficient to explain the observed images at different wavelengths simultaneously, while models with a shallower grain size distribution in the inner disk work well. We found clear evidence that larger grains are trapped in the first bright ring. Our results imply that dust evolution has already taken place in the disk at a relatively young (i.e., ~1Myr) age. We compared the midplane temperature distribution, optical depth, and properties of various dust rings with those reported previously. Using the Toomre parameter, we briefly discussed the gravitational instability as a potential mechanism for the origin of the dust clump detected in the first bright ring via the VLA observations.

We present the first images of the 28SiO v=1, J=2-1 maser emission around the closest known massive young stellar object Orion Source I observed at 86 GHz (3mm) with the VLBA. These images have high spatial (~0.3 mas) and spectral (~0.054 km/s) resolutions. We find that the 3mm masers lie in an X-shaped locus consisting of four arms, with blue-shifted emission in the south and east arms and red-shifted emission in the north and west arms. Comparisons with previous images of the 28SiO v=1,2, J=1-0 transitions at 7mm (observed in 2001-2002) show that the bulk of the J=2-1 transition emission follows the streamlines of the J=1-0 emission and exhibits an overall velocity gradient consistent with the gradient at 7mm. While there is spatial overlap between the 3mm and 7mm transitions, the 3mm emission, on average, lies at larger projected distances from Source I (~44 AU compared with ~35 AU for 7mm). The spatial overlap between the v=1, J=1-0 and J=2-1 transitions is suggestive of a range of temperatures and densities where physical conditions are favorable for both transitions of a same vibrational state. However, the observed spatial offset between the bulk of emission at 3mm and 7mm possibly indicates different ranges of temperatures and densities for optimal excitation of the masers. We discuss different maser pumping models that may explain the observed offset. We interpret the 3mm and 7mm masers as being part of a single wide-angle outflow arising from the surface of an edge-on disk rotating about a northeast-southwest axis, with a continuous velocity gradient indicative of differential rotation consistent with a Keplerian profile in a high-mass proto-binary.

Young stellar objects are known to exhibit strong radio variability on timescales of weeks to months, and a few reports have documented extreme radio flares, with at least an order of magnitude change in flux density on timescales of hours to days. However, there have been few constraints on the occurrence rate of such radio flares or on the correlation with pre-main sequence X-ray flares, although such correlations are known for the Sun and nearby active stars. Here we report simultaneous deep VLA radio and Chandra X-ray observations of the Orion Nebula Cluster, targeting hundreds of sources to look for the occurrence rate of extreme radio variability and potential correlation with the most extreme X-ray variability. We identify 13 radio sources with extreme radio variability, with some showing an order of magnitude change in flux density in less than 30 minutes. All of these sources show X-ray emission and variability, but only on timescales <1h do we find clear correlations with extreme radio flaring. Strong X-ray variability does not predict the extreme radio sources and vice versa. Radio flares thus provide us with a new perspective on high-energy processes in YSOs and the irradiation of their protoplanetary disks. Finally, our results highlight implications for interferometric imaging of sources violating the constant-sky assumption.

We present $\sim1-4"$ resolution Very Large Array (VLA) observations of four CH$_3$OH $J_2-J_1$-$E$ 25~GHz transitions ($J$=3, 5, 8, 10) along with 1.3~cm continuum toward 20 regions of active massive star formation containing Extended Green Objects (EGOs), 14 of which we have previously studied with the VLA in the Class~I 44~GHz and Class~II 6.7~GHz maser lines (Cyganowski et al. 2009). Sixteen regions are detected in at least one 25~GHz line ($J$=5), with 13 of 16 exhibiting maser emission. In total, we report 34 new sites of CH$_3$OH maser emission and ten new sites of thermal CH$_3$OH emission, significantly increasing the number of 25~GHz Class I CH$_3$OH masers observed at high angular resolution. We identify probable or likely maser counterparts at 44~GHz for all 15 of the 25~GHz masers for which we have complementary data, providing further evidence that these masers trace similar physical conditions despite uncorrelated flux densities. The sites of thermal and maser emission of CH$_3$OH are both predominantly associated with the 4.5 $\mu$m emission from the EGO, and the presence of thermal CH$_3$OH emission is accompanied by 1.3~cm continuum emission in 9 out of 10 cases. Of the 19 regions that exhibit 1.3~cm continuum emission, it is associated with the EGO in 16 cases (out of a total of 20 sites), 13 of which are new detections at 1.3~cm. Twelve of the 1.3~cm continuum sources are associated with 6.7~GHz maser emission and likely trace deeply-embedded massive protostars.

Based on sub-arcsecond Atacama Large Millimeter/submillimeter Array (ALMA) and Submillimeter Array (SMA) 1.3 mm continuum images of the massive protocluster NGC 6334I obtained in 2015 and 2008, we find that the dust emission from MM1 has increased by a factor of 4.0$\pm$0.3 during the intervening years, and undergone a significant change in morphology. The continuum emission from the other cluster members (MM2, MM4 and the UCHII region MM3 = NGC 6334F) has remained constant. Long term single-dish maser monitoring at HartRAO finds that multiple maser species toward NGC 6334I flared beginning in early 2015, a few months before our ALMA observation, and some persist in that state. New ALMA images obtained in 2016 July-August at 1.1 and 0.87 mm confirm the changes with respect to SMA 0.87 mm images from 2008, and indicate that the (sub)millimeter flaring has continued for at least a year. The excess continuum emission, centered on the hypercompact HII region MM1B, is extended and elongated ($1.6" \times 1.0" \approx 2100 \times 1300$ au) with multiple peaks, suggestive of general heating of the surrounding subcomponents of MM1, some of which may trace clumps in a fragmented disk rather than separate protostars. In either case, these remarkable increases in maser and dust emission provide direct observational evidence of a sudden accretion event in the growth of a massive protostar yielding a sustained luminosity surge by a factor of $70\pm20$, analogous to the largest events in simulations by Meyer et al. (2017). This target provides an excellent opportunity to assess the impact of such a rare event on a protocluster over many years.

Binary and multiple star systems are a frequent outcome of the star formation process, and as a result, almost half of all sun-like stars have at least one companion star. Theoretical studies indicate that there are two main pathways that can operate concurrently to form binary/multiple star systems: large scale fragmentation of turbulent gas cores and filaments or smaller scale fragmentation of a massive protostellar disk due to gravitational instability. Observational evidence for turbulent fragmentation on scales of $>$1000~AU has recently emerged. Previous evidence for disk fragmentation was limited to inferences based on the separations of more-evolved pre-main sequence and protostellar multiple systems. The triple protostar system L1448 IRS3B is an ideal candidate to search for evidence of disk fragmentation. L1448 IRS3B is in an early phase of the star formation process, likely less than 150,000 years in age, and all protostars in the system are separated by $<$200~AU. Here we report observations of dust and molecular gas emission that reveal a disk with spiral structure surrounding the three protostars. Two protostars near the center of the disk are separated by 61 AU, and a tertiary protostar is coincident with a spiral arm in the outer disk at a 183 AU separation. The inferred mass of the central pair of protostellar objects is $\sim$1 M$_{sun}$, while the disk surrounding the three protostars has a total mass of $\sim$0.30 M$_{\sun}$. The tertiary protostar itself has a minimum mass of $\sim$0.085 M$_{sun}$. We demonstrate that the disk around L1448 IRS3B appears susceptible to disk fragmentation at radii between 150~AU and 320~AU, overlapping with the location of the tertiary protostar. This is consistent with models for a protostellar disk that has recently undergone gravitational instability, spawning one or two companion stars.

Gravitational forces are expected to excite spiral density waves in protoplanetary disks, disks of gas and dust orbiting young stars. However, previous observations that showed spiral structure were not able to probe disk midplanes, where most of the mass is concentrated and where planet formation takes place. Using the Atacama Large Millimeter/submillimeter Array we detected a pair of trailing symmetric spiral arms in the protoplanetary disk surrounding the young star Elias 2-27. The arms extend to the disk outer regions and can be traced down to the midplane. These millimeter-wave observations also reveal an emission gap closer to the star than the spiral arms. We argue that the observed spirals trace shocks of spiral density waves in the midplane of this young disk.

We present VLA and ALMA imaging of the deeply-embedded protostellar cluster NGC6334I from 5 cm to 1.3 mm at angular resolutions as fine as 0.17 (220 AU). The dominant hot core MM1 is resolved into seven components at 1.3 mm, clustered within a radius of 1000 AU. Four of the components have brightness temperatures >200 K, radii ~300 AU, minimum luminosities ~10000 Lsun, and must be centrally heated. We term this new phenomenon a "hot multi-core". Two of these objects also exhibit compact free-free emission at longer wavelengths, consistent with a hypercompact HII region (MM1B) and a jet (MM1D). The spatial kinematics of the water maser emission centered on MM1D are consistent with it being the origin of the high-velocity bipolar molecular outflow seen in CO. The close proximity of MM1B and MM1D (440~AU) suggests a proto-binary or a transient bound system. Several components of MM1 exhibit steep millimeter SEDs indicative of either unusual dust spectral properties or time variability. In addition to resolving MM1 and the other hot core (MM2) into multiple components, we detect five new millimeter and two new centimeter sources. Water masers are detected for the first time toward MM4A, confirming its membership in the protocluster. With a 1.3 mm brightness temperature of 97 K coupled with a lack of thermal molecular line emission, MM4A appears to be a highly optically-thick 240 Lsun dust core, possibly tracing a transient stage of massive protostellar evolution. The nature of the strongest water maser source CM2 remains unclear due to its combination of non-thermal radio continuum and lack of dust emission.

There is substantial evidence for disk formation taking place during the early stages of star formation and for most stars being born in multiple systems; however, protostellar multiplicity and disk searches have been hampered by low resolution, sample bias, and variable sensitivity. We have conducted an unbiased, high-sensitivity Karl G. Jansky Very Large Array (VLA) survey toward all known protostars (n = 94) in the Perseus molecular cloud (d~230 pc), with a resolution of ~15 AU (0.06") at 8 mm. We have detected candidate protostellar disks toward 17 sources (with 12 of those in the Class 0 stage) and we have found substructure on < 50AU scales for three Class 0 disk candidates, possibly evidence for disk fragmentation. We have discovered 16 new multiple systems (or new components) in this survey; the new systems have separations < 500 AU and 3 by < 30 AU. We also found a bi-modal distribution of separations, with peaks at ~75 AU and ~3000 AU, suggestive of formation through two distinct mechanisms: disk and turbulent fragmentation. The results from this survey demonstrate the necessity and utility of uniform, unbiased surveys of protostellar systems at millimeter and centimeter wavelengths.

The first long-baseline ALMA campaign resolved the disk around the young star HL Tau into a number of axisymmetric bright and dark rings. Despite the very young age of HL Tau these structures have been interpreted as signatures for the presence of (proto)planets. The ALMA images triggered numerous theoretical studies based on disk-planet interactions, magnetically driven disk structures, and grain evolution. Of special interest are the inner parts of disks, where terrestrial planets are expected to form. However, the emission from these regions in HL Tau turned out to be optically thick at all ALMA wavelengths, preventing the derivation of surface density profiles and grain size distributions. Here, we present the most sensitive images of HL Tau obtained to date with the Karl G. Jansky Very Large Array at 7.0 mm wavelength with a spatial resolution comparable to the ALMA images. At this long wavelength the dust emission from HL Tau is optically thin, allowing a comprehensive study of the inner disk. We obtain a total disk dust mass of 0.001 - 0.003 Msun, depending on the assumed opacity and disk temperature. Our optically thin data also indicate fast grain growth, fragmentation, and formation of dense clumps in the inner densest parts of the disk. Our results suggest that the HL Tau disk may be actually in a very early stage of planetary formation, with planets not already formed in the gaps but in the process of future formation in the bright rings.

Observations of debris disks offer a window into the physical and dynamical properties of planetesimals in extrasolar systems through the size distribution of dust grains. In particular, the millimeter spectral index of thermal dust emission encodes information on the grain size distribution. We have made new VLA observations of a sample of seven nearby debris disks at 9 mm, with 3" resolution and $\sim5$ $\mu$Jy/beam rms. We combine these with archival ATCA observations of eight additional debris disks observed at 7 mm, together with up-to-date observations of all disks at (sub)millimeter wavelengths from the literature to place tight constraints on the millimeter spectral indices and thus grain size distributions. The analysis gives a weighted mean for the slope of the power law grain size distribution, $n(a)\propto a^{-q}$, of $\langle q \rangle = 3.36\pm0.02$, with a possible trend of decreasing $q$ for later spectral type stars. We compare our results to a range of theoretical models of collisional cascades, from the standard self-similar, steady-state size distribution ($q=3.5$) to solutions that incorporate more realistic physics such as alternative velocity distributions and material strengths, the possibility of a cutoff at small dust sizes from radiation pressure, as well as results from detailed dynamical calculations of specific disks. Such effects can lead to size distributions consistent with the data, and plausibly the observed scatter in spectral indices. For the AU Mic system, the VLA observations show clear evidence of a highly variable stellar emission component; this stellar activity obviates the need to invoke the presence of an asteroid belt to explain the previously reported compact millimeter source in this system.

We present a deep centimeter-wavelength catalog of the Orion Nebula Cluster (ONC), based on a 30h single-pointing observation with the Karl G. Jansky Very Large Array in its high-resolution A-configuration using two 1 GHz bands centered at 4.7 GHz and 7.3 GHz. A total of 556 compact sources were detected in a map with a nominal rms noise of 3 muJy/bm, limited by complex source structure and the primary beam response. Compared to previous catalogs, our detections increase the sample of known compact radio sources in the ONC by more than a factor of seven. The new data show complex emission on a wide range of spatial scales. Following a preliminary correction for the wideband primary-beam response, we determine radio spectral indices for 170 sources whose index uncertainties are less than +/-0.5. We compare the radio to the X-ray and near-infrared point-source populations, noting similarities and differences.

Locating planets in circumstellar Habitable Zones is a priority for many exoplanet surveys. Space-based and ground-based surveys alike require robust toolsets to aid in target selection and mission planning. We present the Catalog of Earth-Like Exoplanet Survey Targets (CELESTA), a database of Habitable Zones around 36,000 nearby stars. We calculated stellar parameters, including effective temperatures, masses, and radii, and we quantified the orbital distances and periods corresponding to the circumstellar Habitable Zones. We gauged the accuracy of our predictions by contrasting CELESTA's computed parameters to observational data. We ascertain a potential return on investment by computing the number of Habitable Zones probed for a given survey duration. A versatile framework for extending the functionality of CELESTA into the future enables ongoing comparisons to new observations, and recalculations when updates to Habitable Zone models, stellar temperatures, or parallax data become available. We expect to upgrade and expand CELESTA using data from the Gaia mission as the data becomes available.

To characterize the mechanisms of planet formation it is crucial to investigate the properties and evolution of protoplanetary disks around young stars, where the initial conditions for the growth of planets are set. Our goal is to study grain growth in the disk of the young, intermediate mass star HD163296 where dust processing has already been observed, and to look for evidence of growth by ice condensation across the CO snowline, already identified in this disk with ALMA. Under the hypothesis of optically thin emission we compare images at different wavelengths from ALMA and VLA to measure the opacity spectral index across the disk and thus the maximum grain size. We also use a Bayesian tool based on a two-layer disk model to fit the observations and constrain the dust surface density. The measurements of the opacity spectral index indicate the presence of large grains and pebbles ($\geq$1 cm) in the inner regions of the disk (inside $\sim$50 AU) and smaller grains, consistent with ISM sizes, in the outer disk (beyond 150 AU). Re-analysing ALMA Band 7 Science Verification data we find (radially) unresolved excess continuum emission centered near the location of the CO snowline at $\sim$90 AU. Our analysis suggests a grain size distribution consistent with an enhanced production of large grains at the CO snowline and consequent transport to the inner regions. Our results combined with the excess in infrared scattered light found by Garufi et al. (2014) suggests the presence of a structure at 90~AU involving the whole vertical extent of the disk. This could be evidence for small scale processing of dust at the CO snowline.

We present a multiplicity study of all known protostars (94) in the Perseus molecular cloud from a Karl G. Jansky Very Large Array (VLA) survey at Ka-band (8 mm and 1 cm) and C-band (4 cm and 6.6 cm). The observed sample has a bolometric luminosity range between 0.1 L$_{\odot}$ and $\sim$33 L$_{\odot}$, with a median of 0.7 L$_{\odot}$. This multiplicity study is based on the Ka-band data, having a best resolution of $\sim$0.065" (15 AU) and separations out to $\sim$43" (10000 AU) can be probed. The overall multiplicity fraction (MF) is found to be of 0.40$\pm$0.06 and the companion star fraction (CSF) is 0.71$\pm$0.06. The MF and CSF of the Class 0 protostars are 0.57$\pm$0.09 and 1.2$\pm$0.2, and the MF and CSF of Class I protostars are both 0.23$\pm$0.08. The distribution of companion separations appears bi-modal, with a peak at $\sim$75 AU and another peak at $\sim$3000 AU. Turbulent fragmentation is likely the dominant mechanism on $>$1000 AU scales and disk fragmentation is likely to be the dominant mechanism on $<$200 AU scales. Toward three Class 0 sources we find companions separated by $<$30 AU. These systems have the smallest separations of currently known Class 0 protostellar binary systems. Moreover, these close systems are embedded within larger (50 AU to 400 AU) structures and may be candidates for ongoing disk fragmentation.

We present the first dust emission results toward a sample of seven protostellar disk candidates around Class 0 and I sources in the Perseus molecular cloud from the VLA Nascent Disk and Multiplicity (VANDAM) survey with ~0.05'' or 12 AU resolution. To examine the surface brightness profiles of these sources, we fit the Ka-band 8 mm dust-continuum data in the u,v-plane to a simple, parametrized model based on the Shakura-Sunyaev disk model. The candidate disks are well-fit by a model with a disk-shaped profile and have masses consistent with known Class 0 and I disks. The inner-disk surface densities of the VANDAM candidate disks have shallower density profiles compared to disks around more evolved Class II systems. The best-fit model radii of the seven early-result candidate disks are R_c > 10 AU; at 8 mm, the radii reflect lower limits on the disk size since dust continuum emission is tied to grain size and large grains radially drift inwards. These relatively large disks, if confirmed kinematically, are inconsistent with theoretical models where the disk size is limited by strong magnetic braking to < 10 AU at early times.

Theoretical models of grain growth predict dust properties to change as a function of protoplanetary disk radius, mass, age and other physical conditions. We lay down the methodology for a multi-wavelength analysis of (sub-)mm and cm continuum interferometric observations to constrain self-consistently the disk structure and the radial variation of the dust properties. The computational architecture is massively parallel and highly modular. The analysis is based on the simultaneous fit in the uv-plane of observations at several wavelengths with a model for the disk thermal emission and for the dust opacity. The observed flux density at the different wavelengths is fitted by posing constraints on the disk structure and on the radial variation of the grain size distribution. We apply the analysis to observations of three protoplanetary disks (AS 209, FT Tau, DR Tau) for which a combination of spatially resolved observations in the range ~0.88mm to ~10mm is available (from SMA, CARMA, and VLA), finding evidence of a decreasing maximum dust grain size (a_max) with radius. We derive large a_max values up to 1 cm in the inner disk between 15 and 30 AU and smaller grains with a_max~1 mm in the outer disk (R > 80AU). In this paper we develop a multi-wavelength analysis that will allow this missing quantity to be constrained for statistically relevant samples of disks and to investigate possible correlations with disk or stellar parameters.

Infrared Dark Clouds (IRDCs) harbor the earliest phases of massive star formation, and many of the compact cores in IRDCs, traced by millimeter continuum or by molecular emission in high critical density lines, host massive young stellar objects (YSOs). We used the Robert C. Byrd Green Bank Telescope (GBT) and the Karl G. Jansky Very Large Array (VLA) to map NH$_{3}$ and CCS in nine IRDCs to reveal the temperature, density, and velocity structures and explore chemical evolution in the dense ($>10^{22}$ cm$^{-2}$) gas. Ammonia is an excellent molecular tracer for these cold, dense environments. The internal structure and kinematics of the IRDCs include velocity gradients, filaments, and possibly colliding clumps that elucidate the formation process of these structures and their YSOs. We find a wide variety of substructure including filaments and globules at distinct velocities, sometimes overlapping at sites of ongoing star formation. It appears that these IRDCs are still being assembled from molecular gas clumps even as star formation has already begun, and at least three of them appear consistent with the morphology of ``hub-filament structures'' discussed in the literature. Furthermore, we find that these clumps are typically near equipartition between gravitational and kinetic energies, so these structures may survive for multiple free-fall times.Keywords: molecular data -- ISM: clouds -- (ISM:) dust, extinction -- ISM: molecules -- Stars: formation -- radio lines: ISM

We present continuum and molecular line observations at 230 GHz and 345 GHz from the Sub-millimeter Array (SMA) toward three protostars in the Perseus L1448N region. The data are from the large project "Mass Assembly of Stellar Systems and their Evolution with the SMA" (MASSES). Three dust continuum sources, Source B, Source NW, and Source A, are detected at both frequencies. These sources have corresponding emission peaks in C18O (J=2-1), 13CO (J=2-1), and HCO+ (J=4-3), and have offsets with N2D+ (J=3-2) peaks. High angular resolution data from a complimentary continuum survey with the Karl G. Jansky Very Large Array show that Source B is associated with three 8 mm continuum objects, Source NW with two, and Source A remains single. These results suggest that multiplicity in L1448N exists at different spatial scales from a few thousand AU to < 100 AU. Velocity gradients in each source obtained from two-dimensional fits to the SMA C18O emission are found to be perpendicular to within 20 degrees of the outflow directions as revealed by 12CO (J=2-1). We have observed that Sources B and NW with multiplicity have higher densities than Source A without multiplicity. This suggests that thermal Jeans fragmentation can be relevant in the fragmentation process. However, we have not observed a difference in the ratio between rotational and gravitational energy between sources with and without multiplicity. We also have not observed a trend between non-thermal velocity dispersions and the level of fragmentation. Our study has provided the first direct and comprehensive comparison between multiplicity and core properties in low-mass protostars, although based on small number statistics.

Magnetic fields can regulate disk formation, accretion and jet launching. Until recently, it has been difficult to obtain high resolution observations of the magnetic fields of the youngest protostars in the critical region near the protostar. The VANDAM survey is observing all known protostars in the Perseus Molecular Cloud. Here we present the polarization data of IRAS 4A. We find that with ~ 0.2'' (50 AU) resolution at {\lambda} = 8.1 and 10.3 mm, the inferred magnetic field is consistent with a circular morphology, in marked contrast with the hourglass morphology seen on larger scales. This morphology is consistent with frozen-in field lines that were dragged in by rotating material entering the infall region. The field morphology is reminiscent of rotating circumstellar material near the protostar. This is the first polarization detection of a protostar at these wavelengths. We conclude from our observations that the dust emission is optically thin with {\beta} ~ 1.3, suggesting that mm/cm-sized grains have grown and survived in the short lifetime of the protostar.

We present new results from the Disks@EVLA program for two young stars: CY Tau and DoAr 25. We trace continuum emission arising from their circusmtellar disks from spatially resolved observations, down to tens of AU scales, at {\lambda} = 0.9, 2.8, 8.0, and 9.8 mm for DoAr25 and at {\lambda} = 1.3, 2.8, and 7.1 mm for CY Tau. Additionally, we constrain the amount of emission whose origin is different from thermal dust emission from 5 cm observations. Directly from interferometric data, we find that observations at 7 mm and 1 cm trace emission from a compact disk while millimeter-wave observations trace an extended disk structure. From a physical disk model, where we characterize the disk structure of CY Tau and DoAr 25 at wavelengths shorter than 5 cm, we find that (1) dust continuum emission is optically thin at the observed wavelengths and over the spatial scales studied, (2) a constant value of the dust opacity is not warranted by our observations, and (3) a high-significance radial gradient of the dust opacity spectral index, {\beta}, is consistent with the observed dust emission in both disks, with low-{\beta} in the inner disk and high-{\beta} in the outer disk. Assuming that changes in dust properties arise solely due to changes in the maximum particle size (amax), we constrain radial variations of amax in both disks, from cm-sized particles in the inner disk (R < 40 AU) to millimeter sizes in the outer disk (R > 80 AU). These observational constraints agree with theoretical predictions of the radial-drift barrier, however, fragmentation of dust grains could explain our amax(R) constraints if these disks have lower turbulence and/or if dust can survive high-velocity collisions.

We report detections of two candidate distant submillimeter galaxies (SMGs), MM J154506.4$-$344318 and MM J154132.7$-$350320, which are discovered in the AzTEC/ASTE 1.1 mm survey toward the Lupus-I star-forming region. The two objects have 1.1 mm flux densities of 43.9 and 27.1 mJy, and have Herschel/SPIRE counterparts as well. The Submillimeter Array counterpart to the former SMG is identified at 890 $\mu$m and 1.3 mm. Photometric redshift estimates using all available data from the mid-infrared to the radio suggest that the redshifts of the two SMGs are $z_{\rm photo} \simeq$ 4-5 and 3, respectively. Near-infrared objects are found very close to the SMGs and they are consistent with low-$z$ ellipticals, suggesting that the high apparent luminosities can be attributed to gravitational magnification. The cumulative number counts at $S_{\rm 1.1mm} \ge 25$ mJy, combined with other two 1.1-mm brightest sources, are $0.70 ^{+0.56}_{-0.34}$ deg$^{-2}$, which is consistent with a model prediction that accounts for flux magnification due to strong gravitational lensing. Unexpectedly, a $z > 3$ SMG and a Galactic dense starless core (e.g., a first hydrostatic core) could be similar in the mid-infrared to millimeter spectral energy distributions and spatial structures at least at $\gtrsim 1"$. This indicates that it is necessary to distinguish the two possibilities by means of broad band photometry from the optical to centimeter and spectroscopy to determine the redshift, when a compact object is identified toward Galactic star-forming regions.

We used the Very Large Array to carry out a multi-epoch radio continuum monitoring of the Orion Nebula Cluster and Orion Molecular Cloud. Our observations reveal the presence of 19 sources. With the exception of the sources BN and C the sources show variability between the different epochs. We have found tentative evidence of variability in the massive object related with source I. Our observations also confirm radio flux density variations of a factor >2 on timescales of hours to days in 5 sources. One of these flaring sources, OHC-E, has been detected for the first time. We conclude that the radio emission arises from: i) highly-variable non-thermal gyrosynchrotron emission produced by electrons accelerated in the magnetospheres of pre-main sequence stars; ii) thermal emission from ionized gas and/or heated dust around massive objects and proplyds. Combining our sample with other radio monitoring and a X-ray catalog, we have studied the properties of 51 radio/X-ray stars. We have found severals hints of a direct relation between the X-ray activity and the mechanisms responsible for (at least some fraction of) the radio emission. We have estimated a radio flaring rate of 0.14 flares day-1 in the densest stellar cluster of the region, suggesting that radio flares from young stars are more common events than previously thought.

We present a CARMA 1.3 mm continuum survey toward 9 Class 0 protostars in the Perseus molecular cloud at $\sim$0.3$^{\prime\prime}$ (70 AU) resolution. This study approximately doubles the number of Class 0 protostars observed with spatial resolutions $<$ 100 AU at millimeter wavelengths, enabling the presence of protostellar disks and proto-binary systems to be probed. We detect flattened structures with radii $>$ 100 AU around 2 sources (L1448 IRS2 and Per-emb-14) and these sources may be strong disk candidates. Marginally-resolved structures with position angles within 30$^{\circ}$ of perpendicular to the outflow are found toward 3 protostars (L1448 IRS3C, IRAS 03282+3035, and L1448C) and are considered disk candidates. Two others (L1448 IRS3B and IRAS 03292+3039) have resolved structure, possibly indicative of massive inner envelopes or disks; L1448 IRS3B also has a companion separated by 0.9$^{\prime\prime}$ ($\sim$210 AU). IC348-MMS does not have well-resolved structure and the candidate first hydrostatic core L1451-MMS is marginally resolved on 1$^{\prime\prime}$ scales. The strong disk candidate sources were followed-up with C$^{18}$O ($J=2\rightarrow1$) observations, detecting velocity gradients consistent with rotation, but it is unclear if the rotation is Keplerian. We compare the observed visibility amplitudes to radiative transfer models, finding that visibility amplitude ratios suggest a compact component (possibly a disk) is necessary for 5 of 9 Class 0 sources; envelopes alone may explain the other 4 systems. We conclude that there is evidence for the formation of large disks in the Class 0 phase with a range of radii and masses dependent upon their initial formation conditions.

We are conducting a Jansky VLA Ka-band (8 mm and 1 cm) and C-band (4 cm and 6.4 cm) survey of all known protostars in the Perseus Molecular Cloud, providing resolution down to $\sim$0.06'' and $\sim$0.35" in Ka-band and C-band, respectively. Here we present first results from this survey that enable us to examine the source NGC 1333 IRAS2A in unprecedented detail and resolve it into a proto-binary system separated by 0.621"$\pm$0.006" ($\sim$143 AU) at 8 mm, 1 cm, and 4 cm. These 2 sources (IRAS2A VLA1 and VLA2) are likely driving the two orthogonal outflows known to originate from IRAS2A. The brighter source IRAS2A VLA1 is extended perpendicular to its outflow in the VLA data, with a deconvolved size of 0.055" ($\sim$13 AU), possibly tracing a protostellar disk. The recently reported candidate companions (IRAS2A MM2 and MM3) are not detected in either our VLA data, CARMA 1.3 mm data, or SMA 850 $\mu$m data. SMA CO ($J=3\rightarrow2$), CARMA CO ($J=2\rightarrow1$), and lower resolution CARMA CO ($J=1\rightarrow0$) observations are used to examine the outflow origins and the nature of the candidate companions to IRAS2A VLA1. The CO ($J=3\rightarrow2$) and ($J=2\rightarrow1$) data show that IRAS2A MM2 is coincident with a bright CO emission spot in the east-west outflow, and IRAS2A MM3 is within the north-south outflow. In contrast, IRAS2A VLA2 lies at the east-west outflow symmetry point. We propose that IRAS2A VLA2 is the driving source of the East-West outflow and a true companion to IRAS2A VLA1, whereas IRAS2A MM2 and MM3 may not be protostellar.

We present sub-arcsecond resolution observations of continuum emission associated with the GG Tau quadruple star system at wavelengths of 1.3, 2.8, 7.3, and 50 mm. These data confirm that the GG Tau A binary is encircled by a circumbinary ring at a radius of 235 AU with a FWHM width of ~60 AU. We find no clear evidence for a radial gradient in the spectral shape of the ring, suggesting that the particle size distribution is spatially homogeneous on angular scales of ~0.1". A central point source, likely associated with the primary component (GG Tau Aa), exhibits a composite spectrum from dust and free-free emission. Faint emission at 7.3 mm is observed toward the low-mass star GG Tau Ba, although its origin remains uncertain. Using these measurements of the resolved, multifrequency emission structure of the GG Tau A system, models of the far-infrared to radio spectrum are developed to place constraints on the grain size distribution and dust mass in the circumbinary ring. The non-negligible curvature present in the ring spectrum implies a maximum particle size of 1-10 mm, although we are unable to place strong constraints on the distribution shape. The corresponding dust mass is 30-300 M_earth, at a temperature of 20-30 K. We discuss how this significant concentration of relatively large particles in a narrow ring at a large radius might be produced in a local region of higher gas pressures (i.e., a particle "trap") located near the inner edge of the circumbinary disk.

We present Karl G. Jansky Very Large Array (VLA) observations of the 7 mm continuum emission from the disk surrounding the young star LkCa 15. The observations achieve an angular resolution of 70 mas and spatially resolve the circumstellar emission on a spatial scale of 9 AU. The continuum emission traces a dusty annulus of 45 AU in radius that is consistent with the dust morphology observed at shorter wavelengths. The VLA observations also reveal a compact source at the center of the disk, possibly due to thermal emission from hot dust or ionized gas located within a few AU from the central star. No emission is observed between the star and the dusty ring, and, in particular, at the position of the candidate protoplanet LkCa 15 b. By comparing the observations with theoretical models for circumplanetary disk emission, we find that if LkCa~15~b is a massive planet (>5 M_J) accreting at a rate greater than 1.e-6 M_J yr^{-1}, then its circumplanetary disk is less massive than 0.1 M_J, or smaller than 0.4 Hill radii. Similar constraints are derived for any possible circumplanetary disk orbiting within 45 AU from the central star. The mass estimate are uncertain by at least one order of magnitude due to the uncertainties on the mass opacity. Future ALMA observations of this system might be able to detect circumplanetary disks down to a mass of 5.e-4 M_J and as small as 0.2 AU, providing crucial constraints on the presence of giant planets in the act of forming around this young star.

For over a decade, the structure of the inner cavity in the transition disk of TW Hydrae has been a subject of debate. Modeling the disk with data obtained at different wavelengths has led to a variety of proposed disk structures. Rather than being inconsistent, the individual models might point to the different faces of physical processes going on in disks, such as dust growth and planet formation. Our aim is to investigate the structure of the transition disk again and to find to what extent we can reconcile apparent model differences. A large set of high-angular-resolution data was collected from near-infrared to centimeter wavelengths. We investigated the existing disk models and established a new self-consistent radiative-transfer model. A genetic fitting algorithm was used to automatize the parameter fitting. Simple disk models with a vertical inner rim and a radially homogeneous dust composition from small to large grains cannot reproduce the combined data set. Two modifications are applied to this simple disk model: (1) the inner rim is smoothed by exponentially decreasing the surface density in the inner ~3 AU, and (2) the largest grains (>100 um) are concentrated towards the inner disk region. Both properties can be linked to fundamental processes that determine the evolution of protoplanetary disks: the shaping by a possible companion and the different regimes of dust-grain growth, respectively. The full interferometric data set from near-infrared to centimeter wavelengths requires a revision of existing models for the TW Hya disk. We present a new model that incorporates the characteristic structures of previous models but deviates in two key aspects: it does not have a sharp edge at 4 AU, and the surface density of large grains differs from that of smaller grains. This is the first successful radiative-transfer-based model for a full set of interferometric data.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations in the dust continuum (690 GHz, 0.45 mm) and 12CO J=6-5 spectral line emission, of the transitional disks surrounding the stars SAO 206462 and SR 21. These ALMA observations resolve the dust-depleted disk cavities and extended gaseous disks, revealing large-scale asymmetries in the dust emission of both disks. We modeled these disks structures with a ring and an azimuthal gaussian, where the azimuthal gaussian is motivated by the steady-state vortex solution from Lyra & Lin (2013). Compared to recent observations of HD 142527, Oph IRS 48, and LkHa 330, these are low-contrast (< 2) asymmetries. Nevertheless, a ring alone is not a good fit, and the addition of a vortex prescription describes these data much better. The asymmetric component encompasses 15% and 28% of the total disk emission in SAO 206462 and SR 21 respectively, which corresponds to a lower limit of 2 MJup of material within the asymmetry for both disks. Although the contrast in the dust asymmetry is low, we find that the turbulent velocity inside it must be large (~20% of the sound speed) in order to drive these azimuthally wide and radially narrow vortex-like structures. We obtain residuals from the ring and vortex fitting that are still significant, tracing non-axisymmetric emission in both disks. We compared these submillimeter observations with recently published H-band scattered light observations. For SR 21 the scattered light emission is distributed quite differently from the submillimeter continuum emission, while for SAO 206462 the submillimeter residuals are suggestive of spiral-like structure similar to the near-IR emission.

We present observations of three Class 0/I protostars (L1157-mm, CB230 IRS1, and L1165-SMM1) using the Karl G. Jansky Very Large Array (VLA) and observations of two (L1165-SMM1 and CB230 IRS1) with the Combined Array for Research in Millimeter-wave Astronomy (CARMA). The VLA observations were taken at wavelengths of $\lambda = 7.3$ mm, 1.4 cm, 3.3 cm, 4.0 cm, and 6.5 cm with a best resolution of $\sim$0\farcs06 (18 AU) at 7.3 mm. The L1165-SMM1 CARMA observations were taken at $\lambda = 1.3$ mm with a best resolution of $\sim0\farcs3$ (100 AU), and the CB230 IRS1 observations were taken at $\lambda = 3.4$ mm with a best resolution of $\sim$3\arcsec\ (900 AU). We find that L1165-SMM1 and CB230 IRS1 have probable binary companions at separations of $\sim$0\farcs3 (100 AU) from detections of secondary peaks at multiple wavelengths. The position angles of these companions are nearly orthogonal to the direction of the observed bipolar outflows, consistent with the expected protostellar disk orientations. We suggest that these companions may have formed from disk fragmentation; turbulent fragmentation would not preferentially arrange the binary companions to be orthogonal to the outflow direction. For L1165-SMM1, both the 7.3 mm and 1.3 mm emission show evidence of a large (R $>$ 100 AU) disk. For the L1165-SMM1 primary protostar and the CB230 IRS1 secondary protostar, the 7.3 mm emission is resolved into structures consistent with $\sim20$ AU radius disks. For the other protostars, including L1157-mm, the emission is unresolved, suggesting disks with radii $< 20$ AU.

In Orion BN/KL, proper motions of 7 mm vibrationally-excited SiO masers trace rotation of a nearly edge-on disk and a bipolar wide-angle outflow 10-100 AU from radio Source I, a binary young stellar object (YSO) of ~20 Msun. Here we map ground-state 7 mm SiO emission with the Very Large Array and track proper motions over 9 years. The innermost and strongest emission lies in two extended arcs bracketing Source I. The proper motions trace a northeast-southwest bipolar outflow 100-1000 AU from Source I with a median 3D motion of ~18 km/s. An overlying distribution of 1.3 cm H2O masers betrays similar flow characteristics. Gas dynamics and emission morphology traced by the masers suggest the presence of a magnetocentrifugal disk-wind. Reinforcing evidence lies in the colinearity of the flow, apparent rotation across the flow parallel to the disk rotation, and recollimation that narrows the flow opening angle ~120 AU downstream. The arcs of ground-state SiO emission may mark the transition point to a shocked super-Alfvenic outflow.

We present a comprehensive view of the W51B HII region complex and the W51C supernova remnant (SNR) using new radio observations from the VLA, VLBA, MERLIN, JCMT, and CSO along with archival data from Spitzer, ROSAT, ASCA, and Chandra. Our VLA data include the first 400 cm (74 MHz) continuum image of W51 at high resolution (88 arcsec). The 400 cm image shows non-thermal emission surrounding the G49.2-0.3 HII region, and a compact source of non-thermal emission (W51B_NT) coincident with the previously-identified OH (1720 MHz) maser spots, non-thermal 21 and 90 cm emission, and a hard X-ray source. W51B_NT falls within the region of high likelihood for the position of TeV gamma-ray emission. Using the VLBA three OH (1720 MHz) maser spots are detected in the vicinity of W51B_NT with sizes of 60 to 300 AU and Zeeman effect magnetic field strengths of 1.5 to 2.2 mG. The multiwavelength data demonstrate that the northern end of the W51B HII region complex has been partly enveloped by the advancing W51C SNR and this interaction explains the presence of W51B_NT and the OH masers. This interaction also appears in the thermal molecular gas which partially encircles W51B_NT and exhibits narrow pre-shock (DeltaV 5 km/s) and broad post-shock (DeltaV 20 km/s) velocity components. RADEX radiative transfer modeling of these two components yield physical conditions consistent with the passage of a non-dissociative C-type shock. Confirmation of the W51B/W51C interaction provides additional evidence in favor of this region being one of the best candidates for hadronic particle acceleration known thus far.

The CORNISH project is the highest resolution radio continuum survey of the Galactic plane to date. It is the 5 GHz radio continuum part of a series of multi-wavelength surveys that focus on the northern GLIMPSE region (10 deg < l < 65 deg), observed by the Spitzer satellite in the mid-infrared. Observations with the Very Large Array in B and BnA configurations have yielded a 1.5" resolution Stokes I map with a root-mean-squared noise level better than 0.4 mJy/beam. Here we describe the data-processing methods and data characteristics, and present a new, uniform catalogue of compact radio-emission. This includes an implementation of automatic deconvolution that provides much more reliable imaging than standard CLEANing. A rigorous investigation of the noise characteristics and reliability of source detection has been carried out. We show that the survey is optimised to detect emission on size scales up to 14" and for unresolved sources the catalogue is more than 90 percent complete at a flux density of 3.9 mJy. We have detected 3,062 sources above a 7-sigma detection limit and present their ensemble properties. The catalogue is highly reliable away from regions containing poorly-sampled extended emission, which comprise less than two percent of the survey area. Imaging problems have been mitigated by down-weighting the shortest spacings and potential artefacts flagged via a rigorous manual inspection with reference to the Spitzer infrared data. We present images of the most common source types found: regions, planetary nebulae and radio-galaxies. The CORNISH data and catalogue are available online at http://cornish.leeds.ac.uk

We present ALMA and VLA observations of the molecular and ionized gas at 0.1-0.3 arcsec resolution in the Class 0 protostellar system IRAS 16293-2422. These data clarify the origins of the protostellar outflows from the deeply embedded sources in this complex region. Source A2 is confirmed to be at the origin of the well known large scale north-east--south-west flow. The most recent VLA observations reveal a new ejection from that protostar, demonstrating that it drives an episodic jet. The central compact part of the other known large scale flow in the system, oriented roughly east-west, is well delineated by the CO(6-5) emission imaged with ALMA and is confirmed to be driven from within component A. Finally, a one-sided blueshifted bubble-like outflow structure is detected here for the first time from source B to the north-west of the system. Its very short dynamical timescale (~ 200 yr), low velocity, and moderate collimation support the idea that source B is the youngest object in the system, and possibly one of the youngest protostars known.

We present high angular resolution Submillimeter Array (SMA) and Karl G. Jansky Very Large Array (VLA) observations of the massive protocluster G18.67+0.03. Previously targeted in maser surveys of GLIMPSE Extended Green Objects (EGOs), this cluster contains three Class I methanol maser sources, providing a unique opportunity to test the proposed role of Class I masers as evolutionary indicators for massive star formation. The millimeter observations reveal bipolar molecular outflows, traced by 13CO(2-1) emission, associated with all three Class I maser sources. Two of these sources (including the EGO) are also associated with 6.7 GHz Class II methanol masers; the Class II masers are coincident with millimeter continuum cores that exhibit hot core line emission and drive active outflows, as indicated by the detection of SiO(5-4). In these cases, the Class I masers are coincident with outflow lobes, and appear as clear cases of excitation by active outflows. In contrast, the third Class I source is associated with an ultracompact HII region, and not with Class II masers. The lack of SiO emission suggests the 13CO outflow is a relic, consistent with its longer dynamical timescale. Our data show that massive young stellar objects associated only with Class I masers are not necessarily young, and provide the first unambiguous evidence that Class I masers may be excited by both young (hot core) and older (UC HII) MYSOs within the same protocluster.

We present dust continuum observations of the protoplanetary disk surrounding the pre-main sequence star AS 209, spanning more than an order of magnitude in wavelength from 0.88 to 9.8 mm. The disk was observed with sub-arcsecond angular resolution (0.2"-0.5") to investigate radial variations in its dust properties. At longer wavelengths, the disk emission structure is notably more compact, providing model-independent evidence for changes in the grain properties across the disk. We find that physical models which reproduce the disk emission require a radial dependence of the dust opacity \kappa_{\nu}. Assuming that the observed wavelength-dependent structure can be attributed to radial variations in the dust opacity spectral index (\beta), we find that \beta(R) increases from \beta<0.5 at \sim20 AU to \beta>1.5 for R>80 AU, inconsistent with a constant value of \beta\ across the disk (at the 10\sigma\ level). Furthermore, if radial variations of \kappa_{\nu} are caused by particle growth, we find that the maximum size of the particle-size distribution (a_{max}) increases from sub-millimeter-sized grains in the outer disk (R>70 AU) to millimeter and centimeter-sized grains in the inner disk regions (R< 70 AU). We compare our observational constraint on a_{max}(R) with predictions from physical models of dust evolution in proto-planetary disks. For the dust composition and particle-size distribution investigated here, our observational constraints on a_{max}(R) are consistent with models where the maximum grain size is limited by radial drift.

We describe the motivation, design and implementation of the CORNISH survey, an arcsecond resolution radio continuum survey of the inner Galactic plane at 5 GHz using the Karl G. Jansky Very Large Array (VLA). It is a blind survey co-ordinated with the northern Spitzer GLIMPSE I region covering 10 deg < l < 65 deg and |b| < 1 deg at similar resolution. We discuss in detail the strategy that we employed to control the shape of the synthesised beam across this survey that covers a wide range of fairly low declinations. Two snapshots separated by 4 hours in hour angle kept the beam elongation to less that 1.5 over 75% of the survey area and less than 2 over 98% of the survey. The prime scientific motivation is to provide an unbiased survey for ultra-compact HII regions to study this key phase in massive star formation. A sensitivity around 2 mJy will allow the automatic distinction between radio loud and quiet mid-IR sources found in the Spitzer surveys. This survey has many legacy applications beyond star formation including evolved stars, active stars and binaries, and extragalactic sources. The CORNISH survey for compact ionized sources complements other Galactic plane surveys that target diffuse and non-thermal sources as well as atomic and molecular phases to build up a complete picture of the ISM in the Galaxy.

We present an interferometric kinematic study of morphologically complex protostellar envelopes based on observations of the dense gas tracers N2H+ and NH3. The strong asymmetric nature of most envelopes in our sample leads us to question the common interpretation of velocity gradients as rotation, given the possibility of projection effects in the observed velocities. Several "idealized" sources with well-ordered velocity fields and envelope structures are now analyzed in more detail. We compare the interferometric data to position-velocity diagrams of kinematic models for spherical rotating collapse and filamentary rotating collapse. For this purpose, we developed a filamentary parametrization of the rotating collapse model to explore the effects of geometric projection on the observed velocity structures. We find that most envelopes in our sample have PV structures that can be reproduced by an infalling filamentary envelope projected at different angles within the plane of the sky. The infalling filament produces velocity shifts across the envelope that can mimic rotation, especially when viewed at single-dish resolutions and the axisymmetric rotating collapse model does not uniquely describe any dataset. Furthermore, if the velocities are assumed to reflect rotation, then the inferred centrifugal radii are quite large in most cases, indicating significant fragmentation potential or more likely another component to the line-center velocity. We conclude that ordered velocity gradients cannot be interpreted as rotation alone when envelopes are non-axisymmetric and that projected infall velocities likely dominate the velocity field on scales larger than 1000 AU.

We present a study of dense molecular gas kinematics in seventeen nearby protostellar systems using single-dish and interferometric molecular line observations. The non-axisymmetric envelopes around a sample of Class 0/I protostars were mapped in the N2H+ (J=1-0) tracer with the IRAM 30m, CARMA and PdBI as well as NH3 (1,1) with the VLA. The molecular line emission is used to construct line-center velocity and linewidth maps for all sources to examine the kinematic structure in the envelopes on spatial scales from 0.1 pc to ~1000 AU. The direction of the large-scale velocity gradients from single-dish mapping is within 45 degrees of normal to the outflow axis in more than half the sample. Furthermore, the velocity gradients are often quite substantial, the average being ~2.3 km\s\pc. The interferometric data often reveal small-scale velocity structure, departing from the more gradual large-scale velocity gradients. In some cases, this likely indicates accelerating infall and/or rotational spin-up in the inner envelope; the median velocity gradient from the interferometric data is ~10.7 km/s/pc. In two systems, we detect high-velocity HCO+ (J=1-0) emission inside the highest-velocity \nthp\ emission. This enables us to study the infall and rotation close to the disk and estimate the central object masses. The velocity fields observed on large and small-scales are more complex than would be expected from rotation alone, suggesting that complex envelope structure enables other dynamical processes (i.e. infall) to affect the velocity field.

Using the Expanded Very Large Array, we have mapped the vicinity of the Orion BN/KL Hot Core with sub-arcsecond angular resolution in seven metastable inversion transitions of ammonia: (J,K)=(6,6) to (12,12). This emission comes from levels up to 1500 K above the ground state, enabling identification of source(s) responsible for heating the region. We used this multi-transition dataset to produce images of the rotational/kinetic temperature and the column density of ammonia for ortho and para species separately and on a position-by-position basis. We find rotational temperature and column density in the range 160-490 K and (1-4)x10^17 cm^-2, respectively. Our spatially-resolved images show that the highest (column) density and hottest gas is found in a northeast-southwest elongated ridge to the southeast of Source I. We have also measured the ortho-para ratio of ammonia, estimated to vary in the range 0.9-1.6. Enhancement of ortho with respect to para and the offset of hot ammonia emission peaks from known (proto)stellar sources provide evidence that the ammonia molecules have been released from dust grains into the gas-phase through the passage of shocks and not by stellar radiation. We propose that the combined effect of Source I's proper motion and its low-velocity outflow impinging on a pre-existing dense medium is responsible for the excitation of ammonia and the Orion Hot Core. Finally, we found for the first time evidence of a slow (5 km/s) and compact (1000 AU) outflow towards IRc7.

We have performed a 1.3 centimeter survey of 24 massive young stellar objects (MYSOs) using the Expanded Very Large Array (EVLA). The sources in the sample exhibit a broad range of massive star formation signposts including Infrared Dark Clouds (IRDCs), UCHII regions, and extended 4.5 micron emission in the form of Extended Green Objects (EGOs). In this work, we present results for G35.03+0.35 which exhibits all of these phenomena. We simultaneously image the 1.3 cm ammonia (1,1) through (6,6) inversion lines, four methanol transitions, two H recombination lines, plus continuum at 0.05 pc resolution. We find three areas of thermal ammonia emission, two within the EGO (designated the NE and SW cores) and one toward an adjacent IRDC. The NE core contains an UCHII region (CM1) and a candidate HCHII region (CM2). A region of non-thermal, likely masing ammonia (3,3) and (6,6) emission is coincident with an arc of 44 GHz methanol masers. We also detect two new 25 GHz Class I methanol masers. A complementary Submillimeter Array 1.3 mm continuum image shows that the distribution of dust emission is similar to the lower-lying ammonia lines, all peaking to the NW of CM2, indicating the likely presence of an additional MYSO in this protocluster. By modeling the ammonia and 1.3 mm continuum data, we obtain gas temperatures of 20-220 K and masses of 20-130 solar. The diversity of continuum emission properties and gas temperatures suggest that objects in a range of evolutionary states exist concurrently in this protocluster.

Since its commissioning in 1980, the Very Large Array (VLA) has consistently demonstrated its scientific productivity. However, its fundamental capabilities have changed little since 1980, particularly in the key areas of sensitivity, frequency coverage, and velocity resolution. These limitations have been addressed by a major upgrade of the array, which began in 2001 and will be completed at the end of 2012. When completed, the Expanded VLA -- the EVLA -- will provide complete frequency coverage from 1 to 50 GHz, a continuum sensitivity of typically 1 microJy/beam (in 9 hours with full bandwidth), and a modern correlator with vastly greater capabilities and flexibility than the VLA's. In this paper we describe the goals of the EVLA project, its current status, and the anticipated expansion of capabilities over the next few years. User access to the array through the OSRO and RSRO programs is described. The following papers in this special issue, derived from observations in its early science period, demonstrate the astonishing breadth of this most flexible and powerful general-purpose telescope.

We present Very Large Array observations of ammonia (NH3) (1,1), (2,2), and CCS (2_1-1_0) emission toward the Infrared Dark Cloud (IRDC) G19.30+0.07 at ~22GHz. The NH3 emission closely follows the 8 micron extinction. The NH3 (1,1) and (2,2) lines provide diagnostics of the temperature and density structure within the IRDC, with typical rotation temperatures of ~10 to 20K and NH3 column densities of ~10^15 cm^-2. The estimated total mass of G19.30+0.07 is ~1130 Msun. The cloud comprises four compact NH3 clumps of mass ~30 to 160 Msun. Two coincide with 24 micron emission, indicating heating by protostars, and show evidence of outflow in the NH3 emission. We report a water maser associated with a third clump; the fourth clump is apparently starless. A non-detection of 8.4GHz emission suggests that the IRDC contains no bright HII regions, and places a limit on the spectral type of an embedded ZAMS star to early-B or later. From the NH3 emission we find G19.30+0.07 is composed of three distinct velocity components, or "subclouds." One velocity component contains the two 24 micron sources and the starless clump, another contains the clump with the water maser, while the third velocity component is diffuse, with no significant high-density peaks. The spatial distribution of NH3 and CCS emission from G19.30+0.07 is highly anti-correlated, with the NH3 predominantly in the high-density clumps, and the CCS tracing lower-density envelopes around those clumps. This spatial distribution is consistent with theories of evolution for chemically young low-mass cores, in which CCS has not yet been processed to other species and/or depleted in high-density regions.

We present new 7mm continuum observations of Orion BN/KL with the VLA. We resolve the emission from the protostar radio Source I and BN at several epochs. Source I is highly elongated NW-SE, and remarkably stable in flux density, position angle, and overall morphology over nearly a decade. This favors the extended emission component arising from an ionized disk rather than a jet. We have measured the proper motions of Source I and BN for the first time at 43 GHz. We confirm that both sources are moving at high speed (12 and 26 km/s, respectively) approximately in opposite directions, as previously inferred from measurements at lower frequencies. We discuss dynamical scenarios that can explain the large motions of both BN and Source I and the presence of disks around both. Our new measurements support the hypothesis that a close (~50 AU) dynamical interaction occurred around 500 years ago between Source I and BN as proposed by Gomez et al. From the dynamics of encounter we argue that Source I today is likely to be a binary with a total mass on the order of 20 Msun, and that it probably existed as a softer binary before the close encounter. This enables preservation of the original accretion disk, though truncated to its present radius of ~50 AU. N-body numerical simulations show that the dynamical interaction between a binary of 20 Msun total mass (I) and a single star of 10 Msun mass (BN) may lead to the ejection of both and binary hardening. The gravitational energy released in the process would be large enough to power the wide-angle flow traced by H2 and CO emission in the BN/KL nebula. Assuming the proposed dynamical history is correct, the smaller mass for Source I recently estimated from SiO maser dynamics (>7 Msun) by Matthews et al., suggests that non-gravitational forces (e.g. magnetic) must play an important role in the circumstellar gas dynamics.

We present and analyze two new high-resolution (approx 0.3 arcsec), high-sensitivity (approx 50 uJy beam-1) Very Large Array 3.6 cm observations of IRAS 16293-2422 obtained in 2007 August and 2008 December. The components A2alpha and A2beta recently detected in this system are still present, and have moved roughly symmetrically away from source A2 at a projected velocity of 30-80 km s-1. This confirms that A2alpha and A2beta were formed as a consequence of a very recent bipolar ejection from A2. Powerful bipolar ejections have long been known to occur in low-mass young stars, but this is -to our knowledge-- the first time that such a dramatic one is observed from its very beginning. Under the reasonable assumption that the flux detected at radio wavelengths is optically thin free-free emission, one can estimate the mass of each ejecta to be of the order of 10^-8 Msun. If the ejecta were created as a consequence of an episode of enhanced mass loss accompanied by an increase in accretion onto the protostar, then the total luminosity of IRAS 16293-2422 ought to have increased by 10-60% over the course of at least several months. Between A2alpha and A2beta, component A2 has reappeared, and the relative position angle between A2 and A1 is found to have increased significantly since 2003-2005. This strongly suggests that A1 is a protostar rather than a shock feature, and that the A1/A2 pair is a tight binary system. Including component B, IRAS 16293-2422 therefore appears to be a very young hierarchical multiple system.

We present multi-epoch VLBA imaging of the 28SiO v=1 & v=2, J=1-0 maser emission toward the massive YSO Orion Source I. Both SiO transitions were observed simultaneously with an angular resolution of ~0.5 mas (~0.2 AU for d=414 pc). Here we explore the global properties and kinematics of the emission through two 19-epoch animated movies spanning 21 months (2001 March 19 to 2002 December 10). These movies provide the most detailed view to date of the dynamics and temporal evolution of molecular material within ~20-100 AU of a massive (~>8M_sun) YSO. The bulk of the SiO masers surrounding Source I lie in an X-shaped locus; emission in the South/East arms is predominantly blueshifted and in the North and West is predominantly redshifted. In addition, bridges of intermediate-velocity emission connect the red and blue sides of the emission distribution. We have measured proper motions of >1000 maser features and find a combination of radially outward migrations along the four arms and motions tangent to the bridges. We interpret the SiO masers as arising from a wide-angle bipolar wind emanating from a rotating, edge-on disk. The detection of maser features along extended, curved filaments suggests that magnetic fields may play a role in launching and/or shaping the wind. Our observations appear to support a picture in which stars with M ~>8 M_sun form via disk-mediated accretion. However, we cannot rule out that the Source I disk may have been formed or altered following a close encounter. (Abridged).

We present high angular resolution PdBI images (beam of ~0.33'') of the J=27-26 line from several vibrational levels (v7=1 and v6=1) of HC3N toward Cepheus A HW2. These images reveal the two main heating sources in the cluster: one centered in the disk collimating the HW2 radio jet (the HW2 disk), and the other associated with a hot core 0.3'' northeast HW2 (the HC). This is the first time that vibrationally excited emission of HC3N is spatially resolved in a disk. The kinematics of this emission shows that the HW2 disk rotates following a Keplerian law. We derive the temperature profiles in the two objects from the excitation of HC3N along the HW2 disk and the HC. These profiles reveal that both objects are centrally heated and show temperature gradients. The inner and hotter regions have temperatures of 350+-30K and 270+-20K for the HW2 disk and the HC, respectively. In the cooler and outer regions, the temperature drops to 250+-30K in the HW2 disk, and to 220+-15K in the HC. The estimated luminosity of the heating source of the HW2 disk is ~2.2E4Lo, and 3000Lo for the HC. The most massive protostar in the HW2 region is the powering source of the HW2 radio jet. We discuss the formation of multiple systems in this cluster. The proximity of the HC to HW2 suggest that these sources likely form a binary system of B stars, explaining the observed precession of the HW2 radio jet.

We have used the Very Large Array (VLA) at 7 mm wavelength to image five rotational transitions (J=1-0) from three SiO isotopologues towards Orion BN/KL: 28SiO v=0,1,2; and 29SiO and 30SiO v=0. For the first time, we have mapped the 29SiO and 30SiO J=1-0 emission, established the maser nature of the emission, and confirmed association with the deeply embedded high-mass young stellar object commonly denoted radio Source I. The 28SiO v=0 maser emission shows a bipolar structure that extends over ~700 AU along a northeast-southwest axis, and we propose that it traces a bipolar outflow driven by Source I. The high-brightness isotopic SiO maser emission imaged with a ~0.2 arcsec resolution has a more compact distribution, generally similar to that of the 28SiO v=1,2 emission, and it probably traces bulk gas flows in a region of diameter <100 AU centered on Source I. On small scales of <10 AU, however, compact 29SiO/30SiO v=0 and 28SiO v=1,2 emission features may be offset from one another in position and line-of-sight velocity. From a radiative transfer analysis based on a large velocity gradient (LVG) pumping model, we derive similar temperatures and densities for the optimum excitation of both 29SiO/30SiO v=0 and 28SiO v=1,2 masers, significantly higher than required for 28SiO v=0 maser excitation. In order to account for the small-scale differences among the isotopologues (v=0) and the main species (v=1,2), follow-up radiative transfer modeling that incorporates non-local line overlap among transitions of all SiO isotopic species may be required.

We have surveyed molecular line emission from Orion BN/KL from 42.3 to 43.6 GHz with the Green Bank Telescope. Sixty-seven lines were identified and ascribed to 13 different molecular species. The spectrum at 7 mm is dominated by SiO, SO2, CH3OCH3, and C2H5CN. Five transitions have been detected from the SiO isotopologues 28SiO, 29SiO, and 30SiO. We report here for the first time the spectra of the 29SiO and 30SiO v=0 J=1-0 emission in Orion BN/KL, and we show that they have double-peaked profiles with velocity extents similar to the main isotopologue. The main motivation for the survey was the search of high-velocity (100-1000 km/s) outflows in the BN/KL region as traced by SiO Doppler components. Some of the unidentified lines in principle could be high-velocity SiO features, but without imaging data their location cannot be established. Wings of emission are present in the v=0 28SiO, 29SiO and 30SiO profiles, and we suggest that the v=0 emission from the three isotopologues might trace a moderately high-velocity (~30-50 km/s) component of the flows around the high-mass protostar Source I in the Orion BN/KL region. We also confirm the 7 mm detection of a complex oxygen-bearing species, acetone (CH3COCH3), which has been recently observed towards the hot core at 3 mm, and we have found further indications of the presence of long cyanopolyynes (HC5N and HC7N) in the quiescent cold gas of the extended ridge.

Kant and Laplace suggested the Solar System formed from a rotating gaseous disk in the 18th century, but convincing evidence that young stars are indeed surrounded by such disks was not presented for another 200 years. As we move into the 21st century the emphasis is now on disk formation, the role of disks in star formation, and on how planets form in those disks. Radio wavelengths play a key role in these studies, currently providing some of the highest spatial resolution images of disks, along with evidence of the growth of dust grains into planetesimals. The future capabilities of EVLA and ALMA provide extremely exciting prospects for resolving disk structure and kinematics, studying disk chemistry, directly detecting proto-planets, and imaging disks in formation.

In this article, we report high-resolution (~ 0.1" -- 0.3"), high-sensitivity (~ 50 -100 uJy beam-1) Very Large Array 0.7 and 1.3 cm observations of the young stellar system IRAS 16293-2422 in rho-Ophiuchus. In the 0.7 cm image, component A to the south-east of the system looks like its usual binary self. In the new 1.3 cm image, however, component A2 appears to have split into two sub-components located roughly symmetrically around the original position of A2. This change of morphology is likely the result of a recent bipolar ejection, one of the very first such events observed in a low-mass source. Also in component A, a marginal detection of 0.7 cm emission associated with the submillimeter component Ab is reported. If confirmed, this detection would imply that Ab is a relatively extended dusty structure, where grain coagulation may already have taken place. With an angular size increasing with frequency, and an overall spectra index of 2, the emission from component B to the north-west of the system is confirmed to be dominated by optically thick thermal dust emission associated with a fairly massive, nearly face-on, circumstellar disk. In the central region, however, we find evidence for a modest free-free contribution that originates in a structure elongated roughly in the east-west direction. We argue that this free-free component traces the base of the jet driving the large-scale bipolar flow at a position angle of about 110 degrees that has long been known to be powered by IRAS 16293-2422.

A comprehensive picture of high-mass star formation has remained elusive, in part because examples of high-mass YSOs tend to be relatively distant, deeply embedded, and confused with other emission sources. These factors have impeded dynamical investigations within tens of AU of high-mass YSOs--scales that are critical for probing the interfaces where outflows from accretion disks are launched and collimated. Using observations of SiO masers obtained with the VLA and the VLBA, the KaLYPSO project is overcoming these limitations by mapping the structure and dynamical/temporal evolution of the material 10-1000 AU from the nearest high-mass YSO: Radio Source I in the Orion BN/KL region. Our data include ~40 epochs of VLBA observations over a several-year period, allowing us to track the proper motions of individual SiO maser spots and to monitor changes in the physical conditions of the emitting material with time. Ultimately these data will provide 3-D maps of the outflow structure over approximately 30% of the outflow crossing time. Here we summarize recent results from the KaLYPSO project, including evidence that high-mass star formation is occurring via disk-mediated accretion.

We present VLA and PdBI subarcsecond images (0.15"-0.6") of the radiocontinuum emission at 7 mm and of the SO2 J=19_{2,18}-18_{3,15} and J=27_{8,20}-28_{7,21} lines toward the Cepheus A HW2 region. The SO2 images reveal the presence of a hot core internally heated by an intermediate mass protostar, and a circumstellar rotating disk around the HW2 radio jet with size 600AUx100AU and mass of 1M_sun. Keplerian rotation for the disk velocity gradient of 5 kms-1 requires a 9 M_sun central star, which cannot explain the total luminosity observed in the region. This may indicate that the disk does not rotate with a Keplerian law due to the extreme youth of this object. Our high sensitivity radiocontinuum image at 7 mm shows in addition to the ionized jet, an extended emission to the west (and marginally to the south) of the HW2 jet, filling the south-west cavity of the HW2 disk. From the morphology and location of this free-free continuum emission at centimeter and millimeter wavelengths (spectral index of 0.4-1.5), we propose that the disk is photoevaporating due to the UV radiation from the central star. All this indicates that the Cepheus A HW2 region harbors a cluster of massive stars. Disk accretion seems to be the most plausible way to form massive stars in moderate density/luminosity clusters.

We have imaged the enigmatic radio source-I (Orion-I) in the Orion-KL nebula with the VLA at 43 GHz with 34 mas angular resolution. The continuum emission is highly elongated and is consistent with that expected from a nearly edge-on disk. The high brightness and lack of strong molecular lines from Orion-I can be used to argue against emission from dust. Collisional ionization and H-minus free-free opacity, as in Mira variables, require a central star with >10^5 Lsun, which is greater than infrared observations allow. However, if significant local heating associated with accretion occurs, lower total luminosities are possible. Alternatively, photo-ionization from an early B-type star and p+/e- bremsstrahlung can explain our observations, and Orion-I may be an example of ionized accretion disk surrounding a forming massive star. Such accretion disks may not be able to form planets efficiently.

We describe phase-retrieval holography measurements of the 100-m diameter Green Bank Telescope using astronomical sources and an astronomical receiver operating at a wavelength of 7 mm. We use the technique with parameterization of the aperture in terms of Zernike polynomials and employing a large defocus, as described by Nikolic, Hills & Richer (2006). Individual measurements take around 25 minutes and from the resulting beam maps (which have peak signal to noise ratios of 200:1) we show that it is possible to produce low-resolution maps of the wavefront errors with accuracy around a hundredth of a wavelength. Using such measurements over a wide range of elevations, we have calculated a model for the wavefront-errors due to the uncompensated gravitational deformation of the telescope. This model produces a significant improvement at low elevations, where these errors are expected to be the largest; after applying the model, the aperture efficiency is largely independent of elevation. We have also demonstrated that the technique can be used to measure and largely correct for thermal deformations of the antenna, which often exceed the uncompensated gravitational deformations during daytime observing. We conclude that the aberrations induced by gravity and thermal effects are large-scale and the technique used here is particularly suitable for measuring such deformations in large millimetre wave radio telescopes.

We present results from SMA observations of the star forming region Cepheus A-East at ~340 GHz (875 micron) with 0.7'' - 2'' resolution. At least four compact submm continuum sources have been detected, as well as a rich forest of hot core line emission. Two kinematically, chemically, and thermally distinct regions of molecular emission are present in the vicinity of the HW2 thermal jet, both spatially distinct from the submm counterpart to HW2. We propose that this emission is indicative of multiple protostars rather than a massive disk as reported by Patel et al. (2005).

9C J1503+4528 is a very young CSS radio galaxy, with an age of order 10^4 years. This source is therefore an ideal laboratory for the study ofthe intrinsic host galaxy/IGM properties, interactions between the radio source and surrounding ISM, links between star formation and AGN activity and the radio source triggering mechanism. Here we present the results of a spectroscopic analysis of this source, considering each of these aspects of radio source physics. We find that shock ionization by the young radio source is important in the central regions of the galaxy on scales similar to that of the radio source itself, whilst evidence for an AGN ionization cone is observed at greater distances. Line and continuum features require the presence of a young stellar population, the best-fit model for which implies an age of 5x10^6 years, significantly older than the radio source. Most interestingly, the relative sizes of radio source and extended emission line region suggest that both AGN and radio source are triggered at approximately the same time. If both the triggering of the radio source activity and the formation of the young stellar population had the same underlying cause, this source provides a sequence for the events surrounding the triggering process. We propose that the AGN activity in 9C J1503+4528 was causedby a relatively minor interaction, and that a super-massive black hole powering the radio jets must have been in place before the AGN was triggered.

We report on the extreme behaviour of the high redshift blazar GB B1428+4217 at z=4.72. A continued programme of radio measurements has revealed an exceptional flare in the lightcurve, with the 15.2 GHz flux density rising by a factor ~3 from ~140 mJy to ~430 mJy in a rest-frame timescale of only ~4 months -- much larger than any previous flares observed in this source. In addition to new measurements of the 1.4-43 GHz radio spectrum we also present the analysis and results of a target-of-opportunity X-ray observation using XMM-Newton, made close to the peak in radio flux. Although the X-ray data do not show a flare in the high energy lightcurve, we are able to confirm the X-ray spectral variability hinted at in previous observations. GB B1428+4217 is one of several high-redshift radio-loud quasars that display a low energy break in the X-ray spectrum, probably due to the presence of excess absorption in the source. X-ray spectral analysis of the latest XMM-Newton data is shown to be consistent with the warm absorption scenario which we have hypothesized previously. Warm absorption is also consistent with the observed X-ray spectral variability of the source, in which the spectral changes can be successfully accounted-for with a fixed column density of material in which the ionization state is correlated with hardness of the underlying power-law emission.

We present multi-frequency, multi-epoch radio imaging of the complex radio source B2151+174 in the core of the cluster, Abell 2390 (z~0.23). From new and literature data we conclude that the FRII-powerful radio source is the combination of a compact, core-dominated `medium-symmetric object' (MSO) with a more extended, steeper spectrum mini-halo. B2151+174 is unusual in a number of important aspects: i) it is one of the most compact and flat spectrum sources in a cluster core known; ii) it shows a complex, compact twin-jet structure in a north-south orientation; iii) the orientation of the jets is 45deg misaligned with apparent structure (ionization cones and dust disk) of the host galaxy on larger scales. Since the twin-jet of the MSO has its northern half with an apparent `twist', it might be that precession of the central supermassive black hole explains this misalignement. B2151+174 may be an example of the early stage (10^3-10^4 yrs duration) of a `bubble' being blown into the ICM where the plasma has yet to expand.

We present 7-mm continuum observations of 14 low-mass pre-main-sequence stars in the Taurus-Auriga star-forming region obtained with the Very Large Array with ~1.5" resolution and ~0.3 mJy rms sensitivity. For 10 objects, the circumstellar emission has been spatially resolved. The large outer disk radii derived suggest that the emission at this wavelength is mostly optically thin. The millimetre spectral energy distributions are characterised by spectral indices alpha = 2.3 to 3.2. After accounting for contribution from free-free emission and corrections for optical depth, we determine dust opacity indices beta in the range 0.5 to 1.6, which suggest that millimetre-sized dust aggregates are present in the circumstellar disks. Four of the sources with beta > 1 may be consistent with submicron-sized dust as found in the interstellar medium. Our findings indicate that dust grain growth to millimetre-sized particles is completed within less than 1 Myr for the majority of circumstellar disks.

We present observations of IRAS 16293-2422 using the Submillimeter Array (SMA) at 300 GHz, and the Very Large Array (VLA) at frequencies from 1.5 to 43 GHz. This source was already known to be a binary system with components A and B separated by ~5". The new SMA data now separate source A into two submm continuum components, denoted Aa and Ab. The strongest of these, Aa, peaks between the cm radio sources A1 and A2. Archival VLA data show proper motion of sources A and B of 17 mas/yr, associated with motion of the rho Ophiuchi cloud. We also find relative motion between A1 and A2 which excludes the possibility that they are gravitationally bound. A2 remains stationary relative to source B, and we identify it as the protostar driving the NE-SW CO outflow. A1 is shock-ionized gas where a precessing jet interacts with dense gas. This jet probably drives the E-W CO outflow. The origin of this jet must be located close to A2, and may be the submm source Aa. Thus source A is now shown to comprise three components within 1". Source B, however, is single, exhibits optically-thick dust emission even at 8 GHz, has a high luminosity, and yet shows no sign of outflow, indicating it is probably very young. The SMA spectrum of IRAS 16293-2422 reports the first astronomical identification of many lines at 300 and 310 GHz. The species detected are typical of hot cores, and are mainly associated with source A. We suggest the shock at A1 may explain the enhanced abundances measured, especially of sulphur-bearing molecules. We use the high brightness temperature of the submm dust emission from source B to demonstrate the unambiguous detection of infall by observing redshifted SO (7_7-6_6) absorption against the dust disk.

We present a complete survey of current star formation in the Perseus molecular cloud, made at 850 and 450 micron with SCUBA at the JCMT. Covering 3 deg^2, this submillimetre continuum survey for protostellar activity is second in size only to that of rho Ophiuchus (Johnstone et al. 2004). Complete above 0.4 msun (5 sigma detection in a 14'' beam), we detect a total of 91 protostars and prestellar cores. Of these, 80% lie in clusters, representative of star formation across the Galaxy. Two of the groups of cores are associated with the young stellar clusters IC348 and NGC1333, and are consistent with a steady or reduced star formation rate in the last 0.5 Myr, but not an increasing one. In Perseus, 40--60% of cores are in small clusters (< 50 msun) and isolated objects, much more than the 10% suggested from infrared studies. Complementing the dust continuum, we present a C^18O map of the whole cloud at 1' resolution. The gas and dust show filamentary structure of the dense gas on large and small scales, with the high column density filaments breaking up into clusters of cores. The filament mass per unit length is 5--11 msun per 0.1 pc. Given these filament masses, there is no requirement for substantial large scale flows along or onto the filaments in order to gather sufficient material for star formation. We find that the probability of finding a submillimetre core is a strongly increasing function of column density, as measured by C^18O integrated intensity, prob(core) proportional to I^3.0. This power law relation holds down to low column density, suggesting that there is no A_v threshold for star formation in Perseus, unless all the low-A_v submm cores can be demonstrated to be older protostars which have begun to lose their natal molecular cloud.

In this work we present VLA observations of the NH3, CCS, and H2O maser emission at 1 cm from the star forming region B1-IRS (IRAS 03301+3057) with ~ 5" (=1750 AU) of angular resolution. CCS emission is distributed in three clumps around the central source. These clumps exhibit a velocity gradient from red- to blueshifted velocities toward B1-IRS, probably due to an interaction with the outflow from an embedded protostar. The outflow and its powering source are traced by a reflection nebula and an associated infrared point source detected in a 2MASS K-band image. We find that this infrared point source is associated with water maser emission distributed in an elongated structure (~ 450 AU size) along the major axis of the reflection nebula and tracing the base of the outflow of the region. Ammonia emission is extended and spatially anticorrelated with CCS. This is the first time that this kind of anticorrelation is observed in a star forming region with such a high angular resolution, and illustrates the importance of time-dependent chemistry on small spatial scales. The relatively large abundance of CCS with respect to ammonia, compared with other star forming regions, suggests an extreme youth for the B1-IRS object (<10E+5 yr). We suggest the possibility that CCS abundance is enhanced via shock-induced chemistry.

We present a search for HCN emission from four high redshift far infrared (IR) luminous galaxies. Current data and models suggest that these high $z$ IR luminous galaxies represent a major starburst phase in the formation of spheroidal galaxies, although many of the sources also host luminous active galactic nuclei (AGN), such that a contribution to the dust heating by the AGN cannot be precluded. HCN emission is a star formation indicator, tracing dense molecular hydrogen gas within star-forming molecular clouds (n(H$_2$) $\sim 10^5$ cm$^{-3}$). HCN luminosity is linearly correlated with IR luminosity for low redshift galaxies, unlike CO emission which can also trace gas at much lower density. We report a marginal detection of HCN (1-0) emission from the $z=2.5832$ QSO J1409+5628, with a velocity integrated line luminosity of $L_{\rm HCN}'=6.7\pm2.2 \times10^{9}$ K km s$^{-1}$ pc$^2$, while we obtain 3$\sigma$ upper limits to the HCN luminosity of the $z=3.200$ QSO J0751+2716 of $L_{\rm HCN}'=1.0\times10^{9}$ K km s$^{-1}$ pc$^2$, $L_{\rm HCN}'=1.6\times10^{9}$ K km s$^{-1}$ pc$^2$ for the $z= 2.565$ starburst galaxy J1401+0252, and $L_{\rm HCN}'=1.0\times10^{10}$ K km s$^{-1}$ pc$^2$ for the $z = 6.42$ QSO J1148+5251. We compare the HCN data on these sources, plus three other high-$z$ IR luminous galaxies, to observations of lower redshift star-forming galaxies. The values of the HCN/far-IR luminosity ratios (or limits) for all the high $z$ sources are within the scatter of the relationship between HCN and far-IR emission for low $z$ star-forming galaxies (truncated).

We present a study of the structure and dynamics of the star forming region B1-IRS (IRAS 03301+3057) using the properties of different molecules at high angular resolution (~4''). We have used VLA observations of NH3, CCS, and H2O masers at 1 cm. CCS emission shows three clumps around the central source, with a velocity gradient from red to blueshifted velocities towards the protostar, probably due to the interaction with outflowing material. Water maser emission is elongated in the same direction as a reflection nebula detected at 2micron by 2MASS, with the maser spots located in a structure of some hundreds of AU from the central source, possibly tracing a jet. We propose a new outflow model to explain all our observations, consisting of a molecular outflow near the plane of the sky. Ammonia emission is extended and anticorrelated with CCS. We have detected for the first time this anticorrelation at small scales (1400 AU) in a star forming region.

We have carried out extensive radio and optical follow-up of 176 sources from the 15 GHz 9th Cambridge survey. Optical identifications have been found for 155 of the radio sources; optical images are given with radio maps overlaid. The continuum radio spectrum of each source spanning the frequency range 1.4 - 43 GHz is also given. Two flux-limited samples are defined, one containing 124 sources complete to 25 mJy and one of 70 sources complete to 60 mJy. Between one fifth and one quarter of sources from these flux-limited samples display convex radio spectra, rising between 1.4 and 4.8 GHz. These rising-spectrum sources make up a much larger fraction of the radio source population at this high selection frequency than in lower frequency surveys. We find that by using non-simultaneous survey flux density measurements at 1.4 and 15 GHz to remove steep spectrum objects, the efficiency of selecting objects with spectra rising between 1.4 and 4.8 GHz (as seen in simultaneous measurements) can be raised to 49 percent without compromising the completeness of the rising spectrum sample.

Using the Very Large Array (VLA) in the A-configuration, we have obtained a high-resolution 3.6 cm map of a hard X-ray source detected by the Chandra X-ray Observatory in a protostellar clump in Orion molecular cloud 3. Two radio continuum sources were detected in the vicinity of the X-ray source, both of which have NIR counterparts. We conclude that these VLA sources are free-free emission produced by shocks in protostellar jets from the NIR class I protostars. Using the centimeter data, we determined the power and orientation of the protostellar jets. The center position of the X-ray emission was found to be ~1--2" offset from the exciting sources of the jets, and the displacement is in the direction of the jets and molecular outflows. We discuss the nature of the X-ray emission as the shock-excited plasma at the shock front where the jet propagates through interstellar medium at a speed of ~1000 km/s.

We report on the results of a search using the VLA for redshifted HCN(1-0) emission from the host galaxy of BR B1202-0725, an optically luminous quasar at z=4.695. The host galaxy emits strongly in the rest-frame far-infrared, and shows characteristics very similar to that of more local, ultraluminous infrared galaxies, in which a significant fraction of the far-infrared emission is powered by star formation. We find a 3-sigma upper limit to the HCN(1-0) emission of 4.9 x 10^10 K kms^-1 pc^2, assuming a lambda-cosmology. This limit is consistent with correlations derived from measurements of the HCN, CO, and far-infrared emission for a sample of more local galaxies, including starbursts (Solomon et al, 1992a).

We have observed the SiO J=5--4 line towards a sample of 25 low-luminosity (L* < 10^3 Lsun) protostellar outflow systems. The line was detected towards 7 of the 25 sources, a detection rate of 28 per cent. The majority (5 out of 7) of sources detected were of class 0 status. We detected a higher fraction of class 0 sources compared with the class I and II sources, although given the small numbers involved the significance of this result should be regarded as tentative. Most of the detected sources showed emission either at or close to the central position, coincident with the protostar. In four cases (HH211, HH25MMS, V-380 OriNE and HH212) emission was also detected at positions away from the center, and was stronger than that observed at the centre position. SiO abundances of 10^(-8) to 8x10^(-7) are derived from LTE analysis. For 2 sources we have additional transitions which we use to conduct statistical equilibrium modeling to estimate the gas density in the SiO-emitting regions. For HH25MMS these results suggest that the SiO emission arises in a higher-density region than the methanol previously observed. We find that the most likely explanation for the preferential detection of SiO emission towards class 0 sources is the greater density of those environments, reinforced by higher shock velocities. We conclude that while not all class 0 sources exhibit SiO emission, SiO emission is a good signpost for the presence of class 0 sources.

High-mass star formation is not well understood chiefly because examples are deeply embedded, relatively distant, and crowded with sources of emission. Using VLA and VLBA observations of water and SiO maser emission, we have mapped in detail the structure and proper motion of material 20-500 AU from the closest high-mass YSO, radio source-I in the Orion KL region. We observe streams of material driven in a rotating, wide angle, bipolar wind from the surface of an edge-on accretion disk. The example of source-I provides strong evidence that high-mass star formation proceeds via accretion

We have carried out a survey of optically-selected dark clouds using the bolometer array SCUBA on the James Clerk Maxwell Telescope, at 850 microns wavelength. The survey covers a total of 0.5 square degrees and is unbiased with reference to cloud size, star formation activity, or the presence of infrared emission. Several new protostars and starless cores have been discovered; the protostars are confirmed through the detection of their accompanying outflows in CO(2-1) emission. The survey is believed to be complete for Class 0 and Class I protostars, and yields two important results regarding the lifetimes of these phases. First, the ratio of Class 0 to Class protostars in the sample is roughly unity, very different from the 1:10 ratio that has previously been observed for the rho Ophiuchi star-forming region. Assuming star formation to be a homogeneous process in the dark clouds, this implies that the Class 0 lifetime is similar to the Class I phase, which from infrared surveys has been established to be approximately 200,000 yr. It also suggests there is no rapid initial accretion phase in Class 0 objects. A burst of triggered star formation some 100,000 yr ago can explain the earlier results for rho Ophiuchus. Second, the number of starless cores is approximately twice that of the total number of protostars, indicating a starless core lifetime of approximately 800,000 yr. These starless cores are therefore very short-lived, surviving only two or three free-fall times. This result suggests that, on size scales of 10,000 AU at least, the dynamical evolution of starless cores is probably not controlled by magnetic processes.

We present images of the HH 211 molecular jet in the SiO v=0, J=1-0 line at 43 GHz made with the Very Large Array at approximately 0.5 arcsec resolution. The SiO emission appears to trace primarily internal bowshocks in the outflow, suggesting that the dust and molecular gas are accelerated via prompt entrainment at internal working surfaces in the jet. There is also some evidence for limb-brightening of the SiO emission, indicating that SiO emission may also arise from entrainment in the jet's boundary layer. Excitation temperatures of >~150-200 K are inferred from the SiO emission. Enhancements in the SiO abundance of ~10^6 over interstellar values are observed, and the possible origin of the SiO is discussed.

High dynamic range imaging of submillimeter dust emission from the envelopes of eight young protostars in the Taurus and Perseus star-forming regions has been carried out using the SCUBA submillimeter camera on the James Clerk Maxwell Telescope. Good correspondence between the spectral classifications of the protostars and the spatial distributions of their dust emission is observed, in the sense that those with cooler spectral energy distributions also have a larger fraction of the submillimeter flux originating in an extended envelope compared with a disk. This results from the cool sources having more massive envelopes rather than warm sources having larger disks. Azimuthally-averaged radial profiles of the dust emission are used to derive the power-law index of the envelope density distributions, p (defined by rho proportional to r^-p), and most of the sources are found to have values of p consistent with those predicted by models of cloud collapse. However, the youngest protostars in our sample, L1527 and HH211-mm, deviate significantly from the theoretical predictions, exhibiting values of p somewhat lower than can be accounted for by existing models. For L1527 heating of the envelope by shocks where the outflow impinges on the surrounding medium may explain our result. For HH211-mm another explanation is needed, and one possibility is that a shallow density profile is being maintained in the outer envelope by magnetic fields and/or turbulence. If this is the case star formation must be determined by the rate at which the support is lost from the cloud, rather than the hydrodynamical properties of the envelope, such as the sound speed.

We have detected the T~Tauri star, DO Tauri, in a 0.6$''$-resolution VLA map of 43.3 GHz ($\lambda$ = 7 mm) continuum emission. The 43 GHz flux density lies on the same power-law slope defined by 89 to 232 GHz measurements, F$_\nu$ $\propto \nu^{\alpha}$ with index $\alpha$ = 2.39$\pm$0.23, confirming that the 43.3 GHz emission is thermal radiation from circumstellar dust. Upper limits to the flux densities at 8.4 and 22.5 GHz constrain the contribution of free-free emission from a compact ionized wind to less than 49\%. The dust emissivity index, $\beta$, is $0.39\pm$0.23, if the emission is optically thin. Fitting a model of a thin circumstellar disk to the observed spectral energy distribution gives $\beta = 0.6\pm0.3$, consistent with the power-law derivation. Both values are substantially lower than is generally accepted for the interstellar medium, suggesting grain growth. Given the youth of DO Tau and the early evolutionary state of its circumstellar disk, this result implies that mm-size grains have already formed by the early T-Tauri phase.