We analyzed high angular resolution 45.5 GHz images of the W49 North massive star forming region obtained in 1998 and 2016 with the Very Large Array. Most of the ultracompact HII regions show no detectable changes over the time interval of the observations. However, subcomponents B1, B2, G2a and G2c have increased its peak flux densities by values in the range of 3.8 to 21.4 \%. Most interestingly, the cometary region C clearly shows proper motions that at the distance of the region are equivalent to a velocity of 76$\pm$6 km s$^{-1}$ in the plane of the sky. We interpret this region as the ionized bowshock produced by a runaway O6 ZAMS star that was ejected from the eastern edge of Welch's ring about 6,400 years ago.

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

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

We describe the remote facilities operated by the Southeastern Association for Research in Astronomy (SARA), a consortium of colleges and universities in the US partnered with Lowell Observatory, the Chilean National Telescope Allocation Committee, and the Instituto de Astrofisica de Canarias. SARA observatories comprise a 0.96m telescope at Kitt Peak, Arizona; a 0.6m instrument on Cerro Tololo, Chile; and the 1m Jacobus Kapteyn Telescope at the Roque de los Muchachos, La Palma, Spain. All are operated using standard VNC or Radmin protocols communicating with on-site PCs. Remote operation offers considerable flexibility in scheduling, allowing long-term observational cadences difficult to achieve with classical observing at remote facilities, as well as obvious travel savings. Multiple observers at different locations can share a telescope for training, educational use, or collaborative research programs. Each telescope has a CCD system for optical imaging, using thermoelectric cooling to avoid the need for frequent local service, and a second CCD for offset guiding. The Arizona and Chile instruments also have fiber-fed echelle spectrographs. Switching between imaging and spectroscopy is very rapid, so a night can easily accommodate mixed observing modes. We present some sample observational programs. For the benefit of other groups organizing similar consortia, we describe the operating structure and principles of SARA, as well as some lessons learned from almost 20 years of remote operations.

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

W49N is a mini-starburst in the Milky Way and thus an ideal laboratory for high-mass star formation studies. Due to its large distance (11.1$_{-0.7}^{+0.9}$ kpc), the kinematics inside and between the dense molecular clumps in W49N are far from well understood. The SMA observations resolved the continuum emission into two clumps. The molecular line observation of SO$_{2}$ (28$_{4,24}$-28$_{3,25}$) suggests that the two clumps have a velocity difference of $\sim$7 km~s$^{-1}$. The eastern clump is very close to two radio sources "G1" and "G2", and the western clump coincides with a radio source "B". The HCN (3-2) line reveals an extremely energetic outflow, which is among the most energetic molecular outflows in the Milky Way. This is the first report of high-velocity molecular outflow detection in W49N. The outflow jet might be in precession, which could account for the distribution, velocity and rotation of water maser spots. Three absorption systems are identified in HCO$^{+}$ (3-2) spectra. The absorption features are blueshifted with respect to the emission of SO$_{2}$ (28$_{4,24}$-28$_{3,25}$) lines, indicating that a cold layer is expanding in front of the warm gas. Further analysis indicates that the expansion is decelerated from the geometric expansion centers.

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

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

The discovery of transiting planets around bright stars holds the potential to greatly enhance our understanding of planetary atmospheres. In this work we present the search for transits of HD 168443b, a massive planet orbiting the bright star HD 168443 V=6.92 with a period of 58.11 days. The high eccentricity of the planetary orbit e=0.53 significantly enhances the a-priori transit probability beyond that expected for a circular orbit, making HD 168443 a candidate for our ongoing Transit Ephemeris Refinement and Monitoring Survey (TERMS). Using additional radial velocities from Keck-HIRES, we refined the orbital parameters of this multi-planet system and derived a new transit ephemeris for HD 168443b. The reduced uncertainties in the transit window make a photometric transit search practicable. Photometric observations acquired during predicted transit windows were obtained on three nights. CTIO 1.0 m photometry acquired on 2010 September 7 had the required precision to detect a transit but fell just outside of our final transit window. Nightly photometry from the T8 0.8 m Automated Photometric Telescope (APT) at Fairborn Observatory, acquired over a span of 109 nights, demonstrates that HD 168443 is constant on a time scale of weeks. Higher-cadence photometry on 2011 April 28 and June 25 shows no evidence of a transit. We are able to rule out a non-grazing transit of HD 168443b.

We have imaged the Sgr B2 Main region with the Very Large Array in the BnA configuration ($\theta_{beam}$ = 0\farcs13) in both the H52$\alpha$ (45.453 GHz) radio recombination line (RRL) and 7 mm continuum emission. At a distance of 8500 pc, this spatial resolution corresponds to a physical scale of 0.005 pc ($\sim$1100 AU). The current observations detect H52$\alpha$ emission in 12 individual ultracompact (UC) and hypercompact (HC) HII regions. Two of the sources with detected H52 $\alpha$ emission have broad ($\Delta$V$_{FWHM}\sim$50 \kms) recombination lines, and two of the sources show lines with peaks at more than one velocity. We use line parameters from the H52$\alpha$ lines and our previous H66$\alpha$ line observations to determine the relative contribution of thermal, pressure and kinematic broadening, and electron density. These new observations suggest that pressure broadening can account for the broad lines in some of the sources, but that gas motions (e.g. turbulence, accretion or outflow) contribute significantly to the broad lines in at least one of the sources (Sgr B2 F3).

We present high spatial resolution mid-IR images of the ring of UCHII regions in W49A obtained at Gemini North, allowing us to identify the driving source of its powerful H2O maser outflow. These data also confirm our previous report that several radio sources in the ring are undetected in the mid-IR because they are embedded deep inside the cloud core. We locate the source of the water maser outflow at the position of the compact mid-IR peak of source G (source G:IRS1). This IR source is not coincident with any identified compact radio continuum source, but is coincident with a hot molecular core, so we propose that G:IRS1 is a hot core driving an outflow analogous to the wide-angle bipolar outflow in OMC-1. G:IRS1 is at the origin of a larger bipolar cavity and CO outflow. The water maser outflow is orthogonal to the bipolar CO cavity, so the masers probably reside near its waist in the cavity walls. Models of the IR emission require a massive protostar of 45Msun, 3e5Lsun, and an effective envelope accretion rate of 1e-3Msun/yr. Feedback from the central star could potentially drive the H2O maser outflow, but it has insufficient radiative momentum to have driven the large-scale CO outflow, requiring that this massive star had an active accretion disk over the past 10^4 yr. Combined with the spatialy resolved morphology in IR images, G:IRS1 in W49 provides compelling evidence for a massive protostar that formed by accreting from a disk, accompanied by a bipolar outflow.

We present Spitzer Space Telescope IRAC and MIPS observations toward a sample of nine high-mass star forming regions at a distance of around 2 kpc. Based on IRAC and MIPS 24 $\mu$m photometric results and 2MASS JHKs data, we carry out a census of young stellar objects (YSOs) in a 5' by 5' field toward each region. Toward seven out of the nine regions, we detect parsec sized clusters with around 20 YSOs surrounded by a more extended and sparse distribution of young stars and protostars. For the other two regions, IRAS 20126+4104 and IRAS 22172+5549, the former has the lowest number of YSOs in the sample and shows no obvious cluster, and the latter appears to be part of a larger, potentially more evolved cluster. The deep IRAC imaging reveals at least twelve outflows in eight out of the nine regions, with nine outflows prominent in the 4.5 $\mu$m band most probably attributed to shocked H$_2$ emission, two outflows dominated by scattered light in the 3.6 and 4.5 $\mu$m bands, and one outflow standing out from its hydrocarbon emission in the 8.0 $\mu$m band. In comparison with previous ground-based observations, our IRAC observations reveal new outflow structures in five regions. The dramatically different morphologies of detected outflows can be tentatively interpreted in terms of possible evolution of massive outflows. The driving sources of these outflows are deeply embedded in dense dusty cores revealed by previous millimeter interferometric observations. We detect infrared counterparts of these dusty cores in the IRAC or MIPS 24 $\mu$m bands. Reflection nebulae dominated by the emission from UV heated hydrocarbons in the 8 $\mu$m band can be found in most regions and they may imply the presence of young B stars.

We have detected C91$\alpha$ (8.5891 GHz) emission toward 4 ultra-compact \HII regions (\UCHII s; W49G, J, L & C) in the W49 North massive star forming region with the Very Large Array (VLA) at 3\arcsec resolution. No carbon line emission was detected toward \UCHII s W49F, A, O, S and Q at this frequency to a 3$\sigma$ level of 2 mJy. We also observed the same region in the C75$\alpha$ line (15.3 GHz) with no detection at a 3$\sigma$ level of 6 mJy with a 1\arcsec.7 beam. Detection of line emission toward these sources add supporting data to the earlier result of \nocite{retal05a}Roshi et al (2005a) that many \UCHII s have an associated photo-dissociation region (PDR). Similarity of the LSR velocities of carbon recombination lines and H$_2$CO absorption toward \UCHII s in W49 North suggests that the PDRs reside in the dense interface zone surrounding these \HII regions. Combining the observed carbon line parameters at 8.6 GHz with the upper limits on line emission at 15.3 GHz, we obtain constraints on the physical properties of the PDRs associated with W49G and J. The upper limit on the number density of hydrogen molecule obtained from carbon line models is $\sim$ $5 \times 10^6$ \cmthree.

We present three epochs of VLBA observations of Galactic HI absorption toward the quasar 3C~138 with resolutions of 20 mas (~ 10 AU). This analysis includes VLBA data from observations in 1999 and 2002 along with a reexamination of 1995 VLBA data. Improved data reduction and imaging techniques have led to an order of magnitude improvement in sensitivity compared to previous work. With these new data we confirm the previously detected milliarcsecond scale spatial variations in the HI opacity at the level of Delta(tau_{max}) =0.50 \pm 0.05. The typical size scale of the optical depth variations is ~ 50 mas or 25 AU. In addition, for the first time we see clear evidence for temporal variations in the HI opacity over the seven year time span of our three epochs of data. We also attempted to detect the magnetic field strength in the HI gas using the Zeeman effect. From this analysis we have been able to place a 3 sigma upper limit on the magnetic field strength per pixel of ~45 muG. We have also been able to calculate for the first time the plane of sky covering fraction of the small scale HI gas of ~10%. We also find that the line widths of the milliarcsecond sizescale HI features are comparable to those determined from previous single dish measurements toward 3C~138, suggesting that the opacity variations cannot be due to changes in the HI spin temperature. From these results we favor a density enhancement interpretation for the small scale HI structures, although these enhancements appear to be of short duration and are unlikely to be in equilibrium.

We have used Very Large Array (VLA) observations of the massive star forming regions W49A and Sgr B2, obtained with resolutions from 2\farcs0 to 0\farcs04, to classify the morphologies of nearly 100 ultracompact HII regions. These high resolution, multi-frequency, multi-configuration VLA observations motivate several modifications of the existing morphological classification scheme for UC HII regions. In this work, we describe the modified morphology scheme and the criteria used in source classification. In particular, we drop the ``core-halo'' classification, add a ``bipolar'' classification, and change the shell classification to ``shell-like''. We tally the percentage of each morphology found in the Sgr B2 and W49A regions and find broad agreement with the Galactic plane surveys in the distribution of morphologies for most types. However, we find that nearly a third of the sources in these regions are shell-like, which is a higher percentage by nearly a factor of ten than found in the surveys of Galactic plane star forming regions by Wood & Churchwell (1989a) and Kurtz et al. (1994). This difference may be due to physical differences in the environments of these two extreme star forming regions. Alternatively, differences in observational technique may be responsible.

We report here on a survey of carbon recombination lines (RLs) near 8.5 GHz toward 17 ultra-compact \HII regions (\UCHII s). Carbon RLs are detected in 11 directions, indicating the presence of dense photodissociation regions (PDRs) associated with the \UCHII s. In this paper, we show that the carbon RLs provide important, complementary information on the kinematics and physical properties of the ambient medium near \UCHII s. Non-LTE models for the carbon line forming region are developed, assuming that the PDRs surround the \UCHII s, and we constrained the model parameters by multi-frequency RL data. Modeling shows that carbon RL emission near 8.5 GHz is dominated by stimulated emission and hence we preferentially observe the PDR material that is in front of the \UCHII continuum. We find that the relative motion between ionized gas and the associated PDR is about half that estimated earlier, and has an RMS velocity difference of 3.3 \kms. Our models also give estimates for the PDR density and pressure. We found that the neutral density of PDRs is typically $>$ 5 $\times$ 10$^5$ \cmthree and \UCHII s can be embedded in regions with high ambient pressure. Our results are consistent with a pressure confined \HII region model where the stars are moving relative to the cloud core. Other models cannot be ruled out, however. Interestingly, in most cases, the PDR pressure is an order of magnitude larger than the pressure of the ionized gas. Further investigation is needed to understand this large pressure difference.

Using the Very Large Array (VLA) of the National Radio Astronomy Observatory in the A and B configurations, we have obtained simultaneous high resolution observations of both the 22 GHz water maser lines as well as the 22 GHz continuum for the H II regions W49N and Sagittarius B2. The angular resolution of both observations is ~0.1", which at the distance of W49N (11.4 kpc; Gwinn, Moran, & Reid 1992) and Sgr B2 (8.5 kpc) corresponds to a physical size of <1000 AU in both sources. The velocity coverage for W49N is +/- 435 km/s; positions for 316 water maser components were obtained. The velocity coverage for Sgr B2 is -40 to +120 km/s; positions for 68 maser components were determined in Sgr B2 Main, 79 in Sgr B2 North, 14 in Sgr B2 Mid-North, and 17 in Sgr B2 South, for a total of 178 water maser positions in Sgr B2. The cross calibration scheme of Reid & Menten (1990, 1997) was used. Using this procedure, high dynamic range continuum images were obtained with accurate registration (rms ~0.01") of the continuum and maser positions. A detailed comparison between H II components and maser positions for both Sgr B2 and W49N is presented. In Sgr B2 Main, the water masers are predominantly located at the outside edge of the high-frequency continuum, lending support to the proposal that entrainment by stellar winds may play an important role in water maser emission.

High resolution 7 mm observations of the W49N massive star forming region have detected recombination line emission from the individual ultracompact (UC) HII regions on 50 milliarcsecond (600 AU) scales. These line observations, combined with multifrequency, high-resolution continuum imaging of the region at 7 mm (VLA) and at 3 mm and 1 mm (BIMA), indicate that five to seven of the eighteen ultracompact sources in W49N are broad recombination line objects (BRLOs) as described by Jaffe & Martin-Pintado (1999). BRLOs have both broad radio recombination lines ($\Delta$V$>$60 \kms) and rising spectra (S$_{\nu}\sim\nu^{\alpha}$), with $\alpha$ values greater than 0.4. The broad line widths of the H52$\alpha$ line are probably related to motions in the ionized gas rather than pressure broadening. A number of models have been proposed to explain the long lifetime of UC HII regions, including the photoevaporated disk model proposed by Hollenbach et al. (1994). This model can also explain the broad lines, rising spectra and bipolar morphologies of some sources. We suggest$-$based on line and continuum observations as well as source morphology$-$that in a subset of the W49N ultracompact sources we may be observing ionized winds that arise from circumstellar disks.

Absorption in the J{K-K+} = 2{11}-2{12} transition of formaldehyde at 2 cm towards the ultracompact HII regions C1 and C2 of W 58 has been observed with the VLA with an angular resolution of ~0.2'' and a velocity resolution of ~1 km/s. The high resolution continuum image of C1 (ON 3) shows a partial shell which opens to the NE. Strong H2CO absorption is observed against W 58 C1. The highest optical depth (tau > 2) occurs in the SW portion of C1 near the edge of the shell, close to the continuum peak. The absorption is weaker towards the nearby, more diffuse compact HII region C2, tau<~0.3. The H2CO velocity (-21.2 km/s) towards C1 is constant and agrees with the velocity of CO emission, mainline OH masers, and the H76 alpha recombination line, but differs from the velocity of the 1720 MHz OH maser emission (~-13 km/s). Observations of the absorption in the J{K-K+} = 1{10}-1{11} transition of formaldehyde at 6 cm towards W 58 C1 and C2 carried out earlier with the WSRT at lower resolution (~4''x7'') show comparable optical depths and velocities to those observed at 2 cm. Based on the mean optical depth profiles at 6 cm and 2 cm, the volume density of molecular hydrogen n(H2) and the formaldehyde column density N(H2CO) were determined. The n(H2) is ~6E4 /cm**3 towards C1. N(H2CO) for C1 is ~8E14 /cm**2 while that towards C2 is ~8E13 /cm**2.

NGC 3576 (G291.28-0.71; l=291.3o, b=-0.7o) and NGC 3603 (G291.58-0.43; l=291.6o, b=-0.5o) are optically visible, luminous HII regions located at distances of 3.0 kpc and 6.1 kpc, respectively. We present 3.4 cm Australian Telescope Compact Array (ATCA) observations of these two sources in the continuum and the H90a, He90a, C90a and H113b recombination lines with an angular resolution of 7" and a velocity resolution of 2.6 km/s. All four recombination lines are detected in the integrated profiles of the two sources. Broad radio recombination lines are detected in both NGC 3576 (DV_{FWHM}>= 50 km/s) and NGC 3603 (DV_{FWHM}>=70 km/s). In NGC 3576 a prominent N-S velocity gradient (~30 km/s/pc) is observed, and a clear temperature gradient (6000 K to 8000 K) is found from east to west, consistent with a known IR color gradient in the source. In NGC 3603, the H90a, He90a and the H113b lines are detected from 13 individual sources. The Y^+ (He/H) ratios in the two sources range from 0.08+/-0.04 to 0.26+/-0.10. We compare the morphology and kinematics of the ionized gas at 3.4 cm with the distribution of stars, 10 micron emission and H_2O, OH, and CH_3OH maser emission. These comparisons suggest that both NGC 3576 and NGC 3603 have undergone sequential star formation.