THE DETECTION OF A HOT MOLECULAR CORE IN THE LARGE MAGELLANIC CLOUD WITH ALMA

, , , and

Published 2016 August 9 © 2016. The American Astronomical Society. All rights reserved.
, , Citation Takashi Shimonishi et al 2016 ApJ 827 72 DOI 10.3847/0004-637X/827/1/72

0004-637X/827/1/72

ABSTRACT

We report the first detection of a hot molecular core outside our Galaxy based on radio observations with ALMA toward a high-mass young stellar object (YSO) in a nearby low metallicity galaxy, the Large Magellanic Cloud (LMC). Molecular emission lines of CO, C17O, HCO+, H13CO+, H2CO, NO, SiO, H2CS, 33SO, 32SO2, 34SO2, and 33SO2 are detected from a compact region (∼0.1 pc) associated with a high-mass YSO, ST11. The temperature of molecular gas is estimated to be higher than 100 K based on rotation diagram analysis of SO2 and 34SO2 lines. The compact source size, warm gas temperature, high density, and rich molecular lines around a high-mass protostar suggest that ST11 is associated with a hot molecular core. We find that the molecular abundances of the LMC hot core are significantly different from those of Galactic hot cores. The abundances of CH3OH, H2CO, and HNCO are remarkably lower compared to Galactic hot cores by at least 1–3 orders of magnitude. We suggest that these abundances are characterized by the deficiency of molecules whose formation requires the hydrogenation of CO on grain surfaces. In contrast, NO shows a high abundance in ST11 despite the notably low abundance of nitrogen in the LMC. A multitude of SO2 and its isotopologue line detections in ST11 imply that SO2 can be a key molecular tracer of hot core chemistry in metal-poor environments. Furthermore, we find molecular outflows around the hot core, which is the second detection of an extragalactic protostellar outflow. In this paper, we discuss the physical and chemical characteristics of a hot molecular core in the low metallicity environment.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

Because cosmic metallicity is increasing in time with the evolution of our universe, interstellar chemistry in low metallicity environments is crucial to understand chemical processes in the past universe. For this purpose, observations of chemically rich objects in nearby low metallicity galaxies and comparative studies with Galactic counterparts play an important role.

Hot molecular cores are one of the early stages of high-mass star formation, and they play a key role in the formation and evolution of complex molecules in space. In terms of physical properties, hot cores are defined as having a small source size (≤0.1 pc), a high density (≥106 cm−3), and a warm gas/dust temperature (≥100 K) (e.g., Kurtz et al. 2000; van der Tak 2004). The chemistry of hot cores is characterized by sublimation of ice mantles, which accumulated in the course of star formation. In cold molecular clouds and prestellar cores, gaseous molecules and atoms are frozen onto dust grains and hydrogenated. As the core is heated by star formation activities, reactions among heavy species become active on the grain surfaces, forming larger molecules. In addition, sublimated molecules, such as CH3OH and NH3, are subject to further gas-phase reactions (e.g., Garrod & Herbst 2006; Garrod et al. 2008; Herbst & van Dishoeck 2009). As a result, hot cores show a wealth of molecular spectral lines in the infrared and radio wavelengths. Thus detailed studies of the chemical properties of hot cores are crucial to understand the complex chemical processes triggered by star formation.

The Large Magellanic Cloud (LMC) is an excellent target to study interstellar and circumstellar chemistry in different metallicity environments owing to its proximity (49.97 ± 1.11 kpc; Pietrzyński et al. 2013) and low metallicity (about one-third of the solar neighborhood; Westerlund 1990). The low dust content in the galaxy results in a harsh radiation environment, and thus the photoprocessing of the interstellar medium (ISM) should be more effective in the LMC than in our Galaxy (Israel et al. 1986). Furthermore, according to gamma-ray observations, the cosmic-ray density in the LMC is estimated to be lower than the Galactic typical values by a factor of four (Abdo et al. 2010). It is therefore highly anticipated that these environmental differences should affect chemical processes, and hot cores in the LMC should provide key information to understand the chemistry, particularly that of complex molecules, in low metallicity environments. However, so far, observations of hot cores have been limited to Galactic sources due to a lack of spatial resolution and sensitivity of radio telescopes.

Most of radio studies on the chemical compositions of molecular gas in the LMC have been performed with single-dish telescopes. Early studies by the SEST 15 m telescope conducted multiline observations toward H ii regions in the LMC and detected molecular species as large as CH3OH, C3H2, and SO2 (Johansson et al. 1994; Chin et al. 1997; Heikkilä et al. 1999; Wang et al. 2009). Submillimeter observations of relatively dense molecular gas in star-forming regions in the LMC are also reported (Paron et al. 2014, 2016). Nishimura et al. (2016) recently conducted deep and unbiased spectral line surveys in the 3 mm window toward a number of molecular clouds in the LMC using the Mopra telescope. They reported low abundances of nitrogen-bearing molecules, a deficiency of CH3OH, and a high abundance of C2H in the LMC compared to Galactic molecular clouds.

The characteristic interstellar chemistry in the LMC is also suggested in previous infrared observations of ices around embedded young stellar objects (YSOs) in the LMC. Shimonishi et al. (2008, 2010) reported that the CO2/H2O ice ratio of high-mass YSOs in the LMC is systematically higher than those of Galactic high-mass YSOs, based on infrared observations with the AKARI satellite. Recently, Shimonishi et al. (2016) reported that the CH3OH ice around high-mass YSOs in the LMC is less abundant compared with Galactic counterparts, based on infrared observations with the Very Large Telescope. The authors suggest that warm ice chemistry (grain surface reactions at a relatively high dust temperature) is responsible for the observed characteristics of ice chemical compositions in the LMC. Furthermore, detailed studies of the 15.2 μm CO2 ice band toward LMC high-mass YSOs with Spitzer suggest a higher degree of thermal processing of ices in the LMC than in our Galaxy (Oliveira et al. 2009; Seale et al. 2011). Since gas–grain chemistry is believed to play an important role in hot cores, the characteristic ice chemistry in the LMC would imply a diverse hot core chemistry in extragalactic environments according to their metallicities.

High spatial resolution interferometry observations toward star-forming regions in the LMC have been reported from the Australia Telescope Compact Array (e.g., Wong et al. 2006; Ott et al. 2008; Seale et al. 2012; Anderson et al. 2014) and recently from the Atacama Large Millimeter/submillimeter Array (ALMA) (e.g., Indebetouw et al. 2013; Fukui et al. 2015). These observations resolve star-forming regions in the LMC down to the parsec- or sub-parsec-scale and investigate the physical properties of dense molecular gas. However, the chemical properties of warm and dense molecular gas associated with a single high-mass YSO in the LMC still remain to be investigated.

In this paper, we report the detection of a hot molecular core in the LMC based on submillimeter interferometric observations with ALMA. In Section 2 we describe the observation details and the data reduction conducted in this work. The obtained images and spectra of the molecular line emission and continuum are presented in Section 3. Analysis of the spectral line data and the derivation of the physical quantities of molecular gas and dust are described in Section 4. The physical and chemical properties of the observed source are discussed in Section 5. Finally, the conclusions of this paper are summarized in Section 6.

2. OBSERVATIONS AND DATA REDUCTION

2.1. Target

The target of the present observations is a high-mass YSO, 2MASS J05264658-6848469 or ST11 (hereafter ST11), located in the LMC. The source has been spectroscopically identified to be a high-mass YSO in previous infrared studies (Seale et al. 2009; Shimonishi et al. 2010). The detailed YSO properties of ST11 are revisited in this work and are discussed in Section 4.5.

2.2. Observations

Observations were carried out with ALMA between 2013 November and 2014 February as a part of the Cycle 1 high priority program 2012.1.01108.S (PI T. Shimonishi). The telescopes were pointed to R.A. = 05h26m46fs63 and decl. = −6848'47farcs10 (J2000), which is the position of ST11 measured by infrared data (Shimonishi et al. 2010). The compact configurations C32-3 and partially C32-2 were used for our observations. The target object was observed in Band 7 with seven spectral bands covering 336.60–337.55, 338.05-339.00, 345.35–346.30, 346.55–347.50, 350.25–351.20, 351.30–352.25, and 356.30–357.25 GHz in the rest-frame frequency. The observed lines and the integration times are presented in Table 1. The primary beam has a full-width at half-maximum (FWHM) of 18''–19'' at these frequencies, which corresponds to the field of view (FOV) of the imaging data. The maximum recoverable angular scale is about 7''. The velocity resolution of the original data is 0.4 km s−1 for the spectral bands and 26 km s−1 for the continuum band.

Table 1.  Line Parameters

Molecule Transition Eu/k Frequency Tmb ΔV $\int {T}_{\mathrm{mb}}{dV}$ VLSR rms Itimea Notes
    (K) (GHz) (K) (km s−1) (K km s−1) (km s−1) (K) (s)  
CO 3–2 33 345.7960 ∼65 ∼16 ∼1150 ∼252 0.42 (1)
C17O 3–2 32 337.0611 2.30 ± 0.08 5.6 13.77 ± 0.98 251.1 0.21
HCO+ 4–3 43 356.7342 29.87 ± 0.08 8.5 270.86 ± 1.73 250.6 0.21
H13CO+ 4–3 42 346.9983 1.59 ± 0.06 7.3 12.34 ± 1.07 250.5 0.20
NO 7/2, 9/2–5/2, 7/2 f 36 350.6895 1.63 ± 0.06 7.2 12.53 ± 0.97 250.6 0.19 (2)
  7/2, 7/2–5/2, 5/2 f 36 350.6908
  7/2, 5/2–5/2, 3/2 f 36 350.6948
NO 7/2, 9/2–5/2, 7/2 e 36 351.0435 0.85 ± 0.08 3.8 3.44 ± 0.69 251.7 0.19
NO 7/2, 7/2–5/2, 5/2 e 36 351.0515 0.71 ± 0.07 6.5 4.86 ± 0.97 251.7 0.19 (3)
  7/2, 5/2–5/2, 3/2 e 36 351.0517
H2CO 51,5–41,4 62 351.7686 2.48 ± 0.07 6.5 17.10 ± 1.02 250.7 0.20
CH3OH 70–60 A+ 65 338.4087 <0.30 <1.9 0.15
HNCO 160,16–150,15 143 351.6333 <0.28 <1.8 0.14
SiO 8–7 75 347.3306 1.64 ± 0.06 7.1 12.32 ± 1.05 250.4 0.20
C34S 7–6 65 337.3965 <0.30 <1.9 0.15
H2CS 101,10–91,9 102 338.0832 0.44 ± 0.07 3.2 1.49 ± 0.50 249.6 0.21
CH3OCH3 84,5–73,4 AE 55 356.5753 <0.30 <1.9 0.15 (3)
  84,5–73,4 EE 55 356.5760
C2H5OH 104,7–93,6 66 357.0674 <0.30 <1.9 0.15
33SO 87–76 81 337.1986 2.84 ± 0.07 5.8 17.63 ± 1.00 250.9 0.21 (4)
HC3N 38–37 323 345.6090 <0.27 <1.7 0.13
HCOOCH3 99,0–88,1 A 80 345.7187 <0.27 <1.7 0.13 (3)
  99,1–88,0 A 80 345.7187
SO2 167,9–176,12 245 336.6696 1.87 ± 0.07 7.8 15.50 ± 1.19 250.6 0.21
SO2 184,14–183,15 197 338.3060 5.25 ± 0.07 6.2 34.54 ± 1.08 250.5 0.21
SO2 201,19–192,18 199 338.6118 6.39 ± 0.07 6.7 45.67 ± 1.02 250.8 0.21
SO2 269,17–278,20 521 345.4490 0.76 ± 0.07 3.4 2.76 ± 0.53 252.0 0.19
SO2 191,19–180,18 168 346.6522 7.99 ± 0.07 6.5 54.93 ± 1.00 250.7 0.20
SO2 106,4–115,7 139 350.8628 2.43 ± 0.06 5.8 15.10 ± 0.85 250.8 0.19
SO2 144,10–143,11 136 351.8739 6.62 ± 0.07 6.6 46.24 ± 1.01 250.7 0.20
SO2 104,6–103,7 90 356.7552 6.19 ± 0.07 7.0 46.21 ± 1.18 250.8 0.21 (5)
SO2 134,10–133,11 123 357.1654 6.03 ± 0.07 6.0 38.25 ± 0.95 250.8 0.21
SO2 ν2 = 1 53,3–42,2 781 357.0872 0.50 ± 0.07 3.0 1.56 ± 0.49 251.5 0.21
33SO2 53,3–42,2 35 346.5901 0.57 ± 0.06 6.2 3.76 ± 1.62 249.7 0.20 (6)
33SO2 134,10–133,11 122 350.7881 0.42 ± 0.05 4.1 1.83 ± 0.55 251.1 0.19 (7)
33SO2 154,12–153,13 149 350.9146 0.40 ± 0.06 4.5 1.89 ± 0.61 249.7 0.19 (7)
33SO2 114,8–113,9 99 350.9951 0.64 ± 0.05 5.3 3.64 ± 0.75 250.9 0.19 (7)
33SO2 54,2–53,3 52 351.6351 0.48 ± 0.06 4.2 2.16 ± 0.59 250.2 0.20 (7)
33SO2 174,14–173,15 179 351.7449 <0.80 <2.4 0.20 (7)
34SO2 132,12–121,11 92 338.3204 0.87 ± 0.06 5.6 5.14 ± 0.83 251.4 0.21
34SO2 144,10–143,11 134 338.7857 0.98 ± 0.06 6.9 7.26 ± 1.07 250.2 0.21
34SO2 74,4–73,5 64 345.5197 1.01 ± 0.07 6.1 6.59 ± 0.93 250.9 0.19
34SO2 64,2–63,3 57 345.5531 0.67 ± 0.06 5.6 4.03 ± 0.77 251.2 0.19
34SO2 54,2–53,3 52 345.6513 0.61 ± 0.05 5.6 3.66 ± 0.84 251.4 0.19
34SO2 44,0–43,1 47 345.6788 0.64 ± 0.08 3.7 2.49 ± 0.72 250.6 0.19
34SO2 174,14–173,15 179 345.9293 0.59 ± 0.05 5.0 3.16 ± 0.67 249.7 0.19
34SO2 282,26–281,27 391 347.4831 0.23 ± 0.06 3.3 0.78 ± 0.54 250.7 0.20
34SO2 214,18–213,19 250 352.0829 0.68 ± 0.07 3.6 2.62 ± 0.57 251.0 0.20
34SO2 200,20–191,19 185 357.1022 1.32 ± 0.07 5.2 7.27 ± 0.80 250.4 0.21

Notes. Uncertainties and upper limits are of 2σ level and do not include systematic errors due to baseline subtraction and adopted spectroscopic constants. (1) Saturated. (2) Blended with three hyperfine components. (3) Blended with two hyperfine components. (4) Blended with four hyperfine components. (5): Partly blended with HCO+(4–3), which is subtracted when fitting this SO2 line. (6) Blended with seven hyperfine components. (7) Blended with ten hyperfine components.

aTotal on-source integration time, where † represents 575 s and ‡ 1452 s.

Download table as:  ASCIITypeset image

2.3. Data Reduction

Raw interferometric data is processed using the Common Astronomy Software Applications (CASA) package. The calibration is performed using CASA 4.1.0, and imaging as well as spectral extraction are performed using CASA 4.3.0. The flux calibrator is J0519-454 and the phase calibrators are J0601-7036 and J0635-7516. The synthesized beam size in the 338 GHz region is approximately 0farcs5 × 0farcs5, which corresponds to 0.12 pc at the assumed distance to the LMC (49.97 kpc, Pietrzyński et al. 2013). The primary beam correction is done by the impbcor task in CASA, but the correction has little effect on the extracted spectra since the target source is located at the center of the FOV and is very compact.

The spectra as well as the continuum flux are extracted from the circular region with a diameter of 0farcs5 centered at R.A. = 05h26m46fs60 and decl. = −6848'47farcs03 (J2000). This region corresponds to the peak position of the 359 GHz (840 μm) continuum emission of ST11 measured in this study, and the diameter corresponds to the beam size. The continuum emission is subtracted from the spectral data using the uvcontsub task in CASA. Several channels are concatenated during the clean process to increase signal-to-noise ratio (S/N) and the channel spacing of the reduced data is 1.5 km s−1 (1.73 MHz) except for CO(3–2). For the CO(3–2) line, the spectral region is not binned since the line is sufficiently strong, and the channel spacing is 0.4 km s−1 (0.46 MHz).

3. RESULTS

3.1. Observed Spectra

Figure 1 shows the spectra extracted from the 0farcs5 diameter region centered at ST11. In the figures the sky frequency is converted to the rest frequency using the LSR velocity of 250.5 km s−1, which is the typical radial velocity of molecular lines detected toward the source. Spectral lines are identified with the aid of the Cologne Database for Molecular Spectroscopy6 (CDMS; Müller et al. 2001, 2005) and the molecular database of the Jet Propulsion Laboratory7 (JPL; Pickett et al. 1998). Molecular emission lines due to CO, C17O, HCO+, H13CO+, H2CO, NO, SiO, H2CS, 33SO, 32SO2, 34SO2, and 33SO2 are detected from a compact region associated with ST11. A number of high excitation lines (Eu > 100 K) are detected for SO2 and its isotopologues. The emission lines of the above molecules in the 345 GHz band are for the first time detected toward the LMC source, except for CO and HCO+. The emission lines of CH3OH, HNCO, CS, HC3N, and complex organic molecules are not detected. Unidentified lines are labeled in the figures, but some of them may be spurious signals.

Figure 1.

Figure 1. ALMA band 7 spectra of ST11 extracted from the central region with a radius of 0farcs5. Detected emission lines are labeled. Tentative detections are indicated by "?" and unidentified lines by "U." The adopted source velocity is 250.5 km s−1.

Standard image High-resolution image

3.2. Synthesized Images

Figures 23 show synthesized images of the continuum and molecular emission lines observed toward ST11. For SO2 and 34SO2, we only show representative lines because the distributions of the emissions are similar in different transitions. The lines of SO2 (ν2 = 1), 33SO2, and H2CS are not included in the figure because they are too weak to visualize the brightness distribution. The images are constructed by integrating each spectrum in the velocity range where the emission line is seen, typically between 240 and 260 km s−1. For CO, the spectrum is integrated between 230 and 285 km s−1 because the line is very broad.

Figure 2.

Figure 2. Flux distribution of the ALMA 840 μm continuum data tracing cold dust (left) and the Gemini/T-ReCS mid-infrared 10 μm image tracing warm dust (right).

Standard image High-resolution image
Figure 3.

Figure 3. Integrated intensity distributions of the CO, C17O, HCO+, H13CO+, H2CO, SiO, NO, 33SO, SO2, and 34SO2 lines. Contours represent the distribution of the 840 μm continuum, and the contour levels are 25%, 50%, and 75% of the peak flux. The synthesized beam size (0farcs5, 0.12 pc at the LMC) is shown by the gray filled circle in each panel.

Standard image High-resolution image

A high spatial resolution mid-infrared image of ST11 is also shown in Figure 2. The image is obtained by T-ReCS at the Gemini South telescope (Program ID: S10B-120, PI: T. Shimonishi) and a broad-band filter in the N-band (7.70–12.97 μm, centered at 10.36 μm) is used for the observation. The ALMA 840 μm continuum traces the distribution of cold dust, while the mid-infrared image traces warm dust.

The source is compact in general and the peak positions of each emission line coincide with the region of dust continuum emission. We estimate a source size (Gaussian FWHM, θ) for the continuum and relatively strong emission lines using a two-dimensional Gaussian fit. The estimated θ in the 840 μm continuum image is about 0farcs6, which is close to the beam size. The mid-infrared emission shows θ = 0farcs5, which is as compact as the 840 μm continuum. The molecular lines of NO, SiO, 33SO, SO2, and 34SO2 typically have θ ∼ 0farcs5, which is indistinguishable from the beam size. Since the distribution is as compact as the beam size, these emissions are considered to be a point source with the present spatial resolution. The H13CO+ and H2CO lines are slightly extended compared with the above lines and have θ ∼ 0farcs8. The lines of CO, C17O, and HCO+ are more extended and θ is 1farcs1–1farcs4. Because the source size is sufficiently smaller than the maximum recoverable angular scale of ∼7'', the emission from ST11 is almost recovered by the present interferometric observations.

4. ANALYSIS

4.1. Spectral Fitting

The line parameters are measured by fitting a single Gaussian profile to the observed lines. For the NO lines at 351.0435 and 351.0515 GHz, we fit a double Gaussian because they are partially blended. In some cases we subtract a local baseline, which is estimated from adjacent line-free regions to correct for weak baseline ripples. For the SO2(104,6–103,7) line, we subtracted the HCO+(4–3) line upon fitting since they are partly blended. In general, good fits are obtained with Gaussian profiles except for the CO(3–2) line, which deviates from a Gaussian. We estimate a peak main-beam brightness temperature, a FWHM, a LSR velocity, and an integrated intensity for each line on the basis of the fitting. For CO, instead, we estimate the peak brightness temperature and the FWHM by visual inspection, and the integrated intensity is estimated by integrating the spectrum in the velocity range between 230 and 285 km s−1. The spectra and the results of the Gaussian fitting are shown in Figures 47. The figures also show the spectral regions of several important non-detection lines. The measured line parameters are summarized in Table 1.

Figure 4.

Figure 4. Spectra of CO, C17O, HCO+, H13CO+, H2CO, SiO, NO, H2CS, and 33SO emission lines extracted from the 0farcs5 diameter region centered at ST11. The blue lines represent Gaussian profiles fitted to the observed spectra. The spectral regions of important non-detection lines including CH3OH, HNCO, C34S, CH3OCH3, C2H5OH, HC3N, and HCOOCH3 are also shown.

Standard image High-resolution image
Figure 5.

Figure 5. Spectra of SO2 emission lines extracted from the 0farcs5 diameter region centered at ST11. The blue lines represent a Gaussian profile fitted to the observed spectra. The spectra are sorted in ascending order of the upper state energy (the emission line with the lowest upper state energy is shown in the upper left panel and that with the highest energy is in the lower right panel). The bottom panel is for the SO2 (ν2 = 1) line.

Standard image High-resolution image
Figure 6.

Figure 6. Spectra of 34SO2 emission lines observed toward ST11 as in Figure 5. The spectra are sorted in ascending order of the upper state energy from the upper left to the lower right.

Standard image High-resolution image
Figure 7.

Figure 7. Spectra of 33SO2 emission lines observed toward ST11 as in Figure 5. The spectra are sorted in ascending order of the upper state energy from the upper left to the lower right. The 174,14–173,15 line is a tentative detection.

Standard image High-resolution image

Multiple hyperfine components sometime exist within the fitted profile according to the spectroscopic catalogs (CDMS or JPL), but these are not resolved due to the low spectral resolution of the present spectra. When the blended lines have comparable upper state energies, we split the measured flux according to the their 2 values and the upper state degeneracy (see Section 4.2 for definitions of S and μ). Using this method, we estimate the line parameters of the strongest hyperfine line, which is used in the subsequent analysis of column densities and rotational temperatures.

4.2. Rotation Diagram Analysis for SO2, 34SO2, and 33SO2

Since we detect multiple emission lines with different excitation energies for SO2, 34SO2, and 33SO2, we perform the rotation diagram analysis assuming an optically thin condition and an local thermodynamic equilibrium (LTE). A column density of molecules in the upper energy level, Nuthin, is derived by the following equation for optically thin lines (e.g., Sutton et al. 1995; Goldsmith & Langer 1999),

Equation (1)

where gu is the degeneracy of the upper level, k is the Boltzmann constant, $\int {T}_{\mathrm{mb}}{dV}$ is the integrated intensity as estimated from the observations, ν is the transition frequency, S is the line strength, and μ is the dipole moment. Under the LTE condition, the total column density, Ntotal, is given by

Equation (2)

where $Q({T}_{\mathrm{rot}})$ is the partition function, Trot is the rotational temperature, and Eu is the upper state energy. This equation is rearranged as follows.

Equation (3)

When ${N}_{u}^{\mathrm{thin}}/{g}_{u}$ is plotted against ${E}_{u}/k$ and data points are fitted by a straight line, the slope and the intercept correspond to Trot and Ntotal, respectively. Thus we can simultaneously determine the rotational temperature and the total column density. All the spectroscopic parameters required in the above analysis are extracted from the CDMS or the JPL database. For the partition function, we interpolate the data given in the databases and estimate the appropriate $Q({T}_{\mathrm{rot}})$ at the derived rotational temperature.

The constructed rotation diagrams for SO2, 34SO2, and 33SO2 are shown in Figure 8 and the derived temperatures and column densities are summarized in Table 2. Uncertainties in the table are of 2σ level and do not include systematic errors due to spectroscopic parameters extracted from the CDMS and JPL databases.

Figure 8.

Figure 8. Rotation diagram analysis of 32SO2 (upper left), 34SO2 (upper right), and 33SO2 (lower left) for ST11. The filled squares (black) are for ST11 and the open squares (green) are for the Orion hot core. The downward triangles represent upper limits. The straight-lines fitted to the ST11 and the Orion data points are shown by the black and green solid lines, respectively. The derived rotational temperatures are shown in Table 3.

Standard image High-resolution image

Table 2.  Estimated Column Densities

Molecule Trot N(X)
  (K) (cm−2)
SO2 190a 8.4 ± 0.3 × 1015a
34SO2 95a 6.2 ± 0.9 × 1014a
33SO2 64a 2.1 ± 0.7 × 1014a
H2 4.5 × 1023
CO 100 ∼8 × 1017
C17O 100 1.2 ± 0.09 × 1016
HCO+ 100 1.4 ± 0.01 × 1014
H13CO+ 100 6.9 ± 0.6 × 1012
H2CO 100 1.0 ± 0.1 × 1014
CH3OH 100 <3.5 × 1014
NO 100 9.1 ± 2.5 × 1015
HNCO 100 <4.3 × 1013
HC3N 100 <1.8 × 1013
SiO 100 1.5 ± 0.1 × 1013
C34S 100 <1.0 × 1013
H2CS 100 2.8 ± 0.9 × 1013
33SO 100 2.7 ± 0.6 × 1014
CH3OCH3 100 <1.3 × 1015
HCOOCH3 100 <7.1 × 1015
C2H5OH 100 <2.2 × 1015

Note. Uncertainties and upper limits are of 2σ level and do not include systematic errors due to adopted spectroscopic constants.

aDerived based on the rotation diagram analysis (see Section 4.2 for details).

Download table as:  ASCIITypeset image

For comparison purposes, we also perform the rotation diagram analysis for the Orion hot core using lines with similar spectroscopic properties to those used in the analysis of ST11. The Orion data are obtained with the 20'' beam size (∼0.04 pc at the Orion region) and are adopted from Schilke et al. (1997). The rotational temperatures of the Orion hot core derived here are somewhat lower than the temperatures derived in Schilke et al. (1997) using complete spectral line survey data in 325–360 GHz. This is because the high excitation lines used in our comparative analysis are fewer than those used in Schilke et al. (1997) to allow for a fair comparison. The above results are shown in Figure 8 and Table 3 together with the results of ST11. The table also shows rotational temperatures of SO2 and 34SO2 measured for W3 (H2O) and G34.3+0.15, which are obtained from the literature (Helmich & van Dishoeck 1997; MacDonald et al. 1996).

Table 3.  Comparison of Rotational Temperatures

Trot (K)
Molecule ST11(1) Orion-KL(2) W3 (H2O)(3) G34.3+0.15(4)
SO2 190 ± 5 124 (121 ± 4)a 184 108, 284b
34SO2 95 ± 7 138 (114 ± 6)a 179 131
33SO2 64 ± 14 104 (73 ± 5)a

Notes. Uncertainties are of 2σ level.

aNumbers in parentheses in the Orion data are derived in this work (see Section 4.2 for details). bTwo-temperature components.

References. Column 1: this work; Column 2: Schilke et al. (1997); Column 3: Helmich & van Dishoeck (1997); Column 4: MacDonald et al. (1996).

Download table as:  ASCIITypeset image

4.3. Column Densities of Other Molecules

The column densities of molecular species other than SO2 and its isotopologues are derived by solving Equation (2) for Ntotal under the assumption of the LTE and optically thin condition. We here assume that the observed molecular species are located in the same region as SO2 and its isotopologues, and thus have similar rotational temperatures; Trot is assumed to be 100 K, which is roughly the average rotational temperature of SO2, 34SO2, and 33SO2.

The similar spatial distributions of emission and the similar radial velocities of SO2 and other molecules support the validity of this assumption, but more rigorous excitation analysis with further spectral data is absolutely necessary in the future.

NO shows multiple transitions with the same upper state energy in the 350–351 GHz region. We estimate the column density for each transition, and the average value is adopted as the final NO column density. The scatter of column densities estimated from different lines is less than 20%.

We also estimate the upper limits of the column densities of important non-detection lines such as CH3OH, HNCO, C34S, and some complex molecules (see Section 5.2 for details of individual lines). Upper limits on peak brightness temperatures are measured using the observed spectra after appropriate channel binning, and then we assume the FWHM of 6 km s−1 to estimate the upper limits on integrated intensities. This FWHM is consistent with the typical velocity width of other molecular lines except for CO and HCO+.

The derived column densities and upper limits (2σ level) are summarized in Table 2.

4.4. Column Density of H2 and Total Gas Mass

The column density of molecular hydrogen, which usually dominates the total mass of embedded sources, is estimated using the dust continuum emission data obtained in our observations. The flux density of the dust continuum, Fν, at the optically thin frequency, ν, is expressed as

Equation (4)

where Ω is a beam solid angle, ${\tau }_{\nu }$ is an optical depth, ${B}_{\nu }({T}_{d})$ is the Planck function, and Td is a dust temperature (Whittet 1992). The optical depth is expressed as

Equation (5)

where ${\rho }_{d}$ is the mass density of dust, κν is the mass absorption coefficient, and L is path length. We use the mass absorption coefficient of dust grains coated by thin ice mantles as presented in Ossenkopf & Henning (1994). Using the dust-to-gas mass ratio, Z, the mass density of dust is expressed as

Equation (6)

where μ is mean atomic mass per hydrogen, ρH is the mass density of hydrogen, NH is the column density of hydrogen, and mH is the hydrogen mass. We here assume μ to be 1.41 according to Cox (2000). We assume that the dust-to-gas mass ratio in the LMC is lower than the typical Galactic value of 0.008 by a factor of three according to Bernard et al. (2008), and we use Z = 0.0027 in this work. By combining Equations (4)–(6), and assuming that all the hydrogen is in the form of H2 (${N}_{{{\rm{H}}}_{2}}={N}_{{\rm{H}}}/2$), the column density of molecular hydrogen is derived by the following equation:

Equation (7)

We measure the flux density per beam solid angle, ${F}_{\nu }/{\rm{\Omega }}$, using the 840 μm dust continuum image of ST11 shown in Figure 2. With the aperture size of 0farcs5 in diameter, the ${F}_{\nu }/{\rm{\Omega }}$ is measured to be 0.34 Jy beam−1. For dust temperature, we assume Td = 40 K, which is the typical dust temperature of high-mass YSOs in the LMC estimated based on far-infrared observations (van Loon et al. 2010). The assumed dust temperature is consistent with the far-infrared SED peak of ST11 (see Section 4.5 for description of the SED). Using Equation (7) and the above parameters, we estimate ${N}_{{{\rm{H}}}_{2}}$ = 4.5 × 1023 cm−2. This H2 column density corresponds to the total gas mass of 115 ${M}_{\odot }$ in the line of sight. The estimated dust opacity at 840 μm is ${\tau }_{840\mu {\rm{m}}}$ = 0.011, which corresponds to the visual extinction ($A{\rm{v}}$) of ∼150 mag.

We also estimate a lower limit on the gas density around ST11. We assume that gas is spherically distributed around a protostar with a radius which is the same as the beam size. With this assumption, the total gas mass around ST11 derived above corresponds to the H2 number density of 2 × 106 cm−3. We emphasize that this is a lower limit because the H2 density increases inversely with the assumed source size.

4.5. Luminosity and Stellar Mass

The spectral energy distribution (SED) of ST11 is shown in Figure 9. Details of the collected data are summarized in Table 4. Most of the energy is emitted in the mid- to far-infrared wavelength regions and the peak of the SED is between 60 μm and 70 μm, which is consistent with the characteristics of high-mass YSOs.

Figure 9.

Figure 9. The SED of ST11. The plotted data are based on IRSF/SIRIUS photometry (pluses, black), AKARI/IRC spectroscopy (solid line, blue), AKARI/IRC photometry (open squares, light blue), Spitzer/MIPS spectroscopy (solid line, green), Spitzer/IRAC and MIPS photometry (open diamonds, light green), Gemini/T-ReCS narrow-band photometry (open circles, light brown), Gemini/T-ReCS N-band spectroscopy (solid line, brown), Herschel/PACS and SPIRE photometry (filled diamonds, orange), the ALMA 840 μm continuum (filled star, red), and the best fitted SED model (dashed line, gray). See Section 4.5 for details.

Standard image High-resolution image

Table 4.  Photometric and Spectroscopic Data of ST11

Instrument Wavelength Flux FWHMa References
  (μm) (mJy) ('')  
IRSF/SIRIUS J 1.23 0.50 ± 0.07 1.3 (1)
IRSF/SIRIUS H 1.66 0.82 ± 0.06 1.2 (1)
IRSF/SIRIUS Ks 2.16 3.98 ± 0.22 1.1 (1)
AKARI /IRC NG 2.5–5 spectroscopy 7.3b (2)
AKARI /IRC N3 3.2 66.0 ± 2.6 4.0 (3)
Spitzer/IRAC Band 1 3.6 88.1 ± 6.6 1.7 (4)
Spitzer/IRAC Band 2 4.5 252.2 ± 10.3 1.7 (4)
Spitzer/IRAC Band 3 5.7 577.8 ± 23.6 1.9 (4)
Spitzer/IRAC Band 4 7.9 1090.0 ± 39.7 2.0 (4)
AKARI /IRC S11 11 885.6 ± 29.7 4.8 (3)
Gemini-S/T-ReCS Si1 7.7 1329 ± 29 0.6 (5)
Gemini-S/T-ReCS Si2 8.7 940 ± 53 0.5 (5)
Gemini-S/T-ReCS Si3 9.7 523 ± 10 0.6 (5)
Gemini-S/T-ReCS Si4 10.4 950 ± 10 0.8 (5)
Gemini-S/T-ReCS Si5 11.7 2297 ± 12 0.8 (5)
Gemini-S/T-ReCS Si6 12.3 3226 ± 28 0.5 (5)
Gemini-S/T-ReCS Lo-Res 8–12 spectroscopy 0.7c (5)
Spitzer/IRS SL, SH, LH 5–36 spectroscopy 4–11c (6)
Spitzer/MIPS 70 70 83810 ± 1313 18 (4)
Herschel/PACS 100 100 47180 ± 3090 8.6 (7)
Herschel/PACS 160 160 26800 ± 1629 13 (7)
Herschel/SPIRE 250 250 9650 ± 589 18 (7)
Herschel/SPIRE 350 350 4359 ± 226 27 (7)
Herschel/SPIRE 500 500 2110 ± 151 41 (7)
ALMA Band 7 837 23.2 ± 2.32 0.5 (5)

Notes.

aFWHM of the point-spread function. bExtraction width of the slitless spectroscopic data. cSlit width.

References. (1) Kato et al. (2007); (2) Shimonishi et al. (2010); (3) Kato et al. (2012); (4) Meixner et al. (2006); (5) This work; (6) Kemper et al. (2010); (7) Meixner et al. (2013).

Download table as:  ASCIITypeset image

The bolometric luminosity of ST11 is estimated to be 5 × 105 L, which is derived by integrating the interpolated SED from 1 μm to 1000 μm. About 60% of the total luminosity is emitted in the 30–100 μm wavelength region. The stellar mass of ST11 is estimated using the Online SED Fitter8 (Robitaille et al. 2007). As the input data of the SED fit, we use 2–840 μm photometric and spectroscopic data, which are obtained from the IRSF/SIRIUS, Spitzer SAGE, AKARI LSLMC, and Herschel HERITAGE databases (Meixner et al. 2006, 2013; Kato et al. 2007, 2012; Kemper et al. 2010; Shimonishi et al. 2010). We also use the results of our mid-infrared narrow-band filter photometry (centered at 7.73, 8.74, 9,69, 10.38, 11.66, and 12.33 μm) and spectroscopy (8–12 μm) conducted with T-ReCS at the Gemini South telescope. The 840 μm flux is estimated using the present ALMA data (see Section 4.4). We exclude the SPIRE 350 μm and 500 μm band data in the fit because they are possibly contaminated by diffuse emission around the YSO due to their large point-spread function (about 27'' and 41'' in FWHM, respectively). The distance to ST11 is assumed to be the same as the distance to the LMC.

The estimated stellar mass of ST11 in the best-fit model is 50 M. The visual extinction derived in the best-fit model is Av ∼80 mag, which differs by a factor of two from the value estimated based on the dust continuum data presented in Section 4.4. The difference may arise from the assumption of optical properties of dust and/or the existence of a dust temperature gradient in the line of sight, although this does not affect the main conclusion of this paper. The result of the SED fit is plotted in Figure 9.

5. DISCUSSION

5.1. Hot Molecular Core Associated with ST11

Hot cores are one of the early stages of high-mass star formation. It is suggested that the hot cores are in a transitional evolutionary stage between a deeply embedded protostellar object and a zero-age main-sequence star with a compact H ii region (e.g., Zinnecker & Yorke 2007). The physical properties of hot cores are characterized by a small source size (≤0.1 pc), a high density (≥106 cm−3), and warm gas/dust temperature (≥100 K) (e.g., Kurtz et al. 2000; van der Tak 2004). In this section, we discuss the presence of a hot core around ST11 on the basis of its physical and spectral properties.

5.1.1. Source Size

A source size of ST11, which is measured by detected emission lines, is consistent with those of Galactic hot cores. The molecular lines of NO, SiO, 33SO, SO2, 34SO2, and 33SO2 as well as the dust continuum show a compact source size, which is equal to or lower than the beam size of ∼0.1 pc. The lines of CO, C17O, HCO+, H13CO+, and H2CO show somewhat extended distributions, but the dominant emission come from the compact region associated with a high-mass YSO. The extended components of these emission lines may arise from their relatively low critical densities and/or efficient formation routes in gas-phase reactions. Protostellar outflows, which are discussed in Section 5.3, could also affect the spatial distribution of these molecular gas.

The spatial association of the molecular emission lines to the infrared source is another key to the nature of the detected warm and dense molecular gas. The separation between a hot core region and an infrared source is reported to be smaller than ∼0.03 pc for typical Galactic sources (De Buizer et al. 2003). With the present spatial resolution, the distribution of molecular lines and the dust continuum should coincide if they arise from a hot core region. This is clearly observed in the present data as shown in Figure 3.

5.1.2. Density

The H2 gas density around ST11 is estimated to be >2 × 106 cm−3 in Section 4.4. The value is a lower limit and actually clear detections of molecular lines with high critical densities such as those of SO2 (ncr ∼ 107 cm−3; Williams & Viti 2013) imply an even higher gas density around ST11. Galactic hot cores typically show gas densities between 106–108 cm−3 (estimated from Table 1 presented in Kurtz et al. 2000). The gas density around ST11 is thus consistent with those of known hot cores.

5.1.3. Gas Temperature

Rotational temperatures of molecular gas around ST11 are compared with those of Galactic hot cores in Table 3. The temperature of SO2, Trot = 195 ± 5 K, is in good agreement with the temperatures of the three Galactic hot cores in the table, which are between 108 K and 284 K. The temperature of 34SO2, Trot = 95 ± 7 K, is somewhat lower than that of SO2, but this may be attributed to the combined effect of the lack of available higher excitation lines of 34SO2 and the slightly larger optical thickness of the SO2 lines. Note that the temperature of 33SO2 is even lower, but this is probably because the detected lines are biased to those with low excitation temperatures (Eu < 200 K) due to the limited frequency coverage of the data. The rotation analysis of 34SO2 and 33SO2 for the Orion data using a similar set of lines to those used in the analysis of ST11 results in Trot = 114 ± 6 K for 34SO2 and Trot = 73 ± 5 K for 33SO2, which are lower than the values derived by Schilke et al. (1997) using a larger number of lines with a wide range of upper state energies (Table 3). Therefore we conclude that the actual temperature of SO2 gas is in the order of ∼100 K, which can trigger hot core chemistry via sublimation of ice mantles. Sublimation of ice mantles around ST11 is also suggested from the infrared observations of ices in Shimonishi et al. (2010), which reported that ST11 shows the second weakest ice absorption bands among 12 high-mass YSOs in the LMC.

5.1.4. Central Protostar

As estimated in Section 4.5, the bolometric luminosity (5 × 105 L) and stellar mass (50 M) of ST11 are consistent with the properties of high-mass YSOs. This indicates that ST11 can form a hot core region with the help of its intense radiation. In addition, the red SED of ST11 (Figure 9) indicates that it is still in an early evolutionary stage, which is consistent with the properties of hot core sources.

5.1.5. Spectral Properties

The velocity widths of emission lines from ST11 are typically 4–7 km s−1 except for CO and HCO+. This is in good agreement with the velocity widths of emission lines from Galactic hot core regions (typically 4–10 km s−1, e.g., Helmich & van Dishoeck 1997). The systemic velocities of the lines are in a very narrow range, typically 250–251 km s−1, which suggests that the detected molecular species spatially coexist in a small region around a high-mass YSO.

SO2 is often detected in hot core sources and is one of the useful tracers of warm and dense gas around a high-mass YSO (e.g., Beuther et al. 2009). However, sometimes warm SO2 gas is also prominently detected in deeply embedded and even younger protostellar objects in which hot cores are not yet formed (e.g., W3 IRS5 in Helmich et al. 1994). The reason for such SO2 enhancement is suggested to be due to shock dominated chemistry triggered by protostellar outflows. These deeply embedded sources often show a deep self-absorption profile in their CO or HCO+ emission lines due to the presence of a significant amount of cold gas in the envelope. As shown in Figure 4, emission from ST11 does not show such a deep self-absorption profile. This suggests that emission from ST11 is mostly dominated by warm gas and that ST11 has already reached the hot core phase.

On the other hand, the relatively weak strengths of the hydrogen recombination lines seen in the near-infrared spectrum of ST11 suggests that a prominent H ii region is not yet formed around the source (see Figure 1 in Shimonishi et al. 2010). Such a transitional evolutionary phase is consistent with the properties of hot cores.

In summary, the compact size of the emitting source, the presence of warm and dense molecular gas around a high-mass protostar, the sublimation of ice mantles, and the detections of rich molecular lines suggest that ST11 is associated with a hot molecular core. This is the first detection of an extragalactic hot molecular core. Note that CH3OH and complex organic molecules, which are often detected in Galactic hot cores, are not detected in ST11. The reason for the lack of these molecular species is discussed in the next section in the context of the characteristic hot core chemistry in low metallicity.

5.2. Molecular Abundances

Hot cores play a key role in the chemical complexity of interstellar and circumstellar molecules. In this section, we compare the chemical compositions of ST11 with those of Galactic hot cores and discuss the impact of metallicity on chemical processes in hot cores.

Fractional abundances of molecules around ST11 are shown in Table 5. We use the following isotope abundances for ST11: 12C/13C = 49, 16O/17O = 3400, and 32S/34S = 15, according to Wang et al. (2009) in which isotope ratios for the star-forming region N113 in the LMC are reported. For 33S, we use 32S/33S = 40, which is estimated in this work (see Section 5.2.12). For comparison purposes, the table also includes the abundances of three Galactic hot cores, Orion, W3 (H2O), and G34.3+0.15. Their abundances are re-estimated using the following isotope ratios of the local ISM and the Sun: 12C/13C = 77, 16O/17O = 1800, 32S/34S = 22, and 32S/33S = 127, according to Wilson & Rood (1994) and Anders & Grevesse (1989). All the abundances are estimated from the 345 GHz region data, except for NO in Orion, which is estimated using the transitions near 150 GHz.

Table 5.  Fractional Abundances

  N(X)/N(H2)  
Molecule ST11a Orionb W3 (H2O)c G34.3+0.15d Notes
CO 9.1 ± 0.7 × 10−5 6.9 × 10−5 1.9 × 10−4 >4.1 × 10−5 (1)
HCO+ 7.5 ± 0.7 × 10−10 1.2 × 10−9 3.1 × 10−9 >1.7 × 10−9 (2)
H2CO 2.2 ± 0.2 × 10−10 2.9 × 10−8 4.2 × 10−9 >3.2 × 10−10
CH3OH <8 × 10−10 1.8 × 10−7 9.2 × 10−8 3.4 × 10−8
NO 2.0 ± 0.6 × 10−8 1.1 × 10−8 >1.1 × 10−8
HNCO <1 × 10−10 1.6 × 10−9 5.0 × 10−9 >4.3 × 10−10
HC3N <4 × 10−11 3.1 × 10−9 1.6 × 10−10 >3.6 × 10−11
SiO 3.3 ± 0.2 × 10−11 4.9 × 10−9 6.1 × 10−10 >2.3 × 10−11
CS <3 × 10−10 8.3 × 10−9 9.9 × 10−9 >2.7 × 10−9 (3)
H2CS 6.2 ± 2.0 × 10−11 9.3 × 10−10 1.6 × 10−9 2.3 × 10−9
SO 2.4 ± 0.8 × 10−8 2.0 × 10−7 1.3 × 10−8 >5.0 × 10−9 (4)
SO2 2.1 ± 0.3 × 10−8 1.2 × 10−7 5.5 × 10−8 3.6 × 10−8 (5)
CH3OCH3 <3 × 10−9 9.1 × 10−9 2.1 × 10−8
HCOOCH3 <2 × 10−8 1.6 × 10−8 7.0 × 10−9 3.0 × 10−8
C2H5OH <5 × 10−9 8.6 × 10−10 6.6 × 10−9

Notes. Uncertainties and upper limits are of 2σ level and do not include systematic errors due to adopted spectroscopic constants. See Section 5.2 for details of the data. (1) Estimated from C17O. (2) Estimated from H13CO+ except for Orion, which is estimated from HC18O+. (3) Estimated from C34S. (4) Estimated from 34SO, except for ST11, which is estimated from 33SO. (5) Estimated from 34SO2.

References.aThis work; bZiurys et al. (1991); Sutton et al. (1995); Schilke et al. (1997), also see Section 5.2; cHelmich & van Dishoeck (1997); dMacDonald et al. (1996).

Download table as:  ASCIITypeset image

Since ST11 is observed with the beam size of 0.12 pc, we here compare the abundances which are derived from the data with similar beam sizes. For Orion, we calculate the average abundance of five regions around the hot core (hot core, compact ridge, extended ridge, northwest plateau, and southeast plateau) using the data presented in Sutton et al. (1995). The abundances of NO and HNCO in Orion are only for the central part of the hot core and the data are adopted from Ziurys et al. (1991) and Schilke et al. (1997), respectively. The averaged Orion abundances roughly reproduce the abundances observed for the ∼0.08 pc region in diameter. The beam sizes of W3 (H2O) and G34.3+0.15 correspond to 0.13 pc and 0.21 pc, respectively. We here assume the distances to Orion, W3 (H2O), and G34.3+0.15 to be 0.41 kpc, 1.95 kpc, and 3.1 kpc (MacDonald et al. 1996; Reid et al. 2009).

Figure 10 compares the fractional abundances of ST11 and Galactic hot cores. In general, most of the molecular species in ST11 show lower abundances compared to Galactic hot cores. Several molecules such as H2CO, CH3OH, HNCO, and CS show significantly lower abundances, which cannot be simply explained by the low abundances of heavy elements in the LMC. On the other hand, NO shows the higher abundance in ST11 than in Galactic sources despite the notably low nitrogen abundance in the LMC. Characteristics of individual molecules are discussed below.

Figure 10.

Figure 10. Molecular abundances of CO, HCO+, H2CO, CH3OH, NO, HNCO, HC3N, SiO, CS, H2CS, SO, and SO2 for ST11 in the LMC (red) and Galactic hot cores: Orion hot core (cyan), W3 H2O (orange), and G34.3+0.15 (gray). The downward and upward triangles represent the upper and lower limits, respectively. The plotted data are summarized in Table 5.

Standard image High-resolution image

5.2.1. CO

Carbon monoxide is detected by the CO(3–2) and C17O(3–2) lines in ST11. The line intensity ratio of CO(3–2)/C17O(3–2) is about 28, and using 16O/17O = 3400 in the LMC (Wang et al. 2009) we estimate optical depths of ∼120 and ∼0.04 for CO(3–2) and C17O(3–2), respectively. The CO(3–2) line is completely optically thick, while the C17O(3–2) line is optically thin as in most of Galactic hot cores (e.g., Helmich & van Dishoeck 1997). The abundance of CO in ST11 estimated from C17O is 9.1 × 10−5, which is very similar to the average CO abundance of 1.0 × 10−4 for three Galactic hot cores in Table 5. Elemental carbon and oxygen are about three times less abundant in the LMC; [C/H]LMC = 1.2 × 10−4 and [O/H]LMC = 2.3 × 10−4 (Korn et al. 2002), while [C/H] = 3.3 × 10−4 and [O/H] = 6.8 × 10−4 (Grevesse & Sauval 1998). Given the low abundances of carbon and oxygen in the LMC, the CO around ST11 is slightly overproduced as compared with Galactic hot cores.

5.2.2. HCO+

Formyl ion is detected via the HCO+(4–3) and H13CO+(4–3) lines. The line intensity ratio of HCO+(4–3)/H13CO+(4–3) is about 19, and using 12C/13C = 49 in the LMC (Wang et al. 2009) we estimate optical depths of ∼2 and ∼0.05 for HCO+(4–3) and H13CO+(4–3), respectively. The HCO+(4–3) line is moderately optically thick, while the H13CO+(4–3) line is optically thin. The abundance of HCO+ in ST11 estimated from H13CO+ is 7.5 × 10−10, which is lower by a factor of ∼3 compared to the average HCO+ abundance of 2.0 × 10−9 for three Galactic hot cores in Table 5.

We here discuss the connection between the observed HCO+ abundance and the cosmic-ray ionization rate in the LMC. One of the possible pathways to form HCO+ in molecular clouds is the gas-phase reaction:

Equation (8)

in which ${{\rm{H}}}_{3}^{+}$ is formed by ionization of H2 via cosmic rays and subsequent reaction with H2 (e.g., Caselli et al. 1998).

We here simply assume that HCO+ is a major cation and the electron is a major anion in dense clouds. Given that the production of cations and electrons by the cosmic-ray ionization balances with their recombination, the equilibrium is formulated as

Equation (9)

where ζ is a cosmic-ray ionization rate, krec = 2.4 × 10−7 (300/T)0.69 cm3 s−1 (5.1 × 10−7 at T = 100 K) is a dissociative recombination rate of HCO+ and electrons (Mitchell 1990), and n indicates the number density of each species. With the assumption that cations and electrons have the the same abundance in neutral clouds, this equation is reformatted as

Equation (10)

where X indicates a fractional abundance relative to H2. This equation indicates that the abundance of HCO+ is proportional to ${\zeta }^{0.5}$ under the above simple assumption.

The cosmic-ray density in the LMC is reported to be 20%–30% of the typical Galactic value of ζ = 3 × 10−17 s−1 based on gamma-ray observations (Abdo et al. 2010). The low cosmic-ray density leads to the low cosmic-ray ionization rate (e.g., Spitzer & Tomasko 1968). According to Equation (10), the HCO+ abundance in the LMC is expected to be lower than the Galactic value by a factor of ∼2 for the same H2 density and gas temperature. This factor is close to the observed difference of the HCO+ abundance between ST11 and Galactic sources. Therefore, we speculate that the low cosmic-ray ionization rate in the LMC is one of the factors that contribute to the slightly low abundance of HCO+ in ST11.

Note that the above assumption is a simplified picture, i.e., the cosmic-ray ionization rate may vary within the LMC and the electron abundance may differ from that of HCO+ due to the presence of other ions. Furthermore, in the above discussion, we consider only one of the possible mechanisms to form HCO+. In the circumstellar environment of hot cores, however, outflows and stellar radiation would affect the ionization condition. Previous studies actually argue that molecular outflows and shocks contribute to the enhancement of HCO+ in star-forming cores (e.g., Girart et al. 1999; Rawlings et al. 2004; Arce & Sargent 2006). Future observations of other molecular ions and shock tracers are required to obtain a comprehensive understanding of ionization conditions in the LMC.

5.2.3. H2CO

Formaldehyde (H2CO) is detected by the ${5}_{\mathrm{1,5}}$–41,4 transition at 351.7686 GHz (Eu = 62 K). The abundance of H2CO in ST11 is estimated to be 2.2 × 10−10, while Galactic hot cores show the abundances between 10−9 and 10−8 as in Table 5, except for G34.3+0.15 which only has a lower limit. The H2CO is less abundant in ST11 than Galactic sources by 1–2 orders of magnitude.

Both gas-phase and grain surface reaction pathways are suggested for the formation of H2CO in dense ISM. The grain surface pathway for the H2CO formation requires the hydrogenation of CO (e.g., Hama & Watanabe 2013, and references therein). However, this pathways is suppressed in the LMC due to higher dust temperatures (see the detailed discussion presented in Section 5.2.4). Therefore we speculate that the suppressed formation of the H2CO ice on grain surfaces contributes to the low abundance of H2CO gas in ST11.

5.2.4. CH3OH

No methanol (CH3OH) lines are detected in ST11. We estimate an upper limit on the fractional abundance to be <8 × 10−10 using the spectrum at the frequency of the CH3OH(70–60 A+) transition at 338.4087 GHz (Eu = 65 K). Abundances of CH3OH in Galactic hot cores are typically between ∼10−8 and ∼10−7 as in Table 5, but sometimes the abundance reaches ∼10−6 (e.g., G5.89-0.39 in Thompson & MacDonald 1999). In ST11, CH3OH is depleted by at least 2–3 orders of magnitude as compared with Galactic hot cores.

The low abundance of CH3OH in the LMC has been suggested in previous studies. Nishimura et al. (2016) reported, based on spectral line surveys toward molecular clouds in the LMC, that thermal emission lines of CH3OH gas are significantly weak in the LMC compared to our Galaxy. Searches for maser emission in the LMC reported the underabundance of CH3OH masers in the LMC (e.g., Ellingsen et al. 2010). Furthermore, Shimonishi et al. (2016) reported that the CH3OH ice is less abundant in the LMC than in our Galaxy based on infrared observations. They reported that all of the observed ten high-mass YSOs in the LMC show CH3OH ice abundances less than 5%–8% relative to water ice, while about one-third of Galactic high-mass YSOs show CH3OH ice abundances between 10% and 40%. The authors suggested that warm ice chemistry is responsible for the low abundance of solid CH3OH in the LMC, i.e., high dust temperatures in the LMC suppress the hydrogenation of CO on the grain surface, which leads to inefficient production of the CH3OH ice. A decrease in the efficiency of CO hydrogenation at an elevated temperature is measured by laboratory experiments (e.g., Watanabe et al. 2003). In addition, the reduced formation of CH3OH in relatively warm molecular clouds is confirmed by numerical simulations of grain surface chemistry dedicated to the LMC environment (Acharyya & Herbst 2015).

The warm gas temperature around ST11 suggests that ice mantles are mostly sublimated. The low elemental abundances of carbon and oxygen in the LMC should partly contribute to the observed low CH3OH abundance. However, we need an additional explanation to account for the CH3OH deficiency by several orders of magnitude, and this should be related to the grain surface chemistry by which CH3OH is mainly formed. We suggest that suppressed production of the CH3OH ice due to warm ice chemistry at the molecular cloud stage or the deeply embedded YSO stage is responsible for the deficiency of CH3OH gas in this LMC hot core.

The above interpretation on the low CH3OH abundance can also be applied to the low abundance of H2CO in ST11 since a major grain surface pathway of both species requires the CO hydrogenation. However, the difference is that H2CO has possible formation routes both in the gas-phase and in the solid-phase, while CH3OH does not have an efficient formation route in the gas-phase under typical molecular cloud conditions. We speculate that the H2CO observed in ST11 is mainly produced by gas-phase reactions. A slightly extended distribution of the H2CO emission (Figure 3) supports this idea, because the species sublimated from ice mantles often show a compact distribution compared to those produced in the gas-phase. Owing to the lack of a possible gas-phase formation route, the degree of depletion around ST11 is actually larger in CH3OH as compared to H2CO.

5.2.5. NO

Nitric oxide (NO) is one of the interesting molecules that show a peculiar abundance in ST11. The abundance of NO in ST11 is estimated to be 2.0 × 10−8. On the other hand, the average abundance and standard deviation of NO in six Galactic sources in Ziurys et al. (1991) including hot cores, high-mass protostellar objects, and Galactic center objects is 8.2 ± 2.9 × 10−9. In the LMC, nitrogen is even more depleted than other major elements such as carbon and oxygen; the elemental abundance of nitrogen in the LMC is [N/H]LMC = 1.0 × 10−5 (Korn et al. 2002), while [N/H] = 8.3 × 10−5 for the Sun (Grevesse & Sauval 1998). Despite the lower nitrogen abundance by a factor of 8 in the LMC, the abundance of NO in ST11 is higher than in Galactic star-forming regions by a factor of 2–3.

Ammonia (NH3) gas, which is often the most abundant nitrogen-bearing species in star-forming regions, is reported to be deficient in the LMC by 1.5–2 orders of magnitude as compared with Galactic star-forming regions (Ott et al. 2010). A previous infrared study suggests that the NH3 ice around an embedded high-mass YSO in the LMC is possibly less abundant as compared with Galactic high-mass YSOs (Shimonishi et al. 2016). Although the NH3 abundance in ST11 is unknown, the relatively high abundance of NO implies that NO could play an important role in nitrogen chemistry around ST11.

It is suggested that NO is formed by a neutral–neutral reaction in the gas-phase (e.g., Herbst & Klemperer 1973; Pineau des Forets et al. 1990):

Equation (11)

and destroyed by

Equation (12)

A numerical simulation of nitrogen chemistry in warm molecular gas suggests that the resultant abundance of NO through the above reactions increases as the temperature of molecular gas increases (Pineau des Forets et al. 1990). However, the enhancement of the NO abundance is suggested to be only a factor of 1.5 as the gas temperature increases from 70 to 250 K. Thus we speculate that pure gas-phase chemistry makes only a limited contribution to the enhanced abundance of NO in ST11.

The compact distribution of the NO emission observed in ST11 may hint at a possible origin in ice sublimation. Laboratory experiments argue that the bombardment of energetic ions on interstellar ice analogs containing H2O, O2, and N2, or CO and N2, produce NO in ice mixtures (Boduch et al. 2012; Sicilia et al. 2012). However, the low cosmic-ray density in the LMC (see Section 5.2.2) implies that the formation of NO by such energetic processing is less efficient around ST11. A reaction pathway though diffusive grain surface chemistry, N + O $\to $ NO, is also suggested in the literature (Tielens & Hagen 1982), but uncertainty remain regarding the efficiency of the reaction in actual grain surfaces.

The reason for the enhanced abundance of NO in ST11 despite the low elemental abundance of nitrogen in the LMC remains to be investigated. Detailed modeling of hot core chemistry in metal-poor environments is obviously necessary to interpret the enhanced abundance of NO in ST11.

5.2.6. HNCO

Isocyanic acid (HNCO) is not detected in ST11. We estimate an upper limit on the fractional abundance of HNCO to be <1 × 10−10 based on the non-detection of the 160,16–150,15 transition at 351.6333 GHz (Eu = 143 K). There seems to be an emission-like feature at the position of the HNCO(163,14–153,13 and 163,13–153,12, 351.4168 GHz) line with VLSR = 250.8 km s−1. However, the high upper state energy of this transition (Eu = 518 K) and the non-detection of the 160,16–150,15 transition suggest that this line is likely a spurious signal. Galactic hot cores and high-mass protostellar objects typically show HNCO abundances between 10−9 and 10−8 (Bisschop et al. 2007), which is consistent with the average abundance of 2.3 × 10−9 for three Galactic hot cores in Table 5. Thus the HNCO abundance in ST11 is lower than Galactic counterparts by at least 1–2 orders of magnitude.

Both gas-phase and grain surface reaction pathways are suggested for the formation of HNCO in dense ISM (e.g., Tielens & Hagen 1982; Turner et al. 1999). The grain surface pathway requires the hydrogenation of CO to form HCO and subsequent reaction of HCO + N $\to $ HNCO. However, as discussed in Section 5.2.4, the hydrogenation of CO is suggested to be less efficient in the LMC due to high dust temperatures, which possibly suppress the grain surface formation of HNCO. This may contribute to the low abundance of HNCO in ST11 along with the low elemental abundance of nitrogen in the LMC.

5.2.7. HC3N

Cyanoacetylene (HC3N), an unsaturated carbon-chain molecule often detected in star-forming regions, is not detected in ST11. We estimate an upper limit on the fractional abundance of HC3N to be <4 × 10−11 based on the non-detection of the 38–37 transition at 345.6090 GHz (Eu = 323 K). Abundances of HC3N in Galactic hot cores vary from ∼10−11 to ∼10−9 as in Table 5. The estimated upper limit on the HC3N abundance in ST11 is thus at the lower end of the Galactic HC3N abundance.

5.2.8. SiO

Silicon monoxide (SiO) is detected by the J = 8–7 transition at 347.3306 GHz (Eu = 75 K). The abundance of SiO for ST11 is 3.3 × 10−11, which seems to be low compared to Galactic sources, but the dispersion in the Galactic SiO abundances is significantly large as shown in Figure 10. The average SiO abundance for the three Galactic hot cores in Table 5 is 1.8 × 10−9, but their abundances range from ∼10−11 to ∼10−9.

The elemental abundance of silicon in the LMC is about three times lower than the solar abundance; [Si/H]LMC = 1.3 × 10−5 (Korn et al. 2002), while [Si/H] = 3.6 × 10−5 (Grevesse & Sauval 1998). Thus the elemental abundance in the LMC partly contributes to the observed low SiO abundance, but additional processes may be necessary to account for the 1–2 orders of magnitude lower SiO abundance.

SiO is often linked with the presence of energetic protostellar outflows since destruction of silicon-bearing dust (e.g., MgSiO3) by shock is believed to be the origin of SiO (e.g., Blake et al. 1996; Wright et al. 1996; Tercero et al. 2011, and references therein). The destruction of dust grains and release of silicon can occur around ST11 since it shows a sign of outflows as described in Section 5.3. On the other hand, the low abundance of SiO would imply that shock chemistry is not very dominant in the present chemical compositions around ST11. Further quantitative discussions require detailed modeling of SiO production in low metallicity environments. Various factors including the degree of shock, elemental abundances, and compositions of dust grains may affect the SiO abundance.

5.2.9. CS

The C34S(7–6) transition at 337.3965 GHz (Eu = 65 K) is not detected toward ST11. We estimate an upper limit on the abundance of carbon monosulfide (CS) to be <3 × 10−10, while Galactic hot cores and high-mass protostellar objects typically show the CS abundances between 10−9 and 10−8 (Table 5; see also van der Tak et al. 2003). CS is significantly less abundant by at least 1–2 orders of magnitude in ST11 compared to Galactic sources, and the low elemental abundance of carbon and sulfur in the LMC cannot by itself explain the depletion of CS. Given the relatively high abundance of SO2 in ST11 as discussed in Section 5.2.12, a significant amount of gas-phase sulfur is possibly incorporated in SO2 due to different circumstellar chemistry around ST11.

5.2.10. H2CS

We detect the emission line of thioformaldehyde (H2CS) at 338.0832 GHz (101,10–91,9, Eu = 65 K), but the S/N of the line is relatively poor. The H2CS abundance is estimated to be 6.2 × 10−11 in ST11, while Galactic hot cores show the abundance of ∼2 × 10−9 on average with a relatively low dispersion (Table 5). H2CS seems to be less abundant in the LMC hot core than in Galactic hot cores as well as CS, and only a small fraction of gas-phase sulfur is incorporated into H2CS in ST11. Since our H2CS abundance is estimated using a single line with poor S/N, further multiline observations of H2CS transitions are necessary for accurate determination of the H2CS abundance in the LMC.

5.2.11. SO

The abundance of sulfur monoxide (SO) is estimated to be 2.4 × 10−8 in ST11, while Galactic hot cores and high-mass protostellar objects typically show SO abundances between 10−9 and 10−8 (Table 5; see also van der Tak et al. 2003). The Orion hot core shows an exceptionally high SO abundance of 2.0 × 10−7. The results suggest that ST11 shows a slightly higher abundance of SO than typical Galactic counterparts despite the low metallicity of the LMC. If we use the solar 32S/33S ratio of 127, then the SO abundance in ST11 is 7.6 × 10−8, which is even higher than typical Galactic abundances.

We emphasize, however, that our SO abundance entails considerable uncertainty. The abundance is estimated using a single 33SO line at 337.1986 GHz (Eu = 81 K), which contains a number of unresolved hyperfine structures, while Galactic SO abundances are estimated from 34SO lines. Multiline observations of SO and its isotopologues are necessary for further discussion.

5.2.12. SO2

Sulfur dioxide (SO2) is in this study a key molecule for which we detect the largest number of transitions; we detect nine SO2 lines, one SO2 (ν2 = 1) line, ten 34SO2 lines, and five 33SO2 lines. The fractional abundance of SO2 is estimated to be 2.1 × 10−8 in ST11 based on rotation diagram analysis of 34SO2 lines. The column density ratio of SO2 and 34SO2 is about 14, suggesting that the optically thin assumption is mostly valid for the observed SO2 lines because the 32S/34S ratio is reported to be 15 in the LMC. If we assume the solar isotope ratio of 32S/34S = 22, the SO2 lines could be moderately optically thick. In either case, we can reasonably assume that 34SO2 and 33SO2 lines are optically thin because these isotopologues are much less abundant than 32SO2. The isotope abundance of 33S in ST11 based on the present results is 32S/33S = 40, while the solar value is 32S/33S = 127.

Galactic hot cores typically show the SO2 abundance from ∼10−8 to ∼10−7, while even younger embedded high-mass protostellar objects show the abundance of ∼10−9 (van der Tak et al. 2003). ST11 shows a factor of ∼3 lower SO2 abundance as compared with the average abundance of 7.0 × 10−8 for the three Galactic hot cores in Table 5. The elemental abundance of sulfur in the LMC is reported to be [S/H]LMC = 5.0 × 10−6, while the solar abundance is [S/H] = 1.9 × 10−5 (Russell & Dopita 1992). Sulfur is less abundant by a factor of ∼4 in the LMC compared with the solar abundance. Thus, the low SO2 abundance in ST11 is well explained by the elemental abundance of sulfur in the LMC. This would suggest that hot core chemistry of SO2 is dependent on elemental abundances of the host galaxy. A multitude of SO2 and its isotopologue line detections in ST11 imply that SO2 can be a key molecular tracer to test hot core chemistry in metal-poor environments.

It should be noted that sulfur chemistry in hot core regions is highly time-dependent (e.g., Charnley 1997). Hence, both the age and the interstellar environment should be taken into account to interpret the chemical compositions of sulfur-bearing species around ST11. Numerical simulations of hot core chemistry dedicated to low metallicity environments are thus highly required.

5.2.13. CH3OCH3, HCOOCH3, C2H5OH

Spectral lines from complex organic molecules are not detected in this work. We estimate upper limits on fractional abundances of three molecules whose relatively strong transitions are covered in the present data: [CH3OCH3/H2] < 3 × 10−9, [HCOOCH3/H2] < 2 × 10−8, and [C2H5OH/H2] < 5 × 10−9. The average abundances of these molecules for Galactic hot cores in Table 5 are [CH3OCH3/H2] = 1.5 × 10−8, [HCOOCH3/H2] = 1.8 × 10−8, and [C2H5OH/H2] = 3.7 × 10−9. These estimates suggest that CH3OCH3 is less abundant by at least an order of magnitude in ST11 compared with Galactic hot cores, while the upper limit abundances of HCOOCH3 and C2H5OH are comparable to the average abundances of Galactic sources. Although the present data do not provide conclusive upper limits on the abundances of the majority of complex organic molecules in the LMC, the possibly lower abundance of CH3OCH3 in ST11 suggests that the formation of large molecules may be less efficient in the LMC.

It is suggested by theoretical studies that CH3OH in ice mantles plays an important role in the formation of complex organic molecules (e.g., Nomura & Millar 2004; Garrod 2008; Herbst & van Dishoeck 2009). Shimonishi et al. (2016) argue that the formation of complex organic molecules from methanol-derived species could be less efficient in the LMC due to the low abundance of the CH3OH ice around high-mass YSOs. We confirm the significant deficiency of CH3OH gas around ST11 in this work. We thus speculate that the low abundance of CH3OH contributes to the low production efficiency of CH3OCH3 and possibly other complex organic molecules in the LMC.

Note that unknown rotational temperatures and limited frequency coverages produce considerable uncertainties on the abundance estimate. Furthermore, the formation of large molecules from other parent species such as H2CO, C2H, c-C3H2, and NO, which are detected in the LMC, should be taken into account for a comprehensive understanding of the complex chemistry around protostars. Hence, further observations with broader frequency coverage and higher sensitivity are critically needed in conjunction with theoretical and experimental efforts to understand complex organic chemistry in metal-poor environments.

5.3. Molecular Outflows

In this section, we discuss the second detection of extragalactic protostellar outflows after Fukui et al. (2015), who reported the detection of protostellar outflows in the star-forming region N159 in the LMC with ALMA. Evidence of protostellar outflows is seen in the CO emission line of ST11. Figure 11(a) shows the spectral profile of the CO(3–2) line extracted from the 0farcs5 diameter region centered at ST11. The profile of the optically thin C17O(3–2) line is also shown for comparison purposes. The CO(3–2) profile shows an apparently broader velocity width compared to other lines detected in ST11. The central velocity of the CO line is nearly consistent with those of other lines (VLSR ∼ 250), but the emission distributes in a wide velocity range from 230 to 280 km s−1. In addition, a prominent redshifted component is seen around VLSR ∼ 270 km s−1. We suggest that these high-velocity wing components are due to protostellar outflows from ST11 because the high-velocity CO gas is well spatially associated with a high-mass YSO as described in the next paragraph. The observed outflow velocity of ∼10–30 km s−1 is consistent with the molecular outflow velocities observed in Galactic high-mass star-forming regions (e.g., Lada 1985). A high-velocity component is not obviously seen in other lines besides CO(3–2).

Figure 11.

Figure 11. (a) Comparison of spectral line profiles of CO(3–2) (thick solid line, black) and C17O(3–2) (thin solid line, gray) observed toward ST11. The spectra are arbitrarily scaled and the horizontal axis is the LSR velocity. The CO line shows a broad line width and a prominent redshifted component, which indicate the presence of outflows in the line of sight. The blue and red vertical dashed lines represent the velocity range of the blueshifted and redshifted high-velocity components, which are visualized in the right panel. (b) Spatial distributions of integrated intensities of high-velocity wings. The blueshifted and redshifted wing components are shown by the blue and red contours. The contour levels are 20%, 40%, 60%, and 80% of the peak intensity. The background is the 840 μm continuum. The synthesized beam size is shown by the gray filled circle at the lower left.

Standard image High-resolution image

Spatial distributions of high-velocity wings are shown in Figure 11(b). Here we define the velocity range of the blueshifted wing to be 230.5 to 240.5 km s−1 and the redshifted wing to be 265.5 to 285.5 km s−1. Both high-velocity wings are spatially associated with a central protostar which is traced by the the continuum emission. The complex structures seen in the distribution of high-velocity gas imply that the circumstellar environment of ST11 is dynamically active. If we assume the spatial extent of high-velocity gas to be 0.24 pc (∼1'') and the outflow velocity to be 20 km s−1, according to the distribution and spectrum of the red wing component, we can roughly estimate an upper limit on the dynamical timescale of the outflow, which is about 104 years. This timescale is lower than a typical formation time of high-mass stars (∼105 years, Zinnecker & Yorke 2007) and thus consistent with high-mass star formation scenarios.

We also make a rough estimate of several outflow parameters, namely, the outflow mass, the mass entrainment rate, the mechanical force, and the energy. The outflow mass is estimated by adding blue and red wing components within the region where outflow gas emission is detected with an S/N higher than six. To convert the integrated intensities of CO(3–2) to the total gas mass, we use a conversion factor of 8.8 M (K km s−1)−1 pc−2, which is derived from the present observations toward the ST11 center. The derived outflow mass is ${M}_{\mathrm{out}}$ = 74 M, where the blueshifted component contains 13 M (Mblue) and the redshifted component contains 61 M (Mred). This outflow mass corresponds to the mass entrainment rate of $\dot{M}={M}_{\mathrm{out}}$/t = 7 × 10−3 M yr−1, where t is the dynamical timescale of outflows discussed above (104 years). The mechanical force (F) and the energy (E) of outflows are derived by the following equations: $F=({M}_{\mathrm{blue}}{V}_{\mathrm{blue}}+{M}_{\mathrm{red}}{V}_{\mathrm{red}})/t$ and $E\,=({M}_{\mathrm{blue}}{V}_{\mathrm{blue}}^{2}+{M}_{\mathrm{red}}{V}_{\mathrm{red}}^{2})/2$. Vblue and Vred are mean velocities of blue and red wing components, and we use Vblue = 13 km s−1 and Vred = 20 km s−1, respectively. Consequently, the derived outflow force is F = 0.14 M km s−1 yr−1 and the outflow energy is E = 3 × 1047 erg. The above outflow parameters for ST11 are roughly consistent with those observed in Galactic high-mass YSOs that have similar luminosities with ST11 (e.g., Beuther et al. 2002).

The present detection increases the number of extragalactic protostellar outflow samples, which should help understand the dynamical processes of high-mass star formation in different metallicity environments. Systematic observations of extragalactic outflow sources are necessary for statistical comparison of Galactic and extragalactic protostellar outflows. Further detailed analysis of outflows around ST11 is beyond the scope of this paper and will be presented in a future work.

6. SUMMARY

We report the first detection of an extragalactic hot molecular core based on radio interferometric observations toward ST11, a high-mass YSO in the LMC, with ALMA. The high spatial resolution (0.12 pc) ALMA Band 7 (345 GHz) spectral and continuum band observation data are presented. We discuss the physical and chemical properties of the source and obtain the following conclusions.

  • 1.  
    Molecular emission lines of CO, C17O, HCO+, H13CO+, H2CO, NO, SiO, H2CS, 33SO, 32SO2, 34SO2, and 33SO2 are detected from a compact region (∼0.1 pc) associated with a high-mass YSO. Furthermore, a number of high excitation lines (Eu > 100 K) of SO2 and its isotopologues are detected. On the other hand, CH3OH, HNCO, CS, HC3N, and complex organic molecules are not detected.
  • 2.  
    The physical properties of ST11 are derived using the obtained data. The H2 gas density around the source is estimated to be at least 2 × 106 cm−3 based on the dust continuum data. The temperature of molecular gas is estimated to be higher than 100 K based on rotation diagram analysis of SO2 and 34SO2 lines. The SED analysis in the 1–1000 μm range suggests that ST11 is a high-mass YSO with the luminosity of 5 × 105 L and the stellar mass of 50 M.
  • 3.  
    The compact size of the emitting source, the warm gas temperature, high density, and rich molecular lines around the high-mass protostar suggest that ST11 is associated with a hot molecular core.
  • 4.  
    We find that the molecular abundances of the hot core in the LMC are significantly different from those of Galactic hot cores. The abundances of CH3OH, H2CO, and HNCO are remarkably lower compared with Galactic sources by at least 1–3 orders of magnitude, although the gas temperature is warm enough for the sublimation of ice mantles. The deficiency of CH3OH gas in a warm and dense region is consistent with the previously reported low abundance of the CH3OH ice in the LMC. We suggest that the chemical compositions of ST11 are characterized by the deficiency of molecules whose formation requires the hydrogenation of CO on grain surfaces.
  • 5.  
    It is interesting that NO shows a higher abundance in ST11 than in Galactic sources despite the notably low abundance of nitrogen in the LMC. This is in contrast to the low abundances of nitrogen-bearing molecules such as NH3, HCN, and HNC in the LMC reported by previous studies. The reason for the enhanced abundance of NO remains to be investigated.
  • 6.  
    The slightly lower abundance of SO2 in ST11 than in Galactic hot cores is well explained by the low abundance of elemental sulfur in the LMC. CS and H2CS are less abundant than Galactic hot cores by at least 1–2 orders of magnitude. The abundance of SO is possibly high in the LMC, but the estimate based on a single 33SO line should be taken with caution. The large number of SO2 and its isotopologue line detections in the LMC hot core imply that SO2 can be a key molecular species to test hot core chemistry in metal-poor environments.
  • 7.  
    We find molecular outflows around ST11, which is the second detection of an extragalactic protostellar outflow. An apparently broad velocity width is seen in the spectral profile of the CO(3–2) line, which ranges in the LSR velocity from 230 to 280 km s−1. A prominent high-velocity component is also seen in the redshifted wing. We estimate an upper limit on the dynamical timescale of outflows to be 104 years, which is consistent with the timescale of high-mass star formation. Several outflow parameters are also estimated based on the present results.

This paper makes use of the following ALMA data: ADS/JAO.ALMA#2012.1.01108.S. ALMA is a partnership of ESO (representing its member states), NSF (USA), and NINS (Japan), together with NRC (Canada), NSC, and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO, and NAOJ. This work has made extensive use of the Cologne Database for Molecular Spectroscopy and the molecular database of the Jet Propulsion Laboratory. Partly based on data obtained at the Gemini Observatory via the time exchange program between Gemini and the Subaru Telescope (Program ID: S10B-120). The Gemini Observatory is operated by the Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini partnership: the National Science Foundation (United States), the National Research Council (Canada), CONICYT (Chile), Ministerio de Ciencia, Tecnología e Innovación Productiva (Argentina), and Ministério da Ciência, Tecnologia e Inovação (Brazil). The authors are grateful to Satoshi Yamamoto for his useful comment on spectral data. T.S. was supported by the ALMA Japan Research Grant of NAOJ Chile Observatory, NAOJ-ALMA-0061. This work is supported by a Grant-in-Aid from the Japan Society for the Promotion of Science (15K17612). Finally, we would like to thank an anonymous referee, whose suggestions greatly improved this paper.

Footnotes

Please wait… references are loading.
10.3847/0004-637X/827/1/72