Free Access
Issue
A&A
Volume 519, September 2010
Article Number A30
Number of page(s) 7
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/200913869
Published online 08 September 2010
A&A 519, A30 (2010)

A new weapon for the interstellar complex organic molecule hunt: the minimum energy principle

M. Lattelais1 - F. Pauzat1 - Y. Ellinger1 - C. Ceccarelli2

1 - UPMC Univ. Paris 06, UMR - CNRS 7616, Laboratoire de Chimie Théorique, 75005 Paris, France
2 - Laboratoire d'Astrophysique, Observatoire de Grenoble, 38041 Grenoble, France

Received 14 December 2009 / Accepted 14 May 2010

Abstract
Context. The hunt for the interstellar complex organic molecules (COMs) supposed to be the building blocks of the molecules at the origin of life is a challenging but very expensive task. It starts with laboratory experiments, associated with theoretical calculations, that give the line frequencies and strengths of the relevant molecules to be identified and finishes with observations at the telescopes.
Aims. The present study aims to suggest possible guidelines to optimize this hunt. Levering on the minimum energy principle (MEP) presented in a previous study, we discuss the link between thermodynamic stability and detectability of a number of structures in the important families of amides, sugars and aminonitriles.
Methods. The question of the relative stability of these different species is addressed by means of quantum density functional theory simulations. The hybrid B3LYP formalism was used throughout. All 72 molecules part of this survey were treated on an equal footing. Each structure, fully optimized, was verified to be a stationary point by vibrational analysis.
Results. A comprehensive screening of 72 isomers of CH3NO, C2H5NO, C3H7NO, C2H4O2, C3H6O3 and C2H4N2 chemical formula has been carried out. We found that formamide, acetamide and propanamide (the first two identified in the Inter-Stellar Medium) are the most stable compounds in their families demonstrating at the same time that the peptide bond >N-C=O at the origin of life is the most stable bond that can be formed. Dihydroxyacetone, whose detection awaits for confirmation, is far from being the most stable isomer of its family while aminoacetonitrile, that has been recently identified, is effectively the most stable species.
Conclusions. The MEP appears to be a useful tool for optimizing the hunt for new species by identifying the potentially more abundant isomers of a given chemical formula.

Key words: ISM: molecules - ISM: abundances

1 Introduction

Beyond all philosophical connotations, the pursuit of the origin of life in extraterrestrial environments has never ceased to attract the attention of the scientific community, fueling the imagination of men for generations. Already proposed by Richter in 1865 and further developed by Arrhenius in 1903, the panspermia hypothesis, which claims that the seeds of life are present all over the Universe, was boosted by the discovery of the first organic molecules in the interstellar medium. However, the role of these molecules in the appearance of life on Earth and, more generally, in the Universe is still far from proven. Nonetheless, understanding what is the level of molecular complexity that can be reached in space is certainly a fundamental information and a major open question for modern astro -physics and -chemistry.

Among the most exciting species discovered, one may consider the first aldehyde H2CO (Snyder et al. 1969), the first alcohol CH3OH (Ball et al. 1970), the first conjugated molecule HC3N (Turner 1971), the first acid HCOOH (Zuckerman et al. 1971) and the first alkyl cyanide CH3CN (Solomon et al. 1971). It is also the time when formamide NH2CHO, the simplest model of the peptide bond was identified in Sgr B2 (Rubin et al. 1971).

In the following thirty years, more and more Complex Organic Molecules (hereafter referred to as COMs) have been shown to survive in different astrophysical environments. Among the extra-solar objects, the Galactic Center Molecular Clouds (Requena-Torres et al. 2006; Hollis et al. 2006a), Hot Cores (Ikeda et al. 2001) and Hot Corinos (Bottinelli et al. 2007) are particularly rich in COMs.

It is worth noticing that, of the $\sim$150 species identified up to today, more than 50 are composed of at least 6 atoms and present the main functional groups necessary to initiate a complex organic chemistry. It is indeed very likely, if not sure, that the COMs so far detected only represent the tip of the iceberg, and that many more complex molecules are synthesized in space. However, the detection of new species is a lengthy process, which involves more than just going to the astronomical observatories. In the vast majority, interstellar species are discovered in space by detection of their rotational lines in the radio to sub-millimeter. The more complex the molecule, the larger the number of the lines, caused by the rotation, bending and torsion of the different pieces of the molecule. The first difficulty is indeed to know what lines to observe to identify the target molecule, namely to know its line spectrum. Actually, in addition to know where the lines lie (in frequency), association with the relevant transition is needed, as much as the knowledge of the line strength. This is a very time-consuming and difficult task, which requires specific skills from spectroscopists and theoreticians. So far, in many cases, it is the lack of spectroscopic information which likely limits the search for new species in space. Published unbiased spectral surveys, in fact, often report a large fraction of unidentified lines.

The scope of this article is to suggest guidelines for ``optimizing'' the hunt of new species, by identifying the potentially more abundant isomers of a given chemical formula. Our guidelines are based on the minimum energy principle (MEP), that we presented in a recent paper (Lattelais et al. 2009). Briefly, the analysis of several species observed in different environments (Molecular clouds, Hot Core/Hot Corinos, Photo-Dissociation Regions and AGB stars) has shown that the most abundant isomers of a same chemical formula are the most stable ones, with only two exceptions over a sample of 14 species under 32 different isomeric forms. In addition, the relative abundance ratio between the most stable and the different isomers appears closely correlated to the energy difference between the two isomers. We called this the MEP. The reason for the existence of MEP in space is not obvious, but this seems to be an observational fact that can be used as such.

Following up, we present here computations of the thermodynamic stability of COMs believed to be potentially important for the emergence of life. We considered several important families: amides, sugars and aminonitriles. In the present article we present a systematic survey of the stability of 72 isomers of these families by means of quantum chemical simulations. Considering that the case of glycine, prototype amino acid whose detection has become a graal for the community, deserves special attention, we chose to present it separately (Lattelais et al. in prep.).

The methodology is summarized in Sect. 2, the results for each of the above families of molecules are reported in Sect. 3, the discussion of the consequences of our computations is presented in Sect. 4.

2 Computational background

The argumentation developed in the following sections relies on the relative positions of isomers on the energy scale, so particular attention has to be given to the methods used for their calculation. In an ideal situation, we would have used the highest possible level of post Hartree-Fock theory with extended basis sets. However, in view of the large number of systems (with chemical size up to 12 atoms) to consider, a compromise had to be made between the level of theory and the computational requirements. This is the reason why density functional theory (DFT) was employed using the hybrid B3LYP formalism (Lee et al. 1988; Becke 1993) coupled with the standard 6-311G(d,p) basis set. This level of theory was designed for the best reproduction of the structural properties, energies and spectroscopic parameters (vibrations, dipole moments, ... ) of ``normal'' organic molecules, that is to say, molecules well represented by a unique Lewis structure.

On the other hand, it is known that it is not appropriate for the study of the species that can be described by more than one mesomeric formula like carbon clusters. In this respect, the case of C3H is illustrative with a swap in the energies of linear and cyclic structures between B3LYP and CCSD(T) calculations (see Appendix A below). This level of theory would not be relevant either for the study of weak interactions, like adsorption at the surface of interstellar grains, because of a wrong behavior of the exchange-correlation functional at large atomic separation and the PW91 (or similar functionals) should be preferred in this case (see for example Wesolowski et al. 1997; Wesolowski 2000; Zhang, et al. 1997).

The extension of the basis set is less important for DFT than for configuration interaction treatments; more than 15 years of experience have shown that a correct description of the electronic densities is already obtained with atomic expansions of triple zeta quality (see for example Martin et al. 1996). It is worth noting that this same B3LYP functional, with equivalent basis set, has been recently used for a study of the photophysics of acetic acid and its isomers (Puletti et al. 2010). For a comprehensive presentation of DFT methods and their applicability, the reader may profitably refer to Lee & Scuseria (1995) and Koch & Holthausen (2001).

In the present work, each structure, fully optimized, was verified to be a stationary point by vibrational analysis. All conformers were studied for each molecule but only the most stable one is reported here. Selected structural parameters are given in Appendix B for the most stable isomer of each family (presenting structural details as well as spectroscopic constants and IR signatures for all the calculated structures is beyond the scope of this article and will be addressed in a forthcoming paper). All the energy differences include the zero point vibrational energy (ZPE). We have verified that including entropic terms to determine Gibbs free energies of formation at the temperature usually considered in the ISM has a marginal impact on the values of the energy differences.

The accuracy of the energetic data has been checked with respect to both the extension of the basis set and the treatment of the electronic correlation on the series of 32 isomers considered in our preceding study (Lattelais et al. 2009). (supporting data are given in Appendix A). Increasing the flexibility of the basis set from triple to quadruple zeta has a marginal effect in our B3LYP calculations. Comparison of the B3LYP/6-311(d,p) energy differences with the results obtained at the ``state of the art'' coupled cluster (CCSD(T)/cc-pVQZ level of theory taken as reference, points to an average discrepancy of 1.7 kcal/mol. The ordering of the isomers on the energy scale is maintained for all systems but one: C3H. In this particular and only case, the CCSD(T) level of theory is necessary to assess the greater stability of the cyclic over the linear structure. At all events, the energy gap is small and concerns non-classical chemical structures. The relative energies reported here for systems with standard bonding should be accurate within $\sim$2 kcal/mol.

All the calculations were performed using the GAUSSIAN03 package (Frisch et al. 2004).

3 Molecules with the peptide bond

Among the possible interstellar molecules, those containing the four most abundant atoms vital for life on Earth, namely H, N, O and C, are of major interest. The interest is even larger for molecules with a peptide bond (the group CO-NH), as the large chains of amino acids forming proteins are glued together by this bond. The question here is: what are the most stable molecules, having the potential to form a precious peptide bond? To answer this question, we explore three groups with an increasing number of carbon atoms.

3.1 The CH3NO isomers

A previous study of the stability of the peptide bond (Chiaramello et al. 2005) has shown that formamide NH2CHO, the first step on the protein route, is the lowest energy isomer of the generic formula CH3NO. Up to now, it is the only CH3NO isomer detected in the interstellar medium (Rubin et al. 1971). We also studied the stability of all other possible isomers. A summary of our computations is given in Table 1.

Table 1:   Computed relative stabilities (kcal/mol) and dipole moments (Debye) of formamide isomers.

There are six stable isomers, five of which have a linear backbone of the three heavy atoms, the last one being a cyclic molecule. Formamide is unambiguously the most stable one. The next isomer on the energy scale is HN=CH-OH; it is the enol form of formamide $\sim$13 kcal/mol higher in energy. The third isomer H2C=NOH is formaldoxime. The next two isomers are nitrogen oxides, i.e. nitrosomethane CH3N=O and the corresponding nitrone H2C=NHO. The last isomer is a three-membered ring, more than 70 kcal/mol above formamide.

It should be noted that the isomers with N-C-O connectivity, i.e. those leading to the prebiotic structures are by far the most stable ones. By contrast, connectivities with a NO bond are more than 50 kcal/mol higher than formamide.

The energy difference between formamide and its enol form is similar to that between acetaldehyde CH3CHO and H2C=CHOH. In addition, the dipole moments of the two enols are $\sim$1 Debye. In these conditions, it is reasonable to think that HN=CH-OH should be the next species to be seen in this family with an abundance ratio NH2CHO/HN=CH-OH probably close to CH3CHO/H2C=CHOH, namely between 2 and 10.

3.2 The C2H5NO isomers

A complete study of 13 isomers, including all those of prebiotic interest, has been carried out that gave the energetic ordering reported in Table 2.

The two isomers of lower energy contain the peptide bond. The most stable one, acetamide CH3CONH2, has been recently observed in Sgr B2 (Hollis et al. 2006b) and remains, to our knowledge, the only isomer observed. N-methyl formamide, CH3NHCHO, is next on the energy scale $\sim$10 kcal/mol higher. In third position, one finds the enolic form of acetamide, CH3COH=NH, $\sim$13 kcal/mol above its parent system. These three low energy isomers have the same N-C-O connectivity.

An interesting result concerning this series is that all species containing the NO bond are $\sim$60 kcal/mol above acetamide. According to the MEP, they can hardly be considered as plausible candidates for detection in the Inter-Stellar Medium. In these conditions, N-methyl formamide, with a dipole moment of $\sim$4 Debye, close to that of acetamide, should be the next molecule that contains the peptide bond to be detected in space. If both amides were to be seen, it would establish that the complexity of the peptide bond can progress on both sides (N and C) of the C-N-O linkage. It is also worth noting that amino acetaldehyde NH2CH2CHO, a precursor of glycine NH2CH2COOH, is more than $\sim$25 kcal/mol above the lowest energy structure.

Table 2:   Computed relative stabilities (kcal/mol) and dipole moments (Debye) of acetamide isomers.

3.3 The C3H7NO isomers

Adding one more carbon and two hydrogen atoms is one step further towards complexity. The same systematic study of 25 isomers of prebiotic interest was carried out, leading to the relative energies reported in Table 3.

As already noted for the smaller systems, the isomers of lower energy contain the peptide bond. All three are within $\sim$9 kcal/mol from the most stable one: propanamide. The presence of alkyl substituents on the carbonyl group (>C=O) has a stabilizing effect whereas it is a destabilizing factor when they are attached to the amino group (>NH). It is illustrated by the higher energy of aldehydes compared to ketones. The enolic structures are more than $\sim$15 kcal/mol above their carbonyl parents. As observed previously, the molecules containing the NO bond are $\sim$50-65 kcal/mol above the most stable species. The first cyclic molecule is a five membered ring, but its high energy ($\sim$40 kcal/mol above propanamide) makes it an unlikely candidate for detection in the ISM. It is also true of the three membered rings that are more than $\sim$25 kcal/mol higher, as in the one and two carbon series.

The three amides have dipole moments close to $\sim$4 Debye. The detection of CH3CONHCH3, although not the most stable isomer, would be an important step in complexity since it is the first symmetric peptide.

Table 3:   Computed relative stabilities (kcal/mol) and dipole moments (Debye) of propanamide isomers.

4 Carbohydrates (or sugars)

Carbohydrates, better known under the name of sugars are among the key molecules, together with phosphorus compounds and nucleic bases, for the emergence of a DNA world. Their general formula Cn(H2O)nshows that they are composed of as many molecules of water as the number of carbon atoms they contain. These molecules belong in two types according to the nature of their carbonyl group: either aldoses, with an aldehyde group (-CHO) at the end of the carbon chain, CH2OH-(CHOH)n-CHO, or ketoses with a ketone (R1R2>CO) group inside the carbon chain, CH2OH-(CHOH)n-CO-(CHOH)m-CH2OH. Finding any molecule of the sugar family is of utmost importance for prebiotic models.

4.1 The C2(H2O)2 isomers

The initial report of glycolaldehyde, CH2OHCHO, the simplest possible aldose, (Hollis et al. 2000) was based on radio observations in the millimeter-wave range in the hot molecular core Sgr B2(N-LMH). It is the last isomer detected with this chemical formula. The two other isomers, acetic acid CH3COOH (Mehringer et al. 1997) and methyl formate HCOOCH3 (Brown et al. 1975) had been observed earlier in the ISM.

Applying an energetic criterion to the successive detections shows that HOCH2CHO is effectively the less stable of the three isomers, more than $\sim$25 kcal/mol above acetic acid (Table 4). We discussed the C2(H2O)2isomers in our previous study (Lattelais et al. 2009), and found that this is one (and the clearest) of the two exceptions to the MEP, as the most abundant isomer of the three is the methyl formate rather than the acetic acid. The reason for this is not clear, but it could be linked to a stronger adsorption on the grains of the acetic acid than of the methyl formate, at the moment of their formation on the grain surfaces. Alternatively, the isomerization of these two species may be particularly difficult because it would require breaking two bonds simultaneously. A recent study of the photophysics of these isomers (Puletti et al. 2010) points also to the possible role of the radiation field. From a larger point of view, it will be interesting to understand whether the sugar group in general share the same problem.

Table 4:   Computed relative stabilities (kcal/mol) and dipole moments (Debye) of C2(H2O)2 isomers.

4.2 The C3(H2O)3 isomers

Detection of dihydroxyacetone (HOCH2COCH2OH), DHA, has been claimed by Widicus Weaver & Blake (2005) but soon after was challenged by Apponi et al. (2006).

What does the energetic criterion say in this case? The relative energies of the eight most stable isomers of DHA are given in Table 5. The two structures with the lowest energy are hydroxy carboxylic acids, both derived from propanoic acid; they are very close to each other. Lactic acid HOCHMeCOOH, an important biological compound, is the most stable, $\sim$1 kcal/mol lower than 1-hydroxypropanoic acid HOCH2CH2COOH. In third position, one has the first example of a new family of molecules, namely dimethyl carbonate only $\sim$6.5 kcal/mol above lactic acid. The next compound on the energy scale is the methyl ester of hydroxy acetic acid. The two sugars, dihydroxyacetone and glyceraldehyde come only in sixth and seventh positions, well above the carboxylic isomers. Finally, the cyclic trimer of formaldehyde is $\sim$12 kcal/mol higher than the sugars, that is $\sim$38 kcal/mol above the most stable species.

The basic trends shown by the C2(H2O)2 isomers are confirmed in the present series. The carboxylic acid is more stable than the methyl ester by more than 10 kcal/mol. The two sugars are $\sim$24-26 kcal/mol above the most stable structure as HOCH2CHO is over CH3COOH. (see Table 5).

Based on the MEP, the two hydroxy carboxylic acids should be the most abundant among the series. However, since their energy difference is at the very limit of precision of the present computational experiment, it is safer to consider that they most probably have equivalent abundances. It is also possible that the same phenomenon that makes the methyl formate more abundant than the acetic acid is also at work here, and that, therefore, the most abundant of the series is the dimethyl carbonate.

Table 5:   Computed relative stabilities (kcal/mol) and dipole moments (Debye) of 1,3-dihydroacetone isomers.

5 Amino acids precursors

Amino acids are molecules containing both amine and carboxyl functional groups. The family of the $\alpha$ amino acids (where the amine and the carboxyl groups are attached by a carbon atom), whose generic formula is H2NCHRCOOH with R a generic organic group, is particularly important because these molecules are the building blocks of proteins. Given their supreme importance in terrestrial life, their presence in the interstellar medium has been long searched for, but without success so far, even for glycine, the simplest member of the family. On the other hand, several complex amino acids are found in the carbonaceous chondrites collected on Earth; for example, more than 70 of those species including 8 out of the 20 present in living organisms (Cooper et al. 1992) were isolated in the Murchison meteorite (fall in 1969). It is however worth noticing that those amino acids are thought to have been formed via aqueous alterations on the parent body.

If one considers that aminoacids are end products in the Inter-Stellar Medium, the presence of any of their precursors is equally important.

In a preceding section we have shown that NH2CH2CHO, that would lead to glycine NH2CH2COOH by oxydation is $\sim$26 kcal/mol above the lowest energy structure of same chemical formula. Its abundance should be very small and its role as intermediate in the formation of glycine very unlikely.

On another hand, the Strecker synthesis has already been proposed to explain the formation of these species in meteorites (Bernstein et al. 2004; Peltzer et al. 1984; Ehrenfreund et al. 2001). A major intermediate step in the process is the formation of an aminonitrile that leads to the desired amino acid by hydrolysis. In this spirit, aminoacetonitrile NH2-CH2-CN, has been proposed as a possible glycine precursor about thirty years ago (Brown et al. 1977). It was detected two years ago in the direction of Sgr B2 (Belloche et al. 2008).

Whether this detection may be consistent with an energetic criterion has been addressed by an extensive survey of the possible isomers with C2H4N2 chemical formula. From the results reported in Table 6, it can be seen that the only isomer observed, NH2-CH2-CN, is effectively the most stable compound that can be formed. The species containing the -CN group are more stable than the other isomers. Isonitrile NH2-CH2-NC is $\sim$19 kcal/mol above its nitrile analogue, which is consistent with the energy differences of $\sim$14 and $\sim$24 kcal/mol found for the HCN/HNC and CH3CN/CH3NC couples respectively (Lattelais et al. 2009). Cyclic molecules are once again high in energy, more than $\sim$25 kcal/mol above the minimum.

Table 6:   Computed relative stabilities (kcal/mol) and dipole moments (Debye) of aminoacetonitrile isomers.

6 Discussion

The MEP says that the most stable isomer of a given atomic distribution has the largest abundance, and for that reason should be observed first unless dipole moment restrictions oppose the observations. The case of the amide molecules is particularly illustrative. These species are the most stable compounds that can be obtained from the CnH2n+1NO set of atoms. This result has important consequences for exobiology since it shows that the peptide bond, one of the few bonds at the origin of life, is the most stable that can be formed whenever possible. In these conditions, N-methyl formamide, CH3-NH-CHO, only $\sim$10 kcal/mol above acetamide with a dipole moment $\sim$4 Debye should be the next molecule of this family to search for. With one more carbon in the chain, propanamide appears as the logical candidate; however, detection of the symmetric peptide model, CH3CONHCH3, that ranks second in the series, would have more implications for exobiology.

The problem with the search for sugars is that these structures are not the most stable that can be formed from the same set of atoms, which is most unfortunate in view of a possible detection. In particular, they are much less stable than the isomers containing a carboxylic function (-COOH). These low-lying hydroxyacids are of utmost interest for future observations, especially lactic acid HOCH(CH3)COOH (the most stable compound) whose detection would be proof of the existence of chiral molecules in space (even if the D or L configuration could not be established). The application of the MEP to the search for chiral molecules in the ISM can be found in another report devoted to these particular molecules (Marloie et al. 2010).

Detection of dihydroxyacetone, that is doubtful because of the confusion of the lines with those of other species and the lack of some specific lines, is equally questionable for energetic reasons. One may also speculate on an additional reason that would explain the difficulties for observing sugars, namely the presence of the OH groups that are able to form hydrogen bonds with the grains and maintain a large part of these molecules attached to the surface. Since the OH groups of the hydroxy carboxylic acids present similar inconvenients, dimethyl carbonate CH3OCOOCH3, although slightly higher on the energy scale, is a candidate to consider since it has no acidic hydrogen for sticking on the grains.

Against this pessimistic view, one may argue that HOCH2CHO and HOCH2CH2OH have been observed. In these molecules, there is a strong intramolecular hydrogen bond via a five membered ring structure that decreases the ability to link to an external environment. The same type of internal hydrogen bonding may occur in lactic acid HOCH(CH3)COOH, making this molecule a most appealing target.

Finally, the detection of aminoacetonitrile is an important step that supports the hypothesis of a route to aminoacids via the Strecker synthesis. This route, however, needs also the presence of water and is most probably linked to grain chemistry where H2O is brought by the icy mantles.

Looking beyond the energetic considerations at the end of this study, one has to point out the role of ice which, in addition to its capability to provide water, could be crucial for the detection or non detection of a large number of molecules, some being largely depleted on the icy mantles of the grains, others not. In any case, hot environments and shock regions remain the best places to consider for the observation of complex organic molecules.

All the structural details, spectroscopic constants and IR signatures of the various isomers/conformers will be collected in a database document (Lattelais et al. in prep.).

Appendix A:

Table A.1 lists a series of 14 systems leading to a body of 32 isomers detected in space. The table reports the relative dipole moment of each isomer and stabilities of the different isomers of the same generic formula, computed according to the methods described in Sect. 2. The energy differences between the isomers are very little affected by the level of basis set, and the steady value of the dipole moment is a consistency check of the calculations. The CCSD(T) values, being the best estimates, have been used as reference for evaluation of the average error bars.

Table A.1:   Relative stability (kcal/mol) and dipole moments (Debye) of isomers with same generic formula.

Appendix B:

Figure B.1 shows selected structural parameters of the most stable species presented in the course of the discussion, as obtained after full geometry optimization at the B3LYP/6-311(d,p) level of theory.

\begin{figure}
\par\includegraphics[scale=0.33, angle=0]{Pauz-Figures-eps/figure_valeur.eps}
\end{figure} Figure B.1:

Stable conformations of the most stable species (distances in Ang.; angles in deg.).

Open with DEXTER

Acknowledgements
This work was supported by CNRS national programs PCMI (Physics and Chemistry of the Interstellar Medium), OPV (Origins of Planets and Life) and ANR (Agence Nationale pour la Recherche, contract ANR-08-BLAN-0225). Part of the calculations were performed using HPC resources from GENCI-CINES (Grant 2009-085128).

References

  1. Apponi, A. J., Halfen, D. T., Ziurys, L. M., et al. 2006, ApJ, 643, L29 [NASA ADS] [CrossRef] [Google Scholar]
  2. Ball, J. A., Gottlieb, C. A., Lilley, A. E., & Radford, H. E. 1970, ApJ, 162, L203 [NASA ADS] [CrossRef] [Google Scholar]
  3. Becke, A. D. 2000, J. Chem. Phys., 98, 5648 [NASA ADS] [CrossRef] [Google Scholar]
  4. Belloche, A., Menten, K. M., Comito, C., et al. 2008, A&A, 482, 179 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Bernstein, M. P., Ashnourn, S. F. M., Sandford, S. A., et al. 2004, ApJ, 601, 365 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bottinelli, S., Ceccarelli, C., Williams, J. P., & Lefloch, B. 2007, A&A, 463, 601 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Brown, R. D., Crofts, J. G., Godfrey, P. D., et al., 1975, ApJ, 197, L29 [NASA ADS] [CrossRef] [Google Scholar]
  8. Brown, R. D., Godfrey, P. D., Ottrey, et al., 1977, J. Mol. Spectr, 68, 359 [NASA ADS] [CrossRef] [Google Scholar]
  9. Chiaramello, J. M., Talbi, D., Berthier, G., & Ellinger, Y. 2005, Int. J. Astrobiol., 4, 125 [CrossRef] [Google Scholar]
  10. Cooper, G. W., Onwo, W. M., & Cronin, J. R. 1992, Geochim. Cosmochim. Acta, 54, 4109 [NASA ADS] [CrossRef] [Google Scholar]
  11. Ehrenfreund, P., Glavin, D. P., Botta, O., et al. 2001, PNAS, 98, 2138 [NASA ADS] [CrossRef] [Google Scholar]
  12. Frisch, M., Trucks, G., & Schlegel, H. 2004, Gaussian 03, Rev.C.02, Gaussian Inc., Wallingford CT., 611, 605 [Google Scholar]
  13. Hollis, J. M., Lovas, F. J., & Jewell, P. R. 2000, ApJ, 540, L107 [NASA ADS] [CrossRef] [Google Scholar]
  14. Hollis, J. M., Remijan, A. J., Jewell, P. R., & Lovas, F. J. 2006a, ApJ, 642, 933 [NASA ADS] [CrossRef] [Google Scholar]
  15. Hollis, J. M., Lovas, F. J., Remijan, A. J., et al. 2006b, ApJ, 643, L25 [NASA ADS] [CrossRef] [Google Scholar]
  16. Ikeda, M., Ohishi, M., Nummelin, A., et al. 2001, ApJ, 560, 792 [NASA ADS] [CrossRef] [Google Scholar]
  17. Koch, W., & Holthausen, M. C. 2001, A Chemist Guide to Density Functional Theory, 2nd edn (Weinheim, Germany: Wiley-VCH) [Google Scholar]
  18. Lattelais, M., Pauzat, F., Ellinger, Y., et al.. 2009, ApJ, 696, L133 [NASA ADS] [CrossRef] [Google Scholar]
  19. Lee, C., Yang, W., & Parr, R. G. 1988, Phys. Rev. B, 37, 785 [NASA ADS] [CrossRef] [Google Scholar]
  20. Lee, T. J., & Scuseria, G. E. 1995, In Quantum mechanical electronic structure calculations with chemical accuracy, ed. Langhoff, S. R. (Dordrecht, The Netherlands: Kluwer Academic Publishers) [Google Scholar]
  21. Marloie, G., Lattelais, M., Pauzat, F., Pilmé, J., & Ellinger, Y. 2010, Interdiscip. Sci. Comput. Life Sci., 2, 1 [CrossRef] [Google Scholar]
  22. Martin, J. M. L., El-Yazal, J., & Francois, J. 1996, J. Phys. Chem., 100, 15358 [CrossRef] [Google Scholar]
  23. Mehringer, D. M., Snyder, L. E., Miao, Y., & Lovas, F. J. 1997, ApJ, 480, L71 [NASA ADS] [CrossRef] [Google Scholar]
  24. Peltzer, E. T., Bada, J. L., Schlesinger, G., et al. 1984, Adv. Space Res., 4, 69 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  25. Puletti, F., Malloci, G., Mulas, G., & Cecchi-Pestellini, C. 2010, MNRAS, 402, 1667 [NASA ADS] [CrossRef] [Google Scholar]
  26. Requena-Torres, M. A., Martín-Pintado, J., Rodríguez-Franco, A., et al. 2006, A&A, 455, 971 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Rubin, R. H., Swenson, Jr., G. W., Benson, R. C., Tigelaar, H. L., & Flygare, W. H. 1971, ApJ, 169, L39 [NASA ADS] [CrossRef] [Google Scholar]
  28. Snyder, L. E., Buhl, D., & Palmer, P. 1969, , 22, 679 [Google Scholar]
  29. Solomon, P. M., Jefferts, K. B., Penzias, A. A., & Wilson, R. W. 1971, ApJ, 168, L107 [NASA ADS] [CrossRef] [Google Scholar]
  30. Turner, B. E. 1971, ApJ, 163, 35 [NASA ADS] [CrossRef] [Google Scholar]
  31. Wesolowski, T. A. 2000, J. Chem. Phys., 113, 1666 [NASA ADS] [CrossRef] [Google Scholar]
  32. Wesolowski, T. A., Parisel, O., Ellinger, Y., & Weber, J. 1997, J. Phys. Chem. B, 101, 7818 [Google Scholar]
  33. Widicus Weaver, S. L., & Blake, G. A. 2005, ApJ, 624, L33 [NASA ADS] [CrossRef] [Google Scholar]
  34. Zhang, Y., Pan, W., & Yang, W. 1997, J. Chem. Phys., 107, 7921 [NASA ADS] [CrossRef] [Google Scholar]
  35. Zuckerman, B., Ball, J. A., & Gottlieb, C. A. 1971, ApJ, 163, L41 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1:   Computed relative stabilities (kcal/mol) and dipole moments (Debye) of formamide isomers.

Table 2:   Computed relative stabilities (kcal/mol) and dipole moments (Debye) of acetamide isomers.

Table 3:   Computed relative stabilities (kcal/mol) and dipole moments (Debye) of propanamide isomers.

Table 4:   Computed relative stabilities (kcal/mol) and dipole moments (Debye) of C2(H2O)2 isomers.

Table 5:   Computed relative stabilities (kcal/mol) and dipole moments (Debye) of 1,3-dihydroacetone isomers.

Table 6:   Computed relative stabilities (kcal/mol) and dipole moments (Debye) of aminoacetonitrile isomers.

Table A.1:   Relative stability (kcal/mol) and dipole moments (Debye) of isomers with same generic formula.

All Figures

  \begin{figure}
\par\includegraphics[scale=0.33, angle=0]{Pauz-Figures-eps/figure_valeur.eps}
\end{figure} Figure B.1:

Stable conformations of the most stable species (distances in Ang.; angles in deg.).

Open with DEXTER
In the text

Copyright ESO 2010

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

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

Initial download of the metrics may take a while.