## Observations of nitrogen isotope fractionation in deeply embedded protostars

S. F. Wampfler 1 2 , , J. K. Jørgensen 2 1 , , M. Bizzarro 1 , and S. E. Bisschop 1 2 ,

- 1 Centre for Star and Planet Formation, Natural History Museum of Denmark, University of Copenhagen, Øster Voldgade 5-7, DK-1350 København K, Denmark
- 2 Niels Bohr Institute, University of Copenhagen, Juliane Maries Vej 30, DK-2100 København Ø, Denmark

Received 7 March 2014 / Accepted 26 July 2014

## ABSTRACT

Context. The terrestrial planets, comets, and meteorites are significantly enriched in 15 N compared to the Sun and Jupiter. While the solar and jovian nitrogen isotope ratio is believed to represent the composition of the protosolar nebula, a still unidentified process has caused 15 N-enrichment in the solids. Several mechanisms have been proposed to explain the variations, including chemical fractionation. However, observational results that constrain the fractionation models are scarce. While there is evidence of 15 Nenrichment in prestellar cores, it is unclear how the signature evolves into the protostellar phases.

Aims.

The aim of this study is to measure the

14

N 15

/

N ratio around three nearby, embedded low- to intermediate-mass protostars.

Methods. Isotopologues of HCN and HNC were used to probe the 14 N 15 / N ratio. A selection of J = 3-2 and 4-3 transitions of H 13 CN, HC 15 N, HN 13 C, and H 15 NC was observed with the Atacama Pathfinder EXperiment telescope (APEX). The 14 N 15 / N ratios were derived from the integrated intensities assuming a standard 12 C 13 / C ratio. The assumption of optically thin emission was verified using radiative transfer modeling and hyperfine structure fitting.

Results. Two sources, IRAS 16293A and R CrA IRS7B, show 15 N-enrichment by a factor of ∼ 1 5 . -2 5 in both HCN and HNC . with respect to the solar composition. IRAS 16293A falls in the range of typical prestellar core values. Solar composition cannot be excluded for the third source, OMC-3 MMS6. Furthermore, there are indications of a trend toward increasing 14 N 15 / N ratios with increasing outer envelope temperature.

Conclusions. The enhanced 15 Nabundances in HCN and HNC found in two Class 0 sources ( 14 N 15 N / ∼ 160 -290) and the tentative trend toward a temperature-dependent 14 N 15 / N ratio are consistent with the chemical fractionation scenario, but 14 N 15 / N ratios from additional tracers are indispensable for testing the models. Spatially resolved observations are needed to distinguish between chemical fractionation and isotope-selective photochemistry.

Key words. Astrochemistry - ISM: abundances - ISM: molecules - Stars: formation - Stars: individual: IRAS 16293-2422, R CrA IRS7B, OMC-3 MMS6

## 1. Introduction

The largest isotopic anomalies in the solar system are found in highly volatile elements like hydrogen, carbon, oxygen, and nitrogen (e.g., Birck 2004; Anders &amp; Grevesse 1989). Thereof, the most extreme variations among solar system objects exist in the deuterium-to-hydrogen ratio, D H, and the ratio of the two / stable nitrogen isotopes, 14 N 15 N / /star . Recent results from the Genesis mission (Marty et al. 2011), which sampled the solar wind, have demonstrated that there is a large discrepancy in the nitrogen isotopic composition between the solar system solids and the gas reservoir: The 14 N 15 N ratios of the terrestrial planets, / comets (e.g., Hutsemékers et al. 2005; McKeegan et al. 2006; Bockelée-Morvan et al. 2008; Manfroid et al. 2009; Jehin et a l. 2009; Mumma &amp; Charnley 2011; Shinnaka et al. 2014), meteorites (e.g., Kerridge 1993; Briani et al. 2009; Bonal et al. 2010), and interplanetary dust particles (IDPs, Messenger 2000; Floss et al. 2006) are significantly enhanced compared to the Sun ( 14 N 15 / N = 441 ± 6, Marty et al. 2011) and Jupiter (435 ± 57,

Owen et al. 2001; Fouchet et al. 2004; Abbas et al. 2004). The derived ratios for the Earth's atmosphere (272, e.g., Junk &amp; Svec 1958; Mariotti 1983), Venus (Ho ff man et al. 1979), and the interior of Mars (Mathew &amp; Marti 2001) are very similar and also lie within the range spanned by most primitive meteorites (e.g., Fig. 2 from Meibom et al. 2007). The nitrogen isotopic composition of the Sun and Jupiter, on the other hand, is substantially 15 N-poorer (Marty et al. 2011; Hashizume et al. 2000) and believed to be representative of the presolar nebula. This idea is supported by the results from an isotopic measurement of a high-temperature condensate - titanium nitride (TiN), named osbornite - in a calcium-aluminum-rich inclusion (CAI) of the 'Isheyevo' meteorite. CAIs are believed to be the first solids that condensated from the gas of the presolar nebula. The 14 N 15 N / ratio measured from this osbornite (424 ± 3, Meibom et al. 2007) is therefore interpreted as reflecting the nitrogen isotope composition of the solar nebula at the time of CAI formation 4567 3 . ± 0 16 Myr ago (Connelly et al. 2012). .

/star In cosmochemistry, nitrogen isotope ratios are generally reported in per mil deviations from the terrestrial atmospheric standard, denoted δ 15 N = [ ( 15 N 14 N / ) sample / ( 15 N 14 N / ) air -1 ] × 10 3 where 15 N 14 Nair / = 3 676 . × 10 -3 .

The aim of this paper is to study the 14 N 15 N ratio around / low-mass protostars to investigate whether fractionation of nitrogen isotopes occurs in the early phases of star formation and to better understand what mechanisms could be responsible for the solar system heterogeneity.

## 1.1. 15 N-enrichments in primitive solar system matter

Meteorites, interstellar dust particles, and comets represent the most pristine solar system material. For meteorites, the nitrogen isotope composition is highly variable among the different classification groups, among meteorites belonging to the same group, and even within individual meteorites (e.g., Alexander et al. 2007; Bonal et al. 2010). The heaviest nitrogen isotope signature in any solar system material was measured to date in the 'Isheyevo' chondrite ( 14 N 15 N / ∼ 50, Briani et al. 2009; Bonal et al. 2010). However, the nitrogen isotope composition in Isheyevo is also extremely variable, ranging from an almost solar-like ( 14 N 15 N / = 394 ± 12) to a very 15 N-heavy isotopic composition. The 15 N-enrichment consists of two components, a di ff use background enrichment ( 14 N 15 N / ∼ 160) and localized micron-sized zones with extremely high 15 N-anomalies, so-called hotspots. Generally, the highest 15 N-enrichments seem to be found in the most pristine meteoritic components (Busemann et al. 2006; Alexander et al. 2007; Briani et al. 2009; Bonal et al. 2010). The carrier of the signature is most likely carbonaceous (Bonal et al. 2010), but the exact functional groups carrying the 15 N-enrichment have not yet been identified.

Interplanetary dust particles (IDPs) are collected in the stratosphere of the Earth and believed to be of cometary or potentially asteroidal origin. They show similar 15 Ncharacteristics as chondrites, with a modest di ff use enrichment and hotspots, but with generally more moderate 15 Nenhancements (Busemann et al. 2006; Floss et al. 2006). The 'Stardust' samples from comet 81P Wild-2 have bulk / nitrogen isotope compositions that cluster around the terrestrial composition, but moderate enhancements are also measured (McKeegan et al. 2006). Hotspot values on submicrometer scale are similar to IDPs. An isotopically light component, consistent with the solar value, is also found in the Wild-2 particles, which is not the case for IDPs.

The average 14 N 15 N ratio from CN spectroscopy in a sam-/ ple of comets is 147 8 . ± 5 6 (Manfroid et al. 2009) with slight . di ff erences for Jupiter family ( 14 N 15 N / = 156 8 . ± 12 2, probable . Kuiper belt origin) and Oort cloud ( 14 N 15 N / = 144 ± 6 5) comets. . These results are consistent with the 14 N 15 N ratio from HCN / measurements by Bockelée-Morvan et al. (2008) in the comets 'Hale-Bopp' and 'Holmes' (see also the review by Jehin et al. 2009) as well as measurements from NH2 in comet C 2012 S1 / (ISON, 14 N / 15 N = 139 ± 38, Shinnaka et al. 2014) and the average spectrum of 12 comets ( ∼ 130, Rousselot et al. 2014).

An overview of 14 N 15 / N ratios in solar system bodies from the literature is provided in Table B.1 in the appendix.

## 1.2. Proposed scenarios: chemical fractionation vs. isotope selective photodissociation

A variety of processes have been proposed to explain the observed nitrogen isotope heterogeneity, including Galactic chemical evolution (GCE, e.g., Adande &amp; Ziurys 2012, and references therein), isotope-selective N2 photodissociation (Clayton 2002; Lyons et al. 2009; Heays et al. 2014), and chemical fractionation (Terzieva &amp; Herbst 2000; Rodgers &amp; Charnley 2008a,b; Wirström et al. 2012; Hily-Blant et al. 2013a).

The term 'Galactic chemical evolution' refers to a scenario in which the observed 14 N 15 N ratio variations are the result of / an enrichment of the interstellar medium (ISM) with the products of stellar nucleosynthesis over time. GCE models predict the 14 N 15 N ratio to vary with distance from the Galactic center / as a result of di ff erent production mechanisms for 14 N and 15 N, the details of which are not yet fully understood (see discussion and references in Adande &amp; Ziurys 2012). Observational support comes from surveys of the nitrogen isotopic composition of the Galaxy (by, e.g., Linke et al. 1977; Wannier et al. 1981; Guesten &amp; Ungerechts 1985; Dahmen et al. 1995; Wilson 1999; Adande &amp; Ziurys 2012). Most of these studies found a radial gradient in the 14 N 15 N ratio of the Galaxy, arguing for Galac-/ tic evolution as the main cause for deviations in the isotopic composition of di ff erent star-forming regions. However, significant source-to-source variations within individual clouds (e.g., Ikeda et al. 2002) as well as evidence for spatial variations in prestellar cores (Hily-Blant et al. 2013a) hint at local processes being more important for the nitrogen isotope composition of the protosolar disk than Galactic inheritance.

Liang et al. (2007) demonstrated that isotopically selective photodissociation of N2 plays an important role in the explanation of the nitrogen isotope anomaly inferred in the upper atmosphere of Titan from HCN observations (Niemann et al. 2005). The radiation field component corresponding to the dissociation energy of 14 N2 cannot penetrate as deep into the atmosphere as the one for 14 N 15 N, because of 14 N2 being more abundant than 14 N 15 N, an e ff ect which is referred to as photochemical selfshielding. To what extent N2 self-shielding could have a ff ected the nitrogen isotope composition of the protosolar nebula is not yet clear, but is currently being investigated (e.g., Lyons et al. 2009; Bertin et al. 2013). Li et al. (2013) and Heays et al. (2014) have recently presented photodissociation rates for 14 N2 and 14 N 15 N. In addition, Heays et al. (2014) have incorporated these rates into chemical models with conditions simulating an interstellar cloud and a protoplanetary disk with grain growth. Their results show 15 N-enhancements by up to a factor of 10 in HCN in the protoplanetary disk model. Based on these results they argue that the observed nitrogen fractionation in protostellar environments could be caused by a combination of chemical and photolytic reactions.

As for chemical fractionation, the astrochemical gas-grain models by Rodgers &amp; Charnley (2008a) predict a considerable 15 N-enhancement for di ff erent nitrogen-bearing molecules in the gas phase of dense cloud cores through isotope exchange reactions in low-temperature chemistry. Their results indicate that the fractionation mechanisms for nitriles, i.e., molecules bearing the nitrile functional group (-CN), and nitrogen hydrides (-NH) are di ff erent. The slow route ( ∼ 10 6 yrs) through the dissociative ionization of N2, leading to N + , is responsible for the fractionation in nitrogen hydrides, including ammonia (NH3) and diazenylium (N2H + ). The more rapid route ( ∼ 10 5 yrs) starts from atomic nitrogen and drives the 15 N-enrichment of (iso-)nitriles and associated species such as cyanide (CN), hydrogen cyanide (HCN), and hydrogen isocyanide (HNC). The largest 15 N-enhancement in this chemical fractionation model is predicted for nitrile species. With substantial 15 N-enrichments predicted by both chemical fractionation and photodissociation models, species like HCN should therefore be particularly good tracers of a potential isotope fractionation in nitrogen. Refined chemical fractionation models by Wirström et al. (2012) include the dependence of the ion-molecule reactions on the H2 spinstate and show that a large 15 N-enhancement is not necessarily expected to correlate with an extreme D-enrichment, depending on the carriers of the two signatures. The lack of a correlation between the 15 N and D anomalies measured in meteorites had previously often been invoked as an argument against an inheritance of the isotopic signature measured in primitive solar system matter from local cloud or disk processes in early phases of solar system formation (e.g., Aléon 2010).

Table 1. Literature values of the 14 N 15 / N ratio from low-mass prestellar cores and protostars (top). For comparison, the most nearby sources from the galactic survey by Adande &amp; Ziurys (2012) are listed as well (bottom).

| Source                                                                 | HCN                                                          | HNC                       | CN                        | NH 3                        | NH 2 D                                                                         | N 2 H +                    | Remarks                                                                         |
|------------------------------------------------------------------------|--------------------------------------------------------------|---------------------------|---------------------------|-----------------------------|--------------------------------------------------------------------------------|----------------------------|---------------------------------------------------------------------------------|
| L1544 L1498 L183 L1521E L1521F                                         | 140 - 360 a 69 - 154 e > 813 f 140 - 250 a 151 ± 16 f > 51 e | > 27 e > 90 e 24 - 31 e   | 510 ± 70 b 476 ± 70 b     |                             | > 700 c                                                                        | 1000 ± 200 d               | prestellar core prestellar core prestellar core prestellar core prestellar core |
| Barnard 1 NGC1333 IRAS4 NGC1333 DCO + NGC1333 L134N(S) L1689N Cha-MMS1 | 330 + 60 - 50 g                                              | 225 + 75 - 45 g 135 i     | 290 + 160 - 80 g          | 300 + 55 - 40 g 344 ± 173 h | 230 + 105 - 55 g > 270 c 360 + 260 - 110 c 530 + 570 - 180 c 810 + 600 - 250 c | 400 + 100 - 65 g , > 600 g | protostar                                                                       |
| IRAS 17233-3606 DR 21 Orion KL Orion Bar                               | 276 ± 62                                                     | 159 ± 31 159 ± 40 338 311 | 243 ± 66 234 ± 47 361 141 |                             |                                                                                |                            | HII region HII region PDR                                                       |

±

±

References. a Hily-Blant et al. (2013a) - the range reflects spatial variations within the core, b Hily-Blant et al. (2013b), c Gerin et al. (2009) based on 15 NH2D observations at the IRAM 30m, d Bizzocchi et al. (2013), revised from the previously published value of 446 ± 71 (Bizzocchi et al. 2010), e Milam &amp; Charnley (2012), f Ikeda et al. (2002), g Daniel et al. (2013), h Lis et al. (2010) based on observations of 15 NH3 inversion lines with the Green Bank Telescope, i reported in Milam &amp; Charnley (2012) based on Tennekes et al. (2006), j Adande &amp; Ziurys (2012)

Table 2. Coordinates and properties of the observed sources.

| Source      | RA [h ms]        | Dec [ ◦ ′ ′′ ]   |   v lsr [km s - 1 ] | Class                     |   d [pc] |   L bol [L /circledot ] |   M env [M /circledot ] |
|-------------|------------------|------------------|---------------------|---------------------------|----------|-------------------------|-------------------------|
| OMC-3       | MMS605:35:23.5 a | - 05:01:32.0 a   |                10   | intermed.-mass (Class 0)  |      450 |                    60   |                   36    |
| Elias 29    | 16:27:09.4 b     | - 24:37:18.6 b   |                 5   | Class I                   |      125 |                    14.1 |                    0.04 |
| IRAS 16293A | 16:32:22.87 c    | - 24:28:36.4 c   |                 3   | Class 0                   |      120 |                    27   |                    5.4  |
| R CrA IRS7B | 19:01:56.4 d     | 36:57:28.3 d     |                 5.6 | transitionary Class 0 / I |      130 |                     4.7 |                    2.2  |

-

References. a Johnstone et al. (2003), b Jørgensen et al. (2008), c Jørgensen et al. (2011), d Lindberg &amp; Jørgensen (2012)

Observations of the 14 N 15 N ratio in individual low-mass pre-/ and protostellar sources in the literature are scarce and have so far led to inconclusive results about whether or not fractionation of nitrogen isotopes occurs in pre- and protostellar environments. Moreover, it still an open question whether chemical fractionation or isotope-selective photodissociation could explain the observed solar system variations. Results from HCN observations of three dense cores by Ikeda et al. (2002) show large variations between the sources and the result for L1521E ( 14 N 15 N / = 151 ± 16) is consistent with the chemical fractionation scenario. Gerin et al. (2009) on the other hand derive 14 N 15 N ratios of / ∼ 360 -810 from 15 NH2D in several interstellar clouds that are consistent with the solar composition and thus argue against nitrogen isotope fractionation. Follow-up observations by Lis et al. (2010) using measurements of 15 NH3 in Barnard 1 and NGC 1333 yield ratios that fall in between the solar and terrestrial values. Bizzocchi et al. (2010) reported a 14 N 15 N ratio of 446 / ± 71 from N 15 NH + observations toward the prestellar core L1544, which was recently revised to be as high as 1000 ± 200 based on new observations of 15 NNH + (Bizzocchi et al. 2013), implying a depletion of 15 N in N2H . + Results from mapping of HCN isotopologues in the two prestel- lar cores L183 and L1544 by Hily-Blant et al. (2013a) show an overall enhancement of 15 N as well as some spatial variations within each source, strongly favoring a fractionation scenario. For L1544, Hily-Blant et al. (2013a) find an enhancement of 15 N in HCN (140 &lt; 14 N 15 N / &lt; 360), but there are no indications for isotope fractionation in CN ( 14 N 15 N / = 500 ± 75, Hily-Blant et al. 2013b). These di ff erences likely arise from the variety of molecular species and isotopologues used to probe the nitrogen isotope ratio in di ff erent sources. For example, the low 14 N 15 N ratios inferred from nitriles like HCN seem to sup-/ port the chemical fractionation scenario, while 14 N 15 N ratios / derived from nitrogen hydrides, like NH3 and N2H , show no + evidence for nitrogen isotope fractionation or even hint at a 15 Ndepletion. Wirström et al. (2012) and Hily-Blant et al. (2013b) proposed refined chemical fractionation scenarios with di ff erent fractionation paths that may be able to explain the discrepancies between the nitrogen isotope composition seen in individual molecular tracers. To constrain the fractionation models, it is thus important to measure isotope ratios from an entire set of species. However, the prestellar core L1544 and the very young Class 0 binary B1b are so far the only sources for which comprehensive observational data sets of 14 N 15 N ratios from both /

nitriles and nitrogen hydrides exist (cf. Table 1). While evidence for a 15 N-enrichment in some molecules is found for L1544 (Hily-Blant et al. 2013a), the results on B1b by Daniel et al. (2013) show no evidence for nitrogen isotope fractionation in any class of molecules. Whether this di ff erence is caused by the more evolved status of B1b, a source-to-source variation, or any other e ff ect, is unclear and needs to be investigated on a larger source sample. A similar data set is not yet available for other sources, in particular not for more evolved evolutionary stages.

In contrast to previous studies, which mainly focused on the prestellar core phase, we present in this paper measurements of the nitrogen isotope composition in the gas around three Class 0 YSOs from HCN and HNC isotopologues. The aim of this work is to investigate whether signatures of nitrogen isotope fractionation also occur in the embedded protostellar stages. The paper is organized as follows: Section 2 describes the observations and the data reduction procedure. The results are presented in Sect. 3. They are compared to previous results and their implications are discussed in Sect. 4. Conclusions are presented in Sect. 5.

## 2. Observations and data reduction

## 2.1. Observations

Observations of HCN and HNC isotopologues were carried out using the Atacama Pathfinder EXperiment telescope 1 (APEX, Güsten et al. 2006) in service mode. Wobbler switching with an amplitude of 150 ′′ was used throughout. A detailed overview on the observations is provided in Table A.1 in the appendix.

1 APEX is a collaboration between the Max-Planck-Institut für Radioastronomie, the European Southern Observatory, and the Onsala Space Observatory.

Table 3. Molecular data of the HCN and HNC isotopologues observed with APEX. HCN isotopologue data are from the CDMS catalog (Müller et al. 2001, 2005) and HNC isotopologues are from the JPL catalogue (Pickett et al. 1998), all obtained through www.splatalogue.net.

| Species     | Transition   | Frequency [GHz]   | E up [K]   | log( A ul ) [s - 1 ]   |
|-------------|--------------|-------------------|------------|------------------------|
|             |              | APEX-1            |            |                        |
| D 13 CN     | J = 3-2      | 213.51987         | 20.49      | -3.36216               |
| DCN         | J = 3-2      | 217.23854         | 20.85      | -3.33964               |
| H 13 C 15 N | J = 3-2      | 251.17523         | 24.11      | -3.15217               |
| HC 15 N     | J = 3-2      | 258.15700         | 24.78      | -3.11641               |
| H 13 CN     | J = 3-2      | 259.01180 /star   | 24.86      | -3.11216               |
| HN 13 C     | J = 3-2      | 261.26331         | 25.08      | -3.18837               |
| HCN         | J = 3-2      | 265.88643         | 25.52      | -3.07662               |
| H 15 NC     | J = 3-2      | 266.58780         | 25.59      | -3.16196               |
| HNC         | J = 3-2      | 271.98114         | 26.11      | -2.94523               |
| DCN         | J = 4-3      | 289.64492         | 34.75      | -2.94907               |
|             |              | APEX-2            |            |                        |
| HC 15 N     | J = 4-3      | 344.20011         | 41.30      | -2.72583               |
| H 13 CN     | J = 4-3      | 345.33977 /star   | 41.43      | -2.72161               |
| HN 13 C     | J = 4-3      | 348.34027         | 41.80      | -2.79774               |
| DCN         | J = 5-4      | 362.04575         | 52.13      | -2.64860               |
| HNC         | J = 4-3      | 362.63030         | 43.51      | -2.55463               |

Notes. ( /star ) Hyperfine structure only partially resolved.

The source sample consists of three Class 0 YSOs, IRAS 16293A, R CrA IRS7B, and OMC-3 MMS6, and the Class I protostar Elias 29. The coordinates and properties of the sources can be found in Table 2. The targets were selected to cover di ff erent source properties: IRAS 16293-2422 is an extremely well-studied (e.g., Blake et al. 1994; van Dishoeck et al. 1995; Schöier et al. 2002; Takakuwa et al. 2007; Caux et al. 2011; Jørgensen et al. 2011; Pineda et al. 2012), isolated pr otobinary in ρ Oph with a hot corino and known deuterium fractionation (e.g., Loinard et al. 2000; Persson et al. 2013). Our observations were pointed at the A component, IRAS 16293A (Jørgensen et al. 2011). R CrA IRS7B is located in the vicinit y of the Herbig star R CrA and therefore exposed to a strong external radiation field (Lindberg &amp; Jørgensen 2012). OMC-3 MMS6 is an intermediate-mass protostar (Johnstone et al. 2003) in the vigorous environment of the Orion star-forming complex, also exposed to enhanced irradiation (Jørgensen et al. 2006). El ias 29 (e.g., Boogert et al. 2002) is one of the most luminous protostars in the ρ Oph cloud and the most evolved source in our sample.

IRAS 16293A was observed in H 13 CN ( J = 3-2 and 4-3), HC 15 N (3-2 and 4-3), H 13 C 15 N(3-2), DCN(3-2, 4-3, and 5-4), and D 13 CN(3-2) on 2012-04-01 and between 2012-06-04 and 2012-06-13 (project 089.F-9316A, PI J. Jørgensen). The observations of HN 13 C(3-2) and H 15 NC(3-2) were carried out on 2013-08-19 20 (project 092.F-9324A, PI S. Wampfler). / OMC3 MMS6, Elias 29, and R CrA IRS7B were observed in HCN(32), H 13 CN (3-2 and 4-3), HC 15 N(3-2), HNC(3-2), HN 13 C (3-2 and 4-3), and H 15 NC(3-2) with APEX between 2012-08-16 and 2012-10-01 (project 090.F-9317A, PI S. Wampfler). The frequencies of the transitions and the molecular data are provided in Table 3.

The Swedish Heterodyne Facility Instrument (SHeFI, Vassilev et al. 2008) single-sideband SIS receivers APEX-1 (211 -275 GHz) and APEX-2 (267 -370 GHz) were used in combination with the eXtended bandwidth Fast Fourier Transform Spectrometer (XFFTS). Two XFFTS units with an instantaneous bandwidth of 2.5 GHz each are available, which partially overlap and o ff er full coverage of the 4 GHz IF bandwidth of the SHeFI receivers. The large bandwidth of this backend type permits several isotopologues to be observed simultaneously, thereby reducing the calibration uncertainty from the absolute value down to the relative calibration uncertainty across the band. For IRAS 16293A, H 13 CN (3-2) and HC 15 N (3-2) were observed in a combined setting, as were H 13 CN (4-3) and HC 15 N (4-3). In the case of OMC-3 MMS6, Elias 29, and R CrA-IRS7B, the observations of H 13 CN (3-2), HC 15 N (3-2), and HN 13 C (3-2) were carried out simultaneously. Moreover, H 15 NC (3-2) and HCN(3-2) were observed together. H 13 CN (43) and HN 13 C (4-3) were also jointly measured.

## 2.2. Data reduction

The data reduction was carried out using the GILDAS-CLASS software package 2 . The data were calibrated to T mb scale using the default forward e ffi ciency of 0.95 (Dumke &amp; Mac-Auli ff e 2010) for APEX, and beam e ffi ciencies 3 of 0.75 for APEX-1 and 0.73 for APEX-2, respectively. The half-power beam width of APEX is 27 ′′ at 230 GHz and 18 ′′ at 345 GHz.

Several observations were split between di ff erent observing dates (cf. Table A.1). Before averaging those data sets, we checked for deviations in the intensity and line profile on strong

2 http://www.iram.fr/IRAMFR/GILDAS/

3 values from http://www.apex-telescope.org/telescope/efficiency

Fig. 1. Spectra of the HCN and HNC isotopologues observed with APEX toward IRAS 16293A. Where needed, the spectra have been scaled by the factors indicated in the individual panels.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0285v1-figure-31bf8247-9893-428e-8a0f-2a16372ff896.jpg)

lines present in the spectrum, and, if the signal-to-noise ratio for the individual observations permits, also on our target lines. The H 15 NC(3-2) observations from R CrA IRS7B carried out on the 09-30-2012 show substantially larger deviations than the calibration uncertainty ( ∼ 30%) when compared to other days and this dataset was consequently discarded. Deviations in the line profiles between di ff erent days also occur for the detected main isotopologues lines HCN and HNC in Elias 29 (rarer isotopologues are not detected). Moreover, the line profiles are a mixture of absorption and emission, which could be an indication for emission in one of the o ff positions (wobbler switching mode). Therefore, we do not report integrated intensities for Elias 29.

Apotential contamination from the o ff positions for the other sources was assessed by re-pipelining some data sets using only one or the other o ff position for H 13 CN and H 13 CO + , a stronger line in the same frequency setup. The resulting spectra, which can be found in appendix A, were then examined for di ff erences. There is no evidence for a substantial emission from the o ff positions.

Table 4. Peak temperatures Tpeak, integrated intensities ∫ T mb dv , line widths dv , and central peak position v of the HCN and HNC isotopologues. The integrated intensites for IRAS 16293A were obtained from integration over the line, while Gaussian fitting was used for the two other sources. Also listed are the obtained noise level (Trms), error on the integrated intensity σ rms , and the total error σ tot (combined error of the integrated intensity and the calibration).

| Transition        | T peak [K]   | ∫ T mb dv [K km s - 1 ]   | dv [km s - 1 ]    | v [km s - 1 ]        | T rms [K]   | σ rms [K km s - 1 ]   | σ tot [K km s - 1 ]   |
|-------------------|--------------|---------------------------|-------------------|----------------------|-------------|-----------------------|-----------------------|
| IRAS 16239A       | IRAS 16239A  | IRAS 16239A               | IRAS 16239A       | IRAS 16239A          | IRAS 16239A | IRAS 16239A           | IRAS 16239A           |
| HCN(3-2)          | 9.66         | 27.52                     | 6 . 81 ± 0 . 02   | 3 . 72 ± 0 . 01      | 0.04        | 0.03                  | 2.3                   |
| H 13 CN(3-2)      | 0.68         | 3.24                      | 5 . 70 ± 0 . 08   | 3 . 75 ± 0 . 03      | 0.03        | 0.02                  | 0.3                   |
| HC 15 N(3-2)      | 0.27         | 1.37                      | 6 . 2 ± 0 . 2     | 3 . 71 ± 0 . 09      | 0.04        | 0.03                  | 0.1                   |
| H 13 CN(4-3)      | 1.08         | 5.56                      | 6 . 11 ± 0 . 09   | 4 . 28 ± 0 . 04      | 0.07        | 0.04                  | 0.8                   |
| HC 15 N(4-3)      | 0.44         | 2.02                      | 6 . 4 0 . 2       | 3 . 11 ± 0 . 09      | 0.06        | 0.04                  | 0.3                   |
| DCN(3-2)          | 0.87         | 2.77                      | ± 3 . 6 ± 0 . 1   | 3 . 87 ± 0 . 04      | 0.05        | 0.03                  | 0.2                   |
| DCN(4-3)          | 0.76         | 2.45                      | 3 . 74 ± 0 . 07   | 3 . 92 ± 0 . 03      | 0.05        | 0.03                  | 0.4                   |
| DCN(5-4)          | 0.80         | 3.04                      | 4 . 47 ± 0 . 14   | 3 . 86 ± 0 . 06      | 0.09        | 0.05                  | 0.4                   |
| D 13 CN(3-2)      | -            | -                         | -                 | -                    | 0.02        | -                     | -                     |
| HNC(3-2)          | 6.97         | 12.93                     | 2 . 934 ± 0 . 002 | 3 . 959 ± 0 . 001    | 0.02        | 0.009                 | 1.1                   |
| HNC(4-3)          | 5.57         | 11.00                     | 3 . 14 ± 0 . 02   | 3 . 671 ± 0 . 001    | 0.11        | 0.05                  | 0.9                   |
| HN 13 C(3-2)      | 0.40         | 0.77                      | 3 . 0 ± 0 . 3     | 3 . 77 ± 0 . 09      | 0.05        | 0.03                  | 0.1                   |
| H 15 NC(3-2) 0.09 |              | 0.22                      | 2 . 5 ± 0 . 4     | 4 . 1 ± 0 . 2        | 0.03        | 0.01                  | 0.03                  |
| R CrA IRS7B       | R CrA IRS7B  | R CrA IRS7B               | R CrA IRS7B       | R CrA IRS7B          | R CrA IRS7B | R CrA IRS7B           | R CrA IRS7B           |
| HCN(3-2) 6.59     |              | 29 . 48 ± 0 . 03          | 4 . 202 ± 0 . 005 | 5 . 987 ± 0 . 002    | 0.02        | 0.01                  | 2.4                   |
| H 13 CN(3-2)      | 0.90         | 1 . 87 ± 0 . 03           | 1 . 96 ± 0 . 03   | 5 . 92 ± 0 . 01      | 0.04        | 0.02                  | 0.2                   |
| HC 15 N(3-2) 0.24 |              | 0 . 45 ± 0 . 02           | 1 . 79 ± 0 . 11   | 5 . 92 ± 0 . 05      | 0.04        | 0.02                  | 0.04                  |
| H 13 CN(4-3)      | 0.46         | 0 . 92 ± 0 . 02           | 1 . 87 ± 0 . 05   | 6 . 05 ± 0 . 02      | 0.04        | 0.01                  | 0.1                   |
| HNC(3-2)          | 6.35         | 14 . 3 ± 0 . 1            | 2 . 12 ± 0 . 02   | 5 . 925 ± 0 . 008    | 0.2         | 0.07                  | 1.2                   |
| HN 13 C(3-2)      | 0.45         | 0 . 75 ± 0 . 02           | 1 . 56 ± 0 . 05   | 5 . 6 ± 0 . 02       | 0.04        | 0.02                  | 0.06                  |
| H 15 NC(3-2)      | 0.14         | 0 . 20 ± 0 . 01           | 1 . 40 ± 0 . 08   | 5 . 82 ± 0 . 03      | 0.02        | 0.006                 | 0.02                  |
| HN 13 C(4-3)      | 0.23         | 0 . 32 ± 0 . 02           | 1 . 32 ± 0 . 08   | 5 . 19 ± 0 . 04      | 0.04        | 0.01                  | 0.05                  |
| OMC-3 MMS6        | OMC-3 MMS6   | OMC-3 MMS6                | OMC-3 MMS6        | OMC-3 MMS6           | OMC-3 MMS6  | OMC-3 MMS6            | OMC-3 MMS6            |
| HCN(3-2) 3.25     |              | 13 . 10 ± 0 . 03          | 3 . 786 ± 0 . 009 | 11 . 119 ± 0 . 002   | 0.01        | 0.01                  | 1.09                  |
| H 13 CN(3-2)      | 0.51         | 0 . 53 ± 0 . 02           | 0 . 98 ± 0 . 04   | 11 . 12 ± 0 . 01     | 0.03        | 0.01                  | 0.05                  |
| HC 15 N(3-2) 0.10 |              | 0 . 10 ± 0 . 02           | 1 . 00 (fixed)    | 11 . 14 ± 0 . 07     | 0.03        | 0.01                  | 0.01                  |
| H 13 CN(4-3)      | 0.18         | 0 . 19 0 . 04             | 1 . 00 (fixed)    | 11 . 09 0 . 10       | 0.06        | 0.02                  | 0.03                  |
| HNC(3-2)          | 5.55         | ± 6 . 86 ± 0 . 09         | 1 . 16 ± 0 . 02   | ± 11 . 177 ± 0 . 007 | 0.19        | 0.06                  | 0.57                  |
| HN 13 C(3-2)      | 0.46         | 0 . 36 ± 0 . 01           | 0 . 74 ± 0 . 03   | 10 . 85 ± 0 . 01     | 0.03        | 0.008                 | 0.03                  |
| H 15 NC(3-2) 0.09 |              | 0 . 054 ± 0 . 004         | 0 . 54 ± 0 . 04   | 11 . 07 ± 0 . 02     | 0.01        | 0.003                 | 0.005                 |
| HN 13 C(4-3)      | -            | -                         | -                 | -                    | 0.01        | -                     | -                     |

First order baselines were applied to the individual scans before averaging them. In particular, the DCN(5-4) spectrum su ff ers from standing waves, which cannot easily be fitted with a sinusoidal baseline, as amplitude and period change. Therefore, the first order polynomial baseline fits were applied locally whenever required. For lines that were covered by both XFFTS units (overlapping spectral range), the processed spectra from the individual parts were finally averaged after visual inspection.

## 3. Results

## 3.1. Line intensities

The line intensities were obtained by Gaussian fitting except for IRAS 16293A, because the line shapes of this source differ substantially from a simple Gaussian profile (cf. Fig. 1). For IRAS 16293A, the intensities were derived from integrating over a velocity range around the source velocity ( -10 km s -1 &lt; /v1 &lt; 16 km s -1 for broad lines and -4 km s -1 &lt; /v1 &lt; 12 km s -1 for narrow lines). Table 4 lists the resulting integrated line intensi- ties, widths, and the obtained noise level (rms). The spectra from R CrA IRS7B and OMC-3 MMS6 are shown in Figs. 2 and 3, respectively.

While the HCN(3-2) line has prominent outflow wings in IRAS 16293A and OMC-3 MMS6, the line width of HNC(3-2) is narrower and does not feature the broad wings at the obtained signal-to-noise ratio.

For IRAS 16293A, H 13 C 15 N(3-2) and D 13 CN(3-2) were targeted as well, but not detected at the achieved noise level. The H 13 C 15 N(3-2) was only observed for a fraction of the time intended, and the sideband was changed during the observation because of a better receiver performance in the upper sideband. Unfortunately, the upper sideband tuning su ff ers from a tuning instability, which went unnoticed during the observations. A reflection from the tuning instability shows up in one channel at the tuning frequency, and a reliable determination of the noise level is therefore not possible. Hence, we do not give an upper limit on the H 13 C 15 N(3-2) intensity.

In the Class I source Elias 29, the only lines detected are HCN(3-2) and HNC(3-2), which are factors of 50-100 times

Fig. 2. Spectra of the HCN and HNC isotopologues observed with APEX toward R CrA IRS7B. Where needed, the spectra have been scaled by the factors indicated in the individual panels

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0285v1-figure-59344d5b-cd5c-48c6-b22d-86ea446ad9af.jpg)

weaker than the same lines from the Class 0 targets. Therefore, none of the rarer isotopologues are detected in Elias 29 at the achieved rms, maybe because of the much smaller envelope mass compared to the Class 0 sources in the sample. Both HCN(3-2) and HNC(3-2) have a line profile that shows emission at velocities higher than the source velocity, and absorption at the lower velocities, which could either be a real source feature or indicate that there was some molecular emission at the o ff position.

## 3.2. 14 N/ 15 N isotopic ratios

Isotope ratios can generally be measured from molecular lines at sub-millimeter wavelengths using three di ff erent approaches, each with its advantages and caveats:

- 1. Singly substituted molecules: The 14 N 15 / N ratio can be measured from, e.g., CN isotopologues (C 15 N Adande &amp; Ziurys 2012; Hily-Blant et al. 2013b), N2H + isotopologues ( 15 NNH , + 15 NNH , + Womack et al. 1992; Bizzocchi et al. 2010, 2013) or ammonia isotopologues ( 15 NH3, Guesten &amp; Ungerechts 1985; Lis et al. 2010). The advantage of this method is that the main isotopologue is usually abundant and thus its lines are strong. On the other hand, the high abundance often implies a high optical depth, which must be accounted for using radiative transfer modeling. If the molecular transition has a resolved hyperfine structure pattern, the intensity ratios of the hyperfine transitions can also be used to infer and correct for the optical depth.
- 2. Doubly substituted molecules: An alternative to the often optically thick main isotopologues are molecules with two atoms replaced by a less abundant isotope, like the ammonia isotopologue 15 NH2D(Gerin et al. 2009). However, the lines

Fig. 3. Spectra of the HCN and HNC isotopologues observed with APEX toward OMC-3 MMS6. Where needed, the spectra have been scaled by the factors indicated in the individual panels

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0285v1-figure-1e954b28-b2e3-4bcd-8f01-663fde2d428e.jpg)

of such doubly substituted molecules are often very weak because of their low abundance.

- 3. Double isotope method: The nitrogen isotope ratio can also be measured from two singly substituted molecules, e.g., H 13 CN and HC 15 N or HN 13 Cand H 15 NC, avoiding the main isotopologue. However, in this example the 12 C 13 / C ratio needs to be known from independent measurement or an assumption has to be made.

Using the double isotope method and assuming that all lines are optically thin, which will be discussed in Sect. 3.3, the 14 N 15 N ratio can be obtained from the H / 13 CN and HC 15 N integrated intensities using the relation

$$\begin{smallmatrix} \text{itio can} & \frac { \frac { 1 4 } { 1 5 } } { 2 5 N } & \frac { 1 4 } { 1 5 N } & = & \frac { \int T d v \, ( H ^ { 1 3 } C N ) } { \int T d v \, ( H C ^ { 1 5 } N ) } \times \frac { 1 2 } { 1 3 C } & & & ( 1 ) \\ \int \frac { 1 } { 2 } \infty. & \dots & \dots &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. & \end{smallmatrix}$$

and from HNC isotopologues accordingly. For IRAS 16293A, the non-detection of D 13 CN hampers a direct determination of the 12 C 13 / Cfrom D 13 CN and DCN, but only provides a 3 σ upper limit of 12 C 13 C / /lessorsimilar 120. Because no 12 C 13 / C ratio measurements for our sources are available from the literature, we assume the local standard ISM value of 69 (Wilson 1999; Milam et al. 2005). The resulting ratios are illustrated in Fig. 4, together with solar system measurements from the literature, and listed in Table 5. The table also includes the ratio from the integrated intensities of H 13 CN(4-3) and HC 15 N(4-3) reported in Watanabe et al. (2012). The H 13 CN(4-3) line may have a contribution of SO2(132 12 , -121 11), but given that the isotopic ratios , derived from H 13 CN(3-2) and (4-3) in IRAS 16293A agree well, the contamination is expected to be minor in this source. For the other sources, H 13 CN(4-3) was not used in the determination of

Fig. 4. Left panel: Illustration of the nitrogen isotope heterogeneity of the solar system, adapted from Meibom et al. (2007); Mumma &amp; Charnley (2011). Comets, meteorites, and interstellar dust particles (IDPs) are located over a range of distances from the Sun and the plotted d Sun values should be regarded as an average value only. Right panel: 14 N 15 / N ratios obtained from APEX measurements (black, this work) and ASTE (grey, Watanabe et al. 2012). Diamonds and triangles indicate values from HCN and HNC isotopologues, respectively. Error bars are 1 σ for combined uncertainties on the calibration and the intensity. The gray bars ('cores') indicate the prestellar core values from HCN by Hily-Blant et al. (2013a).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0285v1-figure-b8ee3279-0dd5-40b8-af16-d68b205ac4f7.jpg)

Table 5. 14 N 15 / N ratios derived from APEX observations of HCN and HNC isotopologues using a standard 12 C 13 / C ratio of 69, except for the H 13 CN HC 15 N(4-3) ratio, where the data are taken from the survey by / Watanabe et al. (2012).

all errors. In the cases where the 14 N 15 N ratios are derived from / HCN isotopologues, the errors listed should therefore be an upper limit on the uncertainty.

| Source      | H 13 CN(3 - 2) HC 15 N(3 - 2)   | HN 13 C(3 - 2) H 15 NC(3 - 2)   | H 13 CN(4 - 3) HC 15 N(4 - 3)   |
|-------------|---------------------------------|---------------------------------|---------------------------------|
| IRAS 16293A | 163 ± 20                        | 242 ± 32                        | 190 ± 38                        |
| R CrA IRS7B | 287 ± 36                        | 259 ± 34                        | (285 ± 139)                     |
| OMC3 MMS6   | 366 86                          | 460 65                          | -                               |

±

±

the 14 N 15 N ratio, just as a constraint for the radiative transfer / models.

The smallest error is expected on the ratio inferred from H 13 CN(3-2) and HC 15 N(3-2) for all sources as well as the H 13 CN(4-3) and HC 15 N(4-3) for IRAS 16293A, because the lines were measured simultaneously. The 14 N 15 / N ratios from the HNC isotopologues are more uncertain, as the involved lines were observed in two individual spectral setups. For the errors, we use a combined error from the intensity measurement and the calibration uncertainty. The error on the integrated intensity is taken from the uncertainty in the fit of the area under the Gaussian profile, provided by the CLASS software, for R CrA IRS7B and OMC-3 MMS6, and for IRAS 16293A calculated as σ rms = √ ∆ /v1 × δ v × T rms where ∆ /v1 is the full width at half maximum of the line (assumed to be the same as that of a Gaussian fit to the line), δ v is the velocity resolution of the spectrum, and T rms is the noise level. For the calibration error σ calib , we used the calibration uncertainties on the integrated intensity provided in the APEX calibration report (Dumke &amp; Mac-Auli ff e 2010) of 8.3% and 14.1% for APEX-1 and APEX-2, respectively. The overall error on a line intensity is then the combination of the two errors, σ tot = √ σ 2 rms + σ 2 calib , and propagated on to the intensity ratio using regular error propagation. Even though some lines were measured in a joint spectral setup, the full calibration uncertainty is included in the determination of

## 3.3. Optical depth considerations

## 3.3.1. Hyperfine structure fitting for IRAS 16293A

The assumption that the lines, in particular H 13 CN(3-2), are optically thin and hence their intensity is a direct measure of the number of molecules, needs to be verified.

For IRAS 16293-2422, we have a H 13 CN(1-0) spectrum (Jørgensen et al. 2004) from the SEST telescope 4 with clearly resolved hyperfine structure (cf. Fig. 5, using a main beam e ffi ciency of 0.75). The pointing is slightly di ff erent ( α = 16:32:22.80, δ = -24:28:33.0) from our recent observations, but given the large beam the spectrum should still provide an estimate of the opacity. The three hyperfine transitions were fitted using the 'HFS method' of the CLASS software to derive the opacity. The HFS routine fits Gaussians to all hyperfine components simultaneously under the assumption of a common excitation temperature, with relative positions and intensities provided by the user (see the CLASS manual or appendix A of Hily-Blant et al. 2013a, for more details). The best fit, illustrated in Fig. 5, yields an opacity for the main component of τ main = 0 274 . ± 0 405 and hence the line is optically thin. .

## 3.3.2. Radiative transfer modeling

The optical depth on higher excited lines can be larger than the opacity of H 13 CN(1-0). Therefore, line radiative transfer models can be used as an alternative method to estimate the optical depth in the isotopologues and check the validity of Eq. 1. The spherically symmetric source model from Schöier et al. (2002) for IRAS 16293-2422, the R CrA IRS7B model with external irradi-

4 The Swedish-ESO 15m Submillimeter Telescope (SEST) was decommissioned in 2003.

Fig. 5. H 13 CN(1-0) spectrum of IRAS 16293-2422 from the SEST telescope. The solid line shows the best fit from the CLASS software.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0285v1-figure-f9e77e9b-0da1-42d2-b403-f4224903a82d.jpg)

ation from Lindberg &amp; Jørgensen (2012), and the best fit model of OMC-3 MMS6 from Jørgensen et al. (2006) are used as an input for the 1D version of the 'RATRAN' radiative transfer code (Hogerheijde &amp; van der Tak 2000). A comparison of the source models, i.e., density and temperature profiles, is shown in Fig. 6. We use molecular data for the HCN isotopologues from the LAMDAdatabase 5 (Schöier et al. 2005). The transition frequencies, energy levels, and Einstein A coe ffi cients in the LAMDA files are based on data from the CDMS catalog (Müller et al. 2001, 2005), and the collisional rate coe ffi cients are scaled and extrapolated HCN-He collision rates from Green &amp; Thaddeus (1974). Individual models were checked against newer collision rates from Dumouchel et al. (2010) and the di ff erences were found to be at the 5-10% level in these cases, i.e., well within the calibration uncertainty.

The abundance profile is chosen according to the results by Jørgensen et al. (2005a), i.e., a drop abundance profile accounting for freeze-out of the molecule at temperatures below a certain threshold and su ffi cient densities. Following Jørgensen et al. (2005a), the freeze-out region for HCN in NGC 1333 IRAS 2A is defined as the part of the envelope with densities above n &gt; 7 × 10 -4 cm -3 and temperatures below T &lt; 90 K, which corresponds to the region where CO and HCO + are also depleted. However, the density at which HCN is frozen out is not a fixed value but time and source dependent (Jørgensen et al. 2005b).

Therefore, models with three di ff erent freezeout evaporation / temperatures (30 K, 60 K, and 90 K) and various abundances were tested for IRAS 16293-2422. The density is higher than n &gt; 7 × 10 -4 cm 3 throughout the envelope -in this model and a depletion zone was therefore not employed. A Doppler b-parameter of db = 2 3 . km s -1 was assumed, corresponding to a line width (FWHM) of ∼ 3 8 km s . -1 . The model FWHM was chosen a bit lower than the observed line

5 http://home.strw.leidenuniv.nl/~moldata/

Fig. 6. Temperature and density as function of radius for IRAS 162932422 (solid line), OMC-3 MMS6 (dotted), and R CrA IRS7B (dashed).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0285v1-figure-087a9d8c-df9f-45ee-a343-cf9e60eb15a0.jpg)

widths, because no additional velocity profile (e.g. infall) was used in the models.

The model with the abundance jump at 90 K and x out ≈ 6 × 10 -12 (for T &lt; 90 K) and x in ≈ 1 × 10 -9 (where T ≥ 90 K) can reproduce the observed H 13 CN(3-2 and 4-3) integrated intensities from IRAS 16293-2422 within a few percent. However, this model underproduces the HC 15 Nlines, in particular the J = 4 -3 transition, by at least 20% for H 13 CN HC 15 N abundance ratios / in the range of 2 6 . -6 3, corresponding to a . 14 N 15 N ratio of / 180 -440 assuming 12 C 13 C / = 69 (Wilson 1999). In addition, an abundance jump of more than two orders of magnitude is almost a factor of five higher than what was found in earlier work (Jørgensen et al. 2005a; Brinch et al. 2009). Therefore, we a lso tested models where the jump occurs at 60 K and 30 K instead of 90 K. For 60 K, the H 13 CN intensities are well modelled using x in ≈ 2 × 10 -10 and x out ≈ 8 × 10 -12 . Using an HC 15 N abundances of x in ≈ 7 7 . × 10 -11 and x out ≈ 3 1 . × 10 -12 , corresponding to 14 N 15 N / = 180 as obtained using Eq. 1, yields a very good fit to the observed integrated intensities of HC 15 N. In contrast, using a solar ratio of 14 N 15 N / = 441 ( x in ≈ 3 1 . × 10 -11 and x out ≈ 1 3 . × 10 -12 ) underproduces the line intensities by ∼ 50%. For 30 K, a good fit to the H 13 CNlines is obtained using x in ≈ 4 × 10 -11 and x out ≈ 5 × 10 -12 . Again, 14 N 15 N / = 180 ( x in ≈ 1 5 . × 10 -11 and x out ≈ 1 9 . × 10 -12 ) gives an excellent fit to the HC 15 N data, whereas a solar abundance ratio, x in ≈ 6 3 . × 10 -12 and x out ≈ 7 8 . × 10 -13 , underestimates the observed integrated intensites by ∼ 50%.

So while the 90 K jump model, which becomes moderately optically thick for H 13 CN ( τ 3 -2 ≈ 1 5 towards line and source . center), can reproduce the H 13 CN lines, it cannot fit HC 15 N for 14 N 15 N ratios ranging from the most extreme values found from / the simple analysis (Eq. 1) to solar abundances. The optically

thin models with the jump occuring at 60 K and 30 K on the other hand are also able to reproduce the HC 15 N integrated intensities for 14 N 15 N / = 180, but not for solar abundances. Therefore, the radiative transfer results support the assumption that the line emission is optically thin, and that Eq. 1 is valid in this case.

For OMC-3 MMS6, a fairly well constrained model is also available. An abundance profile with x ≈ 8 × 10 -13 for T &lt; 90 K and n &gt; 1 × 10 6 cm 3 -and x ≈ 7 × 10 -11 elsewhere yields integrated intensities that are consistent with the observed values within the uncertainties ( ∼ 20%). The corresponding opacities are in the range of τ ∼ 0 5 for . db = 0 6 (FWHM . ∼ 1 km s -1 ) and indicate that the H 13 CN emission from this source is optically thin. For HC 15 N(3-2), the integrated intensity can be reproduced within 20% using models with 14 N 15 N / = 370 (as derived from Eq. 1) as well as 14 N 15 N / = 441 (solar ratio). As already concluded in Sec. 3.2, it is therefore not possible to distinguish between a slightly fractionated and a solar 14 N 15 N ra-/ tio for OMC-3 MMS6. The higher freeze-out density compared to IRAS 16293-2422 could indicate that this source is in a less evolved state.

Because of the complex physical structure of R CrA IRS7B, including a disk (Lindberg et al. 2014), the spherically symmetric model is less reliable for this source than for IRAS 162932422 and OMC-3 MMS6. An agreement within 20% of the modeled and the observed intensities can be obtained with an H 13 CN abundance of x ≈ 2 × 10 -10 for T &lt; 90 K and n &gt; 7 × 10 4 cm -3 and x ≈ 2 × 10 -8 elsewhere. The optical depth in the central line channel and for a ray directly through the source center are τ ∼ 1 -3. The observed HC 15 N emission can be reproduced within the uncertainties with a range of 14 N 15 N ratios / ( ∼ 270 -441) because the model is not well constrained and only one HC 15 N line is available. However, the required high abundances are an indication that the physical model underestimates the observed column density, which results in a higher optical depth. Compared to the IRAS 16293-2422 model, which has a similar outer radius, but higher densities, R CrA IRS7B is not expected to be significantly optically thick. Given the model uncertainties and the degeneracies in the parameter space, the simple analysis (Eq. 1) is likely to provides the most reliable results for this source.

## 4. Discussion

At least two measurements of the 14 N 15 N ratio have been ob-/ tained in all three sources. For IRAS 16239A, the two measurements from the J = 3-2 and 4-3 lines of H 13 CN and HC 15 N, 14 N 15 / N = 163 ± 20 and 14 N 15 / N = 190 ± 38, agree well. The 14 N 15 N / = 242 ± 32 value from the HNC(3-2) isotopologues is di ff erent from the HCN values at the 1 σ level. In the case of R CrA IRS7B, the HCN and HNC values are consistent within the error bars, including also the ratio derived from the data by Watanabe et al. (2012). The 14 N 15 N ratios from HCN and / HNC in OMC-3 MMS6 agree within errors, but the uncertainties are larger in this source than in the two others. Overall, the agreement between di ff erent measurements of the same source is fairly good. The comparison to other measurements and literature values is discussed in the subsections below.

Comparing Table 5 with the temperatures in the outer regions of the spherically symmetric envelope models of the three sources in the sample, it is apparent that the lowest 14 N 15 N ratios are found in the source with the lowest temper-/ ature in the outermost regions, IRAS 16293A. It is followed by R CrA IRS7B, which has a somewhat elevated temperature because it is being irradiated by a nearby Herbig Be star,

Fig. 7. Average 14 N 15 / Nratios of the three young stellar objects plotted versus their outer envelope temperature at the projected beam radius. Error bars reflect the spread in measured 14 N 15 / N ratios.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0285v1-figure-af905ece-fbf4-4d76-96f2-5bdce88cc403.jpg)

R CrA (Lindberg &amp; Jørgensen 2012). The highest, almost solar , 14 N 15 N ratios are from OMC-3 MMS6, located in the active re-/ gion of Orion and thus subject to a strong radiation field, which has the highest temperatures in the outer envelope. More quantitatively, this behavior is illustrated in Fig. 7, showing the average of the 14 N 15 N ratio for each source versus the temperature / at the projected radius of the beam. The APEX measurements listed in Table 5 were averaged for each source and the errors bars reflect the spread of the measured ratios. The temperature was taken from the spherically symmetric source models, evaluated at the radius of the APEX beam ( ∼ 24 ′′ ) at the corresponding source distance ( T = 16 5 K for IRAS 16239A, 25 5 K for . . R CrA IRS7B, and 37 K for OMC-3 MMS6). A larger sample would be needed to draw any conclusions about whether or not there is a direct correlation between temperature and 15 Nenhancement, but our result provides at least an indication for a temperature dependence of the 14 N 15 N ratio in these embedded / protostars.

## 4.1. Comparison with previous HCN and HNC measurements in the target sources

Isotopologues of HCN and HNC in IRAS 16293-2422 were previously measured in the spectral surveys by Caux et al. (2011, IRAM and JCMT) and Blake et al. (1994); van Dishoeck et al. (1995, JCMT and CSO). The integrated intensities from the Caux et al. (2011) and van Dishoeck et al. (1995) surveys are listed in Table B.2 in the appendix. Note that the observations by Caux et al. (2011) were pointed at the B component of the protobinary, whereas our observations are pointed directly at the A component. The two components are separated by ∼ 4 ′′ and thus lie both within the beam. However, in the 230 GHz band observations at the IRAM 30m, the flux of the A component is considerably attenuated (see discussion in Caux et al. 2011). The resulting 14 N 15 N ratios are, assuming a calibration uncer-/ tainty of 15%, 276 ± 74 from the H 13 CN(1-0) and HC 15 N(1-0) intensities, 219 ± 71 for the 3-2 transitions, and 320 ± 106 for the 4-3 lines and thus comparable to our values. No H 15 NC intensities are reported by Caux et al. (2011).

From the integrated intensities of HC 15 N(3-2) and H 13 CN(32) measured by van Dishoeck et al. (1995) we find 14 N 15 N / = 131. From the fact that HC 15 N is detected, they argue - pre-

sumably based on the assumption of a local interstellar medium value for the nitrogen isotope ratio (Wilson &amp; Rood 1994) - that H 13 CN is still slightly optically thick. However, from our model results and the fit to the hyperfine structure of H 13 CN(1-0) we conclude that the lines are optically thin and that the assumption of a standard nitrogen isotope ratio likely does not hold because nitrogen is fractionated.

R CrA IRS7B was surveyed in the 345 GHz spectral window by Watanabe et al. (2012) using the ASTE telescope. The integrated intensities are 0 23 . ± 0 09 K km s . -1 for HC 15 N(4-3), 0 95 . ± 0 10 K km s . -1 for H 13 CN(4-3), 0 32 . ± 0 06 K km s . -1 for HN 13 C(4-3), and 0 05 . ± 0 04 K km s . -1 for H 15 NC(4-3). Based on exactly the same procedure as for our own measurements but with a calibration uncertainty of 19% for ASTE (Watanabe et al. 2012), the inferred 14 N 15 N / ratios are 285 ± 139 from HCN and 442 ± 382 from HNC. The ratios derived from the work by Watanabe et al. (2012) are consistent with our measurements, but the error bars are significantly larger than for our APEX observations.

Schöier et al. (2006) observed HCN(4-3) and HNC(4-3) with APEX in R CrA IRS7B, resulting in integrated intensities of 30 4 K km s . -1 and 11 3 K km s . -1 , respectively. Their intensities are very similar to the values found from our observations of HCN(3-2) and HNC(3-2), which are 30 2 K km s . -1 and 14 2 K km s . -1 . For optically thin HCN lines, similar line strengths for the 4-3 and 3-2 transitions are predicted by radiative transfer models if a constant abundance is assumed (see also Jørgensen et al. 2005a, for a discussion of the HCN abundance structure). Therefore, the similar intensities of the 4-3 and 32 transitions could be an indication for little variations in the HCN abundance with radius in this source, and thus no significant freeze-out of HCN.

No HCN and HNC isotopologue measurements were found in the literature for OMC-3 MMS6.

## 4.2. Comparison with 14 N/ 15 N ratios from the literature

14 N 15 N ratios were previously measured in a selection of low-/ mass pre- and protostellar environments from a variety of molecular tracers. The literature results are summarized in Table 1. The table clearly illustrates the discrepancy of 14 N 15 N ratios / derived from HCN and NH3, which may be caused by different fractionation routes for nitriles and nitrogen hydrides (Wirström et al. 2012; Hily-Blant et al. 2013a,b). It also highlights the lack of a comprehensive set of 14 N 15 N ratios from / di ff erent tracers for a source, which makes the interpretation and comparison of individual measurements very challenging. For instance, Adande &amp; Ziurys (2012) found a reasonable agreement between the 14 N 15 / N ratios inferred from HNC and CN. We find that the ratios derived from HCN and HNC are relatively similar. On the other hand, the prestellar core nitrogen isotope composition varies between HCN and CN (Hily-Blant et al. 2013a,b), so it is not clear if there is any correlation between the 14 N 15 N ratios from HCN, HNC, and CN despite the / fact that they all belong to to the 'nitrile' fractionation route (Wirström et al. 2012). Moreover, the HNC values reported in Milam &amp; Charnley (2012) are substantially lower than any other measurement for reasons which are unclear.

For IRAS 16293A, the inferred ratios from HCN (two transitions) and HNC (one measurement) fall into the range of the typical prestellar core HCN values (Hily-Blant et al. 2013a). The similarity of the nitrogen isotope ratio in IRAS 16293A to the prestellar cores is an indication that either fractionation (not necessarily chemical) continues in the envelope of this YSO or that an enhanced 15 N signature from the prestellar core phase has survived into the Class 0 phase. An enhanced radiation field seems to lower or even suppress nitrogen fractionation, as can be seen for the two sources with an enhanced external radiation field, R CrA IRS7B and OMC-3 MMS6. The 14 N 15 N ratio in / R CrA IRS7B lies at the upper end of the prestellar core values and the 14 N 15 N ratio from OMC-3 MMS6 is even higher, / outside the range spanned by the prestellar cores.

## 4.3. Comparison to the local interstellar medium

To decide whether a source shows an enhancement of 15 N in a molecule, one would like to compare the 14 N 15 N / ratio from the source to a reference value of the interstellar medium (ISM). Di ff erent 14 N 15 N measurements of the local interstel-/ lar medium are available from the literature. Wilson (1999) list 14 N 15 / N = 388 ± 32 for the local ISM (updated from the earlier value of 14 N 15 N / = 450 ± 22, Wilson &amp; Rood 1994) based on the HCN surveys data from Wannier et al. (1981, for the northern sources) and Dahmen et al. (1995, for the southern sky) as well as the NH3 survey by Guesten &amp; Ungerechts (1985). However, Wielen &amp; Wilson (1997) argue that the ISM value should be 14 N 15 N / = 414 ± 32, taking into account that the Sun formed at a smaller galactocentric radius that its current location. A local ISM value of 14 N 15 N / = 290 ± 40 is reported by Adande &amp; Ziurys (2012) from measurements of HNC and CN in warm clouds and the 12 C 13 / C ratios from Milam et al. (2005). While the local value from the work by Adande &amp; Ziurys (2012) is inferred from a linear fit to the data throughout the Galaxy, the local ratio reported in Wilson (1999) is the average of 8 sources close to the solar galactic distance. The Wilson (1999) and Adande &amp; Ziurys (2012) results are not consistent within error bars, and the di ff erence presumably arises from the di ff erent carbon isotope ratios ( 12 C 13 / C) assumed for the double isotope method in the case of HCN and HNC. The 12 C 13 / C ratio from CN and CO isotopologues, used in Adande &amp; Ziurys (2012), are generally lower than those obtained from H2CO isotopologues, included in Wilson (1999). However, the local 12 C 13 / C ratios from Wilson (69 ± 6 from CO and H2CO, 1999) and Milam et al. (2005, 68 ± 15 from CO, CN, and H2CO) agree well, but the values for individual sources close to the solar neighborhood from Milam et al. (2005) are partly lower, which is then transferred onto the 14 N 15 N ratios from HNC (Adande &amp; Ziurys 2012). The / Adande &amp; Ziurys (2012) result agrees with the average value of 14 N 15 N / = 237 + 27 -21 measured in two di ff use clouds along the line of sight toward a strong extragalactic source from absorption spectroscopy of HCN (Lucas &amp; Liszt 1998). It was however noticed that there is a significant scatter of the data points from the Galactic isotope surveys, and thus it is unclear how meaningful the comparison of the obtained 14 N 15 N ratios to an average / local ISM value is to decide whether or not the nitrogen isotopes are fractionated in a specific source. Moreover, the local ISM value is also subject to change, as galactic evolution models by Romano &amp; Matteucci (2003) for instance predict that the 14 N 15 N ratio of the ISM increases with galactocentric distance / and decreases with time.

Our IRAS 16293A results from HCN lie outside the uncertainty range of both the Adande &amp; Ziurys (2012) and Wilson (1999) local ISM values, while the ratio inferred from HNC lies within the uncertainties of Adande &amp; Ziurys (2012), but clearly outside Wilson (1999), as do the R CrA HCN and HNC ratios. The HCN measurement from OMC-3 MMS6 is agreement with both ISM values, while the ratio from HNC, which is the most uncertain of all data points, is consistent with the Wilson (1999)

ISMresult but falls above the one from Adande &amp; Ziurys (2012). Hence, from the comparison with ISM measurements, the nitrogen isotopes in IRAS 16293A show indications for fractionation, and fractionation may possibly also occur in R CrA IRS7B, while OMC-3 MMS6 shows no evidence for fractionation if the comparison with a local ISM value is meaningful. The 12 C 13 C / ratio should also be measured for our sources to check and account for any deviations from the standard.

## 4.4. Comparison to solar system values

A direct comparison of the nitrogen isotope composition with that of solar system bodies is di ffi cult as possible mechanisms for a separation of the protostellar material into a 15 N-poor gas phase and 15 N-rich solids are not yet understood. The freezeout of 15 N-substituted molecular species, in particular ammonia, onto dust grains and their incorporation into the ice mantles is believed to be important (Rodgers &amp; Charnley 2008b; Wirström et al. 2012), but the detailed processes are still unclear. Existing observations of the nitrogen isotope composition in star-forming regions do not point at a large nitrogen isotope fractionation in ammonia (Gerin et al. 2009; Daniel et al. 2013). Gerin et al. (2009) and Bizzocchi et al. (2013) even argue for a depletion of 15 N in NH3 and N2H , respectively, in the prestel-+ lar core L1544. The comparison is further complicated by the fact that the chemical carrier of the elevated 15 N-signature or the functional group carrying the isotopic anomaly in primitive solar system matter has not yet been identified. Furthermore, the cosmomaterials that can nowadays be sampled may represent only a minor reservoir of the presolar nebula. It may have also underdone significant processing during the evolution of the solar system, such as heating or aqueous alteration, which could a ff ect the isotopic signatures (e.g., Alexander et al. 2007).

Bulk 15 N-enrichments in meteorites are comparable to the values inferred from star-forming regions. However, meteoritic 15 N-hotspots feature enrichments exceeding the one most extreme values for prestellar cores (Hily-Blant et al. 2013a) or our protostar sample. Within errors, the cometary ratios are also in agreement with our HCN results for IRAS 16293A, but lower than the OMC-3 MMS6 and R CrA IRS7B results. Overall, the 14 N 15 N ratios in primitive solar system matter, i.e., meteorites, / comets, and IDPs, are comparable to the ones measured around prestellar cores and protostars (cf. comparison in Fig. 4), but the ambiguities on both the meteoritic and the astrochemical side do not yet allow concluding whether or not the 15 N-signature in meteorites was directly inherited from a 15 N-enriched precursor molecule in the presolar nebula or whether it is the result of a secondary process.

## 4.5. Chemical fractionation or photochemistry?

A possible explanation for the isotopic variations among our sample could be the di ff erent amounts of external irradiation the sources are exposed to. IRAS 16293A, being a relatively isolated source, has the lowest temperatures in the outer envelope, whereas OMC-3 MMS6 and R CrA IRS7B are both exposed to an enhanced radiation field, resulting in warmer outer envelope temperatures (Jørgensen et al. 2006; Lindberg &amp; Jørgen sen 2012). There is a tentative trend toward a decreasing 15 Nenhancement with increasing outer envelope temperature in our sample, but the statistics is too small to draw any conclusions. If the trend can be confirmed, it would argue in favor of the chemical fractionation scenario over isotope-selective photochem- istry. If photodissociation of N2 was the process dominating the nitrogen isotope composition, the highest 15 N-enrichments should be present in the two externally irradiated sources, not IRAS 16293A. However, the single-dish observations only provide a spatially averaged value of the 14 N 15 N ratio. / Because fractionation may not occur everywhere in the envelope, and because the mechanism that dominates the nitrogen isotope composition may even vary throughout the envelope depending on physical parameters such as temperature, irradiation, and opacity, spatially resolved observations are clearly needed to study potential spatial variations in 14 N 15 N ratio. / For instance, photochemistry could determine the 14 N 15 N ratio in the outermost / and maybe innermost parts close to the protostar, where the radiation field is intense, whereas chemical fractionation is expected to be important in cold, shielded regions. Our understanding of the evolution of nitrogen isotopic composition during star- and planet formation would also benefit from a more comprehensive data set with nitrogen isotope ratios measured from a census of molecular tracers and sources in di ff erent evolutionary stages. Nonetheless, our results indicate that enhanced 15 N-values in HCN and HNC are present not just in the prestellar stage, but also around at least some Class 0 protostars. Whether fractionation is still ongoing in these sources or just preserved from the precedent prestellar core phase and what mechanism is causing the 15 N-enrichment in HCN and HNC is not yet clear and needs to be further investigated.

## 5. Conclusions

We have measured the 14 N 15 N isotope ratio in the three Class 0 / protostars IRAS 16293A, R CrA IRS7B, and OMC-3 MMS6 from 15 N- and 13 C-substituted isotopologues of HCN and HNC, observed with the APEX telescope. Optically thin emission, the same excitation temperature for both isotopologues, and a standard carbon isotope ratio of 12 C 13 C / = 69 (Wilson 1999; Milam et al. 2005) were assumed to derive the 14 N 15 N ratios di-/ rectly from the integrated line intensities. The assumption of optically thin emission is supported by the results from H 13 CN(10) hyperfine structure fitting for IRAS 16293A and by non-LTE line radiative transfer modeling for all three sources.

The ratios derived from the HCN and HNC isotopologues are relatively similar, but there is significant variation between the tree sources. The values for IRAS 16293A ( 14 N / 15 N ∼ 160 -240) are in the range of typical prestellar core values ( 14 N 15 N / ∼ 140 -360, Hily-Blant et al. 2013a), whereas the results for R CrA IRS7B ( 14 N 15 N / ∼ 260 -290) indicate that the nitrogen isotopes in this source are less fractionated. For OMC-3MMS6( 14 N / 15 N ∼ 370 -460), a solar and hence unfractionated nitrogen isotope composition cannot be excluded. The di ff erences between the sources possibly stem from the di ff erent outer envelope temperatures, presumably caused by external irradiation heating the envelope.

Our results indicate that fractionation of nitrogen isotopes in HCN and HNC either continues in the Class 0 stage or is at least preserved from the prestellar core phase. Additional tracers should also be studied to learn more about what molecules carry enhanced 15 N and could be responsible for the 15 N-enrichment observed in solar system materials.

Chemical fractionation seems a more likely scenario than N2 self-shielding to explain the nitrogen isotope measurements in our data set because of a tentative trend toward a decreasing 14 N 15 N with decreasing envelope temperature. / However, spatially resolved observations are crucial to pinpoint which mech-

anism causes the observed 15 N-enrichment in some Class 0 protostars.

## References

Abbas, M. M., LeClair, A., Owen, T., et al. 2004, ApJ, 602, 1063 Adande, G. R. &amp; Ziurys, L. M. 2012, ApJ, 744, 194 Aléon, J. 2010, ApJ, 722, 1342 Alexander, C. M. O´D., Fogel, M., Yabuta, H., &amp; Cody, G. D.

Geochim. Cosmochim. Acta, 71, 4380

- Anders, E. &amp; Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197
- Bertin, M., Fayolle, E. C., Romanzin, C., et al. 2013, ApJ, 779, 120
- Birck, J. L. 2004, Rev. Mineral. Geochem., 55, 25
- Bizzocchi, L., Caselli, P., &amp; Dore, L. 2010, A&amp;A, 510, L5 +
- Bizzocchi, L., Caselli, P., Leonardo, E., &amp; Dore, L. 2013, A&amp;A, 555, A109
- Blake, G. A., van Dishoeck, E. F., Jansen, D. J., Groesbeck, T. D., &amp; Mundy, L. G. 1994, ApJ, 428, 680

Bockelée-Morvan, D., Biver, N., Jehin, E., et al. 2008, ApJ,679, L49

- Böhlke, J. K., de Laeter, J. R., De Bièvre, P., et al. 2005, J. P hys. Chem. Ref. Data, 34, 57
- Bonal, L., Huss, G. R., Krot, A. N., et al. 2010, Geochim. Cosmochim. Acta, 74, 6590
- Boogert, A. C. A., Hogerheijde, M. R., Ceccarelli, C., et al. 2002, ApJ, 570, 708 Briani, G., Gounelle, M., Marrocchi, Y., et al. 2009, Proceedings of the National Academy of Sciences, 106

Brinch, C., Jørgensen, J. K., &amp; Hogerheijde, M. R. 2009, A&amp;A,

502, 199

Busemann, H., Young, A. F., O'D. Alexander, C. M., et al. 2006, Science, 312,

727

- Caux, E., Kahane, C., Castets, A., et al. 2011, A&amp;A, 532, A23

Clayton, R. N. 2002, Nature, 415, 860

- Connelly, J. N., Bizzarro, M., Krot, A. N., et al. 2012, Science, 338, 651

Coplen, T. B., Krouse, H. R., &amp; Böhlke, J. K. 1992, Pure Appl. Chem, 64, 907

Dahmen, G., Wilson, T. L., &amp; Matteucci, F. 1995, A&amp;A, 295, 194

- Daniel, F., Gérin, M., Roue , E., et al. 2013, A&amp;A, 560, A3 ff
- Dumke, M. &amp; Mac-Auli ff e, F. 2010, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7737

Dumouchel, F., Faure, A., &amp; Lique, F. 2010, MNRAS, 406, 2488

Fletcher, L. N., Greathouse, T. K., Orton, G. S., et al. 2014, Icarus, 238, 170

Floss,

C.,

Stadermann,

F.

J.,

Bradley,

J.

P., et

al.

2006,

Geochim.

Cos- mochim. Acta, 70, 2371

Fouchet, T., Irwin, P. G. J., Parrish, P., et al. 2004, Icarus, 172, 50

Fouchet, T., Lellouch, E., Bézard, B., et al. 2000, Icarus, 1 43, 223

- Gerin, M., Marcelino, N., Biver, N., et al. 2009, A&amp;A, 498, L9

Green, S. &amp; Thaddeus, P. 1974, ApJ, 191, 653

- Guesten, R. &amp; Ungerechts, H. 1985, A&amp;A, 145, 241

Güsten, R., Nyman, L. Å., Schilke, P., et al. 2006, A&amp;A, 454, L 13

- Hashizume, K., Chaussidon, M., Marty, B., &amp; Robert, F. 2000, Science, 290, 1142
- Heays, A. N., Visser, R., Gredel, R., et al. 2014, A&amp;A, 562, A61
- Hily-Blant, P., Bonal, L., Faure, A., &amp; Quirico, E. 2013a, Icarus, 223, 582
- Hily-Blant, P., Pineau des Forêts, G., Faure, A., Le Gal, R.,&amp; Padovani, M. 2013b, A&amp;A, 557, A65
- Ho man, J. H., Hodges, R. R., McElroy, M. B., Donahue, T. M., &amp; Kolpin, M. ff 1979, Science, 205, 49

Hogerheijde, M. R. &amp; van der Tak, F. F. S. 2000, A&amp;A, 362, 697

Hutsemékers, D., Manfroid, J., Jehin, E., et al. 2005, A&amp;A, 4 40, L21

Ikeda, M., Hirota, T., &amp; Yamamoto, S. 2002, ApJ, 575, 250

- Jehin, E., Manfroid, J., Hutsemékers, D., Arpigny, C., &amp; Zuc coni, J.-M. 2009, Earth Moon and Planets, 105, 167
- Johnstone, D., Boonman, A. M. S., &amp; van Dishoeck, E. F. 2003, A&amp;A, 412, 157 Jørgensen, J. K., Bourke, T. L., Myers, P. C., et al. 2005a, Ap J, 632, 973

Jørgensen, J. K., Bourke, T. L., Nguyen Luong, Q., &amp; Takakuwa , S. 2011, A&amp;A,

534, A100

- Jørgensen, J. K., Johnstone, D., Kirk, H., et al. 2008, ApJ, 6 83, 822
- Jørgensen, J. K., Johnstone, D., van Dishoeck, E. F., &amp; Doty, S. D. 2006, A&amp;A, 449, 609

Jørgensen, J. K., Schöier, F. L., &amp; van Dishoeck, E. F. 2004, A &amp;A, 416, 603 Jørgensen, J. K., Schöier, F. L., &amp; van Dishoeck, E. F. 2005b, A&amp;A, 435, 177 Junk, G. &amp; Svec, H. J. 1958, Geochim. Cosmochim. Acta, 14, 234 Kerridge, J. F. 1993, Reviews of Geophysics, 31, 423 Li, X., Heays, A. N., Visser, R., et al. 2013, A&amp;A, 555, A14 Liang, M.-C., Heays, A. N., Lewis, B. R., Gibson, S. T., &amp; Yung, Y. L. 2007,

ApJ, 664, L115 Lindberg, J. E. &amp; Jørgensen, J. K. 2012, A&amp;A, 548, A24 Lindberg, J. E., Jørgensen, J. K., Brinch, C., et al. 2014, A&amp; A, 566, A74 Linke, R. A., Goldsmith, P. F., Wannier, P. G., Wilson, R. W., &amp; Penzias, A. A. 1977, ApJ, 214, 50

Lis, D. C., Wootten, A., Gerin, M., &amp; Roue ff , E. 2010, ApJ, 710, L49 Loinard, L., Castets, A., Ceccarelli, C., et al. 2000, A&amp;A, 359, 1169

- 2007,
- Lucas, R. &amp; Liszt, H. 1998, A&amp;A, 337, 246
- Lyons, J. R., Bergin, E. A., Ciesla, F. J., et al. 2009, Geochim. Cosmochim. Acta, 73, 4998
- Manfroid, J., Jehin, E., Hutsemékers, D., et al. 2009, A&amp;A, 5 03, 613
- Mariotti, A. 1983, Nature, 303, 685
- Marty, B., Chaussidon, M., Wiens, R. C., Jurewicz, A. J. G., &amp; Burnett, D. S. 2011, Science, 332, 1533
- Mathew, K. J. &amp; Marti, K. 2001, J. Geophys. Res., 106, 1401
- McKeegan, K. D., Aléon, J., Bradley, J., et al. 2006, Science , 314, 1724

Meibom, A., Krot, A. N., Robert, F., et al. 2007, ApJ, 656, L33

- Messenger, S. 2000, Nature, 404, 968
- Milam, S. N. &amp; Charnley, S. B. 2012, in Lunar and Planetary Inst. Technical Report, Vol. 43, Lunar and Planetary Institute Science Conference Abstracts, 2618
- Milam, S. N., Savage, C., Brewster, M. A., Ziurys, L. M., &amp; Wycko ff , S. 2005, ApJ, 634, 1126
- Müller, H. S. P., Schlöder, F., Stutzki, J., &amp; Winnewisser, G. 2005, Journal of Molecular Structure, 742, 215
- Müller, H. S. P., Thorwirth, S., Roth, D. A., &amp; Winnewisser, G. 2001, A&amp;A, 370, L49
- Mumma, M. J. &amp; Charnley, S. B. 2011, ARA&amp;A, 49, 471
- Niemann, H. B., Atreya, S. K., Bauer, S. J., et al. 2005, Nature, 438, 779
- Owen, T., Maha ff y, P. R., Niemann, H. B., Atreya, S., &amp; Wong, M. 2001, ApJ, 553, L77
- Persson, M. V., Jørgensen, J. K., &amp; van Dishoeck, E. F. 2013, A &amp;A, 549, L3

Pickett,

H. M., Poynter, R. L., Cohen, E. A., et al. 1998, J. Quant. Spec. Ra- diat. Transf., 60, 883

- Pineda, J. E., Maury, A. J., Fuller, G. A., et al. 2012, A&amp;A, 544, L7
- Rodgers, S. D. &amp; Charnley, S. B. 2008a, ApJ, 689, 1448
- Rodgers, S. D. &amp; Charnley, S. B. 2008b, MNRAS, 385, L48
- Romano, D. &amp; Matteucci, F. 2003, MNRAS, 342, 185
- Rousselot, P., Pirali, O., Jehin, E., et al. 2014, ApJ, 780, L17
- Schöier, F. L., Jørgensen, J. K., Pontoppidan, K. M., &amp; Lundg ren, A. A. 2006, A&amp;A, 454, L67
- Schöier, F. L., Jørgensen, J. K., van Dishoeck, E. F., &amp; Blake , G. A. 2002, A&amp;A, 390, 1001
- Schöier, F. L., van der Tak, F. F. S., van Dishoeck, E. F., &amp; Black, J. H. 2005, A&amp;A, 432, 369
- Shinnaka, Y., Kawakita, H., Kobayashi, H., Nagashima, M., &amp; Boice, D. C. 2014, ApJ, 782, L16
- Takakuwa, S., Ohashi, N., Bourke, T. L., et al. 2007, ApJ, 662, 431
- Tennekes, P. P., Harju, J., Juvela, M., &amp; Tóth, L. V . 2006, A&amp;A, 456, 1037

Terzieva, R. &amp; Herbst, E. 2000, MNRAS, 317, 563

- van Dishoeck, E. F., Blake, G. A., Jansen, D. J., &amp; Groesbeck, T. D. 1995, ApJ, 447, 760
- Vassilev, V., Meledin, D., Lapkin, I., et al. 2008, A&amp;A, 490, 1157
- Wannier, P. G., Linke, R. A., &amp; Penzias, A. A. 1981, ApJ, 247, 522
- Watanabe, Y., Sakai, N., Lindberg, J. E., et al. 2012, ApJ, 745, 126
- Wielen, R. &amp; Wilson, T. L. 1997, A&amp;A, 326, 139
- Wilson, T. L. 1999, Reports on Progress in Physics, 62, 143
- Wilson, T. L. &amp; Rood, R. 1994, ARA&amp;A, 32, 191
- Wirström, E. S., Charnley, S. B., Cordiner, M. A., &amp; Milam, S. N. 2012, ApJ, 757, L11
- Womack, M., Ziurys, L. M., &amp; Wycko ff , S. 1992, ApJ, 387, 417
- Wong, M. H., Atreya, S. K., Maha ff y, P. N., et al. 2013, Geophys. Res. Lett., 40, 6033

Acknowledgements. The authors thank the referee for constructive comments that helped to improve the article. Research for this project was supported by a Danish National Research Foundation Centre of Excellence grant to M.B. (grant number DNRF97) and from a Lundbeck Foundation Junior Group Leader Fellowship to J.K.J. M.B. further acknowledges funding from the European Research Council under ERC Consolidator grant agreement 616027-STARDUST2ASTEROIDS. The research leading to these results has received funding from the European Commission Seventh Framework Programme (FP 2007-2013) under grant agreement No 283393 (RadioNet3). / The authors would like to thank the Onsala APEX sta / ff for carrying out the observations. Per Bergman and the Onsala APEX sta / ff are also thanked for help assessing potential contamination from the o ff position.

## Appendix A: Observational details

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0285v1-figure-231c3785-3efd-46c5-8591-4949159549bf.jpg)

Fig. A.1. Comparison of the H 13 CO (3-2) line from IRAS 16293A ob-+ served with APEX. The spectra resulting from phase 1 (corresponding to o ff position on one side of the source in symmetric wobbler switching mode) and phase 2 (o ff position on the other side) are shown in blue and red, respectively. The spectra have been baseline subtracted and normalized to the average peak intensity of the two spectra. Data from 2012-04-01.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0285v1-figure-7ea21c86-8f6e-4540-b746-66aa4fac6519.jpg)

Fig. A.2. Like Fig A.1 but for H 13 CN(3-2). Includes a sinusoidal baseline for standing wave removal.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0285v1-figure-53288479-84d9-48bc-96ac-87efc8987b20.jpg)

Fig. A.3. Like Fig A.1 but for R CrA IRS7B. Data from 2012-09-27.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0285v1-figure-a06dbcbb-fcc6-4ae4-873a-5ee4a72625d9.jpg)

Fig. A.4. Like Fig A.3 but for H 13 CN(3-2).

Fig. A.5. Like Fig A.1 but for OMC-3 MMS6. Data from 2012-11-22.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0285v1-figure-9fd82699-9d53-4c2f-a494-c55ec9cccf9f.jpg)

Fig. A.6. Like Fig A.5 but for H 13 CN(3-2).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0285v1-figure-886e2204-1cab-4cad-996f-e795c4270fe3.jpg)

Table A.1. Observational details.

| Transition(s)                            | Date         | # scans      | time (on) [min]   | Remarks                                                                          |
|------------------------------------------|--------------|--------------|-------------------|----------------------------------------------------------------------------------|
| IRAS 16293A                              | IRAS 16293A  | IRAS 16293A  | IRAS 16293A       | IRAS 16293A                                                                      |
| H 13 CN(3-2), HC 15 N(3-2)               | 2012-04-01   | 67           | 20                |                                                                                  |
| H 13 CN(4-3), HC 15 N(4-3)               | 2012-06-04   | 41           | 12                |                                                                                  |
| H 13 C 15 N(3-2)                         | 2012-06-13   | 150          | 44                | USB (30 scans), LSB (120 scans) LSB su ff ers from strong spur                   |
| DCN(3-2)                                 | 2012-06-08   | 29           | 9                 | scans 47-48 ignored (LO locked out)                                              |
| DCN(4-3)                                 | 2012-06-07   | 30           | 9                 |                                                                                  |
| DCN(5-4), HNC(4-3)                       | 2012-06-04   | 40           | 12                | baseline ripples                                                                 |
| D 13 CN(3-2)                             | 2012-06-08   | 120          | 35                | scans 541-542 ignored (LO locked out)                                            |
|                                          | 2012-06-09   | 41           | 12                | scans 215-218 ignored (LO locked out)                                            |
|                                          | 2012-06-10   | 50           | 15                |                                                                                  |
|                                          | 2012-06-11   | 60           | 18                |                                                                                  |
| HN 13 C(3-2), H 13 CN(3-2), HC 15 N(3-2) | 2013-08-20   | 36           | 21                | H 13 CN and HC 15 N intensities higher than those observed on 2012-04-01         |
| H 15 NC(3-2), HCN(3-2)                   | 2013-08-19   | 36           | 21                |                                                                                  |
| R CrA IRS 7B                             | R CrA IRS 7B | R CrA IRS 7B | R CrA IRS 7B      | R CrA IRS 7B                                                                     |
| HN 13 C(3-2), H 13 CN(3-2), HC 15 N(3-2) | 2012-08-16   | 30           | 9                 |                                                                                  |
|                                          | 2012-09-27   | 50           | 15                |                                                                                  |
| HN 13 C(4-3), H 13 CN(4-3)               | 2012-08-22   | 50           | 15                |                                                                                  |
|                                          | 2012-09-25   | 30           | 9                 |                                                                                  |
| H 15 N(3-2)C, HCN(3-2)                   | 2012-08-16   | 108          | 32                |                                                                                  |
|                                          | 2012-09-27   | 140          | 41                |                                                                                  |
|                                          | 2012-09-29   | 90           | 26                |                                                                                  |
|                                          | 2012-09-30   | 40           | 12                | intensities &line profiles o ff , data discarded                                 |
|                                          | 2012-11-20   | 30           | 9                 |                                                                                  |
|                                          | 2012-11-21   | 40           | 12                |                                                                                  |
|                                          | 2012-11-23   | 100          | 29                |                                                                                  |
|                                          | 2012-11-24   | 10           | 3                 |                                                                                  |
|                                          | 2012-11-28   | 50           | 15                |                                                                                  |
| HNC(3-2)C                                | 2012-09-27   | 20           | 6                 |                                                                                  |
| OMC3 MMS6                                | OMC3 MMS6    | OMC3 MMS6    | OMC3 MMS6         | OMC3 MMS6                                                                        |
| HN 13 C(3-2), H 13 CN(3-2), HC 15 N(3-2) | 2012-08-16   | 28           | 8                 |                                                                                  |
|                                          | 2012-11-22   | 120          | 35                |                                                                                  |
| HN 13 C(4-3), H 13 CN(4-3)               | 2012-09-26   | 30           | 9                 |                                                                                  |
| H 15 NC(3-2), HCN(3-2)                   | 2012-08-16   | 70           | 21                |                                                                                  |
|                                          | 2012-09-25   | 170          | 50                |                                                                                  |
|                                          | 2012-11-24   | 90           | 26                |                                                                                  |
|                                          | 2012-11-25   | 70           | 21                |                                                                                  |
|                                          | 2012-11-26   | 190          | 56                |                                                                                  |
| HNC(3-2)C                                | 2012-09-27   | 10           | 3                 |                                                                                  |
| Elias 29                                 | Elias 29     | Elias 29     | Elias 29          | Elias 29                                                                         |
| HN 13 C(3-2), H 13 CN(3-2), HC 15 N(3-2) | 2012-08-16   | 30           | 9                 | intensity on the two days is very di ff erent (factor of ∼ 1 . 5), data not used |
| HN 13 C(4-3), H 13 CN(4-3)               | 2012-09-30   | 80           | 24                | intensity on the two days is very di ff erent (factor of ∼ 1 . 5), data not used |
| H 15 N(3-2)C, HCN(3-2)                   | 2012-08-16   | 50           | 15                | intensity on the two days is very di ff erent (factor of ∼ 1 . 5), data not used |
|                                          | 2012-10-01   | 76           | 23                | intensity on the two days is very di ff erent (factor of ∼ 1 . 5), data not used |
| HNC(3-2)C                                | 2012-09-30   | 160          | 47                | intensity on the two days is very di ff erent (factor of ∼ 1 . 5), data not used |
|                                          | 2012-10-01   | 50           | 15                |                                                                                  |

Appendix B: Literature values

Table B.1. Measurements of the 14 N 15 / N ratio in solar system materials from the literature. (B) stands for 'bulk' measurements and (HS) for 'hotspot'. Boldface indicates the unit(s) reported in the original publication (if none of the values is in boldface, the original publication reports 15 N 14 / N). Values marked with a star were used for Fig. 4.

| Object          | 14 N / 15 N                      | δ 15 N [ /permil ]                  | Ref.      | Remarks                                                  |
|-----------------|----------------------------------|-------------------------------------|-----------|----------------------------------------------------------|
| Sun             | 441 ± 6                          | - 383 ± 8                           | a , /star | Genesis                                                  |
| Lunar soils     | ≥ 358                            | - 240                               | b         | trapped solar wind                                       |
| Meteorite (CAI) | 424 ± 3                          | ≤ - 359 ± 6                         | c , /star | Isheyevo, osbornite (TiN)                                |
| Jupiter         | 526                              | 423                                 | d         | ISO, NH 3 , 15 N / 14 N = 1 . 9 + 0 . 9 - 1 . 0 × 10 - 3 |
| Jupiter         | 435 ± 57                         | - - 375 ± 82                        | e , /star | Galileo, NH 3                                            |
| Jupiter         | 448 62                           | 393 84                              | f         | Cassini-CIRS, NH                                         |
| Jupiter         | ± 450 ± 106                      | - ± - 396 ± 142                     | g         | 3 Cassini-CIRS, NH 3                                     |
| Saturn          | ≥ 357                            | - 238                               | h         | ground-based, NH 3                                       |
| Titan           | 183 ± 5                          | ≤ 486 ± 41                          | i         | Huygens, N 2 , secondary origin                          |
| Earth (atmo.)   | 272                              | 0                                   | j , /star |                                                          |
| Venus (atmo.)   | 272 ± 54                         | 0 ± 200                             | k , /star |                                                          |
| Mars (rock)     | 280                              | 30                                  | l , /star | ALH84001 meteorite                                       |
| Mars (atmo.)    | 173 ± 11                         | - 572 ± 82                          | m         | Curiosity, secondary origin                              |
| Comets          | 130 to 220                       | 230 to 1100                         | n , /star | approx. range, CN                                        |
| Comet           | 139 ± 26                         | 957 ± 366                           | o         | 17P / Holmes, HCN                                        |
| Comet           | 165 ± 40                         | 648 ± 400                           | o         | 17P / Holmes, CN                                         |
| Comet           | 139 ± 38                         | 957 ± 535                           | p         | Comet C / 2012 S1 (ISON), NH 3                           |
| Comets          | 127 ( 80 - 190 )                 | 1142                                | q         | avg. spectrum (12 comets), NH 3                          |
| Comet (B)       | 181 to 272                       | 0 to 500                            | r         | Comet 81P / Wild2, Stardust                              |
| Comet (HS)      | 118 to 368                       | - 260 to 1300 ± 400                 | r         | Comet 81P / Wild2, Stardust                              |
| IDPs (B)        | 180 . 7 ± 2 . 3 to 305 . 1 ± 3 . | - 108 ± 9 to 500 ± 20               | s         | IDPs 'Kipling', 'Eliot'                                  |
| IDPs (HS)       | 118 ± 6 to 171 ± 4               | 590 ± 40 to 1300 ± 120              | s         | IDPs 'Kipling', 'Pirandello'                             |
| Cluster IDPs    | 184 to 300                       | - 93 ± 4 to 480 ± 40                | t         | IDPs 'Porky', 'Dragonfly'                                |
| Meteorites (B)  | 192 to 237                       | 146 to 415                          | u , /star |                                                          |
| Meteorites (HS) | 65 to 193                        | 410 ± 130 to 3200 ± 700             | u         |                                                          |
| Meteorites      | 192 to 291                       | - 65 . 6 ± 2 . 9 to 415 . 3 ± 1 . 6 | v , /star | insoluble organic matter (IOM)                           |
| Meteorite (HS)  | 50                               | 4450                                | w , /star | Isheyevo clast                                           |
| Meteorite (HS)  | 46                               | 4900 300                            | x         | Isheyevo clast                                           |

±

References. a Marty et al. (2011), b Hashizume et al. (2000), c Meibom et al. (2007), d Fouchet et al. (2000), e Owen et al. (2001), f Abbas et al. (2004), g Fouchet et al. (2004), h Fletcher et al. (2014), i Niemann et al. (2005) j Junk &amp; Svec (1958); Coplen et al. (1992); Böhlke et al. (2005), k Ho man et al. (1979), ff l Mathew &amp; Marti (2001), m Wong et al. (2013), n Manfroid et al. (2009); Jehin et al. (2009); Mumma &amp; Charnley (2011) and references therein, o Bockelée-Morvan et al. (2008), p Shinnaka et al. (2014), q Rousselot et al. (2014), r McKeegan et al. (2006), s Floss et al. (2006), t Messenger (2000), u Busemann et al. (2006), v Alexander et al. (2007), w Bonal et al. (2010), x Briani et al. (2009)

Table B.2. Integrated intensities from the surveys by Caux et al. (2011) and Blake et al. (1994); van Dishoeck et al. (1995).

| Transition                                                                                                                        | ∫ Tdv [K km / s]   | Telescope   |
|-----------------------------------------------------------------------------------------------------------------------------------|--------------------|-------------|
| H 13 CN (1 2 -0 1 ) H 13 CN (1 1 -0 1 ) H 13 CN (1 0 -0 1 ) H 13 CN (3-2) H 13 CN (4-3) HC 15 N (1-0) HC 15 N (2-1) HC 15 N (3-2) | 1 . 04 ± 0 . 12    | IRAM        |
|                                                                                                                                   | 0 . 75 ± 0 . 09    | IRAM        |
|                                                                                                                                   | 0 . 25 ± 0 . 03    | IRAM        |
|                                                                                                                                   | 8 . 16 ± 1 . 39    | IRAM        |
|                                                                                                                                   | 6 . 90 ± 1 . 24    | JCMT        |
|                                                                                                                                   | 0 . 51 ± 0 . 06    | IRAM        |
|                                                                                                                                   | 1 . 88 ± 0 . 32    | IRAM        |
|                                                                                                                                   | 2 . 57 ± 0 . 45    | IRAM        |
| HC 15 N (4-3)                                                                                                                     | 1 . 49 ± 0 . 27    | JCMT        |
| H 13 CN (3-2)                                                                                                                     | 3 . 45             | JCMT        |
| H 13 CN (4-3)                                                                                                                     | 8 . 13             | CSO         |
| HC 15 N (3-2)                                                                                                                     | 1 . 82             | JCMT        |
| HN 13 C (3-2)                                                                                                                     | 1 . 33             | JCMT        |
S. F. Wampfler J. K. Jorgensen M. Bizzarro S. E. Bisschop
Observations of nitrogen isotope fractionation in deeply embedded protostars
(Abridged) The terrestrial planets, comets, and meteorites are significantly
enriched in 15N compared to the Sun and Jupiter. While the solar and jovian
nitrogen isotope ratio is believed to represent the composition of the
protosolar nebula, a still unidentified process has caused 15N-enrichment in
the solids. Several mechanisms have been proposed to explain the variations,
including chemical fractionation. However, observational results that constrain
the fractionation models are scarce. While there is evidence of 15N-enrichment
in prestellar cores, it is unclear how the signature evolves into the
protostellar phases. Our aim is to measure the 14N/15N ratio around three
nearby, embedded low-to-intermediate-mass protostars. Isotopologues of HCN and
HNC were used to probe the 14N/15N ratio. A selection of H13CN, HC15N, HN13C,
and H15NC transitions was observed with the APEX telescope. The 14N/15N ratios
were derived from the integrated intensities assuming a standard 12C/13C ratio.
The assumption of optically thin emission was verified using radiative transfer
modeling and hyperfine structure fitting. Two sources, IRAS 16293A and R CrA
IRS7B, show 15N-enrichment by a factor of around 1.5-2.5 in both HCN and HNC
with respect to the solar composition. Solar composition cannot be excluded for
the third source, OMC-3 MMS6. Furthermore, there are indications of a trend
toward increasing 14N/15N ratios with increasing outer envelope temperature.
The enhanced 15N abundances in HCN and HNC found in two Class~0 sources
(14N/15N of 160-290) and the tentative trend toward a temperature-dependent
14N/15N ratio are consistent with the chemical fractionation scenario, but
14N/15N ratios from additional tracers are indispensable for testing the
models. Spatially resolved observations are needed to distinguish between
chemical fractionation and isotope-selective photochemistry.