A&A 433, 583-590 (2005)
DOI: 10.1051/0004-6361:20042443
Z. Peeters 1 - O. Botta 1,2 - S. B. Charnley 3 - Z. Kisiel 4 - Y.-J. Kuan 5,6 - P. Ehrenfreund 1
1 - Astrobiology Laboratory, Leiden Institute of Chemistry, PO Box 9502, 2300 RA Leiden, The Netherlands
2 - International Space Science Institute, Hallerstrasse 6, 3012 Bern, Switzerland
3 - Space Science Division, NASA Ames Research Center, MS 245-3, Moffett Field, CA 94035, USA
4 - Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
5 - Department of Earth Sciences, National Taiwan Normal University, 88 Sec. 4 Ting-Chou Road, Taipei 116, Taiwan
6 - Institute of Astronomy and Astrophysics, Academia Sinica, PO box 23-141, Taipei 106, Taiwan
Received 26 November 2004 / Accepted 14 December 2004
Abstract
Nitrogen-containing cyclic organic molecules (N-heterocycles) play important roles in terrestrial biology, for example as the nucleobases in genetic material. It has previously been shown that nucleobases are unlikely to form and survive in interstellar and circumstellar environments. Also, they were found to be unstable against ultraviolet (UV) radiation. However, nucleobases were detected in carbonaceous meteorites, suggesting their formation and survival is possible outside the Earth. In this study, the nucleobase precursor pyrimidine and the related N-heterocycles pyridine and s-triazine were tested for UV stability. All three N-heterocycles were found to photolyse rapidly and their stability decreased with an increasing number of nitrogen atoms in the ring. The laboratory results were extrapolated to astronomically relevant environments. In the diffuse interstellar medium (ISM) these N-heterocycles in the gas phase would be destroyed in
10-100 years, while in the Solar System at 1 AU distance from the Sun their lifetime would not extend beyond several hours. The only environment where small N-heterocycles could survive, is in dense clouds. Pyridine and pyrimidine, but not s-triazine, could survive the average lifetime of such a cloud. The regions of circumstellar envelopes where dust attenuates the UV flux, may provide a source for the detection of N-heterocycles. We conclude that these results have important consequences for the detectability of N-heterocycles in astronomical environments.
Key words: astrobiology - methods: laboratory - stars: individual: CRL 618
In the extended gas and dust envelopes of late-type stars, polymerisation of acetylene can produce benzene, larger polycyclic aromatic hydrocarbons (PAHs) and large carbonaceous dust particles (e.g. Cherchneff et al. 1992). N-heterocycles may form as by-products of acetylene polymerisation. Intermediate replacement of C2H2 by HCN could lead to direct incorporation of a nitrogen atom into a growing aromatic ring (Ricca et al. 2001), forming pyridine. Due to a small HCN/C2H2 ratio in these envelopes, this scenario predicts decreasing abundances of aromatic molecules containing more than one nitrogen atom per ring, such as pyrimidine (two nitrogens) and s-triazine (three nitrogens), see Fig. 1.
Inclusion of nitrogen into the symmetric skeleton of aromatic hydrocarbon compounds induces a permanent dipole moment, which allows the molecule to be observed in a pure rotational spectrum. Pyridine has been the subject of more than 20 published studies on rotational spectroscopy. The millimeter wave rotational spectrum of the parent isotopomer was investigated up to J = 60 and 205 GHz by Wlodarczak et al. (1988) and accurate spectroscopic constants were reported. The rotational spectrum of pyrimidine has recently been thoroughly reinvestigated with the use of several different spectrometers (Kisiel et al. 1999) and precise spectroscopic constants of the ground state and several excited vibrational states, as well as the molecular geometry have been determined. s-Triazine is a symmetric top molecule with zero permanent dipole moment and no results from rotational spectroscopy are available.
Of the four molecules in Fig. 1 only benzene has been detected, by ISO observations in post-asymptotic giant branch (AGB) object CRL 618 (Cernicharo et al. 2001). In parallel with the experimental work reported here and by Peeters et al. (2003), we have undertaken an observational program to search for various N-heterocycles in interstellar and circumstellar environments (Kuan et al. 2004,2003c). Searches for pyridine towards carbon-rich AGBs did not yield any positive results. The upper limits of total molecular column densities that were derived for pyridine, are
cm-2 in IRC+10216 and
cm-2 in CRL 618 (Kuan et al. 2003a,2004). For nucleic acid building block pyrimidine upper limits of
,
and
cm-2 were found in the objects Sgr B2(N), Orion KL and W51 e1/e2, respectively (Kuan et al. 2003c). Obviously, no searches have been conducted for s-triazine.
![]() |
Figure 1: Four six-membered cyclic aromatic molecules: 1. benzene; 2. pyridine; 3. pyrimidine; and 4. 1,3,5-triazine (or s-triazine), having zero, one, two, and three nitrogens per ring, respectively. |
Open with DEXTER |
It is crucial to investigate the photostability of small N-heterocycles in order to understand their role in interstellar and circumstellar gas phase chemistry. Their destruction rates determine the lifetime and abundance of these molecules in interstellar regions and, consequently, their detectability. In the present paper, we measured the UV photostability of the nucleobase precursor molecule pyrimidine, along with the N-heterocycles pyridine and s-triazine, to investigate the effect of the number of nitrogen atoms in the ring on the photostability of these molecules. UV photolysis experiments on benzene have been reported by Ruiterkamp et al. (2005) and the results for this hydrocarbon are compared with the results for the nitrogen-containing heterocycles. In the following sections we will describe the photostability experiments and their results, and discuss the implications of these results for astronomical observations of small N-heterocycles in different space environments.
For the matrix material 99.999% pure argon (Praxair) was used. Pyridine (99.8%, Fluka) and pyrimidine (99%, Aldrich) were prepared by four freeze-pump-thaw cycles to remove dissolved gases. The hygroscopic s-triazine (98%, Merck) was handled exclusively under a nitrogen atmosphere in a glove-box, followed by four freeze-pump-thaw cycles. The vapour of each of these molecules was diluted in argon to a 1:750 ratio, using a glass gas-mixing line with a background pressure of 10-6 mbar. The absolute pressures of the gases were measured with a 0.1-200 mbar and a 1-2000 mbar manometer (Leybold). The resulting N-heterocycle/argon-mixtures were deposited onto the 12 K CsI window to a total thickness of 1
m, controlled by standard laser interferometry. A
1
m thick sample layer allowed for full penetration of the UV light through the sample.
For the photolysis experiments a microwave-excited hydrogen flow lamp with a flux of
photons cm-2 s-1 was used. The UV lamp was calibrated by actinometry as described by Cottin et al. (2003) before and after each series of experiments. UV irradiation was performed in
10 s intervals to a total of 2 min of irradiation (5 s intervals to a total of 1 min irradiation for s-triazine) followed by 1 min intervals until the sample was completely destroyed. Photodestruction of the N-heterocycles was monitored with in situ Fourier-transform infra-red (FT-IR) spectroscopy using an Excalibur FTS-4000 spectrometer (BioRad) in the range 4000-500 cm-1 at 1 cm-1 resolution. In these experiments, the combination of an accurate gas mixing ratio, a laser-interferometry controlled sample thickness, and a calibrated lamp flux allowed for accurate and reproducible photodestruction rate measurements, which will be shown in the next section.
![]() |
Figure 2:
The FT-IR spectra of pyridine, pyrimidine and s-triazine
( top to bottom) isolated in an argon matrix at 12 K in the ranges 3500-3000 cm-1
( left panels) and 1800-600 cm-1 ( right panels) with a resolution of 1 cm-1. The matrix isolated molecules were irradiated with UV for 10 s, 1 and 10 min. The vertical scale bars indicate the infra-red absorption in absorption units. Band assignments are presented in Table 1. New bands appearing after photolysis are designated by asterisks (![]() |
Open with DEXTER |
Table 1:
The wavenumbers (
in cm-1) of the most prominent peaks of the heterocycles with their assignments, used to monitor the UV photodestruction.
denotes a stretching mode,
is a bending mode and
refers to a torsion mode. Taken from a Destexhe et al. (1994) and b Morrison et al. (1997).
The UV photolysis of the heterocyclic molecules was measured by recording the infra-red spectrum after short periods of exposure to UV radiation. In Fig. 2 the spectra after 10 s and 1 and 10 min of photolysis are shown for each N-heterocyclic molecule. The diminishing of the peaks in Fig. 2 after UV irradiation shows that all three molecules were destroyed by the UV. The rate of destruction was calculated by plotting the natural log of the normalised integrated peak areas against time and fitting the data points with a linear regression line. Only the first 120 s of photolysis time (the first 35 s for s-triazine) were considered, because at longer photolysis times the destruction rate showed a deviation from first order behaviour.
From the slopes of the linear fits, the destruction cross-sections and half-lives were calculated, which are compiled into Table 2. Half-lives were calculated from the laboratory data and extrapolated to the diffuse ISM, a dense cloud and the Solar System environment at 1 AU from the Sun, using the appropriate UV fluxes:
photons cm-2 s-1 for the diffuse ISM (Mathis et al. 1983),
photons cm-2 s-1 for a dense cloud (Prasad & Tarafdar 1983) and
photons cm-2 s-1 for the solar UV flux >6 eV at 1 AU distance from the Sun. The resulting half-lives for those astronomical environments were added to Table 2. Finally, the half-lives for the N-heterocycles were compared to benzene (Ruiterkamp et al. 2005) and plotted against the number of nitrogen atoms in the ring, resulting in Fig. 3.
Table 2:
The UV destruction cross-sections (
)
and half-lives in the experimental set-up in the laboratory and in the diffuse interstellar environment (DISM), dense clouds (DC) and the Solar System at 1 AU from the Sun (SoSy). Also included in this table is the cyclic hydrogen-bonded trimer of HCN, which is the first photoproduct of s-triazine. Upon continued irradiation, c-(HCN)3 is also photolysed, see Fig. 4.
The photodissociation of s-triazine in an argon matrix has been previously described by Satoshi et al. (1997). New bands appearing after 10 s of UV irradiation,
![]() |
Figure 3: The half-lives of the three N-heterocycles plotted against the number of nitrogen atoms in the ring. Also included in this graph is the half-life of benzene, taken from Ruiterkamp et al. (2005). This plot shows that the photostability of small heterocyclic molecules decreases when an increasing number of nitrogen atoms is incorporated in the ring. The error bars indicate the standard deviation of the linear fit used in the calculation of the half-lives of those molecules. |
Open with DEXTER |
Table 3: Wavenumbers (in cm-1) of the photoproducts that appeared after UV irradiation of pyridine, pyrimidine or s-triazine. The band assignments were taken from a Milligan & Jacox (1973); b Ruiterkamp et al. (2005); c Satoshi et al. (1997); and d Stroh & Winnewisser (1989). The bands for HCN and C2H2 in the pyrimidine photoproducts are shifted with respect to the bands given in the references, so a complex between HCN and C2H2 is inferred (see text).
The change of the relative abundances of s-triazine and its photoproducts during irradiation - as measured by the integrated peak areas, relative to the highest value per peak - is shown in
Fig. 4. The relative abundances in Fig. 4
were determined from the peaks at 1552 cm-1 for s-triazine, 3264 cm-1 for c-(HCN)3 and 2295 cm-1 for CNCN. From the results described above, the following s-triazine photolysis pathway is suggested:
,
see also the inset in
Fig. 4. After 30 s of irradiation only
of the initial amount of s-triazine was remaining. Assuming the suggested photolysis pathway, 90% of the s-triazine was converted into c-(HCN)3 in the first 30 s. After this period c-(HCN)3 degraded with only little replenishment by further degradation of s-triazine. From the c-(HCN)3 peaks in the infra-red spectra between 30 and 100 s of irradiation, the UV destruction cross-section and half-lives for c-(HCN)3 could be calculated in the same way as for pyridine, pyrimidine and s-triazine. The results were added to Table 2.
![]() |
Figure 4:
The relative column densities of s-triazine and its photoproducts during photolysis. Shown in the graph are s-triazine (graph: ![]() ![]() ![]() |
Open with DEXTER |
In the photolysis experiments of pyrimidine only two new bands were found after 10 min of irradiation with UV. These bands were located at 3297.8 and 3314.7 cm-1. Since s-triazine has been shown to photolyse into three HCN molecules (Satoshi et al. 1997) and benzene into acetylene molecules (Ruiterkamp et al. 2005), it can be expected that pyrimidine would degrade into two HCN molecules and one acetylene molecule. These photoproducts are likely to form a complex in some way, since they are enclosed by an argon matrix. The peaks for the HCN monomer in an argon matrix are found at 3306 and 3202 cm-1, while for a linear dimer of HCN these peaks shift to 3304 and 3210 cm-1 (Satoshi et al. 1997). For acetylene isolated in an argon matrix a peak is found at 3253 cm-1 while a peak at 3284 cm-1 might be assigned to a complex of multiple acetylene molecules (Ruiterkamp et al. 2005). Using these data, the new peak at 3297.8 cm-1 is tentatively assigned to an acetylene mode and the peak at 3314.7 cm-1 to an HCN mode in a complex between acetylene and HCN.
The only new band appearing in the infra-red spectra of pyridine after UV photolysis, is found at 904.1 cm-1. This band is found in every experiment, including pyrimidine and s-triazine and is assigned to a proton-argon complex, [HAr2]+ (Milligan & Jacox 1973). All photoproduct bands and their assignments are listed in Table 3.
Upon irradiation of N-heterocycles with UV, these molecules are destroyed as can be seen by the extinction of their peaks in the infra-red spectra in Fig. 2. The UV destruction cross-section and half-lives in Table 2 show that this destruction proceeds rapidly under laboratory conditions. When these laboratory results are extrapolated to astrophysical conditions, pyridine, pyrimidine and s-triazine are expected to have very short lifetimes in the diffuse interstellar and Solar System environments. Only in dense interstellar clouds, where these molecules are shielded from the interstellar UV field, would pyridine and pyrimidine be able to survive the average life-span of such a cloud (106 years, Hartmann et al. 2001; Elmegreen 2000). s-Triazine, with a half-life of 250 000 years in a dense cloud environment, would have only 6% of its initial amount remaining at the end of the cloud's lifetime. The c-(HCN)3 complex that forms upon photolysis of s-triazine, however, has a half-life that is comparable to the half-life of pyrimidine and might survive the lifetime of a dense cloud.
Around 20 molecules have been detected in CRL 618. According to models by Herpin & Cernicharo (2000) and Woods et al. (2002) benzene and, possibly, N-heterocycles are formed in the circumstellar envelope of CRL 618. The physical structure of this circumstellar envelope is very complex. The strong UV field from the central star has produced a compact H II region that is confined by a dense circumstellar disk (Cernicharo 2004). The benzene absorption arises from the innermost regions of the dense circumstellar disk, where dust absorption attenuates the strong stellar UV flux. Models reported by Woods et al. (2002) indicate that when the extinction decreases to 10 mag at a distance of
1016 cm from the central star, the abundances of molecules such as benzene rapidly drop, indicating destruction by interstellar UV photons. In such environments N-heterocycles would be rapidly destroyed as well.
The experiments described in this paper are performed employing the low temperature matrix isolation technique, while the results are extrapolated to the gas phase in different interstellar and circumstellar environments. Although the argon matrix provides a good spectroscopic comparison to the gas phase, there are some small differences between argon matrix and gas phase in terms of photochemistry. In the matrix, energy absorbed by the molecule can be dissipated in the matrix, thereby lowering the chance of photodestruction. Also, after photolysis photoproducts remain in close contact with each other in the matrix, because they are enclosed by solid argon. These two factors will increase the real gas phase destruction rate (lower the gas phase half-life) and the matrix isolated destruction cross-sections presented in this paper can be regarded as lower limits (upper limits for the half-lives). Ruiterkamp et al. (2005) showed a factor 10 difference between gas phase and low temperature matrix isolation data for benzene.
Woods et al. (2002) have shown that circumstellar ion-molecule chemistry, based on the findings of McEwan et al. (1999), could in fact produce benzene in the abundance observed in CRL 618. In this scenario, the reaction sequence is initiated by HCO+ reacting with C2H2. A very high ionization rate (about 100 times the canonical interstellar value due to cosmic ray particles) is required to produce c-C6H6 efficiently. Woods et al. (2002) attributed soft X-rays or shock waves as the source of enhanced ionization. In situ erosion of carbon dust to PAHs and simpler molecules, either by sputtering in shocks or by photolysis (Wooden et al. 2004), could provide a viable source of benzene to compete with its photodestruction in the diffuse ISM.
N-heterocycles are not formed as efficiently as benzene in circumstellar environments, as was discussed in Sect. 1. A low formation rate combined with a low stability against UV radiation, makes the detection of N-heterocycles in circumstellar environments difficult. Nevertheless, sources with higher HCN/C2H2 ratios in the dust formation zone could have higher production rates of N-heterocycles and thus serve as better targets than CRL 618 to search for these molecules. For example, the HCN/C2H2 ratio in CRL 618 is observed to be about 0.16 (Woods et al. 2002), whereas the estimated HCN/C2H2 ratio around the carbon-rich AGB star IRC+10216 is about 0.6 (Cherchneff et al. 1993). This elevated HCN/C2H2 ratio, along with a much slower outflow velocity and lower UV flux, suggests that the circumstellar envelope of IRC+10216 would be a better source for finding N-heterocycles.
The high photodestruction rates make it unlikely that small N-heterocycles could have survived the conditions during Solar System formation, unless they were well protected against radiation in the disk mid-plane at all times.
It is generally accepted that the amino acids found in carbonaceous meteorites, are produced during the phase of aqueous alteration in the parent body. These compounds can be synthesized from interstellar precursor molecules such as HCN, ammonia and carbonyl compounds, including formaldehyde, acetaldehyde and acetone. This formation route, the so-called Strecker-cyanohydrin reaction, can take place inside the meteorite parent body during aqueous alteration (Cronin & Chang 1993). These conditions could also be favourable for the formation of purines and pyrimidines.
The first synthetic step towards the abiotic synthesis of purines could be HCN polymerisation, which would lead to the formation of he HCN-tetramer diaminomaleonitrile (DAMN) as an intermediate. This compound, which is obtained in yields of 0.3 to 0.5% in aqueous 0.1-1.0 M HCN solutions, is also a possible intermediate for the formation of certain amino acids (Ferris & Hagan 1984). The formation of DAMN requires basic pH values (9.2 in the laboratory experiment), while the liquid inside the meteorite parent bodies is thought to be buffered to pH 6 (slightly acidic) by the mineral composition present in the meteorite (DuFresne & Anders 1962).
Several possible pathways exist for the reaction of DAMN with other compounds such as ammonia and/or formamide to lead to the formation of adenine, the simplest purine. The reaction of DAMN with urea leads to the formation of guanine and xanthine in yields of 5-10%. Another prebiotic pathway for the formation of the adenine derivative 8-hydroxymethyladenine (HMA) as well as traces of adenine from HCN and formaldehyde, was experimentally demonstrated by Schwartz & Bakker (1989). Finally, it has been demonstrated recently that guanine can be produced from ammonium cyanide polymerisations (Levy et al. 1999).
The abiotic synthesis of pyrimidines is more difficult to achieve. The first proposal for the formation of a pyrimidine from a cyano-compound was the synthesis of cytosine from cyanoacetylene and cyanate, both possibly prebiotic compounds, followed by hydrolysis, leading to uracil (Ferris et al. 1968). Starting from DAMN, a variety of pathways for the formation of pyrimidine derivatives, for example the reaction of DAMN with guanidine that yields 5-hydroxyuracil, were proposed by Ferris et al. (1978). Only uracil was detected in carbonaceous chondrites (Stoks & Schwartz 1979), a fact that may reflect the apparent difficulty in the abiotic synthesis of these compounds.
In summary, carbonaceous chondrites appear to contain several classes of N-heterocyclic compounds, including purines, pyrimidines, quinolines/isoquinolines and pyridines. Stable isotope data of these compounds, which would be very helpful in determining their extraterrestrial origin, have yet to be reported.
Acknowledgements
The authors wish to thank T. Millar, P. Woods and E. Herbst for discussions. P.E. and Z.P. are supported by grant NWO-VI 016.023.003, Y.J.K. is supported by grant NSC 93-2112-M-003, and S.B.C. is supported by NASA's Exobiology program through NASA Ames Cooperative Agreement No. NCC2-1412 to the SETI Institute. O.B. was an ESA post-doctoral fellow.