Search

We present deep, sub-arcsecond ($\sim$2000 AU) resolution ALMA 0.82 mm observations of the former high-mass prestellar core candidate G11.92-0.61 MM2, recently shown to be an $\sim$500 AU-separation protobinary. Our observations show that G11.92-0.61 MM2, located in the G11.92-0.61 protocluster, lies on a filamentary structure traced by 0.82 mm continuum and N$_2$H$^+$(4-3) emission. The N$_2$H$^+$(4-3) spectra are multi-peaked, indicative of multiple velocity components along the line of sight. To analyse the gas kinematics, we performed pixel-by-pixel Gaussian decomposition of the N$_2$H$^+$ spectra using SCOUSEPY and hierarchical clustering of the extracted velocity components using ACORNS. Seventy velocity- and position-coherent clusters (called "trees") are identified in the N$_2$H$^+$-emitting gas, with the 8 largest trees accounting for >60% of the fitted velocity components. The primary tree, with $\sim$20% of the fitted velocity components, displays a roughly north-south velocity gradient along the filamentary structure traced by the 0.82 mm continuum. Analysing a $\sim$0.17 pc-long substructure, we interpret its velocity gradient of $\sim$10.5 km s$^{-1}$pc$^{-1}$ as tracing filamentary accretion towards MM2 and estimate a mass inflow rate of $\sim$1.8$\times10^{-4}$ to 1.2$\times10^{-3}$ M$_\odot$ yr$^{-1}$. Based on the recent detection of a bipolar molecular outflow associated with MM2, accretion onto the protobinary is ongoing, likely fed by the larger-scale filamentary accretion flows. If 50% of the filamentary inflow reaches the protostars, each member of the protobinary would attain a mass of 8 M$_\odot$ within $\sim1.6\times$10$^5$ yr, comparable to the combined timescale of the 70 $\mu$m- and MIR-weak phases derived for ATLASGAL-TOP100 massive clumps using chemical clocks., Comment: 21 pages, 15 figures, 2 tables, Accepted for publication in MNRAS

Complex organic molecules (COMs) have been observed to be abundant in the gas phase toward protostars. Deep line surveys have been carried out only for a limited number of well-known high-mass star forming regions using the Atacama Large Millimeter/submillimeter Array (ALMA), which has unprecedented resolution and sensitivity. Statistical studies on oxygen-bearing COMs (O-COMs) in high-mass protostars using ALMA are still lacking. With the recent CoCCoA survey, we are able to determine the column density ratios of six O-COMs with respect to methanol (CH$_3$OH) in a sample of 14 high-mass protostellar sources to investigate their origin through ice and/or gas-phase chemistry. The selected species are: acetaldehyde (CH$_3$CHO), ethanol (C$_2$H$_5$OH), dimethyl ether (DME, CH$_3$OCH$_3$), methyl formate (MF, CH$_3$OCHO), glycolaldehyde (GA, CH$_2$OHCHO), and ethylene glycol (EG, (CH$_2$OH)$_2$). DME and MF have the highest and most constant ratios within one order of magnitude, while the other four species have lower ratios and exhibit larger scatter by 1-2 orders of magnitude. We compare the O-COM ratios of high-mass CoCCoA sources with those of 5 low-mass protostars available from the literature, along with the results from experiments and simulations. We find that the O-COM ratios with respect to methanol are on the same level in both the high- and low-mass samples, which suggests that these species are mainly formed in similar environments during star formation, probably in ice mantles on dust grains during early pre-stellar stages. Current simulations and experiments can reproduce most observational trends with a few exceptions, and hypotheses exist to explain the differences between observations and simulations/experiments, such as the involvement of gas-phase chemistry and different emitting areas of molecules., Comment: Accepted for publication in Astronomy & Astrophysics. 17 pages and 7 figures in the main body, 15 pages and 8 figures in the appendix

High-mass protostars (M$_{\star} >$ 8 M$_{\odot}$) are thought to gain the majority of their mass via short, intense bursts of growth. This episodic accretion is thought to be facilitated by gravitationally unstable and subsequently inhomogeneous accretion disks. Limitations of observational capabilities, paired with a lack of observed accretion burst events has withheld affirmative confirmation of the association between disk accretion, instability and the accretion burst phenomenon in high-mass protostars. Following its 2019 accretion burst, a heat-wave driven by a burst of radiation propagated outward from the high-mass protostar G358.93-0.03-MM1. Six VLBI (very long baseline interferometry) observations of the raditively pumped 6.7 GHz methanol maser were conducted during this period, tracing ever increasing disk radii as the heat-wave propagated outward. Concatenating the VLBI maps provided a sparsely sampled, milliarcsecond view of the spatio-kinematics of the accretion disk covering a physical range of $\sim$ 50 - 900 AU. We term this observational approach `heat-wave mapping'. We report the discovery of a Keplerian accretion disk with a spatially resolved four-arm spiral pattern around G358.93-0.03-MM1. This result positively implicates disk accretion and spiral arm instabilities into the episodic accretion high-mass star formation paradigm., Comment: Published in Nature Astronomy in 2023

High-mass stars are thought to accumulate much of their mass via short, infrequent bursts of disk-aided accretion. Such accretion events are rare and difficult to observe directly but are known to drive enhanced maser emission. In this Letter we report high-resolution, multi-epoch methanol maser observations toward G358.93-0.03 which reveal an interesting phenomenon; the sub-luminal propagation of a thermal radiation "heat-wave" emanating from an accreting high-mass proto-star. The extreme transformation of the maser emission implies a sudden intensification of thermal infrared radiation from within the inner (40 mas, 270 au) region. Subsequently, methanol masers trace the radial passage of thermal radiation through the environment at $\geq$ 4-8\% the speed of light. Such a high translocation rate contrasts with the $\leq$ 10 km s$^{-1}$ physical gas motions of methanol masers typically observed using very long baseline interferometry (VLBI). The observed scenario can readily be attributed to an accretion event in the high-mass proto-star G358.93-0.03-MM1. While being the third case in its class, G358.93-0.03-MM1 exhibits unique attributes hinting at a possible `zoo' of accretion burst types. These results promote the advantages of maser observations in understanding high-mass star formation, both through single-dish maser monitoring campaigns and via their international cooperation as VLBI arrays., Comment: Published in Nature Astronomy in 2020

Noise-induced gate errors remain one of the main obstacles to realizing a broad range of quantum information technologies. Dynamical error suppression using carefully designed control schemes is critical for overcoming this challenge. Such schemes must be able to correct against multiple noise sources simultaneously afflicting a qubit in order to reach error correction thresholds. Here, we present a general framework for designing control fields that simultaneous suppress both noise in the fields themselves as well as transverse dephasing noise. Using the recently developed Space Curve Quantum Control formalism, in which robust quantum evolution is mapped to closed geometric curves in a multidimensional Euclidean space, we derive necessary and sufficient conditions that guarantee the cancellation of both types of noise to leading order. We present several techniques for solving these conditions and provide explicit examples of error-resistant control fields. Our work also sheds light on the relation between holonomic evolution and the suppression of control field errors., Comment: 12 pages, 4 figures

We present high-resolution (< ~160 au) Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm observations of the high-mass prestellar core candidate G11.92-0.61 MM2, which reveal that this source is in fact a protobinary system with a projected separation of 505 au. The binary components, MM2E and MM2W, are compact (radii < 140 au) sources within the partially optically thick dust emission with alpha$_{0.9 cm-1.3 mm}$=2.47-2.94. The 1.3 mm brightness temperatures, T$_b$=68.4/64.6 K for MM2E/MM2W, imply internal heating and minimum luminosities L$_*$ > 24.7 Lsun for MM2E and L$_*$ > 12.6 Lsun for MM2W. The compact sources are connected by a "bridge" of lower-surface-brightness dust emission and lie within more extended emission that may correspond to a circumbinary disk. The circumprotostellar gas mass, estimated from ~0.2" resolution VLA 0.9 cm observations assuming optically thin emission, is 6.8+/-0.9 Msun. No line emission is detected towards MM2E and MM2W in our high-resolution 1.3 mm ALMA observations. The only line detected is 13CO J=2-1, in absorption against the 1.3 mm continuum, which likely traces a layer of cooler molecular material surrounding the protostars. We also report the discovery of a highly asymmetric bipolar molecular outflow that appears to be driven by MM2E and/or MM2W in new deep, ~0.5" resolution (1680 au) ALMA 0.82 mm observations. This outflow, traced by low-excitation CH3OH emission, indicates ongoing accretion onto the protobinary system. Overall, the super-Alfvenic models of Mignon-Risse et al. (2021) agree well with the observed properties of the MM2E/MM2W protobinary, suggesting that this system may be forming in an environment with a weak magnetic field., Comment: 11 pages, 3 figures, 2 tables; accepted for publication in ApJL, no changes to manuscript

We present the most complete to date interferometric study of the centimeter wavelength methanol masers detected in G358.93-0.03 at the burst and post-burst epochs. A unique, NIR/(sub)mm-dark and FIR-loud MYSO accretion burst was recently discovered in G358.93-0.03. The event was accompanied by flares of an unprecedented number of rare methanol maser transitions. The first images of three of the newly-discovered methanol masers at 6.18, 12.23, and 20.97 GHz are presented in this work. The spatial structure evolution of the methanol masers at 6.67, 12.18, and 23.12 GHz is studied at two epochs. The maser emission in all detected transitions resides in a region of $\sim$0.2$^{\prime\prime}$ around the bursting source and shows a clear velocity gradient in the north-south direction, with red-shifted features to the north and blue-shifted features to the south. A drastic change in the spatial morphology of the masing region is found: a dense and compact "spiral" cluster detected at epoch I evolved into a disperse, "round" structure at epoch II. During the transition from the first epoch to the second, the region traced by masers expanded. The comparison of our results with the complementary VLA, VLBI, SMA, and ALMA maser data is conducted. The obtained methanol maser data support the hypothesis of the presence of spiral-arm structures within the accretion disk, which was suggested in previous studies of the source.

Thioacetamide (CH3CSNH2) is the sulfur analog to acetamide (CH3CONH2) and it is a viable candidate to search for in astronomical environments specifically toward regions where other S-bearing molecules have been found and, if possible, that also contain a detection of CH3CONH2. If detected, it would not only continue to expand the view of molecular complexity in astronomical environments, but also help to better elucidate the possible formation pathways of these types of species in these environments. The rotational spectrum of CH3CSNH2 was investigated up to 650 GHz. Using the newly refined spectrum of CH3CSNH2, as well as additional spectroscopic data on the chemically related species CH3CONH2, a variety of astronomical sources were searched including data from the following large surveys: The PRIMOS conducted with the Green Bank Telescope (GBT); Exploring molecular complexity with ALMA (EMoCA) conducted with ALMA; and Astrochemical Surveys at IRAM (ASAI) conducted with the Institut de Radioastronomie Millimetrique (IRAM) 30m Telescope. A total of 1428 transitions from the vt=0 state with maximum values J=47 and Ka=20 in the range up to 330 GHz, and J=95 and Ka=20 in the range from 400 - 660 GHz were assigned. We also assigned 321 transitions from the vt=1 state with the maximum values J=35 and Ka=9 up to 330 GHz. The final fit is based on the rho-axis-method (RAM) Hamiltonian model that includes 40 parameters. An astronomical search for CH3CSNH2 was conducted based on all the new spectroscopic data. No transitions of CH3CSNH2 were detected toward any of the sources contained in our survey. Using the appropriate telescope and physical parameters for each astronomical source, upper limits to the column densities were found for CH3CSNH2 toward each source., Comment: 16 pages, 8 figures, accepted for publication in Astronomy & Astrophysics

We have observed a sample of 9 Extended Green Objects (EGOs) at 1.3 and 5 cm with the VLA with sub-arcsecond resolution and ~7-14 uJy/beam sensitivities in order to characterize centimeter continuum emission as it first appears in these massive protoclusters. We find EGO-associated continuum emission - within 1 arcsec of the extended 4.5 um emission - in every field, which is typically faint (order 10^1-10^2 uJy) and compact (unresolved at 0.3-0.5 arcsec). The derived spectral indices of our 36 total detections are consistent with a wide array of physical processes, including both non-thermal (19% of detections) and thermal free-free processes (e.g. ionized jets and compact HII regions, 78% of sample), and warm dust (1 source). We also find EGO-associated 6.7 GHz CH$_3$OH and 22 GHz H$_2$O maser emission in 100% of the sample, and NH$_3$ (3,3) masers in ~45%; we do not detect any NH$_3$ (6,6) masers at ~5.6 mJy/beam sensitivity. We find statistically-significant correlations between radio-distance and bolometric luminosities at two physical scales and three frequencies, consistent with thermal emission from ionized jets, but no correlation between water-maser and radio-distance luminosities for our sample. From these data, we conclude that EGOs likely host multiple different centimeter continuum-producing processes simultaneously. Additionally, at our ~1000 au resolution, we find that all EGOs except G18.89$-$0.47 contain 1 to 2 massive sources based on the presence of methanol maser groups, which is consistent with our previous work suggesting that these are typical massive protoclusters in which only one to a few of the YSOs are massive., Comment: 37 pages, 13 figures (incl. 2 figure sets), 8 tables (incl. 4 MRTs)

In recent years, dramatic outbursts have been identified toward massive protostars via infrared and millimeter dust continuum and molecular maser emission. The longest lived outburst ($>6$ yr) persists in NGC6334I-MM1, a deeply-embedded object with no near-IR counterpart. Using FORCAST and HAWC+ on SOFIA, we have obtained the first mid-infrared images of this field since the outburst began. Despite being undetected in pre-outburst ground-based 18 $\mu$m images, MM1 is now the brightest region at all three wavelengths (25, 37, and 53 $\mu$m), exceeding the ultracompact HII region MM3 (NGC6334F). Combining the SOFIA data with ALMA imaging at four wavelengths, we construct a spectral energy distribution of the combination of MM1 and the nearby hot core MM2. The best-fit Robitaille radiative transfer model yields a luminosity of $(4.9\pm0.8)\times10^4 L_\odot$. Accounting for an estimated pre-outburst luminosity ratio MM1:MM2 = $2.1\pm0.4$, the luminosity of MM1 has increased by a factor of $16.3\pm4.4$. The pre-outburst luminosity implies a protostar of mass 6.7 $M_\odot$, which can produce the ionizing photon rate required to power the pre-outburst hypercompact HII region surrounding the likely outbursting protostar MM1B. The total energy and duration of the outburst exceed the S255IR-NIRS3 outburst by a factor of $\gtrsim3$, suggesting a different scale of event involving expansion of the protostellar photosphere (to $\gtrsim$ 20 $R_\odot$), thereby supporting a higher accretion rate ($\gtrsim$0.0023 $M_\odot$ yr$^{-1}$) and reducing the ionizing photon rate. In the grid of hydrodynamic models of Meyer et al. 2021, the combination of outburst luminosity and magnitude (3) places the NGC6334I-MM1 event in the region of moderate total accretion ($\sim$0.1-0.3 $M_\odot$) and hence long duration ($\sim$40-130 yr)., Comment: 11 pages, 3 figures, 2 tables; accepted for publication in The Astrophysical Journal Letters

Class II methanol masers are signs of massive young stellar objects (MYSOs). Recent findings show that MYSO accretion bursts cause flares of these masers. Thus, maser monitoring can be used to identify such bursts. Burst-induced SED changes provide valuable information on a very intense phase of high-mass star formation. In mid-January 2019, a maser flare of the MYSO G358.93-0.03 was reported. ALMA and SMA imaging resolved the core of the star forming region and proved the association of the masers with the brightest continuum source MM1. However, no significant flux rise of the (sub)mm dust continuum was found. Thus, we performed NIR imaging with GROND and IFU spectroscopy with FIFI-LS aboard SOFIA to detect possible counterparts to the (sub)mm sources, and compare their photometry to archival measurements. The comparison of pre-burst and burst SEDs is of crucial importance to judge whether a luminosity increase due to the burst is present and if it triggered the maser flare. The FIR fluxes of MM1 measured with FIFI-LS exceed those from Herschel significantly, which clearly confirms the presence of an accretion burst. The second epoch data, taken about 16 months later, still show increased fluxes. Our RT modeling yielded major burst parameters and suggests that the MYSO features a circumstellar disk which might be transient. From the multi-epoch SEDs, conclusions on heating and cooling time-scales could be drawn. Circumstances of the burst-induced maser relocation have been explored. The verification of the accretion burst from G358 is another confirmation that Class II methanol maser flares represent an alert for such events. The few events known to date already indicate that there is a broad range in burst strength and duration as well as environmental characteristics. The G358 event is the shortest and least luminous MYSO accretion burst so far., Comment: 22 pages, 11 figures, accepted by A&A. Abstract abridged for arxiv submission

This paper reports observations of a 22 GHz water maser `superburst' in the G25.65+1.05 massive star forming region, conducted in response to an alert from the Maser Monitoring Organisation (M2O). Very long baseline interferometry (VLBI) observations using the European VLBI Network (EVN) recorded a maser flux density of $1.2 \times 10^{4}$ Jy. The superburst was investigated in the spectral, structural and temporal domains and its cause was determined to be an increase in maser path length generated by the superposition of multiple maser emitting regions aligning in the line of sight to the observer. This conclusion was based on the location of the bursting maser in the context of the star forming region, its complex structure, and its rapid onset and decay.

The Pr-rich end of the alloy series Pr$_{1-x}$Nd$_x$Os$_4$Sb$_{12}$ has been studied using muon spin rotation and relaxation. The end compound PrOs$_4$Sb$_{12}$ is an unconventional heavy-fermion superconductor, which exhibits a spontaneous magnetic field in the superconducting phase associated with broken time-reversal symmetry. No spontaneous field is observed in the Nd-doped alloys for x $>$ 0.05. The superfluid density is insensitive to Nd concentration, and no Nd$^{3+}$ static magnetism is found down to the lowest temperatures of measurement. Together with the slow suppression of the superconducting transition temperature with Nd doping, these results suggest anomalously weak coupling between Nd spins and conduction-band states., Comment: 17 pages, 6 figures

We report the detection of new 12.178, 12.229, 20.347, and 23.121 GHz methanol masers in the massive star-forming region G358.93-0.03, which are flaring on similarly short timescales (days) as the 6.668 GHz methanol masers also associated with this source. The brightest 12.178 GHz channel increased by a factor of over 700 in just 50 d. The masers found in the 12.229 and 20.347 GHz methanol transitions are the first ever reported and this is only the fourth object to exhibit associated 23.121 GHz methanol masers. The 12.178 GHz methanol maser emission appears to have a higher flux density than that of the 6.668 GHz emission, which is unusual. No associated near-infrared flare counterpart was found, suggesting that the energy source of the flare is deeply embedded., Comment: 9 pages, 3 figures, 2 tables, Published in MNRAS

We present (sub)millimeter imaging at 0.5'' resolution of the massive star-forming region G358.93-0.03 acquired in multiple epochs at 2 and 3 months following the recent flaring of its 6.7 GHz methanol maser emission. Using SMA and ALMA, we have discovered 14 new Class II methanol maser lines ranging in frequency from 199 GHz to 361 GHz, which originate mostly from vt=1 torsionally-excited transitions and include one vt=2 transition. The latter detection provides the first observational evidence that Class II maser pumping involves levels in the vt=2 state. The masers are associated with the brightest continuum source (MM1), which hosts a line-rich hot core. The masers present a consistent curvilinear spatial velocity pattern that wraps around MM1, suggestive of a coherent physical structure 1200 au in extent. In contrast, the thermal lines exhibit a linear pattern that crosses MM1 but at progressive position angles that appear to be a function of either increasing temperature or decreasing optical depth. The maser spectral profiles evolved significantly over one month, and the intensities dropped by factors of 3.0 to 7.2, with the vt=2 line showing the largest change. A small area of maser emission from only the highest excitation lines closest to MM1 has disappeared. There are seven additional dust continuum sources in the protocluster, including another hot core (MM3). We do not find evidence for a significant change in (sub)millimeter continuum emission from any of the sources during the one month interval, and the total protocluster emission remains comparable to prior single dish measurements., Comment: 10 pages, 2 figures, 3 tables, accepted for publication in The Astrophysical Journal Letters. All maser transitions detected in this object can be displayed via an option on the advanced interface of http://splatalogue.net

Centimeter continuum observations of protostellar jets have revealed the presence of knots of shocked gas where the flux density decreases with frequency. This spectrum is characteristic of nonthermal synchrotron radiation and implies the presence of both magnetic fields and relativistic electrons in protostellar jets. Here, we report on one of the few detections of nonthermal jet driven by a young massive star in the star-forming region G035.02$+$0.35. We made use of the NSF's Karl G. Jansky Very Large Array (VLA) to observe this region at C, Ku, and K bands with the A- and B-array configurations, and obtained sensitive radio continuum maps down to a rms of 10 $\mu$Jy beam$^{-1}$. These observations allow for a detailed spectral index analysis of the radio continuum emission in the region, which we interpret as a protostellar jet with a number of knots aligned with extended 4.5 $\mu$m emission. Two knots clearly emit nonthermal radiation and are found at similar distances, of approximately 10,000 au, each side of the central young star, from which they expand at velocities of hundreds km s$^{-1}$. We estimate both the mechanical force and the magnetic field associated with the radio jet, and infer a lower limit of $0.4\times10^{-4} $M$_{\odot}$ yr$^{-1}$ km s$^{-1}$ and values in the range $0.7-1.3 $mG, respectively., Comment: 10 pages, 4 figures, 5 tables, accepted by Astronomy & Astrophysics

We present high resolution ($\sim$300 au) Atacama Large Millimeter/submillimeter Array (ALMA) observations of the massive young stellar object G11.92-0.61 MM 1. We resolve the immediate circumstellar environment of MM 1 in 1.3 mm continuum emission and CH$_{3}$CN emission for the first time. The object divides into two main sources - MM 1a, which is the source of a bipolar molecular outflow, and MM 1b, located 0.57'' (1920 au) to the South-East. The main component of MM 1a is an elongated continuum structure, perpendicular to the bipolar outflow, with a size of $0.141'' \times 0.050''$ ($480\times170$ au). The gas kinematics toward MM 1a probed via CH$_{3}$CN trace a variety of scales. The lower energy $J=12-11$ $K=3$ line traces extended, rotating gas within the outflow cavity, while the $v$8=1 line shows a clearly-resolved Keplerian rotation signature. Analysis of the gas kinematics and dust emission shows that the total enclosed mass in MM 1a is $40\pm5$ M$_{\odot}$ (where between 2.2-5.8 M$_{\odot}$ is attributed to the disk), while MM 1b is $<0.6$ M$_{\odot}$. The extreme mass ratio and orbital properties of MM 1a and MM 1b suggest that MM 1b is one of the first observed examples of the formation of a binary star via disk fragmentation around a massive young (proto)star., Comment: 8 pages, 3 figures, accepted for publication in ApJ Letters

We compare multi-epoch sub-arcsecond VLA imaging of the 22 GHz water masers toward the massive protocluster NGC6334I observed before and after the recent outburst of MM1B in (sub)millimeter continuum. Since the outburst, the water maser emission toward MM1 has substantially weakened. Simultaneously, the strong water masers associated with the synchrotron continuum point source CM2 have flared by a mean factor of 6.5 (to 4.2 kJy) with highly-blueshifted features (up to 70 km/s from LSR) becoming more prominent. The strongest flaring water masers reside 3000 au north of MM1B and form a remarkable bow shock pattern whose vertex coincides with CM2 and tail points back to MM1B. Excited OH masers trace a secondary bow shock located ~120 au downstream. ALMA images of CS (6-5) reveal a highly-collimated north-south structure encompassing the flaring masers to the north and the non-flaring masers to the south seen in projection toward the MM3-UCHII region. Proper motions of the southern water masers over 5.3 years indicate a bulk projected motion of 117 km/s southward from MM1B with a dynamical time of 170 yr. We conclude that CM2, the water masers, and many of the excited OH masers trace the interaction of the high velocity bipolar outflow from MM1B with ambient molecular gas. The previously-excavated outflow cavity has apparently allowed the radiative energy of the current outburst to propagate freely until terminating at the northern bow shock where it strengthened the masers. Additionally, water masers have been detected toward MM7 for the first time, and a highly-collimated CS (6-5) outflow has been detected toward MM4., Comment: 19 pages, 9 figures, accepted for publication by The Astrophysical Journal

As a product of the maser monitoring program with the 26m telescope of the Hartebeesthoek Radio Astronomy Observatory (HartRAO), we present an unprecedented, contemporaneous flaring event of 10 maser transitions in hydroxyl, methanol, and water that began in 2015 January in the massive star-forming region NGC6334I in the velocity range -10 to -2 km/s. The 6.7 GHz methanol and 22.2 GHz water masers began flaring within 22 days of each other, while the 12.2 GHz methanol and 1665 MHz hydroxyl masers flared 80 and 113 days later respectively. The 1665 MHz, 6.7 GHz, and 22.2 GHz masers have all remained in their flared state for nearly 3 years. The brightest flaring components increased by factors of 66, 21, 26, and 20 in the 12.2 and 6.7 GHz methanol, 1665 MHz hydroxyl and 22.2 GHz water maser transitions respectively; some weaker components increased by up to a factor of 145. We also report new maser emission in the 1720, 6031, and 6035 MHz OH lines and the 23.1 GHz methanol line, along with the detection of only the fifth 4660 MHz OH maser. We note the correlation of this event with the extraordinary (sub)millimeter continuum outburst from the massive protostellar system NGC6334I-MM1 and discuss the implications of the observed time lags between different maser velocity components on the nature of the outburst. Finally, we identify two earlier epoch maser flaring events likely associated with this object, which suggest a recurring accretive phenomenon that generates powerful radiative outbursts., Comment: 18 pages, 14 figures, accepted for publication in MNRAS

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., Comment: Accepted for publication in The Astrophysical Journal; 17 pages, 9 figures, 8 tables

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., Comment: 6 pages, 4 figures, accepted for publication in the Astrophysical Journal Letters; corrections from proof applied

We present 1.05 mm ALMA observations of the deeply embedded high-mass protocluster G11.92-0.61, designed to search for low-mass cores within the accretion reservoir of the massive protostars. Our ALMA mosaic, which covers an extent of ~0.7 pc at sub-arcsecond (~1400 au) resolution, reveals a rich population of 16 new millimetre continuum sources surrounding the three previously-known millimetre cores. Most of the new sources are located in the outer reaches of the accretion reservoir: the median projected separation from the central, massive (proto)star MM1 is ~0.17 pc. The derived physical properties of the new millimetre continuum sources are consistent with those of low-mass prestellar and protostellar cores in nearby star-forming regions: the median mass, radius, and density of the new sources are 1.3 Msun, 1600 au, and n(H2)~10^7 cm^-3. At least three of the low-mass cores in G11.92-0.61 drive molecular outflows, traced by high-velocity 12CO(3-2) (observed with the SMA) and/or by H2CO and CH3OH emission (observed with ALMA). This finding, combined with the known outflow/accretion activity of MM1, indicates that high- and low-mass stars are forming (accreting) simultaneously within this protocluster. Our ALMA results are consistent with the predictions of competitive-accretion-type models in which high-mass stars form along with their surrounding clusters., Comment: 17 pages, 5 figures, 5 tables, accepted for publication in MNRAS; v2: proper acknowledgement to NRAO added

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., Comment: 18 pages, 7 figures, 5 tables, accepted for publication in The Astrophysical Journal

The formation process of massive stars is not well understood, and advancement in our understanding benefits from high resolution observations and modelling of the gas and dust surrounding individual high-mass (proto)stars. Here we report sub-arcsecond (<1550 au) resolution observations of the young massive star G11.92-0.61 MM1 with the SMA and VLA. Our 1.3 mm SMA observations reveal consistent velocity gradients in compact molecular line emission from species such as CH$_3$CN, CH$_3$OH, OCS, HNCO, H$_2$CO, DCN and CH$_3$CH$_2$CN, oriented perpendicular to the previously reported bipolar molecular outflow from MM1. Modelling of the compact gas kinematics suggests a structure undergoing rotation around the peak of the dust continuum emission. The rotational profile can be well fit by a model of a Keplerian disc, including infall, surrounding an enclosed mass of 30-60M$_{\odot}$, of which 2-3M$_{\odot}$ is attributed to the disc. From modelling the CH$_3$CN emission, we determine that two temperature components, of 150 K and 230 K, are required to adequately reproduce the spectra. Our 0.9 and 3.0cm VLA continuum data exhibit an excess above the level expected from dust emission; the full centimetre-submillimetre wavelength spectral energy distribution of MM1 is well reproduced by a model including dust emission, an unresolved hypercompact H{\i}{\i} region, and a compact ionised jet. In combination, our results suggest that MM1 is an example of a massive proto-O star forming via disc accretion, in a similar way to that of lower mass stars., Comment: 18 pages, 9 figures. Accepted to MNRAS

There have been several recent detections of candidate Keplerian discs around massive young protostars. Given the relatively large disc-to-star mass ratios in these systems, and their young ages, it is worth investigating their propensity to becoming self-gravitating. To this end, we compute self-consistent, semi-analytic models of putative self-gravitating discs for five candidate disc systems. Our aim is not to fit exactly the observations, but to demonstrate that the expected dust continuum emission from marginally unstable self-gravitating discs can be quite weak, due to high optical depth at the midplane even at millimetre wavelengths. In the best cases, the models produce "observable" disc masses within a factor of <1.5 of those observed, with midplane dust temperatures comparable to measured temperatures from molecular line emission. We find in two cases that a self-gravitating disc model compares well with observations. If these discs are self-gravitating, they satisfy the conditions for disc fragmentation in their outer regions. These systems may hence have as-yet-unresolved low mass stellar companions, and are thus promising targets for future high angular resolution observations., Comment: 9 pages, 20 figures, accepted for publication in MNRAS

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy., Comment: 11 pages, 7 figures, 2 tables; accepted for publication in the Astrophysical Journal Letters; this version with small changes to affiliations

We have modelled Atacama Large Millimeter/sub-millimeter Array (ALMA) long baseline imaging of the strong gravitational lens system H-ATLAS J090311.6+003906 (SDP.81). We have reconstructed the distribution of band 6 and 7 continuum emission in the z=3.042 source and we have determined its kinematic properties by reconstructing CO(5-4) and CO(8-7) line emission in bands 4 and 6. The continuum imaging reveals a highly non-uniform distribution of dust with clumps on scales of 200 pc. In contrast, the CO line emission shows a relatively smooth, disk-like velocity field which is well fit by a rotating disk model with an inclination angle of 40+/-5 degrees and an asymptotic rotation velocity of 320 km/s. The inferred dynamical mass within 1.5kpc is (3.5+/-0.5)x10^{10} M_sol which is comparable to the total molecular gas masses of (2.7+/-0.5)x10^{10} M_sol and (3.5+/-0.6)x10^{10} M_sol from the dust continuum emission and CO emission respectively. Our new reconstruction of the lensed HST near-infrared emission shows two objects which appear to be interacting, with the rotating disk of gas and dust revealed by ALMA distinctly offset from the near-infrared emission. The clumpy nature of the dust and a low value of the Toomre parameter of Q=0.3 suggest that the disk is in a state of collapse. We estimate a star formation rate in the disk of 470+/-80 M_sol/yr with an efficiency 65 times greater than typical low-redshift galaxies. Our findings add to the growing body of evidence that the most infra-red luminous, dust obscured galaxies in the high redshift Universe represent a population of merger induced starbursts., Comment: 11 pages, 7 figures, resubmitted to MNRAS following referee's comments. Includes a correction to the molecular gas mass from CO