We report results of 0.05"-resolution observations toward the O-type proto-binary system IRAS 16547-4247 with the Atacama Large Millimeter/submillimeter Array (ALMA). We present dynamical and chemical structures of the circumbinary disk, circumstellar disks, outflows and jets, illustrated by multi-wavelength continuum and various molecular lines. In particular, we detect sodium chloride, silicon compounds, and vibrationally-excited water lines as probes of the individual protostellar disks at a scale of 100 au. These are complementary to typical hot-core molecules tracing the circumbinary structures on a 1000-au scale. The H2O line tracing inner-disks has an upper-state energy of Eu/k>3000K, indicating a high temperature of the disks. On the other hand, despite the detected transitions of NaCl, SiO, and SiS not necessarily having high upper-state energies, they are enhanced only in the vicinity of the protostars. We interpret that these molecules are the products of dust destruction, which only happens in the inner disks. This is the second detection of alkali metal halide in protostellar systems after the case of the disk of Orion Source I, and also one of few massive protostellar disks associated with high-energy transition water and silicon compounds. These new results suggest these "hot-disk" lines may be common in innermost disks around massive protostars, and have great potential for future research of massive star formation. We also tentatively find that the twin disks are counter-rotating, which might give a hint of the origin of the massive proto-binary system IRAS 16547-4247.

We investigate the imaging performance of an interferometric array in the case of wide field, high resolution, narrow band, snapshot imaging. We find that, when uv-cell sizes are sufficiently small (ie. image sizes are sufficiently large), each instantaneous visibility record is gridded into its own uv-cell. This holds even for dense arrays, like the core of the next generation VLA. In this particular, application, Uniform weighting of the gridded visibilities approaches Natural weighting, with its often deleterious consequences on the resulting synthesized beam. For a core-dominated array, we show that the resulting image noise is highly correlated on scales comparable to the spatial frequencies of the core baselines. In general, this study accentuates the fact that, for imaging applications that require high resolution (Plains array and greater), many of the core antennas can be employed as a separate subarray for low resolution science, without sacrificing the quality of the high resolution science.

We present multi-wavelength images observed with SOFIA-FORCAST from ~10 to 40 $\mu$m of 14 protostars, selected as intermediate-mass protostar candidates, as part of the SOFIA Massive (SOMA) Star Formation Survey. We build protostellar spectral energy distributions (SEDs) with the SOFIA observations, together with archival data from Spitzer, Herschel and IRAS. We then fit the SEDs with radiative transfer (RT) models of Zhang & Tan (2018), based on Turbulent Core Accretion theory, to estimate key properties of the protostars. The SEDs generally indicate the validity of these RT models down to intermediate-mass and/or early-stage protostars. The protostars analyzed so far in the SOMA survey span a range of luminosities from ~$10^{2}$ to ~$10^{6} L_{\odot}$, current protostellar masses from ~0.5 to ~30 $M_{\odot}$, and ambient clump mass surface densities, $\Sigma_{\rm cl}$ of 0.1 - 3 g cm$^{-2}$. A wide range of evolutionary states of the individual protostars and of the protocluster environments are also probed. The 19 to 37 $\mu$m spectral index of the sources correlates with outflow cavity opening angle, ratio of this angle to viewing angle, and evolutionary stage. We have also added a sample of ~50 protostellar sources identified from within Infrared Dark Clouds and expected to be at the earliest stages of their evolution. With this global sample, most of the evolutionary stages of high- and intermediate-mass protostars are probed. From the best fitting models of the protostars, there is no evidence of a threshold value of protocluster clump mass surface density being needed to form protostars up to ~25 $M_\odot$. However, to form more massive protostars, there is tentative evidence that $\Sigma_{\rm cl}$ needs to be at least 1 g cm$^{-2}$. We discuss how this is consistent with expectations from core accretion models that include internal feedback from the forming massive star.

Massive protostars generate strong radiation feedback, which may help set the mass they achieve by the end of the accretion process. Studying such feedback is therefore crucial for understanding the formation of massive stars. We report the discovery of a photoionized bipolar outflow towards the massive protostar G45.47+0.05 using high-resolution observations at 1.3 mm with the Atacama Large Millimeter/Submillimeter Array (ALMA) and at 7 mm with the Karl G. Jansky Very Large Array (VLA). By modeling the free-free continuum, the ionized outflow is found to be a photoevaporation flow with an electron temperature of 10,000 K and an electron number density of ~1.5x10^7 cm^-3 at the center, launched from a disk of radius of 110 au. H30alpha hydrogen recombination line emission shows strong maser amplification, with G45 being one of very few sources to show such millimeter recombination line masers. The mass of the driving source is estimated to be 30-50 Msun based on the derived ionizing photon rate, or 30-40 Msun based on the H30alpha kinematics. The kinematics of the photoevaporated material is dominated by rotation close to the disk plane, while accelerated to outflowing motion above the disk plane. The mass loss rate of the photoevaporation outflow is estimated to be ~(2-3.5)x10^-5 Msun/yr. We also found hints of a possible jet embedded inside the wide-angle ionized outflow with non-thermal emissions. The possible co-existence of a jet and a massive photoevaporation outflow suggests that, in spite of the strong photoionization feedback, accretion is still on-going.

In this study we analyze 70 radio continuum sources associated with dust clumps and considered to be candidates for the earliest stages of high-mass star formation. The detection of these sources was reported by Rosero et al. (2016), who found most of them to show weak (${\scriptstyle <}$1 mJy) and compact (${\scriptstyle <}\,$0.6$^{\prime \prime}$) radio emission. Herein, we used the observed parameters of these sources to investigate the origin of the radio continuum emission. We found that at least $\sim 30\%$ of these radio detections are most likely ionized jets associated with high-mass protostars, but for the most compact sources we cannot discard the scenario that they represent pressure-confined HII regions. This result is highly relevant for recent theoretical models based on core accretion that predict the first stages of ionization from high-mass stars to be in the form of jets. Additionally, we found that properties such as the radio luminosity as a function of the bolometric luminosity of ionized jets from low and high-mass stars are extremely well-correlated. Our data improve upon previous studies by providing further evidence of a common origin for jets independently of luminosity.

We present multi-wavelength images observed with SOFIA-FORCAST from $\sim$10 to 40 $\mu$m of seven high luminosity massive protostars, as part of the SOFIA Massive (SOMA) Star Formation Survey. Source morphologies at these wavelengths appear to be influenced by outflow cavities and extinction from dense gas surrounding the protostars. Using these images, we build spectral energy distributions (SEDs) of the protostars, also including archival data from Spitzer, Herschel and other facilities. Radiative transfer (RT) models of Zhang & Tan (2018), based on Turbulent Core Accretion theory, are then fit to the SEDs to estimate key properties of the protostars. Considering the best five models fit to each source, the protostars have masses $m_{*} \sim 12-64 \: M_{\odot}$ accreting at rates of $\dot{m}_{*} \sim 10^{-4}-10^{-3} \: M_{\odot} \: \rm yr^{-1}$ inside cores of initial masses $M_{c} \sim 100-500 \: M_{\odot}$ embedded in clumps with mass surface densities $\Sigma_{\rm cl} \sim 0.1-3 \: \rm g \: cm^{-2}$ and span a luminosity range of $10^{4} -10^{6} \: L_{\odot}$. Compared with the first eight protostars in Paper I, the sources analyzed here are more luminous, and thus likely to be more massive protostars. They are often in a clustered environment or have a companion protostar relatively nearby. From the range of parameter space of the models, we do not see any evidence that $\Sigma_{\rm cl}$ needs to be high to form these massive stars. For most sources the RT models provide reasonable fits to the SEDs, though the cold clump material often influences the long wavelength fitting. However, for sources in very clustered environments, the model SEDs may not be such a good description of the data, indicating potential limitations of the models for these regions.

We study centimeter continuum emission of eight high- and intermediate-mass protostars that are part of the SOFIA Massive (SOMA) Star Formation Survey, thus building extended spectral energy distributions (SEDs) from the radio to the infrared. We discuss the morphology seen in the centimeter continuum images, which are mostly derived from archival VLA data, and the relation to infrared morphology. We use the SEDs to test new models of high-mass star formation including radiative and disk-wind feedback and associated free-free and dust continuum emission (Tanaka, Tan, & Zhang 2016). We show that interferometric data of the centimeter continuum flux densities provide additional, stringent tests of the models by constraining the ionizing luminosity of the source and help to break degeneracies encountered when modeling the infrared-only SEDs, especially for the protostellar mass. Our derived parameters are consistent with physical parameters estimated by other methods such as dynamical protostellar masses. We find a few examples of additional stellar sources in the vicinity of the high-mass protostars, which may be low-mass young stellar objects. However, the stellar multiplicity of the regions, at least as traced by radio continuum emission, appears to be relatively low.

The next-generation Very Large Array (ngVLA) is an astronomical observatory planned to operate at centimeter wavelengths (25 to 0.26 centimeters, corresponding to a frequency range extending from 1.2 to 116 GHz). The observatory will be a synthesis radio telescope constituted of approximately 244 reflector antennas each of 18 meters diameter, and 19 reflector antennas each of 6 meters diameter, operating in a phased or interferometric mode. We provide a technical overview of the Reference Design of the ngVLA. This Reference Design forms a baseline for a technical readiness assessment and the construction and operations cost estimate of the ngVLA. The concepts for major system elements such as the antenna, receiving electronics, and central signal processing are presented.

We report 3 mm continuum, CH3CN(5-4) and 13CS(2-1) line observations with CARMA, in conjunction with 6 and 1.3 cm continuum VLA data, and 12 and 25 micron broadband data from the Subaru Telescope toward the massive proto-star IRAS18566+0408. The VLA data resolve the ionized jet into 4 components aligned in the E-W direction. Radio components A, C, and D have flat cm SEDs indicative of optically thin emission from ionized gas, and component B has a spectral index alpha = 1.0, and a decreasing size with frequency proportional to frequency to the -0.5 power. Emission from the CARMA 3 mm continuum, and from the 13CS(2-1), and CH3CN(5-4) spectral lines is compact (i.e. < 6700 AU), and peaks near the position of VLA cm source, component B. Analysis of these lines indicates hot, and dense molecular gas, typical for HMCs. Our Subaru telescope observations detect a single compact source, coincident with radio component B, demonstrating that most of the energy in IRAS18566+0408 originates from a region of size < 2400 AU. We also present UKIRT near-infrared archival data for IRAS18566+0408 which show extended K-band emission along the jet direction. We detect an E-W velocity shift of about 10 km/sec over the HMC in the CH3CN lines possibly tracing the interface of the ionized jet with the surrounding core gas. Our data demonstrate the presence of an ionized jet at the base of the molecular outflow, and support the hypothesis that massive protostars with O-type luminosity form with a mechanism similar to lower mass stars.

We present a high sensitivity radio continuum survey at 6 and 1.3$\,$cm using the Karl G. Jansky Very Large Array towards a sample of 58 high-mass star forming regions. Our sample was chosen from dust clumps within infrared dark clouds with and without IR sources (CMC-IRs, CMCs, respectively), and hot molecular cores (HMCs), with no previous, or relatively weak radio continuum detection at the $1\,$mJy level. Due to the improvement in the continuum sensitivity of the VLA, this survey achieved map rms levels of $\sim$ 3-10 $\mu$Jy beam$^{-1}$ at sub-arcsecond angular resolution. We extracted 70 centimeter continuum sources associated with 1.2$\,$mm dust clumps. Most sources are weak, compact, and are prime candidates for high-mass protostars. Detection rates of radio sources associated with the mm dust clumps for CMCs, CMC-IRs and HMCs are 6$\%$, 53$\%$ and 100$\%$, respectively. This result is consistent with increasing high-mass star formation activity from CMCs to HMCs. The radio sources located within HMCs and CMC-IRs occur close to the dust clump centers with a median offset from it of 12,000$\,$AU and 4,000$\,$AU, respectively. We calculated 5 - 25$\,$GHz spectral indices using power law fits and obtain a median value of 0.5 (i.e., flux increasing with frequency), suggestive of thermal emission from ionized jets. In this paper we describe the sample, observations, and detections. The analysis and discussion will be presented in Paper II.

We present results of Chandra ACIS-I and Karl G. Jansky Very Large Array (VLA) 6 cm continuum observations of the IRAS 20126+4104 massive star forming region. We detect 150 X-ray sources within the 17 arcmin x 17 arcmin ACIS-I field, and a total of 13 radio sources within the 9'.2 primary beam at 4.9 GHz. Among these are the first 6 cm detections of the central sources reported by Hofner et al. (2007), namely I20N1, I20S, and I20var. A new variable radio sources is also reported. Searching the 2MASS archive we identified 88 NIR counterparts to the X-ray sources. Only 4 of the X-ray sources had 6 cm counterparts. Based on an NIR color-color analysis, and on the Besancon simulation of Galactic stellar populations (Robin et al. 2003), we estimate that about 90 X-ray sources are associated with this massive star forming region. We detect an increasing surface density of X-ray sources toward the massive protostar and infer the presence of a cluster of at least 46 YSOs within a distance of 1.2 pc from the massive protostar.

We report 1.3 cm and 6 cm continuum observations toward the massive proto-stellar candidate G11.11-0.12P1 using the Karl G. Jansky Very Large Array (VLA). We detect a string of four unresolved radio continuum sources coincident with the mid-IR source in G11P1. The continuum sources have positive spectral indices consistent with a thermal (free-free) ionized jet. The most likely origin of the ionized gas are shocks due to the interaction of a stellar wind with the surrounding high-density material. We also present NIR United Kingdom Infrared Telescope (UKIRT) archival data which show an extended structure detected only at K-band (2.2 $\mu$m), which is oriented perpendicular to the jet, and that may be scattered light from a circumstellar disk around the massive protostar. Our observations plus the UKIRT archival data thus provide new evidence that a disk/jet system is present in the massive protostellar candidate located in the G11.11-0.12P1 core.

We report the results of a survey in the CH3CN J= 12-11 transition toward a sample of massive proto-stellar candidates. The observations were carried out with the 10 m Submillimeter telescope on Mount Graham, AZ. We detected this molecular line in 9 out of 21 observed sources. In six cases this is the first detection of this transition. We also obtained full beam sampled cross-scans for five sources which show that the lower K-components can be extended on the arcminute angular scale. The higher K-components however are always found to be compact with respect to our 36" beam. A Boltzmann population diagram analysis of the central spectra indicates CH3CN column densities of about 10^14 cm^(-2), and rotational temperatures above 50 K, which confirms these sources as hot molecular cores. Independent fits to line velocity and width for the individual K-components resulted in the detection of an increasing blue shift with increasing line excitation for four sources. Comparison with mid-infrared images from the SPITZER GLIMPSE/IRAC archive for six sources show that the CH3CN emission is generally coincident with a bright mid-IR source. Our data clearly show that the CH3CN J= 12-11 transition is a good probe of the hot molecular gas near massive protostars, and provide the basis for future interferometric studies.

The Wilkinson Microwave Anisotropy Probe (WMAP) free-free foreground emission map is used to identify diffuse ionized regions (DIR) in the Galaxy (Rahman & Murray 2010). It has been found that the 18 most luminous WMAP sources produce more than half of the total ionizing luminosity of the Galaxy. We observed radio recombination lines (RRLs) toward the luminous WMAP source G49.75-0.45 with the Green Bank Telescope near 1.4 GHz. Hydrogen RRL is detected toward the source but no helium line is detected, implying that n_He+/n_H+ < 0.024. This limit puts severe constraint on the ionizing spectrum. The total ionizing luminosity of G49 (3.05 x 10^51 s^-1) is ~ 2.8 times the luminosity of all radio HII regions within this DIR and this is generally the case for other WMAP sources. Murray & Rahman (2010) propose that the additional ionization is due to massive clusters (~ 7.5 x10^3 Msun for G49) embedded in the WMAP sources. Such clusters should produce enough photons with energy \geq 24.6 eV to fully ionize helium in the DIR. Our observations rule out a simple model with G49 ionized by a massive cluster. We also considered 'leaky' HII region models for the ionization of the DIR, suggested by Lockman and Anantharamaiah, but these models also cannot explain our observations. We estimate that the helium ionizing photons need to be attenuated by > ~10 times to explain the observations. If selective absorption of He- ionizing photons by dust is causing this additional attenuation, then the ratio of dust absorption cross sections for He- and H- ionizing photons should be > ~6.

We report the orbital parameters for ROXR1 14 and RX J1622.7-2325Nw, two young, low-mass, and double-lined spectroscopic binaries recently discovered in the Ophiuchus star forming region. Accurate orbital solutions were determined from over a dozen high-resolution spectra taken with the Keck II and Gemini South telescopes. These objects are T Tauri stars with mass ratios close to unity and periods of ~5 and ~3 days, respectively. In particular, RX J1622.7-2325Nw shows a non-circularized orbit with an eccentricity of 0.30, higher than any other short-period pre-main sequence spectroscopic binary known to date. We speculate that orbit of RX J1622.7-2325Nw has not yet circularized because of the perturbing action of a ~1" companion, itself a close visual pair. A comparison of known young spectroscopic binaries and main sequence (MS) spectroscopic binaries in the eccentricity-period plane shows an indistinguishable distribution of the two populations, implying that orbital circularization occurs in the first 1 Myr of a star's lifetime. With the results presented in this paper we are increasing by ~4% the small sample of PMS spectroscopic binary stars with known orbital elements.