Deuterated species in extragalactic star-forming regions
We present a theoretical study of the deuterated species detectability in various types of extragalactic star-forming regions based on our predictions of chemical abundances. This work is motivated by the past and current attempts at observing deuterated species in external galaxies such as NGC 253, IC 342 and the LMC. Here, we investigate the influence of the density, the temperature, the FUV radiation field, the cosmic ray ionisation, and the metallicity on the fractional abundances and D/H abundance ratios of about 20 deuterated species. Without modelling any particular source, we determined how the deuterium chemistry behaves in different physical environments such as starburst, cosmic-rays enhanced environments, low metallicity and high redshift galaxies. In general, our predicted column densities seem in good agreement with those derived from the current limited dataset of observations in external galaxies. We provide, for the first time, a list of key deuterated species whose abundances are high enough to be possibly detectable by the Atacama Large Millimeter Array (ALMA) and Herschel, as a function of galactic nuclear activity and redshift.
The past few years have seen a boom in studies of molecular deuteration in galactic low-mass and high-mass star forming regions, triggered by the discovery of a large fraction of mono-, doubly- and even multiply-deuterated molecules (e.g. Loren & Wootten 1985; Vrtilek et al. 1985; Walmsley et al. 1987; Mauersberger et al. 1988; Jacq et al. 1990; Gerin et al. 1992; Schilke et al. 1992; Jacq et al. 1993; van Dishoeck et al. 1995; Ceccarelli et al. 1998; Stark et al. 1999; Roueff et al. 2000; Parise et al. 2002; Bacmann et al. 2003; Caselli et al. 2003; Vastel et al. 2004; Roueff et al. 2005; Miettinen et al. 2009). Studying deuterium chemistry, especially in highly embedded regions, is important because in such environments, CO is depleted and deuterated species are thus driving the chemistry. In embedded regions in our Galaxy, a large deuterium fractionation is observed, well above the elemental abundance ratio D/H of (Oliveira, 2003). High molecular deuteration can occur by two different mechanisms. One involves gas phase chemistry and ion-molecule deuterium exchange reactions taking place at low temperatures (Watson, 1980). The second mechanism is based on grain chemistry (Tielens, 1983).
In this work, we investigate deuteration in extragalactic star-forming environments, where the physical conditions can be drastically different from what we see in our own Galaxy. To date, very few detections of deuterated species have been attempted in extragalactic environments (e.g. Mauersberger et al. 1995; Martín et al. 2006). In a previous paper, Bayet et al. (2008) investigated the non-deuterated chemistry in a variety of extragalactic environments and obtained theoretical predictions for the warmer and dense gas, traced by non-deuterated species. Here we implement the mono-deuterated species (hereafter called D-species) network to our existing one used in Bayet et al. (2008), in order to investigate the differences in the D-fractionation and the D/H ratios in regions of star-formation (i.e. gas with n(H)cm) in various types of galaxies (i.e starburst, cosmic-rays enhanced environments, low metallicity and high redshift). We also aim at guiding observers with potentially detectable deuterated tracers of so called dense star-forming molecular gas in external galaxies when using Herschel or the forthcoming ALMA, as based on the strength of their chemical abundances.
The paper is divided in five sections; Section 2 describes the model used and presents the choice of model parameters used for mimicking various extragalactic environments. Section 3 presents the D-species fractional abundances and D/H abundance ratios obtained in the various gas components and environments investigated. Finally, in Section 4, we discuss the results and conclude.
2 Model descriptions
We have used the same chemical model as the one used in Bayet et al. (2008), but with an extended chemical network to include mono-deuterated species. The chemical code, UCL_Chem, is presented in Viti & Williams (1999) and Viti et al. (2004). We refer the reader to these papers for a full description of the code.
To summarise, UCL_Chem model is a time and depth dependent chemical model. As in Bayet et al. (2008), we first model the collapse of a 10K core (Phase I); we then follow the chemical evolution of the region once the stars are born (Phase II). The temperature in both Phases is the typical average for starless and hot cores, respectively, since here we are considering an ensemble of cores within a large region of space. Bayet et al. (2010b) do in fact explore in great detail the influence of changes in temperature and density within a single hot core when subject to various environmental conditions. They find that thermal variations due to temperature structure would flatten out when averaged over a large number of cores.
The presence of an infrared source in the centre or in the vicinity of the region is simulated by subjecting the gas and the dust to an increase in temperature. In our models, we have assumed equal dust and gas temperatures. This is justified since in Phase II the opacity is very high (see values listed in Table 4) and thus the thermal coupling between gas and dust is valid.
In both phases, the chemical network, now improved with D-species, is based on more than 4000 chemical reactions adapted from the UMIST 2006 database (Millar et al., 1997; Le Teuff et al., 2000; Woodall et al., 2007) involving about 300 species (92 of them being deuterated) of which 62 are surface mantle species. The relevant surface reactions included in this work involve mainly simple hydrogenation, apart from the formation of the CHCN which is believed to be formed via reactions involving HCN on grains (Garrod & Herbst, 2006). We note that our approach for surface reactions is rather simplified with respect to that adopted by Garrod & Herbst (2006) and Cazaux et al. (2010). Nevertheless, a qualitative comparison of the mantle composition at the end of Phase I with the results from Cazaux et al. (2010) mantle show good agreement (c.f. HO and HDO abundances on grains as a function of the radiation field for example). As a first approximation and to keep the chemical network to a reasonable size, we include only the mono-D species. DCO and its corresponding ion are the only doubly-D-species in our network.
The deuterium chemical network is similar to that of Roberts & Millar (2000a, b) and Roberts et al. (2004). We assume that D atoms react in the same way as H atoms when reacting with other atoms and molecules.
One of the outputs of the UCL_Chem is the fractional abundance (with respect to the total number of hydrogen nuclei and as a function of time) of gas and surface species from which we calculate the D/H ratios.
The above assumptions (i.e. no doubly-D-species except DCO formed on grain nor triply-D-species included in the code) do affect the reliability of the predictions for a couple of species namely DCO, HD and ND (see Roberts et al. 2003). On the other hand, only mono-D-species have been observed so far in extragalactic environments, and the relative doubly-D/mono-D or triply-D/mono-D species abundance ratios seen in galactic sources are small (see e.g. CDOH/CHDOH = 0.46 and CHDOH/CHDOH = 0.2 from Parise et al. 2004, NDH/NHD = 0.05 from Roueff et al. 2005, and, DCO/HDCO=0.02-0.33 from Roueff & Gerin 2003).
To investigate the deuterium chemistry in high-mass star-forming regions in various extragalactic environments, we have run a grid of 15 chemical models (see Tables 1, 2 and 3). For each type of galaxy (see below), we have studied the deuterium chemistry coming from a moderately dense (cm - index 1 used in tables and figures), a dense (cm - index 2) and a very dense gas (cm - index 3) component. This has allowed us to potentially distinguish which gas component produces the highest fractional abundances of D-molecules, in each type of environment. This also allowed us to determine which gas component is likely to be the best match to the observed extragalactic D/H ratios (see Sect.3).
For reproducing the various types of galaxies, we have assumed the following :
For the Starburst (SB) case, we have used a FUV radiation field of 1000 times the standard value, and a temperature of 500 K, in Phase II, instead of 300 K. These models aim at mimicking environments such as NGC 253 or M 82.
The Cosmic-ray enhanced environments (SB) case is represented in our models by a cosmic ray ionisation rate of 100 times the standard value. This type of models aims at reproducing supernovae and massive star formation as seen in galaxies such as Arp 220 (supernovae rate of 4 year, see e.g. Parra et al. 2007).
The Low metallicity (Low-met) case is represented, in our models, by a decrease in the metallicity to a fifth of the solar value coupled with an increase of the gas-to-dust mass ratio by a factor of 5. Indeed, in the code, quantities such as metallicity, optical depth () and the H formation rate coefficient are coupled (see values presented in Table 4 for each model). In this galaxy type, we have also assumed an increase of FUV radiation field of 1000 times the standard value, and a temperature of 500 K since we are particularly interested in sources such as IC 10, which also host a strong star-formation activity.
Finally, the High redshift (High-z) case is represented by models with a metallicity (and its coupled quantities) of a fifth solar as measured in IC10 (Zaritsky et al., 1994), which is a local dwarf galaxy often used as a good template for archetypical higher-z galaxy populations (e.g. Madden et al. 1997; Leroy et al. 2006; Yin et al. 2010); an increase of FUV radiation field by 1000 times the standard value; a temperature of 500 K and an increase of the cosmic ray ionisation rate by a factor of 100 times the standard value. Indeed, in distant objects (at redshifts greater than 0.1), it is expected to often have a combination of several nuclear activities such as AGN, SB and supernovae (see for instance Seymour 2009 or the studies on the quasar APM08+279 located at z from Weiß et al. 2007 and Riechers et al. 2009). We note that we may indeed underestimate the metallicity at high redshift in our study since optical observations of quasars at z (e.g. Jiang et al. 2007; Juarez et al. 2009) recently revealed that a solar metallicity value is more likely than a sub-solar one for reproducing their observations. However, we emphasize that we have not tried in our present models to match the physical conditions for any high-z source in particular but that we model here what one can consider as an archetypical high redshift environment. Since we know that the early Universe is more populated by dwarf galaxies than any other galaxy type (e.g. White & Frenk 1991), we adopt the values for IC10 as our high redshift standard.
3 Sensitivity of deuterated chemical abundances and D/H abundance ratios to variations in the physical and chemical parameters.
The main aim of this work is to study the sensitivity of deuterated species in models of star-forming regions to variations in physical and chemical parameters that may be characteristics of various galaxy types. A collection of likely detectable tracers as seen by the strength of their chemical abundances will be given in Subsect. 3.3, but first we simply analyse the trends the deuterium chemistry follows as the environment (i.e. galaxy type) and the gas density vary. These trends are summarized in Figs. 1 to 6.
3.1 Deuterated fractional abundances sensitivity
Influence of the gas density
When the gas density increases from n(Hcm to n(Hcm, regardless of the galaxy type, we can divide deuterated species into three categories: those which are insensitive to density changes (e.g. HDO, NHD and HDCS), those that increase with density (e.g. HDS, CD and CHOD) and those which show a decrease in their fractional abundances (e.g. DCN, DNC and CHDOH). The fractional abundances of the D-species belonging to the first category actually tend to converge after yrs, making these species good D-tracers of dense gas (n(Hcm) whatever the environments since they show abundances above 10. On the contrary, we show in Figs. 1 to 3 that especially DCN and DNC decrease with an increase of density (e.g. compare Models NS and NS). This is no surprise since the same behavior is seen for both HCN and HNC when the density increases (Awad et al., 2010). This phenomenon reinforces the conclusions of Bayet et al. (2009) that HCN, HNC and their deuterated counterparts are not particularly good tracers of the very dense gas in galaxies. HDS, CD and CHOD increase with density. The enhancements in their fractional abundance, by more than a factor of 5 (Figs. 1 to 3) are particularly interesting since observing their molecular emissions may allow us to discriminate between the different gas components in a galaxy.
Influence of the FUV radiation field and the temperature (SB case)
Variation in FUV or temperature do not lead to significant changes in abundances for species such as HDO, NHD, DCN, DCO, HDCS, CHDOH and CHOD (i.e. variations in chemical abundances less than a factor of 2-3). However, HDCO and HDS abundances decrease between a factor of 5 and 15 when the FUV radiation field increases. This is due to dissociations by FUV photons. At early times ( yrs), and at low density (n(Hcm) the abundance of CD increases by three orders of magnitude with the increase of both FUV radiation field and temperature. Hence, this species is likely to be a good D-tracer of excitation in galaxies since very sensitive to changes in both temperature and FUV radiation field. At later times, the increase of the CD abundance is less pronounced, but still significant (i.e. about two orders of magnitude). On the other hand, DNC shows the largest decrease i.e. by more than one order of magnitude, when the FUV radiation field and the temperature are increased (linked as seen in Table 2). This may be due to cosmic-ray-induced photodissociation, which has a significant temperature dependence.
Influence of the cosmic ray ionisation rate (SB case)
Increasing the cosmic ray ionisation rate () by a factor of 100 (Models SB to SB - See Table 3 and Fig. 2) leads to most abundances been significantly reduced by several orders of magnitude (e.g. HDCS, NHD, DCO, HDCO, CHDOH, HDO, DNC). As already found by Bayet et al. (2008), by Meijerink & Spaans (2005), and Meijerink et al. (2006), the effects of increasing the cosmic ray ionisation rate are complex, and they generally lead to a chemistry that approaches steady-state more quickly (see Fig.2). The general decrease of abundances seen when increases is probably due to a release of reactive ions such as ionized carbon, reacting quickly with many oxygen-bearing species such as HDCO, DCO, etc. However, some species such as DCO, ND (both at low density only), DCN and DCN (whatever the gas component) show abundances increased by several orders of magnitude with . In the case of nitrogen-bearing species, we suspect that the main reservoir of nitrogen, N, becomes ionized and gives rise to N, a reactive species that promotes the formation of nitrogen-bearing species. This behaviour is also found for the hydrogen counterparts (e.g. Bayet et al. 2008). What is interesting here is that D-species seem to be affected in the same way as their hydrogen counterpart. Species such as ND and DCN are thus likely to be good D-tracers of cosmic ray ionisation-excited environments since they show high fractional abundances and are sensitive to their change. In principle, we thus can use deuteration to constrain although the effects of high cosmic rays ionisation rates on the chemistry are very complex (Bayet et al. 2010, in prep) and highly dependent on initial conditions as well as the available gas coolants (Bayet et al., 2010a).
Influence of the metallicity (Low-met case)
Low metallicity environments such as those described in Sect. 2 (see Models Low-met to Low-met in Table 3), lead to few surprises for most species in that their abundance is simply reduced accordingly (see Fig. 3 for HDCO, DCO, DNC, NHD and HDS). However, some species, such as CD or DCN, show an increase in their abundances with a decrease of the metallicity whereas some other species such as HDCS, CHDOH, CHOD do not seem sensitive to metallicity changes. Therefore these latest species (potentially) are good D-tracers of dense gas (since they have high abundances) whatever the metallicity (especially interesting at high redshift), similarly to their hydrogen counterpart (such as CHOH) already noted by Bayet et al. (2008) and Röllig et al. (2007).
Influence of the above parameters coupled together (High-z case)
In extreme conditions such as those used for mimicking high redshift environments (see Models High-z to High-z in Table 3 and Fig. 3), most of the chemical abundances of D-species decrease: CHOD, HDS disappear completely from Fig. 3 as compared to the normal spiral case (see Fig. 1); NHD, DCO, DNC, and CHDOH have chemical abundances above the limit of detectability only at early time (i.e. before 1.5-3 yrs) then they drop severely at later time; HDCS, HDO, HDCO, even if still above the limit of detectability at all times, decrease by more than two orders of magnitude in high redshift environments as compared to the normal spiral case. These species nonetheless remain the best D-tracers of dense gas components in high redshift environments since very sensitive to the parameters changes. Only DCN and CD show fractional abundance either boosted by up to four orders of magnitude as compared to the Models NS-NS or stay unchanged, respectively. In the case of the high-z, the formation of DCN is dominated by the reaction D + HCN which is inefficient in our models for normal spirals because of the deficiency of atomic deuterium. In the models for normal spirals, the fractional abundance of D is indeed five orders of magnitude lower than in the models for high-z cores.
3.2 D/H abundance ratios sensitivity
As seen in Figs. 4 to 6, the D-fractionation calculated from the chemical abundances of all species do not show similar D/H ratio values, regardless of the environment. In protostellar galactic cores, such as IRAS 16293-2422, such behaviour is already well-known (e.g. van Dishoeck et al. 1995; Ceccarelli et al. 1998; Parise et al. 2004). Authors indeed observed NHD/NH=0.1 whereas they obtained DCO/HCO=0.009. They interpret the higher D/H ratios they obtained as coming from a cooler and more extended gas located in the envelope around the source and not coming from the hot gas in the core. In our study, we can see (Fig. 4) that the highest ratios are indeed obtained, on average, for the galactic cases (i.e. Models NS), for Model NS which has got the smallest density (see Table 5).
In roughly all the models, there are three groups of molecules which converge to similar D/H values: the group including DCO/HCO, ND/NH and HDCO/HCO; the group containing NHD/NH and DCO/HCO and the group containing CHDOH/CHOD, CHDCN/CHCN, CHOD/CHOH. D/H however varies by several orders of magnitude amongst those groups. With the current limited D/H dataset on galactic and extragalactic sources, confirmation of such results is not possible.
The HD/H does not vary as much as the other D/H ratios from a galaxy type to another. On average the largest discrepancy between D/H ratios group estimates is shown to be for the SB case where the spread between all the D/H values is the largest (see Fig.5).
Despite these uncertainties and for providing future observers with D/H values, steady-state D/H ratios for various molecules are listed in Table 5.
3.3 Predictions for observers
In Table 5, we list the D-species that should be likely detectable in external galaxies as based on the values obtained for their chemical abundance. We define the limit of detectability as [n(X)/n] = , as is typical for dense gas molecules in the Milky Way. Of course, detectability will also depend on the excitation conditions as well as the geometry of the source. According to the estimates of Lintott et al. (2005), however these abundances may be sufficiently large that unresolved active galaxies should be likely detectable in these species even at high redshift.
Of particular interest are HDO, DCN and HDCO which are very abundant ([n(X)/n]) in all the investigated environments. These species are ideal molecules to be used for a first observational campaign aiming at validating the model predictions described here.
A second category of interesting species includes DCO, DCN, DNC and ND as they seem to arise only from one gas component in a single galaxy type: SB (for DCO), SB (for DCN), NS (for DNC) and SB (for ND). In a similar vein, CD and HDCS are predicted to be excessively reactive to the presence of cosmic rays ionisation whereas they survive easily in FUV-enhanced environments. HDS is more sensitive to FUV changes than CD and HDCS and it survives to both FUV and cosmic ray ionisation as long as the gas is dense. Here, HDS may thus be a good discriminant of the different gas components in SB and SB environments since it shows high abundances for the densest gas component. However, as soon as the metallicity drops, it becomes undetectable whereas HDCS and CD abundances are still high.
The third category of species contains DCO, CHDOH and CHOD. These molecules do not survive in high-redshift nor in low-metallicity galaxies. The fact that these species survive nonetheless in SB sources may reflect the poor influence of FUV photons on highly embedded material.
Species such as CHDCN and HD are likely undetectable, whatever the conditions since their chemical abundances are low.
Finally, DCN, HDO, HDCO and to a lesser extent DCO, DCO and NHD seem also to be good candidates i.e. abundant enough at various redshifts.
Of course, our predictions are qualitative in nature since based
only on an abundance criterion. For more quantitative predictions,
radiative transfer models, as well as the knowledge of the source
geometry, are needed. DCN, HDCO, DCO, DCO and
NHD expected line intensities can not be derived due to the
lack of collisional rates. For HDO however, we have been able to
derive very rough estimates of line intensities using
4 Discussion and conclusions
Although we do not aim at modelling any specific source, we have
qualitatively compared our models with extragalactic D-species
observations published so far. Mauersberger et al. (1995) obtained upper
limits for the emission of the DCN J=2-1 line toward NGC 253 and
IC 342. More recently, Martín et al. (2006) observed DCO, DCN,
DNC and ND in NGC 253. When we convert these
observations and our predicted fractional abundances into
approximate column densities
Due to the lack of additional observations in extragalactic environments, we are not able to further refine our predictions. Nevertheless, our simple approach offers some significant insights on the major trends of the deuterium chemistry when subject to various physical conditions. Our models can also potentially help discriminate cold and dense star-forming gas between SB, SB, low-metallicity, etc activity. These predictions may be of particular interest when looking at extragalactic very dense high-mass star-forming regions where CO is expected to be depleted onto grains (similarly to what it is seen in our own Galaxy), and where thus D-species drive the chemistry.
EB acknowledges financial support from STFC. ZA would like to thank both the ORS and the Perren studentship schemes for funding. The authors thank the referee for helpful comments on the original version of this paper.
|Initial number density (phase 1)||300 H cm|
|Final number density (phase 1 and 2)||1 cm|
|Temperature (phase 1)||10 K|
|Temperature (phase 2)||300 K|
|External UV radiation intensity||1 Habing|
|Cosmic ray ionization rate||1.3 s|
|Visual extinction||580.6 mag|
|H formation rate coefficient||1.0 cms|
|Metallicity||z||solar values, see Table 3|
: The collapse is treated as a ‘modified free-fall’, as defined by Rawlings et al. (1992), and the parameter is introduced to allow for collapse in free-fall 1.0, or somewhat more slowly (1.0) if gas or magnetic pressure resists the collapse.
|(cm)||(pc)||(z)||mass ratio||Abund. ratios||()||(K)||()||parameters|
: z = 1 corresponds to solar values of the elemental abundances ratios; : expressed in units of = 1.3 s.
|ratio of the number|
|densities of grains||1.010||1.010||1.010||1.010||1.010||1.010||2.010||2.010||2.010|
|to hydrogen nuclei|
: The H formation rate coefficient is computed at the beginning of the phase 2 of the UCL hot core model and it is expressed in cms.
|Normal Spiral||Starburst||Cosmic- Rays||Low||High|
- affiliation: Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, UK.
- affiliation: Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, UK.
- affiliation: Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, UK.
- See http://www.sron.rug.nl/ vdtak/radex/radex.php
- See http://www.strw.leidenuniv.nl/ moldata/
- See primary beam size values in http://www.eso.org/sci/facilities/alma/observing/specifications/
- We performed such conversion by simply multiplying our predicted fractional abundance by the column density of hydrogen at the relevant visual extinction.
- Aladro, R., Martin-Pintado, J., Martin, S., Mauersberger, R., & Bayet, E. 2010, ArXiv e-prints 2010arXiv1009.1831A
- Awad, Z., Bayet, E., & Viti, S. 2010, MNRAS, in prep.
- Bacmann, A., Lefloch, B., Ceccarelli, C., Steinacker, J., Castets, A., & Loinard, L. 2003, ApJ, 585, L55
- Bayet, E., Aladro, R., Martín, S., Viti, S., & Martín-Pintado, J. 2009, ApJ, 707, 126
- Bayet, E., Hartquist, T. W., Viti, S., Williams, D. A., & Bell, T. A. 2010a, ArXiv e-prints 2010arXiv1007.1533B
- Bayet, E., Viti, S., Williams, D. A., & Rawlings, J. M. C. 2008, ApJ, 676, 978
- Bayet, E., Yates, J., & Viti, S. 2010b, MNRAS, submitted
- Carilli, C. L., Kohno, K., Kawabe, R., Ohta, K., Henkel, C., Menten, K. M., Yun, M. S., Petric, A., & Tutui, Y. 2002, AJ, 123, 1838
- Caselli, P., van der Tak, F. F. S., Ceccarelli, C., & Bacmann, A. 2003, A&A, 403, L37
- Cazaux, S., Cobut, V., Marseille, M., Spaans, M., & Caselli, P. 2010, ArXiv e-prints 2010arXiv1007.1061C
- Ceccarelli, C., Castets, A., Loinard, L., Caux, E., & Tielens, A. G. G. M. 1998, A&A, 338, L43
- Garrod, R. T. & Herbst, E. 2006, A&A, 457, 927
- Gerin, M., Combes, F., Wlodarczak, G., Jacq, T., Guelin, M., Encrenaz, P., & Laurent, C. 1992, A&A, 259, L35
- Israel, F. P. & Baas, F. 2003, A&A, 404, 495
- Jacq, T., Walmsley, C. M., Henkel, C., Baudry, A., Mauersberger, R., & Jewell, P. R. 1990, A&A, 228, 447
- Jacq, T., Walmsley, C. M., Mauersberger, R., Anderson, T., Herbst, E., & De Lucia, F. C. 1993, A&A, 271, 276
- Jiang, L., Fan, X., Vestergaard, M., Kurk, J. D., Walter, F., Kelly, B. C., & Strauss, M. A. 2007, AJ, 134, 1150
- Juarez, Y., Maiolino, R., Mujica, R., Pedani, M., Marinoni, S., Nagao, T., Marconi, A., & Oliva, E. 2009, A&A, 494, L25
- Knauth, D. C., Andersson, B.-G., McCandliss, S. R., & Moos, H. W. 2003, ApJ, 596, L51
- Le Teuff, Y. H., Millar, T. J., & Markwick, A. J. 2000, A&AS, 146, 157
- Leroy, A., Bolatto, A., Walter, F., & Blitz, L. 2006, ApJ, 643, 825
- Lintott, C. J., Viti, S., Williams, D. A., Rawlings, J. M. C., & Ferreras, I. 2005, MNRAS, 360, 1527
- Loren, R. B. & Wootten, A. 1985, ApJ, 299, 947
- Madden, S. C., Poglitsch, A., Geis, N., Stacey, G. J., & Townes, C. H. 1997, ApJ, 483, 200
- Martín, S., Mauersberger, R., Martín-Pintado, J., Henkel, C., & García-Burillo, S. 2006, ApJS, 164, 450
- Mauersberger, R., Henkel, C., & Chin, Y. 1995, A&A, 294, 23
- Mauersberger, R., Henkel, C., Jacq, T., & Walmsley, C. M. 1988, A&A, 194, L1
- Meijerink, R. & Spaans, M. 2005, A&A, 436, 397
- Meijerink, R., Spaans, M., & Israel, F. P. 2006, ApJ, 650, L103
- Meyer, D. M., Jura, M., & Cardelli, J. A. 1998, ApJ, 493, 222
- Miettinen, O., Harju, J., Haikala, L. K., Kainulainen, J., & Johansson, L. E. B. 2009, A&A, 500, 845
- Millar, T. J., MacDonald, G. H., & Gibb, A. G. 1997, A&A, 325, 1163
- Oliveira, C. M. 2003, in Bulletin of the American Astronomical Society, Vol. 35, Bulletin of the American Astronomical Society, 1317–+
- Parise, B., Castets, A., Herbst, E., Caux, E., Ceccarelli, C., Mukhopadhyay, I., & Tielens, A. G. G. M. 2004, A&A, 416, 159
- Parise, B., Ceccarelli, C., Tielens, A. G. G. M., Herbst, E., Lefloch, B., Caux, E., Castets, A., Mukhopadhyay, I., Pagani, L., & Loinard, L. 2002, A&A, 393, L49
- Parra, R., Conway, J. E., Diamond, P. J., Thrall, H., Lonsdale, C. J., Lonsdale, C. J., & Smith, H. E. 2007, ApJ, 659, 314
- Rawlings, J. M. C., Hartquist, T. W., Menten, K. M., & Williams, D. A. 1992, MNRAS, 255, 471
- Riechers, D. A., Walter, F., Carilli, C. L., & Lewis, G. F. 2009, ApJ, 690, 463
- Roberts, H., Herbst, E., & Millar, T. J. 2003, ApJ, 591, L41
- —. 2004, A&A, 424, 905
- Roberts, H. & Millar, T. J. 2000a, A&A, 364, 780
- —. 2000b, A&A, 361, 388
- Röllig, M., Abel, N. P., Bell, T., Bensch, F., Black, J., Ferland, G. J., Jonkheid, B., Kamp, I., Kaufman, M. J., Le Bourlot, J., Le Petit, F., & etc, . 2007, A&A, 467, 187
- Roueff, E. & Gerin, M. 2003, Space Science Reviews, 106, 61
- Roueff, E., Lis, D. C., van der Tak, F. F. S., Gerin, M., & Goldsmith, P. F. 2005, A&A, 438, 585
- Roueff, E., Tiné, S., Coudert, L. H., Pineau des Forêts, G., Falgarone, E., & Gerin, M. 2000, A&A, 354, L63
- Schilke, P., Walmsley, C. M., Pineau Des Forets, G., Roueff, E., Flower, D. R., & Guilloteau, S. 1992, A&A, 256, 595
- Sembach, K. R. & Savage, B. D. 1996, ApJ, 457, 211
- Seymour, N. 2009, in Panoramic Radio Astronomy: Wide-field 1-2 GHz Research on Galaxy Evolution
- Snow, T. P., Rachford, B. L., & Figoski, L. 2002, ApJ, 573, 662
- Sofia, U. J., Cardelli, J. A., Guerin, K. P., & Meyer, D. M. 1997, ApJ, 482, L105+
- Stark, R., van der Tak, F. F. S., & van Dishoeck, E. F. 1999, ApJ, 521, L67
- Tielens, A. G. G. M. 1983, A&A, 119, 177
- van der Tak, F. F. S., Black, J. H., Schöier, F. L., Jansen, D. J., & van Dishoeck, E. F. 2007, A&A, 468, 627
- van Dishoeck, E. F., Blake, G. A., Jansen, D. J., & Groesbeck, T. D. 1995, ApJ, 447, 760
- Vastel, C., Phillips, T. G., & Yoshida, H. 2004, ApJ, 606, L127
- Viti, S., Collings, M. P., Dever, J. W., McCoustra, M. R. S., & Williams, D. A. 2004, MNRAS, 354, 1141
- Viti, S. & Williams, D. A. 1999, MNRAS, 305, 755
- Vrtilek, J. M., Gottlieb, C. A., Langer, W. D., Thaddeus, P., & Wilson, R. W. 1985, ApJ, 296, L35
- Walmsley, C. M., Hermsen, W., Henkel, C., Mauersberger, R., & Wilson, T. L. 1987, A&A, 172, 311
- Watson, W. D. 1980, in Les Spectres des Molécules Simples au Laboratoire et en Astrophysique, 526–544
- Weiß, A., Downes, D., Neri, R., Walter, F., Henkel, C., Wilner, D. J., Wagg, J., & Wiklind, T. 2007, A&A, 467, 955
- White, S. D. M. & Frenk, C. S. 1991, ApJ, 379, 52
- Woodall, J., Agúndez, M., Markwick-Kemper, A. J., & Millar, T. J. 2007, A&A, 466, 1197
- Yin, J., Magrini, L., Matteucci, F., Lanfranchi, G. A., Gonçalves, D. R., & Costa, R. D. D. 2010, ArXiv e-prints 2010arXiv1005.3500Y
- Zaritsky, D., Kennicutt, R. C., & Huchra, J. P. 1994, ApJ, 420, 87