Free Access
Issue
A&A
Volume 600, April 2017
Article Number A16
Number of page(s) 10
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201630183
Published online 21 March 2017

© ESO, 2017

1. Introduction

The Large Magellanic Cloud (LMC) at a distance of ~50 kpc (Pietrzyński et al. 2013) is the nearest low-metallicity (Rolleston 2002) star-forming galaxy to the Milky Way. The relatively face-on view offered by the LMC provides an ideal perspective for studying star formation, particularly massive star formation associated with its numerous stellar clusters. In the LMC, the far-ultraviolet (FUV) field is stronger than in the Milky Way (Westerlund 1990). It is thus well suited for studying the properties of the interstellar medium, the evolution of molecular clouds, and star formation in an active low-metallicity galaxy, and it also provides a link to galaxies at high redshift.

The physical properties of the molecular gas in the LMC, in particular the kinetic temperature, are not well constrained. The easily thermalized and optically thick rotational CO transitions are good temperature tracers of the cold and dense gas in local clouds. Generally, however, they often suffer from a lack of information on the beam-filling factor in extragalactic clouds. Multi-level observations of suitable molecules deliver the most reliable temperature determinations. The metastable lines of ammonia (NH3) are frequently used as a standard molecular cloud thermometer in molecular clouds within our Galaxy and also in external galaxies (Ho & Townes 1983; Walmsley & Ungerechts 1983; Danby et al. 1988; Henkel et al. 2000, 2008; Weiß et al. 2001; Mauersberger et al. 2003; Ao et al. 2011; Lebrón et al. 2011; Ott et al. 2011; Wienen et al. 2012; Mangum et al. 2013a). However, the ammonia abundance can vary strongly in different molecular environments (e.g., 10-5 in dense molecular “hot cores” around newly formed massive stars, Mauersberger et al. 1987; 10-8 in dark clouds, Benson & Myers 1983; Chira et al. 2013; 10-10 in a massive star-forming cloud of the LMC, Ott et al. 2010) and is extremely affected by a high UV flux. Thus, in a low-metallicity environment with high UV flux and a lack of shielding dust grains, ammonia is of limited use as a reliable probe to trace the gas kinetic temperature. To exacerbate this problem, the LMC is also a galaxy with particulary low nitrogen abundance (~10% solar, Wang et al. 2009; Ott et al. 2010).

Formaldehyde (H2CO) is a ubiquitous molecule in the Galactic interstellar medium (ISM) of our galaxy and external galaxies (Downes et al. 1980; Cohen & Few 1981; Bieging et al. 1982; Cohen et al. 1983; Baan et al. 1986, 1990, 1993; Henkel et al. 1991; Zylka et al. 1992; Hüttemeister et al. 1997; Heikkilä et al. 1999; Wang et al. 2004, 2009; Mangum et al. 2008, 2013b; Zhang et al. 2012; Ao et al. 2013; Tang et al. 2013; Ginsburg et al. 2015, 2016; Guo et al. 2016). H2CO is thought to be formed on the surface of dust grains by successive hydrogenation of CO (Watanabe & Kouchi 2002; Woon 2002; Hidaka et al. 2004): CO HCO H2CO. Variations in the fractional abundance of H2CO do not exceed one order of magnitude. For example, the fractional abundance of H2CO is similar across various subregions of the well-studied Orion-KL nebula, that is, in the “hot core” and the “compact ridge” (Mangum et al. 1990, 1993; Caselli et al. 1993; Johnstone et al. 2003).

Para-H2CO has a rich variety of millimeter and submillimeter transitions. Line ratios of para-H2CO involving different Ka ladders are good tracers of the kinetic temperature, such as para-H2CO JKaKc = 322–221/303–202, 423–322/404–303, and 523–422/505–404, since the relative populations of the Ka ladders of para-H2CO are predominantly governed by collisions (Mangum & Wootten 1993; Mühle et al. 2007). In these para-H2CO lines, the above three transitions with rest frequencies of 218.222, 218.475, and 218.760 GHz, respectively, are particularly suitable for use as a thermometer because they are strong enough for extragalactic observations and because they can be measured simultaneously within a bandwidth of 1 GHz. Temperatures determined from these ratios are free from uncertainties related to pointing accuracy, calibration, or different beam sizes. Since the line emission is optically thin and the levels are located at up to 68 K above the ground state, the line ratios are sensitive to gas kinetic temperatures of up to 50 K with relatively small uncertainties (Mangum & Wootten 1993; Ao et al. 2013). Para-H2CO 3–2 line ratios have been used to measure the molecular gas kinetic temperatures in our Galactic center (Qin et al. 2008; Ao et al. 2013; Johnston et al. 2014; Ginsburg et al. 2016; Immer et al. 2016), star formation regions (Mangum & Wootten 1993; Mitchell et al. 2001; Watanabe & Mitchell 2008; Tang et al. 2017), and in external galaxies (e.g., Mühle et al. 2007).

Multitransition observations of molecular clouds in the LMC (Johansson et al. 1998; Heikkilä et al. 1999; Israel et al. 2003; Kim et al. 2004; Bolatto et al. 2005; Pineda et al. 2008; Mizuno et al. 2010; Minamidani et al. 2008, 2011; Fujii et al. 2014; Paron et al. 2016) suggest that the molecular gas traced by CO may be warmer and/or denser than in our Galaxy. NH3(1, 1) and (2, 2) lines have been surveyed toward seven star-forming regions in the LMC by Ott et al. (2010) using the Australia Telescope Compact Array (ATCA). Emission is only detected in the massive star-forming region N159W. This represents the only detection of NH3 in the Magellanic Clouds to date. The gas kinetic temperature derived from NH3 (2, 2)/(1, 1) is cold (~16 K), which is twice lower than the derived dust temperature of 30–40 K (Heikkilä et al. 1999; Bolatto et al. 2000; Gordon et al. 2014). Ott et al. (2010) also found a low fractional NH3 abundance of ~4 × 10-10, which is lower by 1.5–5 orders of magnitude than those observed in Galactic star-forming regions. Previous observations of formaldehyde have been made in the LMC (e.g., Whiteoak & Gardner 1976; Johansson et al. 1994; Heikkilä et al. 1999; Wang et al. 2009), and H2CO has been detected in many dense clumps. These observations show that the fractional abundance of para-H2CO ranges from 1 to 6 × 10-10, which agrees with the values found in our Galactic molecular clouds (e.g., Güsten & Henkel 1983; Zylka et al. 1992; Johnstone et al. 2003; Ao et al. 2013; Tang et al. 2017).

Table 1

Source coordinates, integration times, and epochs.

For this paper, we have carried out deep observations of the six star-forming regions 30 Dor, N44BC, N113, N159E, N159S, and N159W in the LMC. Targeting three transitions of para-H2CO (JKAKC = 303–202, 322–221, and 321–220) as well as C18O 2–1, we simultaneously determine kinetic temperatures and spatial densities at high precision. In Sects. 2 and 3 we introduce our observations of the para-H2CO triplet and the data reduction, and describe the main results. These are then discussed in Sect. 4. Our main conclusions are summarized in Sect. 5.

thumbnail Fig. 1

Para-H2CO (303–202, 322–221, and 321–220) and C18O (2–1) spectra.

2. Observations and data reduction

Our observations were carried out in 2008 and 2014 (summarized in Table 1) with the Atacama Pathfinder EXperiment (APEX) 12 m telescope located on Chajnantor (Chile) using the APEX-1 (SHeFI) receiver. The beam size is ~30′′ (~7 pc at 50 kpc distance) at 218 GHz. The main beam efficiency and the forward efficiency were 0.75 and 0.97, respectively. N113 and N159W were observed in 2008 with an old fast Fourier transform spectrometer (FFTS), which consists of two units with a bandwidth of 1 GHz each and a velocity resolution of 0.1675 km s-1. The frequency is centered at 218.480 GHz. Our data include all three of the 218 GHz para-H2CO lines. 30 Dor, N159E, N159S, and N44BC were observed in 2014. Here we used the new extended bandwidth fast Fourier transform spectrometer (XFFTS) backend with two spectral windows of 2.5 GHz bandwidth and a velocity resolution of 0.1047 km s-1. The central frequency is set at 218.550 GHz. These data do not only include the three para-H2CO lines, but also the 219.560 GHz C18O (2–1) transition.

Toward each of the six star-forming regions in the LMC Wong et al. (2011) took a single-pointing high-sensitivity (~5 mK rms, Tmb scale, beam size ~30′′) spectrum centered on its CO emission peak, which is used to estimate N(H2) in Sect. 4.1. Ancillary C18O 2–1 data have been published by Heikkilä et al. (1998) and Wang et al. (2009). In addition, ammonia (NH3) data from the LMC (Ott et al. 2010) and Herschel infrared data also play an important role.

Data reduction for spectral lines was performed using CLASS from the GILDAS package1. To enhance the signal-to-noise ratios (S/N) in individual channels, we smoothed contiguous channels to a velocity resolution of ~1.0 km s-1. The observed sources are listed in Table 1.

Table 2

Para-H2CO and C18O spectral parameters.

3. Results

The three para-H2CO lines are detected in all sources except in N159S. There, only the strongest para-H2CO line, the 303–202 transition, is detected. C18O (2–1), measured in 30 Dor, N159E, N159S, and N44BC, is detected in all sources. The para-H2CO and C18O line spectra are presented in Fig. 1. Line parameters are listed in Table 2, where velocity-integrated intensity, Tmb\hbox{$\int T_{\rm mb}$}dv, local standard of rest velocity, Vlsr, full width at half maximum line width (FWHM), peak main beam brightness temperature, Tmb, and rms noise, were obtained from Gaussian fits.

3.1. Kinetic temperature and spatial density

To determine gas kinetic temperatures and spatial densities, we used the RADEX non-local thermal equilibrium (non-LTE) model (van der Tak et al. 2007) offline code2 with H2CO collision rates from Wiesenfeld & Faure (2013) and C18O collision rates from Yang et al. (2010). The RADEX code needs five input parameters: background temperature, kinetic temperature, H2 density, column density, and line width. For the background temperature, we adopted 2.73 K. Model grids for the para-H2CO and C18O lines encompass 30 densities (n(H2) = 103–107 cm-3) and 30 temperatures ranging from 10 to 110 K. For the line width, we used the observed line width value (Table 2). The total beam-averaged column density of C18O can be obtained from the J = 2–1 integrated intensity following Batrla & Wilson (2003), N(C18O)=5.3×1014Tmb(C18O,J=21)dv,\begin{equation} N({\rm C^{18}O}) = 5.3 \times 10^{14} \int T_{\rm mb} ({\rm C^{18}O}, J=2{-}1){\rm d}v, \end{equation}(1)where Tmb(C18O,J=21)dv\hbox{$\int T_{\rm mb} ({\rm C^{18}O}, J=2-1){\rm d}v$} is the C18O (2–1) integrated intensity. The results are listed in Table 3. In the LMC, the gas is well mixed so that isotope ratios are about the same throughout the galaxy (Johansson et al. 1994; Chin et al. 1997; Heikkilä et al. 1998, 1999; Wang et al. 2009). In the local ISM, the 16O/18O ratio is approximately 500 (Wilson & Rood 1994). We have to keep in mind, however, that fractional abundances in the LMC differ from those in the local ISM. In the LMC, 18O is underabundant with respect to 16O by about a factor of ~4 (16O/18O ~2000 in the LMC, Chin 1999; Wang et al. 2009), while 18O is underabundant by a factor of ~2.4 with respect to 17O (locally, 18O/17O ~4.1, Zhang et al. 2007; Wouterloot et al. 2008; in the LMC ~1.7, Heikkilä et al. 1998; Wang et al. 2009). Therefore, N(H2)/N(C18O) ratios should be higher by a factor of ~3 in the LMC with respect to the local ISM. This correlation is used to derive H2 column densities. The results for the H2 column density are listed in Table 3. Assuming a CO to H2 conversion factor of XCO ~ 4 × 1020 cm-2 (K km s-1)-1 for the LMC (Pineda et al. 2009; Wang et al. 2009), we can also derive a column density of N(H2) from the CO(1–0) integrated intensity moment map reported by Wong et al. (2011). The results of N(H2) are listed in Table 3. The H2 column densities derived from C18O and CO are consistent within the estimated uncertainty by a factor of two.

Previous observational results on transitions of C18O (1–0, 2–1) and para-H2CO (202–101, 303–202, 322–221 and 321–220) using the 15 m Swedish-ESO Submillimetre Telescope (SEST) show that the column density ratio of N(C18O)/N(para-H2CO) is ~100 in dense clumps of the LMC (Heikkilä et al. 1998, 1999; Wang et al. 2009). Assuming that the N(C18O)/N(para-H2CO) ratio is the same in our sample, we estimate the column density of para-H2CO from the N(C18O) column density. The results are listed in Table 3.

Table 3

Integrated intensity ratios, column densities, spatial densities, and temperatures.

The para-H2CO (303–202) line is the strongest of the three 218 GHz para-H2CO transitions observed by us. In order to avoid low uncertain values in the denominator, the para-H2CO 322–221/303–202 and 321–220/303–202 ratios are most suitable to derive the kinetic temperature. The para-H2CO 322–221 and 321–220 transitions have similar energy above the ground state, Eu ≃ 68 K, similar line brightness, and are often detected at the same time (e.g., Mühle et al. 2007; Bergman et al. 2011; Wang et al. 2012; Lindberg & Jørgensen 2012; Ao et al. 2013; Immer et al. 2014, 2016; Treviño-Morales et al. 2014; Ginsburg et al. 2016; Tang et al. 2017); therefore, para-H2CO 322–221/303–202 and 321–220/303–202 ratios are both good thermometers to determine the gas temperature. The kinetic temperature is traced by these two ratios with an uncertainty of 25% below 50 K (Mangum & Wootten 1993). While at n(H2) 105 cm-3 both ratios are similarly suitable, the para-H2CO 322–221/303–202 line ratio is affected by gas density at n(H2) < 105 cm-3 (Tang et al. 2017). For our sample, spatial density measurements probed by molecular tracers like CS, SO, CO, CI, H2CO, HCO+, and HCN (Heikkilä et al. 1998, 1999; Kim et al. 2004; Wang et al. 2009) show a range of 0.3–10 × 105 cm-3. Therefore, we use the para-H2CO 321–220/303–202 integrated intensity ratio to derive the kinetic temperature (Ginsburg et al. 2016; Immer et al. 2016).

thumbnail Fig. 2

RADEX non-LTE modeling of the gas kinetic temperature and spatial density, using the para-H2CO 321–220/303–202 (red solid and dashed lines are observed values and uncertainties) and para-H2CO 303–202/C18O 2–1 (blue solid and dashed lines) integrated intensity ratios. For details, see the notes to Table 3.

In our Galaxy, spatial densities, n(H2), derived from para-H2CO (303–202) are higher than from C18O (2–1). However, in the LMC, where lower density regions are often photoionized, this may be different. Therefore, here the para-H2CO 303–202/C18O 2–1 integrated intensity ratio is used to constrain the spatial density assuming the two tracers have a similar spatial extent and sample the same region. For both lines we find similar line parameters (e.g., Vlsr, FWHM, see Table 2) in our sample, so our assumption is reasonable. Figure 2 shows how the parameters are constrained by the line ratio distribution of para-H2CO 321–220/303–202 and para-H2CO 303–202/C18O 2–1 in the Tkinn(H2) parameter space. The determined results are listed in Table 3. The spatial density of our sample derived from the para-H2CO 303–202/C18O 2–1 ratio shows a relatively narrow range of 0.4–2.9 × 105 cm-3 with an average of 1.5 ± 0.4 × 105 cm-3 (see Fig. 2 and Table 3; errors given here and elsewhere are unweighted standard deviations of the mean), which is consistent with the results for the same dense clumps found from observations of SO, HCO+, c-C3H2, CH3OH, and H2CO (Heikkilä et al. 1999; Wang et al. 2009), for example.

As we described, our three para-H2CO (3–2) transitions are sensitive to gas at density ~105 cm-3 (Immer et al. 2016). To highlight how much the derived kinetic temperature depends on the derived density based on the admittedly uncertain assumption that para-H2CO 303–202 and C18O 2–1 trace the same gas, we have plotted in Fig. 3 the relation between kinetic temperature and the para-H2CO 321–220/303–202 ratio at spatial density n(H2) = 105 cm-3 with column density of N(para-H2CO) = 2.1 × 1012 cm-2 and line width of 6 km s-1 (these are rough average values for our sample) using RADEX. When this is compared with the actually obtained values for individual sources derived from our para-H2CO 321–220/303–202 and para-H2CO 303–202/C18O 2–1 ratios, the plot shows that the temperatures derived in the two different ways agree well.

LTE is a good approximation for the H2CO level populations under optically thin high-density conditions (Mangum & Wootten 1993; Watanabe & Mitchell 2008). Although LTE and RADEX non-LTE models use different approximations and assumptions, it is useful to check how the temperatures derived by the two methods differ. The para-H2CO line intensity ratios 322–221/303–202 and 321–220/303–202 can be used to measure the LTE kinetic temperature because (see Sect. 1) the Ka = 0 and 2 ladders of para-H2CO are mainly connected by collisions. The LTE kinetic temperature can be calculated assuming that the lines are optically thin and originate from a high-density region (Mangum & Wootten 1993), Tkin=47.1ln(0.556I(303202)I(322221)),\begin{equation} T_{\rm kin} = \frac{47.1}{\ln\left(0.556\frac{I(3_{03}-2_{02})}{I(3_{22}-2_{21})}\right)}, \end{equation}(2)where I(303–202)/I(322–221) is the para-H2CO integrated intensity ratio. The LTE kinetic temperatures are listed in Table 3. The uncertainty is 30% for this method of temperature measurement (Mangum & Wootten 1993). Considering this uncertainty, the temperatures derived from LTE and the RADEX non-LTE model are consistent with each other.

thumbnail Fig. 3

RADEX non-LTE modeling of the relation between the gas kinetic temperature and the para-H2CO 321–220/303–202 integrated intensity ratio (solid line), assuming n(H2) = 105 cm-3, N(para-H2CO) = 2.1 × 1012 cm-2, and a line width of 6 km s-1. Blue points are the results derived from the measured para-H2CO 321–220/303–202 and para-H2CO 303–202/C18O 2–1 ratios. For N159S, 3σ upper limits for the line ratio and the kinetic temperature have been taken.

3.2. Individual sources

Below we discuss the individual results for the six sources covered by this study.

3.2.1. 30 Dor

Our three para-H2CO and C18O 2–1 transitions have been observed toward 30 Dor by Heikkilä et al. (1999) with the SEST (beam size ~23′′), but only para-H2CO 303–202 was detected. With our higher sensitivity we detect all four lines and confirm their para-H2CO results. The kinetic temperature derived from para-H2CO (321–220/303–202) is ~63 K, which is the highest value determined in our sample. 30 Dor, hosting a cluster of O3 stars that rivals super star clusters, is the most spectacular star-forming region in the Local Group (Walborn & Blades 1997; Massey & Hunter 1998), which makes this high Tkin value comprehensible. The spatial density derived from the para-H2CO 303–202/C18O 2–1 ratio is ~1.5 × 105 cm-3 at this temperature, which agrees with results derived from other species (e.g., CS, SO, HCO+, Heikkilä et al. 1999).

3.2.2. N113

All three 218 GHz para-H2CO lines as well as C18O 21 transition have been observed by Wang et al. (2009) and Heikkilä et al. (1998) with the SEST. Para-H2CO 303–202 and C18O 2–1 were detected. We detect the three transitions of para-H2CO for the first time. The kinetic temperature derived from para-H2CO (321–220/303–202) is ~54 K. N113 is the LMC star-forming region with the most luminous 22 GHz H2O maser (Whiteoak & Gardner 1986; Lazendic et al. 2002; Oliveira et al. 2006), not quite as spectacular as 30 Dor, but it hosts a few bright HII regions whose stellar energy feedback is likely to have elevated its temperature, thus leading to the second highest Tkin value. The spatial density derived from para-H2CO 303–202/C18O 2–1 (ratio data from Wang et al. (2009) and Heikkilä et al. (1998), assuming that the para-H2CO 303–202/C18O 2–1 ratio at the SEST beam size, ~23′′, is similar to that in the APEX beam size, ~30′′, see Sect. 3.2.4) is ~8.9 × 104 cm-3 at temperature ~54 K, which agrees with results from other molecules (e.g., CS, HCO+, HCN, Wang et al. 2009).

3.2.3. N44BC

C18O 2–1 has been detected in N44BC by Heikkilä et al. (1998) with the SEST. We detect the three 218 GHz para-H2CO transitions and the C18O line. Our observations confirm their C18O 2–1 result. The kinetic temperature derived by para-H2CO (321–220/303–202) is ~47 K. Such a high temperature probably results from the stellar energy feedback from the adjacent superbubble on the molecular cloud, which also shows bright mid-IR emission (Chen et al. 2009). The spatial density derived from the para-H2CO 303–202/C18O 2–1 ratio is ~1.1 × 105 cm-3 at this temperature, which agrees with results from other molecules (e.g., CO, CI, Heikkilä et al. 1998; Kim et al. 2004).

3.2.4. N159W

The three para-H2CO 218 GHz transitions as well as C18O 2–1 have been detected in N159W by Heikkilä et al. (1998, 1999) with the SEST. Our observations confirm their para-H2CO results. The kinetic temperature derived by para-H2CO (321–220/303–202) is ~35 K. The spatial density derived from the para-H2CO 303–202/C18O 2–1 ratio (data from Heikkilä et al. 1998, 1999) is ~1.0 × 105 cm-3 at this temperature, which agrees with results measured by other species (e.g., CS, SO, Heikkilä et al. 1999). To quantify potential differences in temperature and density derived from APEX and SEST data, we determined the temperature and the density with the same method using the SEST para-H2CO and C18O data from Heikkilä et al. (1998, 1999). The derived temperature and density are 30.3+5.7-5.4\hbox{$^{+5.7}_{-5.4}$} K and 1.14+0.35-0.23×105\hbox{$^{+0.35}_{-0.23} \times10^5$} cm-3, respectively. This indicates nearly the same spatial density and a difference of a few Kelvin for the kinetic temperature. This temperature difference is similar to its 1σ uncertainty (~6 K). Therefore, the temperature gradient and density gradient is small when moving from a beam size of 7.3 pc (APEX) to 5.6 pc (SEST).

3.2.5. N159E

The three para-H2CO 218 GHz transitions as well as C18O 2–1 are detected. To our knowledge, this is the first detection of para-H2CO in the N159E region. The source shows similar properties as N159W. The kinetic temperature derived from para-H2CO (321–220/303–202) is ~37 K. The spatial density derived from the para-H2CO 303–202/C18O 2–1 ratio is ~2.9 × 105 cm-3 at this temperature.

3.2.6. N159S

Two velocity components are detected by C18O, at 232.4 and 236.4 km s-1. Para-H2CO 303–202 is detected at a velocity of 235.9 km s-1 (see Table 2). However, the para-H2CO 322–221 and 321–220 lines are not detected. The dust temperature ranges from 20 to 30 K (Gordon et al. 2014). N159S appears to be a cold cloud (Heikkilä et al. 1999) and has been shown to host no massive star formation at present and during the last 10 Myr (Chen et al. 2010). The upper limit to the kinetic temperature derived from para-H2CO (321–220/303–202) based on our observational 3σ limit for the 321–220 line is ~27 K. The spatial density derived from the para-H2CO 303–202/C18O 2–1 ratio is >4.3 × 104 cm-3, which is consistent with results measured by other species (e.g., CS, SO, Heikkilä et al. 1999).

4. Discussion

4.1. Comparison of temperatures derived from H2CO, CO, NH3, and the dust

The kinetic temperatures of molecular clumps in the LMC have been calculated by multitransition data of CO (Johansson et al. 1998; Heikkilä et al. 1999; Israel et al. 2003; Kim et al. 2004; Bolatto et al. 2005; Pineda et al. 2008; Mizuno et al. 2010; Minamidani et al. 2008, 2011; Fujii et al. 2014; Paron et al. 2016). These observations show that the higher temperatures (Tkin ≳ 100 K) occur in cloud regions that are of lower density (103 cm-3) and that the gas is colder (Tkin = 10–80 K) in regions of higher density (104–105 cm-3). As described in Sect. 3.1, the spatial density range of our sample derived from the para-H2CO 303–202/C18O 2–1 ratio with respect to the Galaxy is 0.4–2.9 × 105 cm-3 with an average of 1.5 ± 0.4 × 105 cm-3. Excluding the quiescent cloud N159S, where only one para-H2CO line could be detected, the gas kinetic temperatures derived from para-H2CO (321–220/303–202), range from 35 to 63 K with an average of 47 ± 5 K. Temperatures and densities derived from CO are for 30 Dor Tkin ~ 40–80 K and n(H2) ~ 3 × 1033 × 105 cm-3 (Johansson et al. 1998; Heikkilä et al. 1999; Minamidani et al. 2008), N159W Tkin ~ 16− > 30 K and n(H2) ~ 3 × 1038 × 105 cm-3 (Johansson et al. 1998; Heikkilä et al. 1999; Bolatto et al. 2005; Minamidani et al. 2008), N159E Tkin> 40 K and n(H2) ~ 1 × 1033 × 105 cm-3 (Minamidani et al. 2008), and N159S Tkin ~ 10–60 K and n(H2) ~ 1 × 1031 × 105 cm-3 (Heikkilä et al. 1999; Minamidani et al. 2008). Temperatures derived from para-H2CO are consistent with but much more precise than the results derived from CO in the dense regions (>103 cm-3).

Except for N159E, all our sources have been surveyed in NH3 (1, 1) and (2, 2) by Ott et al. (2010). These lines are only detected in the massive star-forming region N159W. We compare fractional abundances of N(para-NH3)/N(H2) and N(para-H2CO)/N(H2) to column density N(H2) obtained from 12CO(1–0) in Fig. 4. These show that formaldehyde has a stable fractional abundance ranging from 0.6 to 5.7 × 10-10 cm-2 with an average of 2.7 ± 1.8 × 10-10 cm-2 in molecular clouds of the LMC with N(H2) column densities ranging from 0.4 to 2 × 1022 cm-2. Ammonia only survives in a high column density environment with N(H2) ~ 2 × 1022 cm-2. The fractional abundance of ammonia is ~10-10–10-9 in N159W and M 82, which is similar to that of formaldehyde in the LMC. The kinetic temperature derived from the NH3 (2, 2)/(1, 1) line ratio in N159W is ~16 K (Ott et al. 2010), which is twice lower than that derived from para-H2CO. Previous para-H2CO (322–221/303–202) and NH3 (2, 2)/(1, 1) observations toward the starburst galaxy M 82 also show significantly different gas kinetic temperatures (Weiß et al. 2001; Mühle et al. 2007). M 82, a satellite galaxy like the LMC, shows a similar environment, involving low metallicity combined with a high UV flux. This only leaves NH3 surviving in the most shielded pockets of molecular gas, resulting in a low fractional abundance (see Fig. 4) and a low kinetic temperature. Furthermore, this abundance demonstrates in an exemplary way that H2CO is less affected by photodissociation and samples a more extended region. Therefore, para-H2CO traces in these instances a higher temperature than NH3 (2, 2)/(1, 1). We conclude that para-H2CO line ratios are a superior thermometer for tracing dense gas temperatures in low-metallicity galaxies with strong UV radiation. Nevertheless, more detailed spatially resolved comparisons of H2CO with NH3 temperatures would be very interesting because this could provide more information on temperature gradients and the location of different kinetic temperature layers.

thumbnail Fig. 4

Abundances of N(para-NH3)/N(H2), red squares, and N(para-H2CO)/N(H2), blue points, vs. column density N(H2) obtained from 12CO. 30 Dor-27 and N160 data taken from Heikkilä et al. (1999). N159W data taken from Ott et al. (2010). M 82 data taken from Weiß et al. (2001) and Mühle et al. (2007).

The temperatures derived from dust and gas are often in agreement in the active and dense clumps of Galactic disk clouds (Dunham et al. 2010; Giannetti et al. 2013; Battersby et al. 2014). However, observed gas and dust temperatures do not agree with each other in the Galactic Central Molecular Zone (CMZ) and external galaxies (Güsten et al. 1981; Ao et al. 2013; Mangum et al. 2013a; Ott et al. 2014; Ginsburg et al. 2016; Immer et al. 2016). Dust temperatures in the LMC have been obtained by Gordon et al. (2014) using Herschel 100 to 500 μm dust continuum emission data. They range approximately from 13 to 73 K. For our sample, the dust temperatures range from 30 to 75 K (Werner et al. 1978; Heikkilä et al. 1999; Bolatto et al. 2000; Wang et al. 2009; Seale et al. 2014; Gordon et al. 2014), while the para-H2CO derived temperature ranges from 35 to 63 K. The temperatures derived from para-H2CO ratios and dust emission are therefore in good agreement (see Table 3). This indicates that the dust and H2CO kinetic temperatures are equivalent in the star-forming regions of the LMC. Assuming that H2CO traces the bulk of the dense molecular gas and that ammonia shows very low abundances, this can be generalized in the sense that dense gas and dust temperatures are generally equivalent.

4.2. Star-forming regions in the Galaxy, the LMC, and other galaxies

The gas temperatures of APEX Telescope Large Area Survey of the GALaxy (ATLASGAL) massive star-forming clumps have been measured by para-H2CO (303–202, 322–221, and 321–220) line ratios (Tang et al. 2017). The gas kinetic temperatures derived thus at density n(H2) = 105 cm-3 with size of 0.3–0.7 pc range from 30 to 61 K with an average of 46 ± 9 K, which agrees remarkably well with the results in the LMC with a beam size of ~7 pc. Large area mapping measurements of kinetic temperatures in the Galactic CMZ with the same transitions of para-H2CO (Ao et al. 2013; Ginsburg et al. 2016) suggest that the mean gas temperature is ~48 K at n(H2) = 105 cm-3 or ~65 K at n(H2) = 104 cm-3 in the whole ~300 pc surveyed region. It shows a higher Tkin value than our observed results. The spatial densities derived from the para-H2CO 303–202/C18O 2–1 ratio in our sample agree with the observed results in the Galactic clumps (Beuther et al. 2002; Motte et al. 2003; Wienen et al. 2012, 2015), HII regions (Henkel et al. 1983; Ginsburg et al. 2011), and giant molecular clouds (GMCs; Wadiak et al. 1988; Ginsburg et al. 2015; Immer et al. 2016). This agreement indicates that the physical conditions of the star-forming regions should be similar in both the LMC and our Galactic disk.

When the three transitions of para-H2CO at ~218 GHz are used to measure the kinetic temperature of the starburst galaxy M 82, we see that the derived kinetic temperature (Tkin(H2CO) ~ 200 K; Mühle et al. 2007) is significantly higher than in the LMC. Kinetic temperatures of starburst galaxies measured with multi-inversion transitions of NH3 show a range from 24 to 250 K (Henkel et al. 2000, 2008; Mauersberger et al. 2003; Ao et al. 2011; Lebrón et al. 2011; Mangum et al. 2013a). The temperatures derived from para-H2CO line ratios in the LMC overlap with the values found for a sample at lower temperature (e.g., M83, NGC 6946). This is most likely due to the inclusion of higher excited ammonia lines, which are probably difficult to detect in the LMC, however, because particularly warm regions irradiated by enhanced UV radiation should be almost devoid of NH3 there. The spatial densities in starburst galaxies derived from ortho-H2CO (211–212/110–111) line ratios (Mangum et al. 2008, 2013b) show a similar range (104.5–105.5 cm-3) to that derived from para-H2CO 303–202/C18O 2–1 ratios in the LMC. This indicates that star formation in the LMC and external galaxies may arise from dense molecular gas (>104 cm-3), but gas-heating rates may be quite different.

4.3. Correlation of gas temperature with star formation

In star-forming galaxies, a lack of correlation between the gas kinetic temperatures traced by NH3 and star formation rate indicated by infrared luminosity was found by Mangum et al. (2013a). To investigate how the kinetic temperatures traced by para-H2CO correlate with star formation in the LMC, we compared the gas kinetic temperature to the Herschel 250 μm dust emission, which indicates the infrared luminosity (Seale et al. 2014). We averaged the Herschel 250 μm data over a 30′′ beam corresponding to our para-H2CO data. A comparison between gas kinetic temperatures derived from para-H2CO and the Herschel 250 μm flux is shown in Fig. 5. It reveals a correlation of gas temperature and 250 μm flux for the five sources with all 218 GHz para-H2CO lines detected (slope = 0.97 ± 0.49, correlation coefficient R ~ 0.75). According to the relation between the infrared luminosity and 250 μm flux LFIRFS2501.2\hbox{$F_{\rm S250}^{1.2}$} in the LMC (Seale et al. 2014), the infrared luminosity and the gas temperature derived from para-H2CO are related by a power law of the form LFIRTkin1.2±0.59\hbox{$T_{\rm kin}^{1.2\pm0.59}$} , where the power-law index is lower than that of the Stefan-Boltzmann law (LTkin4\hbox{$T_{\rm kin}^4$}). This suggests that this picture is an oversimplification, and that star formation occurs in extended regions, which leads to the formation of stellar clusters with multiple far-infrared (FIR) sources (e.g., Chen et al. 2009, 2010). Assuming that the 250 μm flux is mostly emitted from clusters of FIR sources that are distributed across the regions from which the H2CO emission arises, this can be generalized in the sense that the gas-heating mechanism must be related to the formation of young massive stars. To determine whether the bulk of the H2CO emission originates from within individual clusters of FIR sources, from in between adjacent FIR clusters, or from both, higher resolution H2CO observations are mandatory.

thumbnail Fig. 5

Comparison of gas kinetic temperatures derived from para-H2CO 321–220/303–202 against the Herschel 250 μm flux. The straight line is the result from a linear fit for the five sources with all 218 GHz para-H2CO lines detected. N159S provides with respect to Tkin a 3σ upper limit. Thus this point may be located much farther to the left, closer to the linear fit obtained from the other five sources.

We need more data points in the LMC to understand the relationship between Tkin and star formation and to then apply this relationship to more distant galaxies with ALMA. Para-H2CO 321–220/303–202 and para-H2CO 303–202/C18O 2–1 line ratios provide a direct estimate of the gas kinetic temperatures and spatial densities for molecular gas on a scale of ~7 pc in the star-forming regions of the LMC. It would be meaningful to use these line ratios to measure the physical properties of the dense molecular gas at smaller linear scales with ALMA and to start systematic investigations in more distant star-forming galaxies.

5. Summary

We have measured the kinetic temperature and spatial density with para-H2CO (JKAKC = 303–202, 322–221, and 321–220) line ratios and the C18O 2–1 line in massive star-forming regions of the LMC. Kinetic temperatures derived from the formaldehyde 218 GHz line triplet were compared with those obtained from the dust, and in one case, also with those from ammonia using the 12 m APEX telescope. The main results are the following:

  • 1.

    Using the RADEX non-LTE program, we derived gas kinetictemperatures and spatial densities modeling the measuredpara-H2CO 321–220/303–202 and para-H2CO 303–202/C18O 2–1 line ratios.

  • 2.

    The gas kinetic temperatures derived from para-H2CO (321–220/303–202) line ratios of the star-forming regions in the LMC are warm, ranging from 35 to 63 K with an average of 47 ± 5 K, which is similar to the temperature obtained from Galactic disk massive star-forming clumps.

  • 3.

    The spatial density derived from the para-H2CO 303–202/C18O 2–1 ratio shows a range of 0.4–2.9 × 105 cm-3 with an average of 1.5 ± 0.4 × 105 cm-3. It agrees with results measured by ortho-H2CO (211–212/110–111) line ratios in Galactic regions of massive star formation.

  • 4.

    Temperatures derived from the para-H2CO line ratios agree with those derived from CO in dense regions (n(H2) > 103 cm-3). The gas temperature derived from the NH3 (2, 2)/(1, 1) line ratio is ~16 K in N159W (Ott et al. 2010), which is twice lower than the temperature derived from the para-H2CO line ratio and the dust. Ammonia only survives in the most shielded pockets of molecular gas in the low-metallicity environment of the LMC, which is affected by a high UV flux. Formaldehyde is less affected by photodissociation and traces a more extended region.

  • 5.

    A comparison of the gas kinetic temperature derived from para-H2CO and the temperature obtained from dust emission shows good agreement. This indicates that the bulk of the dense gas and dust are in approximate thermal equilibrium in the dense star formation regions of the LMC.

  • 6.

    A correlation between the gas kinetic temperatures derived from para-H2CO and infrared luminosity indicated by the 250 μm flux suggests that our kinetic temperatures traced by para-H2CO are closely associated with extended star formation in the LMC.


Acknowledgments

We thank the staff of the APEX telescope for their assistance in observations. The authors are also thankful for the helpful comments of the anonymous referee. This work was funded by The National Natural Science Foundation of China under grant 11433008 and The Program of the Light in Chinas Western Region (LCRW) under grant Nos. XBBS201424 and The National Natural Science Foundation of China under grant 11373062. C.H. acknowledges support by a visiting professorship for senior international scientists of the Chinese Academy of Sciences (2013T2J0057). This research has used NASA’s Astrophysical Data System (ADS).

References

  1. Ao, Y., Henkel, C., Braatz, J. A., et al. 2011, A&A, 529, A154 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Ao, Y., Henkel, C., Menten, K. M., et al. 2013, A&A, 550, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Baan, W. A., Güsten, R., & Haschick, A. D. 1986, ApJ, 305, 830 [NASA ADS] [CrossRef] [Google Scholar]
  4. Baan, W. A., Henkel, C., Schilke, P., et al. 1990, ApJ, 353, 132 [NASA ADS] [CrossRef] [Google Scholar]
  5. Baan, W. A., Haschick, A. D., & Uglesich, R. 1993, ApJ, 415, 140 [NASA ADS] [CrossRef] [Google Scholar]
  6. Batrla, W., & Wilson, T. L. 2003, A&A, 408, 231 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Battersby, C., Bally, J., Dunham, M., et al. 2014, ApJ, 786, 116 [NASA ADS] [CrossRef] [Google Scholar]
  8. Benson, P. J., & Myers, P. C. 1983, ApJ, 270, 589 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bergman, P., Parise, B., Liseau, R., & Larsson, B. 2011, A&A, 527, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Beuther, H., Schilke, P., Menten, K. M., et al. 2002, ApJ, 566, 945 [NASA ADS] [CrossRef] [Google Scholar]
  11. Bieging, J. H., Wilson, T. L., & Downes, D. 1982, A&AS, 49, 607 [NASA ADS] [Google Scholar]
  12. Bolatto, A. D., Jackson, J. M., Israel, F. P., et al. 2000, ApJ, 545, 234 [NASA ADS] [CrossRef] [Google Scholar]
  13. Bolatto, A. D., Israel, F. P., & Martin, C. L. 2005, ApJ, 633, 210 [NASA ADS] [CrossRef] [Google Scholar]
  14. Caselli, P., Hasegawa, T. I., & Herbst, E. 1993, ApJ, 408, 548 [NASA ADS] [CrossRef] [Google Scholar]
  15. Chen, C.-H. R., Chu, Y. H., Gruendl, R. A., et al. 2009, ApJ, 695, 511 [NASA ADS] [CrossRef] [Google Scholar]
  16. Chen, C.-H. R., Indebetouw, R., Chu, Y. H., et al. 2010, ApJ, 721, 1206 [NASA ADS] [CrossRef] [Google Scholar]
  17. Chin, Y. 1999, IAU Symp., 190, 279 [NASA ADS] [Google Scholar]
  18. Chin, Y. N., Henkel, C., Whiteoak, J. B., et al. 1997, A&A, 317, 548 [NASA ADS] [Google Scholar]
  19. Chira, R. A., Beuther, H., Linz, H., et al. 2013, A&A, 552, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Cohen, R. J., & Few, R. W. 1981, MNRAS, 194, 711 [NASA ADS] [CrossRef] [Google Scholar]
  21. Cohen, R. J., Matthews, N., Few, R. W., et al. 1983, MNRAS, 203, 1123 [NASA ADS] [CrossRef] [Google Scholar]
  22. Danby, G., Flower, D. R., Valiron, P., et al. 1988, MNRAS, 235, 229 [NASA ADS] [CrossRef] [Google Scholar]
  23. Downes, D., Wilson, T. L., Bieging, J., & Wink, J. 1980, A&AS, 40, 379 [NASA ADS] [Google Scholar]
  24. Dunham, M. K., Rosolowsky, E., Evans, N. J. II., et al. 2010, ApJ, 717, 1157 [NASA ADS] [CrossRef] [Google Scholar]
  25. Fujii, K., Minamidani, T., Mizuno, N., Onishi, T., et al. 2014, ApJ, 796, 123 [NASA ADS] [CrossRef] [Google Scholar]
  26. Frerking, M. A., Langer, W. D., & Wilson, R. W. 1982, ApJ, 262, 590 [NASA ADS] [CrossRef] [Google Scholar]
  27. Giannetti, A., Brand, J., Sánchez-Monge, Á., et al. 2013, A&A, 556, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Ginsburg, A., Darling, J., Battersby, C., et al. 2011, ApJ, 736, 149 [NASA ADS] [CrossRef] [Google Scholar]
  29. Ginsburg, A., Bally, J., Battersby, C., et al. 2015, A&A, 573, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Ginsburg, A., Henkel, C., Ao, Y., et al. 2016, A&A, 586, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Gordon, K. D., Roman-Duval, J., Bot, C., et al. 2014, ApJ, 797, 85 [NASA ADS] [CrossRef] [Google Scholar]
  32. Guo, W. H., Esimbek, J., Tang, X. D., et al. 2016, Ap&SS, 361, 264 [NASA ADS] [CrossRef] [Google Scholar]
  33. Güsten, R., & Henkel, C. 1983, A&A, 125, 136 [NASA ADS] [Google Scholar]
  34. Güsten, R., Walmsley, C. M., & Pauls, T. 1981, A&A, 103, 197 [NASA ADS] [Google Scholar]
  35. Heikkilä, A., Johansson, L. E. B., & Olofsson, H. 1998, A&A, 332, 493 [NASA ADS] [Google Scholar]
  36. Heikkilä, A., Johansson, L. E. B., & Olofsson, H. 1999, A&A, 344, 817 [NASA ADS] [Google Scholar]
  37. Henkel, C., Wilson, T. L., Walmsley, C. M., & Pauls, T. 1983, A&A, 127, 388 [NASA ADS] [Google Scholar]
  38. Henkel, C., Baan, W. A., & Mauersberger, R. 1991, A&ARv, 3, 47 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  39. Henkel, C., Mauersberger, R., Peck, A. B., et al. 2000, A&A, 361, L45 [NASA ADS] [Google Scholar]
  40. Henkel, C., Braatz, J. A., Menten, K. M., & Ott, J. 2008, A&A, 485, 451 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Hidaka, H., Watanabe, N., Shiraki, T., et al. 2004, ApJ, 614, 1124 [NASA ADS] [CrossRef] [Google Scholar]
  42. Ho, P. T. P., & Townes, C. H. 1983, ARA&A, 21, 239 [NASA ADS] [CrossRef] [Google Scholar]
  43. Hüttemeister, S., Mauersberger, R., & Henkel, C. 1997, A&A, 326, 59 [NASA ADS] [Google Scholar]
  44. Immer, K., Galván-Madrid, R., König, C., et al. 2014, A&A, 572, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Immer, K, Kauffmann, J., Pillai, T., et al. 2016, A&A, 595, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Israel, F. P., Johansson, L. E. B., Rubio, M., et al. 2003, A&A, 406, 817 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Johansson, L. E. B., Olofsson, H., Hjalmarson, A., et al. 1994, A&A, 291, 89 [NASA ADS] [Google Scholar]
  48. Johansson, L. E. B., Greve, A., Booth, R. S., et al. 1998, A&A, 331, 857 [NASA ADS] [Google Scholar]
  49. Johnston, K. G., Beuther, H., Linz, H., et al. 2014, A&A, 568, A56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Johnstone, D., Boonman, A. M. S., & van Dishoeck, E. F., 2003, A&A, 412, 157 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Kim, S., Walsh, W., & Xiao, K. 2004, ApJ, 616, 865 [NASA ADS] [CrossRef] [Google Scholar]
  52. Lazendic, J. S., Whiteoak, J. B., Klamer, I., et al. 2002, MNRAS, 331, 969 [NASA ADS] [CrossRef] [Google Scholar]
  53. Lebrón, M., Mangum, J. G., Mauersberger, R., et al. 2011, A&A, 534, A56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Lindberg, J. E., & Jørgensen, J. K. 2012, A&A, 548, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Mangum, J. G., & Wootten, A. 1993, ApJS, 89, 123 [NASA ADS] [CrossRef] [Google Scholar]
  56. Mangum, J. G., Wootten, A., Wadiak, E. J., & Loren, R. B. 1990, ApJ, 348, 542 [NASA ADS] [CrossRef] [Google Scholar]
  57. Mangum, J. G., Wootten, A., & Plambeck, R. L. 1993, ApJ, 409, 282 [NASA ADS] [CrossRef] [Google Scholar]
  58. Mangum, J. G., Darling, J., Menten, K. M., & Henkel, C. 2008, ApJ, 673, 832 [NASA ADS] [CrossRef] [Google Scholar]
  59. Mangum, J. G., Darling, J., Henkel, C., et al. 2013a, ApJ, 779, 33 [NASA ADS] [CrossRef] [Google Scholar]
  60. Mangum, J. G., Darling, J., Henkel, C., & Menten, K. M. 2013b, ApJ, 766, 108 [NASA ADS] [CrossRef] [Google Scholar]
  61. Massey, P., & Hunter, D. A. 1998, ApJ, 493, 180 [NASA ADS] [CrossRef] [Google Scholar]
  62. Mauersberger, R., Henkel, C., & Wilson, T. L. 1987, A&A, 173, 352 [NASA ADS] [Google Scholar]
  63. Mauersberger, R., Henkel, C., Weiß, A., et al. 2003, A&A, 403, 561 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Minamidani, T., Mizuno, N., Mizuno, Y., et al. 2008, ApJS, 175, 485 [NASA ADS] [CrossRef] [Google Scholar]
  65. Minamidani, T., Tanaka, T., Mizuno, Y., et al. 2011, AJ, 141, 73 [NASA ADS] [CrossRef] [Google Scholar]
  66. Mitchell, G. F., Johnstone, D., Moriarty-Schieven, G., et al. 2001, ApJ, 556, 215 [NASA ADS] [CrossRef] [Google Scholar]
  67. Mizuno, Y., Kawamura, A., Onishi, T., et al. 2010, PASJ, 62, 51 [NASA ADS] [Google Scholar]
  68. Motte, F., Schilke, P., & Lis, D. C. 2003, ApJ, 582, 277 [NASA ADS] [CrossRef] [Google Scholar]
  69. Mühle, S., Seaquist, E. R., & Henkel, C. 2007, ApJ, 671, 1579 [NASA ADS] [CrossRef] [Google Scholar]
  70. Oliveira, J. M., van Loon, J. Th., Stanimirović, S., & Zijlstra, A. A. 2006, MNRAS, 372, 1509 [NASA ADS] [CrossRef] [Google Scholar]
  71. Ott, J., Henkel, C., Staveley-Smith, L., & Weiß, A. 2010, ApJ, 710, 105 [NASA ADS] [CrossRef] [Google Scholar]
  72. Ott, J., Henkel, C., Braatz, J. A., & Weiß, A. 2011, ApJ, 742, 95 [NASA ADS] [CrossRef] [Google Scholar]
  73. Ott, J., Weiß, A., Staveley-Smith, L., et al. 2014, ApJ, 785, 16 [Google Scholar]
  74. Paron, S., Ortega, M. E., Fariña, C., et al. 2016, MNRAS, 455, 518 [NASA ADS] [CrossRef] [Google Scholar]
  75. Pineda, J. L., Mizuno, N., Stutzki, J., et al. 2008, A&A, 482, 197 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  76. Pineda, J. L., Ott, J., Klein, U., et al. 2009, ApJ, 703, 736 [NASA ADS] [CrossRef] [Google Scholar]
  77. Pietrzyński, G., Graczyk, D., Gieren, W., et al. 2013, Nature, 495, 76 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  78. Qin, S. L., Zhao, J. H., Moran, J. M., et al. 2008, ApJ, 677, 353 [NASA ADS] [CrossRef] [Google Scholar]
  79. Rolleston, W. R. J., Trundle, C., & Dufton, P. L. 2002, A&A, 396, 53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  80. Seale, J. P., Meixner, M., Sewiło, M., et al. 2014, AJ, 148, 124 [NASA ADS] [CrossRef] [Google Scholar]
  81. Tang, X. D., Esimbek, J., Zhou, J. J., et al. 2013, A&A, 551, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  82. Tang, X. D., Henkel, C., Menten, K. M., et al. 2017, A&A, 598, A30 (Paper I) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  83. Treviño-Morales, S. P., Pilleri, P., Fuente, A., et al. 2014, A&A, 569, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  84. van der Tak, F. F. S., Black, J. H., Schöier, F. L., et al. 2007, A&A, 468, 627 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  85. Wadiak, E. J., Rood, R. T., & Wilson, T. L. 1988, ApJ, 324, 931 [NASA ADS] [CrossRef] [Google Scholar]
  86. Walborn, N. R., & Blades, J. C. 1997, ApJS, 112, 457 [NASA ADS] [CrossRef] [Google Scholar]
  87. Walmsley, C. M., & Ungerechts, H. 1983, A&A, 122, 164 [NASA ADS] [Google Scholar]
  88. Wang, M., Henkel, C., Chin, Y.-N., et al. 2004, A&A, 422, 883 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  89. Wang, M., Chin, Y.-N., Henkel, C., et al. 2009, ApJ, 690, 580 [NASA ADS] [CrossRef] [Google Scholar]
  90. Wang, Y., Beuther, H., Zhang, Q., et al. 2012, ApJ, 754, 87 [NASA ADS] [CrossRef] [Google Scholar]
  91. Watanabe, N., & Kouchi, A. 2002, ApJ, 571, 173 [Google Scholar]
  92. Watanabe, T., & Mitchell, G. 2008, AJ, 136, 1947 [NASA ADS] [CrossRef] [Google Scholar]
  93. Weiß, A., Neininger, N., Henkel, C., Stutzki, J., & Klein, U. 2001, ApJ, 554, 143 [Google Scholar]
  94. Werner, M.W., Becklin, E.E., Gatley, I., et al. 1978, MNRAS, 184, 365 [NASA ADS] [CrossRef] [Google Scholar]
  95. Westerlund, B.E., 1990, A&ARv, 2, 29 [NASA ADS] [CrossRef] [Google Scholar]
  96. Whiteoak, J. B., & Gardner, F. F. 1976, MNRAS, 174, 51 [NASA ADS] [CrossRef] [Google Scholar]
  97. Whiteoak, J. B., & Gardner, F. F. 1986, MNRAS, 222, 513 [NASA ADS] [CrossRef] [Google Scholar]
  98. Wienen, M., Wyrowski, F., Schuller, F., et al. 2012, A&A, 544, A146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  99. Wienen, M., Wyrowski, F., Menten, K. M., et al. 2015, A&A, 579, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  100. Wiesenfeld, L., & Faure, A. 2013, MNRAS, 432, 2573 [NASA ADS] [CrossRef] [Google Scholar]
  101. Wilson, T. L., & Rood, R. 1994, ARA&A, 32, 191 [NASA ADS] [CrossRef] [Google Scholar]
  102. Wong, T., Hughes, A., Ott, J., et al. 2011, ApJS, 197, 16 [NASA ADS] [CrossRef] [Google Scholar]
  103. Woon, D. 2002, ApJ, 569, 541 [NASA ADS] [CrossRef] [Google Scholar]
  104. Wouterloot, J. G. A., Henkel, C., Brand, J., Davis, G. R. 2008, A&A, 487, 237 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  105. Yang, B., Stancil, P.C., Balakrishnan, N., & Forrey, R. C. 2010, ApJ, 718, 1062 [NASA ADS] [CrossRef] [Google Scholar]
  106. Zhang, J. S., Henkel, C., Mauersberger, R., et al. 2007, A&A, 465, 887 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  107. Zhang, C. P., Esimbek, J., Zhou, J. J., et al. 2012, Ap&SS, 337, 283 [NASA ADS] [CrossRef] [Google Scholar]
  108. Zylka, R., Güsten, R., Henkel, C., & Batrla, W. 1992, A&AS, 96, 525 [NASA ADS] [Google Scholar]

All Tables

Table 1

Source coordinates, integration times, and epochs.

Table 2

Para-H2CO and C18O spectral parameters.

Table 3

Integrated intensity ratios, column densities, spatial densities, and temperatures.

All Figures

thumbnail Fig. 1

Para-H2CO (303–202, 322–221, and 321–220) and C18O (2–1) spectra.

In the text
thumbnail Fig. 2

RADEX non-LTE modeling of the gas kinetic temperature and spatial density, using the para-H2CO 321–220/303–202 (red solid and dashed lines are observed values and uncertainties) and para-H2CO 303–202/C18O 2–1 (blue solid and dashed lines) integrated intensity ratios. For details, see the notes to Table 3.

In the text
thumbnail Fig. 3

RADEX non-LTE modeling of the relation between the gas kinetic temperature and the para-H2CO 321–220/303–202 integrated intensity ratio (solid line), assuming n(H2) = 105 cm-3, N(para-H2CO) = 2.1 × 1012 cm-2, and a line width of 6 km s-1. Blue points are the results derived from the measured para-H2CO 321–220/303–202 and para-H2CO 303–202/C18O 2–1 ratios. For N159S, 3σ upper limits for the line ratio and the kinetic temperature have been taken.

In the text
thumbnail Fig. 4

Abundances of N(para-NH3)/N(H2), red squares, and N(para-H2CO)/N(H2), blue points, vs. column density N(H2) obtained from 12CO. 30 Dor-27 and N160 data taken from Heikkilä et al. (1999). N159W data taken from Ott et al. (2010). M 82 data taken from Weiß et al. (2001) and Mühle et al. (2007).

In the text
thumbnail Fig. 5

Comparison of gas kinetic temperatures derived from para-H2CO 321–220/303–202 against the Herschel 250 μm flux. The straight line is the result from a linear fit for the five sources with all 218 GHz para-H2CO lines detected. N159S provides with respect to Tkin a 3σ upper limit. Thus this point may be located much farther to the left, closer to the linear fit obtained from the other five sources.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.