The American Astronomical Society (AAS), established in 1899 and based in Washington, DC, is the major organization of professional astronomers in North America. Its membership of about 7,000 individuals also includes physicists, mathematicians, geologists, engineers, and others whose research and educational interests lie within the broad spectrum of subjects comprising contemporary astronomy. The mission of the AAS is to enhance and share humanity's scientific understanding of the universe.
The Institute of Physics (IOP) is a leading scientific society promoting physics and bringing physicists together for the benefit of all. It has a worldwide membership of around 50 000 comprising physicists from all sectors, as well as those with an interest in physics. It works to advance physics research, application and education; and engages with policy makers and the public to develop awareness and understanding of physics. Its publishing company, IOP Publishing, is a world leader in professional scientific communications.
A publishing partnership
Letters
A NEW RADIO RECOMBINATION LINE MASER OBJECT TOWARD THE MonR2 H ii REGION
I. Jiménez-Serra1, A. Báez-Rubio2, V. M. Rivilla2, J. Martín-Pintado2, Q. Zhang1, M. Dierickx1, and N. Patel1
1 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA; ijimenez-serra@cfa.harvard.edu, qzhang@cfa.harvard.edu, mdierickx@cfa.harvard.edu, npatel@cfa.harvard.edu
2 Centro de Astrobiología (CSIC/INTA), Ctra. de Torrejón a Ajalvir km 4, E-28850 Torrejón de Ardoz, Madrid, Spain; ryvendel@gmail.com, jmartin@cab.inta-csic.es, baezra@cab.inta-csic.es
We report the detection of a new radio recombination line (RRL) maser object toward the IRS2 source in the MonR2 ultracompact H ii region. The continuum emission at 1.3 mm and 0.85 mm and the H30α and H26α lines were observed with the Submillimeter Array (SMA) at angular resolutions of ∼05–3''. The SMA observations show that the MonR2-IRS2 source is very compact and remains unresolved at spatial scales ⩽400 AU. Its continuum power spectrum at millimeter wavelengths is almost flat (α = −0.16, with Sν∝να), indicating that this source is dominated by optically thin free–free emission. The H30α and H26α RRL emission is also compact and peaks toward the position of the MonR2-IRS2 source. The measured RRL profiles are double peaked with the H26α line showing a clear asymmetry in its spectrum. Since the derived line-to-continuum flux ratios (∼80 and 180 km s−1 for H30α and H26α, respectively) exceed the LTE predictions, the RRLs toward MonR2-IRS2 are affected by maser amplification. The amplification factors are, however, smaller than those found toward the emission-line star MWC349A, indicating that MonR2-IRS2 is a weakly amplified maser. Radiative transfer modeling of the RRL emission toward this source shows that the RRL masers arise from a dense and collimated jet embedded in a cylindrical ionized wind, oriented nearly along the direction of the line of sight. High-angular resolution observations at submillimeter wavelengths are needed to unveil weakly amplified RRL masers in very young massive stars.
Calculations of the level populations of atomic hydrogen in H ii regions have shown that global population inversions can exist across the Rydberg levels of the hydrogen atom3 (Cillié 1936; Baker & Menzel 1938). Population inversions lead to the formation of recombination lines under non-LTE conditions so that stimulated emission can occur. The general understanding is that most H ii regions in our Galaxy are optically thin at wavelengths ⩾1 cm and stimulated emission is practically negligible (Walmsley 1990). However, toward some ultracompact (UC) H ii regions, the presence of ionized stellar winds modifies the internal electron density structure of these sources (with ne ⩾ 107 cm−3), leading to the formation of optically thick cores where radio recombination line (RRL) masers can form (Martín-Pintado et al. 1989, 1994). RRL maser amplification is therefore expected to be a common phenomenon in UC H ii regions with evidence of stellar winds (Martín-Pintado 2002).
The first RRL maser object was discovered by Martín-Pintado et al. (1989) toward the emission-line star MWC349A. This source has an edge-on disk and a bipolar ionized flow (Cohen et al. 1985; Planesas et al. 1992; Weintroub et al. 2008; Martín-Pintado et al. 2011), whose continuum power-law spectrum (α ∼ 0.6, with Sν∝να; Altenhoff et al. 1981) is characteristic of a constant velocity stellar wind (Olnon 1975). Martín-Pintado et al. (1989) showed that while the RRLs at 3 mm are faint with single Gaussian profiles, the RRLs at 1.3 mm are double peaked with intensity factors of ⩾50 brighter than the lines at 3 mm. Since all RRLs with λ ⩽ 2 mm are largely amplified (Thum et al. 1998), Martín-Pintado et al. (1989) concluded that the RRLs toward MWC349A are masers. It has been proposed that η Carinae and Cepheus A HW2 also show RRL maser emission (Cox et al. 1995; Jiménez-Serra et al. 2011). However, MWC349A still remains as the only RRL maser object firmly detected to date.
In this Letter, we report the detection of a new RRL maser object toward the Monoceros R2 (MonR2) UC H ii region with the SMA (d ∼ 830 pc; Herbst & Racine 1976). MonR2 is a blister-type H ii region (diameter of ∼27''; Massi et al. 1985; Wood & Churchwell 1989) that hosts a cluster of IR sources (e.g., Carpenter et al. 1997). Among these sources, IRS2 is a compact young stellar object (YSO; Álvarez et al. 2004) with a luminosity of ∼5000 L☉ (Howard et al. 1994). This source (hereafter MonR2-IRS2) is responsible for the spherical reflection nebula reported by Aspin & Walther (1990) in the MonR2 UC H ii region. The detection of strong blueshifted asymmetries in the H26α RRL toward this source by the SMA, reveals that RRL maser amplification can form in dense ionized winds toward very young massive stars.
2. OBSERVATIONS
Observations of the H30α line (231.9 GHz) toward the MonR2-IRS2 source were carried out with the SMA4 in the very extended (VEX) configuration in two tracks in single-receiver mode (4 GHz bandwidth per sideband). In addition, the H30α and the H26α (353.6 GHz) RRLs were simultaneously observed in a third track in compact (COM) configuration in dual-receiver mode, which provided 2 GHz of total bandwidth per sideband and receiver. The instrumental parameters of the SMA observations are reported in Table 1. The phase center of the observations was set at α(J2000) = 06h07m4583, δ(J2000) = −06°22'5350. We used a uniform spectral resolution of 0.8 MHz, which provided a velocity resolution of ∼1.1 km s−1 at 231.9 GHz, and of ∼0.7 km s−1 at 353.6 GHz. Data calibration was carried out within the IDL MIR software package, and continuum subtraction, imaging, and deconvolution were done within MIRIAD. The uncertainty in the flux calibration was within 20%.
Table 1. Instrumental Parameters of the SMA Observations
In Figure 1 (left panels), we present the continuum images at 1.3 mm measured toward MonR2-IRS2 in COM and VEX configurations (gray scale: upper and middle panels), and at 0.85 mm in COM (lower panel). The derived parameters of this emission, obtained by performing two-dimensional Gaussian fits, are given in Table 2. The position of the continuum peak is very close to that reported for the IRS2 source in the near-IR and X-rays (Carpenter et al. 1997; Nakajima et al. 2003). The derived angular sizes for the MonR2-IRS2 source indicate that this object is very compact and unresolved in the SMA images (Table 2). Since the continuum peak intensity at 1.3 mm is ∼0.15 Jy beam−1 in both VEX and COM, MonR2-IRS2 does not show any structure at subarcsecond scales (04–05 or 320–400 AU). This is consistent with the results by Álvarez et al. (2004) from near-IR speckle imaging. By smoothing the 1.3 mm and 0.85 mm images in COM to the same angular resolution (∼28), we derive a decreasing spectral index of α = −0.16 (with Sν∝να) for MonR2-IRS2. This spectral index is consistent with optically thin free–free continuum emission. We note that our estimate of the spectral index is not affected by missing flux because the MonR2-IRS2 source is very compact. Indeed, its measured flux at 1.3 mm is the same in the VEX and COM beams (see Table 2), confirming the lack of large-scale structures in this source.
Figure 1. Left panels: integrated intensity images of the H30α line from −35.5 to 47.0 km s−1 (thick contours) observed with the SMA in COM and VEX (upper and middle panels), and of the H26α RRL measured toward MonR2-IRS2 in COM (lower panel). These images are superimposed on the continuum emission observed at 1.3 mm in COM and VEX (gray scale: upper and middle panels) and at 0.85 mm in COM (lower panel). Offsets are given in arcseconds with respect to the phase center of the observations. Note that the spatial scales shown for the VEX images (middle panel) are ∼5 times smaller than those of the COM data (upper and lower panels). The first contour and step level for the 1.3 mm and H30α maps in COM (upper panel) are, respectively, 13.5 (3σ) and 27 mJy beam−1, and 2.4 (3σ) and 2.4 Jy beam−1 km s−1. For the 1.3 mm and H30α maps in VEX (middle panel), the first contour and step level are 4.5 (3σ) and 18 mJy beam−1, and 0.75 (3σ) and 2.25 Jy beam−1 km s−1, respectively. For the 0.85 mm and H26α images in COM (lower panel), these levels are 30 (3σ) and 30 mJy beam−1, and 4.5 (3σ) and 4.5 Jy beam−1 km s−1, respectively. Right panels: spectra of the H30α and H26α RRLs measured with the SMA in COM and VEX toward the MonR2-IRS2 continuum peak. The dashed curves show the individual Gaussian fits of the different velocity components detected toward MonR2-IRS2. The solid lines show the total Gaussian fit of the H30α and H26α RRLs. The vertical dotted line indicates the radial velocity of the source (vLSR = 10.5 km s−1; Torrelles et al. 1983).
Notes. aThe error in the peak flux corresponds to the 1σ noise level in the continuum images.
bThe error in the integrated intensity flux of the RRLs is calculated as σArea = 1, with δv being the velocity resolution in the RRL spectra (see Section 3) and Δv being the line width derived from the Gaussian fit of the RRL emission.
cThe error in the peak intensity of the H30α and H26α lines is given by the 1σ noise level in the RRL images.
Superimposed on the continuum images, Figure 1 reports the integrated intensity emission of the H30α and H26α lines measured toward MonR2-IRS2 from −35.5 km s−1 to 47.0 km s−1 (thick contours). The RRL maps show a compact structure centered at the continuum peak, which is almost unresolved even in the VEX images at angular scales of ∼05. The spectra of the H30α and H26α RRLs extracted from the position of the continuum peak are shown in the right panels of Figure 1. These spectra have been smoothed to a velocity resolution of 2.2 km s−1 for H30α, and of 2.8 km s−1 for H26α. The H30α and H26α lines have bright double-peaked line profiles whose components are redshifted and blueshifted by ∼20–25 km s−1 with respect to the radial velocity of the source (vLSR ∼ 10.5 km s−1; Torrelles et al. 1983). In addition, some faint emission at vLSR ∼ 10.5 km s−1 is also detected, likely associated with the extended MonR2 UC H ii region. This would explain why the double-peaked line profile of H30α is not as clearly seen in COM as in VEX, since the contribution from the extended emission would be larger within the COM beam.
From Figure 1, we also find that while the double-peaked H30α RRLs have similar peak intensities of ∼0.2–0.3 Jy beam−1, the H26α line profile shows a clear asymmetry with the blueshifted gas being factors of ∼2 brighter than the redshifted emission (∼0.9 Jy beam−1 versus ∼0.5 Jy beam−1, respectively). This behavior can only be explained if stimulated emission plays a key role in the formation of the H26α RRL (see Martín-Pintado et al. 1993; Jiménez-Serra et al. 2011; Báez-Rubio et al. 2012). As shown in Section 5, the H26α RRL blueshifted asymmetry is a consequence of the presence of more background emission to be amplified by the blueshifted foreground ionized gas than for the redshifted material, assuming that the RRL masers are formed in an expanding ionized jet.
The observed parameters of the three velocity components of the H30α and H26α RRLs are reported in Table 2. These parameters were derived by fixing the peak velocity and line width of the ionized component at ambient velocities to, respectively, 10.5 km s−1 (the radial velocity of the MonR2-IRS2 source) and 30 km s−1 (the expected thermal line width for RRLs in H ii regions with electron temperatures of Te = 1–2 × 104 K; see, e.g., Keto et al. 2008). The derived peak velocities of the H30α and H26α components are 29–33 km s−1 for the redshifted gas, and −11 to −14 km s−1 for the blueshifted emission. The measured line widths are ∼10–15 km s−1 (Table 2), except for the H30α COM data whose line widths are broader due to the contamination from the extended MonR2 UC H ii region.
By integrating the emission of the H30α RRL from −18 to −4.8 km s−1 and from 27.1 to 38.1 km s−1, the VEX images of the blueshifted and redshifted H30α gas reveal a spatial shift of 0045 (or 360 AU) in the southeast–northwest direction (blueshifted emission toward the southeast, redshifted gas toward the northwest). Although this shift is small compared to the angular resolution of the VEX data (05 × 04; Table 1), higher positional accuracy of (θbeam/(2 × S/N)) (with S/N being the signal-to-noise ratio of the H30α emission) can be achieved thanks to the bright integrated intensities of the blueshifted and redshifted H30α RRL components. By fitting two-dimensional Gaussians to the integrated intensity images of these velocity components, we obtain an accuracy in the H30α Gaussian centroid position of ∼0008. This implies that the detected H30α spatial shift of 0045 ± 0008 is at the 5.5σ confidence level. This shift could be associated with a collimated ionized jet propagating at a velocity of 20 km s−1 (see Section 5).
4. RRL MASERS IN MonR2-IRS2
The double-peaked RRL profiles with clear asymmetries in the low-n RRL transitions toward MonR2-IRS2 (i.e., in the H26 α line; Section 3) resemble the behavior of the RRLs detected toward MWC349A and suggest that the RRL emission toward MonR2-IRS2 at λ ⩽ 1.3 mm is affected by maser amplification. The integrated line-to-continuum flux ratios, ILTRs,5 derived toward this source support this idea (Table 3). As in MWC349A (Martín-Pintado et al. 1989; Thum et al. 1994), the ILTRs derived toward MonR2-IRS2 (80 km s−1 for H30α and 180 km s−1 for H26α) exceed the values predicted under LTE conditions (65 km s−1 and 103 km s−1, respectively). However, unlike MWC349A, the RRL maser amplification factors toward MonR2-IRS2 with respect to LTE are significantly smaller than those found in MWC349A (∼1.2–1.7 for MonR2-IRS2 versus ⩾5–7 for MWC349A; Table 3 and Strelnitski et al. 1996; Thum et al. 1998). This indicates that MonR2-IRS2 is a weakly amplified RRL maser. This is due to the fact that MWC349A has a larger optically thick core than MonR2-IRS2, allowing a much larger maser amplification of the RRLs (see Ponomarev et al. 1994 and Section 5).
Notes. aCalculated for optically thin continuum emission, T*e = 104 K and N(He+)/N(H+) = 0.08.
bEstimated from the ratio between the total area of the H30α and H26α lines and the continuum peak fluxes at 1.3 mm and 0.85 mm (Table 2).
cFrom Martín-Pintado et al. (1989).
dDerived from the average H26α line flux measured by Thum et al. (1994) and from the expected free–free continuum flux at 0.85 mm assuming a spectral index of α ∼ 0.6 (Altenhoff et al. 1981).
By comparing the H30α and H26α RRLs measured toward MonR2-IRS2 with radiative transfer modeling of the RRL emission, we can constrain the physical structure and kinematics of the ionized gas toward this YSO. To do this, we have used the three-dimensional RRL radiative transfer code MORELI developed by Martín-Pintado et al. (2011) and Báez-Rubio et al. (2012). In our calculations, we have used the departure coefficients, bn and βn, derived by Walmsley (1990). For the geometry of the source, we have assumed that MonR2-IRS2 has a central mass of 13 M☉, which is consistent with a B1 star on the zero-age main sequence (luminosity of ∼5000 L☉; Panagia 1973) and with the Lyα photon flux (∼6 × 1045 photons s−1) derived from the 1.3 mm and 0.85 mm continuum data (Table 2). In addition, we consider that MonR2-IRS2 powers a cylindrical, isothermal (Te = 104 K) ionized outflow, which expands at a velocity of 20 km s−1. The cylindrical outflow is defined so that its outer radius is 100 AU and its electron density distribution ne decreases as r−2, with r being the distance to the star and with ne = 108 cm−3 at a distance r ∼7 AU. The total length of the cylinder is L = 520 AU and its axis is oriented along the direction of the line of sight. The assumed radial velocity for the MonR2-IRS2 source is vLSR = 10 km s−1.
To reproduce the optically thin free–free continuum emission of MonR2-IRS2 at 1.3 mm and 0.85 mm, we consider a cavity around the star with radius r = 13 AU (see Báez-Rubio et al. 2012). The predicted continuum fluxes at 1.3 mm and 0.85 mm are 0.158 Jy beam−1 and 0.151 Jy beam−1, respectively, which differ from the observed values by less than 5%.
In Figure 2 (left panels), we compare the RRL profiles of the H30α and H26α lines measured with the SMA toward the MonR2-IRS2 source (in VEX and COM, respectively), with the synthetic line profiles predicted by the MORELI code for the cylindrical, ionized outflow assumed for MonR2-IRS2 (Case a). Although our model reproduces the double-peaked profile of the H30α RRL relatively well, it clearly fails to predict the asymmetry detected in the H26α line (Figure 2). Indeed, the predicted line profiles for both the H30α and H26α RRLs are symmetric because the continuum emission of MonR2-IRS2 is optically thin throughout the ionized outflow, which prevents the strong amplification and the asymmetries of the RRLs.
Figure 2. Comparison of the H30α (upper panels) and H26α (lower panels) RRL spectra measured with the SMA toward MonR2-IRS2 (in VEX and COM, respectively), with those predicted by the MORELI code for two geometries of the MonR2-IRS2 source: (a) a cylindrical, ionized outflow oriented along the line of sight (left panels) and (b) a cylindrical, ionized outflow with an elongated and dense inhomogeneity resembling a collimated jet (right panels). The vertical dotted lines indicate the assumed radial velocity of the source vLSR = 10 km s−1.
In order to explain the strongly asymmetric H26α line profile, larger electron densities are required for the inner regions of the ionized wind so that the continuum emission becomes locally optically thick. To simulate this, we have assumed two elongated inhomogeneities with ne = 108 cm−3 located at the axis of the blueshifted and redshifted lobes of the cylindrical ionized outflow, resembling a collimated ionized jet. The size of the inhomogeneities is 2 AU× 2 AU× 60 AU, and they are placed on the axis of the cylindrical outflow at a distance of 43 AU from the central star.
In Figure 2 (right panels), we report the H30α and H26α line profiles predicted by our model with a cylindrical ionized outflow powered by a dense, collimated jet (Case b). While the predicted H30α emission shows a similar line profile to that obtained in Case (a) for only a cylindrical ionized outflow, the model with the ionized outflow+jet perfectly matches the asymmetry observed in the H26α line profile. This is due to the fact that the inner regions of the jet are optically thick, allowing the maser amplification of the H26α line. Since the optically thick region in the plane of the sky of the ionized jet is very small (only 2 AU× 2 AU), its contribution to the total optical depth of the MonR2-IRS2 continuum source is practically negligible. We note that only the H26α RRL gets amplified because the departure coefficient βn of the H26α transition at an electron density of ne = 108 cm−3 is significantly larger than that of the H30α line (see Figure 8 in Strelnitski et al. 1996). The electron densities ne assumed for the elongated inhomogeneities do not likely exceed 108 cm−3, because the βn coefficients for both the H30α and H26α RRLs get close to zero, hindering the RRL maser amplification (see Figure 5 in Strelnitski et al. 1996).
In our model for MonR2-IRS2 with an ionized outflow+jet (Case b), the elongation of the dense inhomogeneities is required to sufficiently amplify the blueshifted peak of the H26α emission. In contrast with the blueshifted peak, the RRL emission from the redshifted lobe of the ionized jet does not get amplified because the optically thick continuum acts as a screen for the redshifted H26α emission.
Finally, we stress that other geometries for the MonR2-IRS2 source have been explored in our study. In the case of an edge-on Keplerian rotating disk with a biconical ionized flow (similar to that assumed for MWC349A; see Martín-Pintado et al. 2011), the model also predicts double-peaked line profiles for large opening angles of the ionized wind (θa ⩾ 60°; see Báez-Rubio et al. 2012, for the definition of this angle). However, the model not only fails to reproduce the blueshifted asymmetry of the H26α line profile but also predicts an RRL emission excess at velocities close to the ambient cloud velocity that is not observed. The presence of high-density inhomogeneities could partially alleviate the discrepancies between the model predictions and the observations, but the best fit is obtained with the cylindrical ionized outflow powered by the elongated jet.
In summary, we report the detection of a new RRL maser object toward the IRS2 source in the MonR2 UC H ii region. Our SMA images reveal that MonR2-IRS2 is a compact object with an spectral index α = −0.16, characteristic of optically thin free–free emission. The line profiles of the H30α and H26α RRLs are double peaked, and their derived ILTRs clearly exceed those predicted under LTE conditions. The RRL emission at λ ⩽ 1.3 mm toward the MonR2-IRS2 source are weakly amplified masers. Our radiative transfer modeling of the RRLs toward the MonR2-IRS2 source suggests that the RRL masers arise from a dense and highly collimated jet embedded in a ionized, cylindrical outflow, nearly oriented along the direction of the line of sight. Interferometric observations at submillimeter wavelengths have the potential to unveil a population of weakly amplified RRL maser objects in UC H ii regions with high-velocity ionized winds.
We are grateful to the SMA director and staff for letting us carry out part of the SMA observations under the Harvard Astronomy Ay191 course, and to the Ay191 students for their enthusiasm during this course. We also acknowledge the anonymous referee for the constructive comments to the manuscript. I.J.-S. acknowledges the Smithsonian Astrophysical Observatory for the support provided through an SMA fellowship. J.M.-P. and I.J.-S. have been partially funded by MICINN grants ESP2007-65812-C02-C01 and AYA2010-21697-C05-01 and AstroMadrid (CAM S2009/ESP-1496).
Facility: SMA - SubMillimeter Array
Footnotes
3
The spontaneous decay rate decreases with quantum number n as n−5, leading to an underpopulation of the lower-n levels of the hydrogen atom (see, e.g., Strelnitski et al. 1996).
4
The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and the Academia Sinica.
5
ILTR is defined as ΔvTL/TC, and depends on the RRL frequency as ν1.1 for LTE and optically thin emission (see, e.g., Martín-Pintado et al. 1989).