No Arabic abstract
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.9pm0.8)times10^4 L_odot$. Accounting for an estimated pre-outburst luminosity ratio MM1:MM2 = $2.1pm0.4$, the luminosity of MM1 has increased by a factor of $16.3pm4.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).
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 $70pm20$, 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.
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.
Following an eruptive accretion event in NGC6334I-MM1, flares in the various maser species, including water masers, were triggered. We report the observed relative proper motion of the highly variable water masers associated with the massive star-forming region, NGC6334I. High velocity H$_2$O maser proper motions were detected in 5 maser clusters, CM2-W2 (bow-shock structure), MM1-W1, MM1-W3, UCHII-W1 and UCHII-W3. The overall average of the derived relative proper motion is 85 km s$^{-1}$. This mean proper motion is in agreement with the previous results from VLA multi-epoch observations. Our position and velocity variance and co-variance matrix analyses of the maser proper motions show its major axis to have a position angle of $-$79.4$^circ$, cutting through the dust cavity around MM1B and aligned in the northwest-southeast direction. We interpret this as the axis of the jet driving the CM2 shock and the maser motion. The complicated proper motions in MM1-W1 can be explained by the combined influence of the MM1 northeast-southwest bipolar outflow, CS(6-5) north-south collimated bipolar outflow, and the radio jet. The relative proper motions of the H$_2$O masers in UCHII-W1 are likely not driven by the jets of MM1B protostar but by MM3-UCHII. Overall, the post-accretion burst relative proper motions of the H$_2$O masers trace shocks of jet motion.
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.
Stochastic heating of small grains is often mentioned as a primary cause of large infrared (IR) fluxes from star-forming galaxies, e.g. at 24mu m. If the mechanism does work at a galaxy-wide scale, it should show up at smaller scales as well. We calculate temperature probability density distributions within a model protostellar core for four dust components: large silicate and graphite grains, small graphite grains, and polycyclic aromatic hydrocarbon particles. The corresponding spectral energy distributions are calculated and compared with observations of a representative infrared dark cloud core. We show that stochastic heating, induced by the standard interstellar radiation field, cannot explain high mid-IR emission toward the centre of the core. In order to reproduce the observed emission from the core projected centre, in particular, at 24mu m, we need to increase the ambient radiation field by a factor of about 70. However, the model with enhanced radiation field predicts even higher intensities at the core periphery, giving it a ring-like appearance, that is not observed. We discuss possible implications of this finding and also discuss a role of other non-radiative dust heating processes.