Abstract
Introduction:
The 21 m emission feature discovered in a small sample of carbon-rich protoplanetary nebulae has remained unidentified for over 30 years. A dozen of different molecular species (both organics and inorganics) have been proposed as the carrier candidates of this important feature, among which polycyclic aromatic hydrocarbons (PAHs) and fulleranes (hydrogenated fullerenes) have yet to be sufficiently examined.
Method:
In this study, we attempt to fit the 21 m features in observed spectra of the above-mentioned astronomical sources via theoretically simulated spectra of various PAHs and fulleranes, aiming to investigate whether the two hydrocarbon families can reproduce the 21 m feature.
Results and Discussion:
Based on the fitting outcomes we conclude that fulleranes can provide a more plausible explanation for the origin of 21 m feature than PAHs.
1 Introduction
As a low- or intermediate-mass star progresses into the final stage of the asymptotic giant branch (AGB), its circumstellar envelope is detached from the stellar surface, and the effective temperature rises (e.g., Kwok, 1993; 2024). When the center star is sufficiently hot, the envelope is ionized, initiating the formation of a planetary nebula (PN). During the transition from the AGB to PN phases, there is a brief evolutionary phase ( yr), called protoplanetary nebula (PPN, see, e.g., Volk Kwok, 1989). Observing PPNe is challenging because their main radiation bands are in the infrared wavelength region and almost without emission lines in the optical region. Infrared instruments such as the Infrared Astronomical Satellite (IRAS, Neugebauer et al., 1984), the Infrared Space Observatory (ISO, Van Winckel, 2003), and the Spitzer Space Telescope (Gehrz et al., 2007) have significantly enhanced our ability to observe and study PPNe.
Despite evolving from an AGB star’s envelope, a PPN’s infrared spectrum exhibits substantial variations with the appearance of a group of Unidentified Infrared Emission (UIE) bands at 3.3, 3.4, 6.2, 6.9, 7.7, 8.6, 11.3, and 12.7 m (Duley and Williams, 1981; Buss Jr et al., 1990). suggesting the formation of complex hydrocarbons with aromatic and aliphatic functional groups during the PPN stage (Kwok, 2004). The chemical processes that lead to the production of these organic molecules in PPN have yet to be uncovered.
In addition, an UIE feature at 21 m was first discovered by Kwok et al. (1989) in four carbon-rich PPNe from the IRAS database. Over the subsequent three decades, researchers have detected a total of 31 sources exhibiting this feature from the IRAS, ISO, and Spitzer Space Telescope databases (García-Lario et al., 1999; Hrivnak et al., 2000; 2009; Cerrigone et al., 2011; Volk et al., 2011; Matsuura et al., 2014; Gładkowski et al., 2019). Among them, 20 are within the Milky Way galaxy, nine are within the Large Magellanic Cloud (LMC), and two are within the Small Magellanic Cloud (SMC). The 21 m emission feature has an asymmetric profile that expands over the wavelength range of 17.4–23.4 m with a maximum flux centered around 20.1 m. The feature is characterized by a rapid rise in the short and gradually declining at the long wavelength side (Volk et al., 2020). The 21 m feature is mostly accompanied by a broad emission at 30 m and a plateau emission at 11–17 m. The exclusive presence of the 21 m feature in PPNe may hold a vital clue in understanding PPN chemistry. However, the identification of its carrier remains highly controversial. So far different carrier candidates have been proposed (see, Volk et al., 2020, for a list). Considering that the carrier must meet two criteria: 1) the involved elements must be abundant enough to account for the intense 21 m emission; 2) other emission bands produced by the substance must resemble the observed spectrum, Zhang et al. (2009) examined nine Si-, Fe-, and Ti-bearing carrier candidates, and found that most of them except FeO nanoparticles do not satisfy the two criteria. Subsequently, Li et al. (2013) found evidence that FeO cannot be responsible for the 21 m feature.
Apart from inorganic compounds, complex organic molecules such as polycyclic aromatic hydrocarbons (PAHs) and hydrogenated fullerenes (fulleranes), could also emit at around 21 m. Distinguishing from a single inorganic species, PAHs and fulleranes are presenting hydrocarbon families with specific molecular structures. PAHs contain a large cyclic conjugate structure, allowing them to be stable in interstellar space, and are being considered as possible carriers of some groups of UIEs (3.3, 6.2, 7.7, 8.6, 11.3, and 12.7 m) (Tielens, 2008; Duley and Williams, 1981). Papoular (2011) found that PAH-like molecules can reproduce the 21 m feature. Fullerene () was first discovered in the PN Tc1 (Cami et al., 2010) and subsequently observed in different circumstellar environments including the PPN (García-Hernández et al., 2010; 2011a; b; Sellgren et al., 2010; Gielen et al., 2011; Zhang and Kwok, 2011; Evans et al., 2012; Roberts et al., 2012; Otsuka et al., 2013). In laboratory conditions, can be quickly hydrogenated into by atomic hydrogen (Cataldo and Iglesias-Groth, 2009) and would undergo rapid dehydrogenation when heated to around 550 (Rüchardt et al., 1993). So far there is no compelling evidence for the presence of circumstellar fulleranes, although tentative detection has been anticipated (Zhang et al., 2017; Palotás et al., 2020). Fullerene cage structure maintains in harsh astronomical environments (Sadjadi and Parker, 2021) to form fulleranes and emit UIEs. The possibility of fulleranes as the carrier of the 21 m feature was raised by Webster (1995) and revisited by Zhang et al. (2020).
In comparison to other astrochemically relevant species, PAHs and fulleranes have certain advantages to account for the 21 m emission. The cosmic abundances of carbon and hydrogen are high. The 21 m features are revealed only in carbon-rich sources, suggesting its carrier is more likely to be carbon-bearing. Moreover, considering the different sizes, structures, charge states, and impurities, the numbers of PAHs and fulleranes could be large, providing a very flexible way to fit the observed spectrum. However, their likehoods as the 21 m feature carrier have not been sufficiently evaluated. To this end, we perform spectral fittings of the observations utilizing the theoretical spectra of PAHs and fulleranes, aiming to investigate how well the two hydrocarbon families reproduce the 21 m feature.
The paper is structured as follows: Section 2 presents the methodology and the data used. In Section 3 we discuss the goodness of different materials as the carrier of the 21m emission feature. Section 4 presents our conclusions.
2 Methodology
2.1 Observational spectra
The Spitzer Space Telescope observed 11 sources exhibiting the 21 m emission feature between 2004 and 2008 as part of the programs No. 20208 (PI: B. Hrivnak) and 93 (PI: D. Cruikshank), including IRAS 04296 + 3,429, IRAS 05113 + 1,347, IRAS 05341 + 0,852, IRAS 06530-0,213, IRAS 07134 + 1,005, IRAS 07430 + 1,115, IRAS 19477 + 2,401, IRAS 20000 + 3,239, IRAS 22223 + 4,327, IRAS 22574 + 6,609, and IRAS 23304 + 6,147. Among them, IRAS 07134 + 1,005 and IRAS 20000 + 3,239 exhibit exceptionally strong 21 m emissions, which are saturated in the observations. The Spitzer spectra of the other nine PPNe were extracted from the Spitzer archive for our analysis (see Zhang et al., 2010, for the details of the spectra and the data processing), covering a wavelength range from 9.9 to 37.2 m. IRAS 04296 + 3,429 and IRAS 22223 + 4,327 were also observed using the short-low (SL) module in programs 30036 (PI: G. Fazio) and 45 (PI: T. Roellig), thus their spectra have a larger wavelength coverage (5.5–37.2 m.)
To subtract the continuum underlying the features, we selected the concave points in the spectra as anchors. A spline interpolation was utilized to construct a curve that passes through all selected anchors. Then constructed curve was subtracted from the observed spectrum for the subsequent fitting. The continuum-subtracted spectra are shown in Figure 1.
FIGURE 1
2.2 Theoretical spectra of PAHs
Theoretical Infrared (IR) spectra of PAH molecules calculated at density functional theory (DFT) (Parr, 1985; Mattioda et al., 2020) were obtained from the NASA Ames PAH IR Spectroscopic Database1 (PAHdb Bauschlicher et al., 2010; 2018). Among the 4000 PAH spectra in PAHdb, we chose 288 for the fitting based on the following criteria: 1) the molecules contain C, H, and N atoms only; 2) the molecules include aromatic C-H bonds that are responsible for UIE; (3) the numbers of selected small-, medium-, and large-sized PAHs (with carbon-atom number , 50–100, and ) should be roughly the same to avoid a bias toward the PAHs with a specific size; 4) for each selected PAH, the spectrum of its cation should be available. The molecular formula and UID of the selected PAHs are presented in Supplementary Appendix S1 of the Supplementary Material. Following previous studies (Rosenberg et al., 2014; Zhang and Kwok, 2015), we performed a line broadening by convoluting with a Gaussian profile of a width of 15 . Using the online tools provided by PAHdb, the normalized fitting spectra were finally deduced under an assumption of thermal excitation with temperatures ranging from 100 to 150 K.
2.3 Theoretical spectra of fulleranes
Theoretically calculated IR spectra of 55 fulleranes (, 2, 4, 6, 8, 10, 12, 14, 16, 18, 20, and 36) were taken from our previous works (Zhang et al., 2017).Two DFT hybrid functionals B3LYP and BH&HLYP in combination with PC1 basis set (Jensen, 2001; 2002) have been applied to the computations. Gaussian09 quantum chemistry package (Frisch, 2009) was used for this purpose. Five structural isomers have been considered for each fulleranes. Fulleranes with an odd number of hydrogen atoms are thermodynamically less favorable than those with an even number and thus were not considered in this study (Kabo et al., 2010). The line broadening was performed by convoluting a Gaussian profile with a width of 15 . To simulate the thermal excitation model, we multiplied the theoretical spectra by a temperature-dependent Planck function. All of the assumptions for constructing fullerane spectra are the same as those for PAH spectra, minimizing potential biases in comparison of their fitting results.
2.4 Spectra fitting
The nine continuum-subtracted spectra were fitted by synthesizing separately the theoretical spectra of PAHs and those of fulleranes using the Markov chain Monte Carlo algorithm. The steps involved were as follows.
1. Import the observational spectrum and the spectrum data matrix of PAHs/fulleranes.
2. Generate a series of random numerical sequences, each containing the same number of digits as the number of molecules in the spectrum data matrix.
3. Multiplying the spectrum data matrix with the numerical sequence, we get a series of synthesized spectra.
4. Compare the observed spectrum and the synthesized spectra, and pick up that with the minimum error.
5. Repeat the above steps until no synthesized spectrum with lower error can be found.
6. Output the synthesized spectrum and the numerical sequence given by the optimal fitting.
The goodness of the fits is evaluated quantitatively by the reduced chi-square , with a smaller indicating a better fit. The calculations of were mainly based on the spectra within the wavelength range from 17.4 to 23.4 m. However, multiple spectral features could be seen outside this wavelength coverage, which might originate from different species. At these wavelength ranges (m and m), the errors were accounted in the calculations only when the synthesized spectrum was more intense than the observed one.
After the optimal fitting was obtained, we can investigate the types of PAHs and fulleranes that are mostly responsible for the 21 m feature. For that purpose, the PAHs were classified according to their sizes, charge states, number of nitrogen atoms, and C/H ratios; the fulleranes were classified according to their hydrogenation degrees.
3 Result and discussion
As thermal excitation has been assumed, the synthesized spectra depend on the preset temperatures. If the preset temperature was too high, in the short-wavelength regions, the synthesized spectra would exhibit too intense features to be compatible with the observations, and vice versa. To optimize the temperature adopted, we performed the fitting using a few different temperature values for PAH and fullerane spectra, as shown in Figures 2, 3 respectively. As shown in the figures, when the preset temperature is 150 K, the synthesized spectra show too strong emission bands near 12.7 m. When it is 100 K, the band around 30 m appears to be too intense. Therefore, we adopted a temperature of 125 K for the fitting.
FIGURE 2
FIGURE 3
The values of the fits are listed in Table 1. It is shown that adopting a temperature of 125 K gives the optimal choice for the PAH spectral fitting. Notably, (fullerane) are generally much smaller than (PAH), suggesting that fulleranes are more plausible to reproduce the 21 m feature than PAHs. Therefore, in the following, we will focus on the discussion of the fitting results of fulleranes. For the PAH spectral fitting, the contributions from different PAH groups are illustrated by Supplementary Appendix S2 in the Supplementary Material.
TABLE 1
| Source | Temperature | (PAH) | (fullerane) |
|---|---|---|---|
| IRAS 06530–0213 | 100 | 6.35 | 0.407 |
| 125 | 2.69 | 0.298 | |
| 150 | 2.89 | 0.263 | |
| IRAS 05341 + 0852 | 100 | 0.753 | 0.0836 |
| 125 | 0.361 | 0.0851 | |
| 150 | 1.91 | 0.0880 | |
| IRAS 07430 + 1115 | 100 | 1.15 | 0.236 |
| 125 | 0.272 | 0.244 | |
| 150 | 0.293 | 0.249 | |
| IRAS 05113 + 1347 | 100 | 0.377 | 0.125 |
| 125 | 0.242 | 0.221 | |
| 150 | 0.401 | 0.262 | |
| IRAS 19477 + 2401 | 100 | 1.35 | 0.0621 |
| 125 | 0.840 | 0.105 | |
| 150 | 3.92 | 0.502 | |
| IRAS 04296 + 3429 | 100 | 0.255 | 0.280 |
| 125 | 0.297 | 0.260 | |
| 150 | 0.394 | 0.326 | |
| IRAS 22223 + 4327 | 100 | 0.187 | 0.0893 |
| 125 | 0.104 | 0.114 | |
| 150 | 0.229 | 0.140 | |
| IRAS 23304 + 6147 | 100 | 4.48 | 0.758 |
| 125 | 3.54 | 0.758 | |
| 150 | 3.92 | 0.835 | |
| IRAS 22574 + 6609 | 100 | 1.72 | 0.121 |
| 125 | 0.697 | 0.133 | |
| 150 | 0.283 | 0.142 |
of the fits.
The fulleranes are divided into two groups according to their hydrogenation degree ( with and ). Figure 4 shows the contributions of the two groups to the best fits (Figure 3). It is clear that slightly hydrogenated cannot contribute to the emission band around 21 m, which are mostly from the fulleranes with moderate hydrogen content (–36). This supports the finding of Zhang (2020) that moderately hydrogenated are a promising carrier material producing the 21 m feature, and is consistent with the experimental results (Cataldo, 2003; Iglesias-Groth et al., 2012).
FIGURE 4
can be formed in the PPN phase (Zhang and Kwok, 2011). However, it seems that the and 21 m PPNe are mutually exclusive. No 21 m source exhibits the bands. The PPN has a hotter central star than the 21 m PPNe, and thus is more evolved (Zhang et al., 2010). The deep spectroscopy of two PNe exhibiting strong bands do not detect the C-H stretching bands of fulleranes around 3.4–3.6 m (Díaz-Luis et al., 2016). Theoretically, The emission bands from carbon-cage vibrations are still visible for slightly hydrogenated , and gradually fade with increasing hydrogenation (Zhang et al., 2017). Consequently, a reasonable hypothesis is that the moderately hydrogenated responsible for the 21 m the feature could be formed in the early PPN stage, and then, with further increasing temperature of the central, is rapidly dehydrogenated by the intense ultraviolet radiation. In the PN stage, hydrogen atoms are completely removed from surface so that no C-H band could be observed.
The binding energy of hydrogen atoms linking on fullerene surface is 3.3 eV and 1.9 eV for fulleranes with even and odd hydrogen-atom numbers, respectively (Abbink et al., 2024). The values are lower than that of ordinary H-H bonds (4.5 eV). Therefore, cannot be effectively hydrogenated if hydrogen is mainly in molecular state. Glassgold and Huggins (1983) found that the ejected hydrogen from an AGB star remains in molecular state until the stellar temperature exceeds 2500 K. Therefore, atomic hydrogen dominates the stellar wind when the AGB envelope enter the PPN stage, and efficiently hydrogenates fullerenes. As a result of the dramatic changes of the physical conditions of PPNe, fulleranes are rapidly formed and then are rapidly dehydrogenated. If fullerenes are heavily hydrogenated, the carbon cage is unstable, and may be destructed. These can explain the transient nature of the 21 m feature.
The C-H streching vibration of moderately hydrogenated fullerenes may provide observable emission features around 3.4 m (Iglesias-Groth et al., 2012). Unfortunately, this feature lies outside the wavelength coverage of Spitzer spectra. Previous observations have revealed that the 21 m sources are enriched with aliphatic features (Kwok et al., 2001), which may partly origin from fulleranes. It is highly desirable to investigate the correlation between the intensities of the 3.4 m and 21 m features. The James Webb Space Telescope could add light on this topic.
4 Conclusion
To examine the possibility of PAHs and fulleranes as the carrier of the 21 m emission feature, we fit the infrared spectra of nine PNe exhibiting the 21 m feature. The result shows that fulleranes can provide a better match than PAHs. Further analysis suggests that if fulleranes are responsible for the 21 m feature, their hydrogenation degree must be moderate. The intense mass loss and rapidly increasing ultraviolet radiation of PPNe provide favorable environments for the formation of moderately hydrogenated fullerenes. During further evolution, fulleranes are readily dehydrogenated or destructed, providing a plausible interpretation for the rarity of the 21 m feature.
Nevertheless, we have no means to draw firm conclusions at this moment. We hope that this work could attract research interests in fulleranes in space as collaborative efforts in observations, theories, and experiments are required.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
X-XL: Writing–original draft, Writing–review and editing, Investigation, Formal Analysis, Software, Visualization. YZ: Investigation, Writing–original draft, Writing–review and editing, Conceptualization, Data curation, Funding acquisition, Methodology, Project administration, Resources, Supervision, Validation. SS: Data curation, Formal Analysis, Validation, Writing–original draft, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. The authors declare that financial support was received for the research, authorship, and/or publication of this article. The financial supports of this work are from the National Natural Science Foundation of China (NSFC, No. 12473027, 12473027, and 12333005), the Guangdong Basic and Applied Basic Research Funding (No. 2024A1515010798), and the science research grants from the China Manned Space Project (NO. CMS-CSST-2021-A09, CMS-CSST-2021-A10, etc.). This article is based upon work from COST Action CA21126 - Carbon molecular nanostructures in space (NanoSpace), supported by COST (European Cooperation in Science and Technology).
Acknowledgments
The NASA Ames PAH IR Spectroscopic Database is gratefully acknowledged.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fspas.2024.1489982/full#supplementary-material
Footnotes
References
1
AbbinkD.FoingB.EhrenfreundP. (2024). A model for the hydrogenation and charge states of fullerene C60: implications for diffuse interstellar band research. Astronomy and Astrophysics684, A165–A175. 10.1051/0004-6361/202347478
2
BauschlicherC. W.BoersmaC.RiccaA.MattiodaA. L.CamiJ.PeetersE.et al (2010). The NASA Ames polycyclic aromatic hydrocarbon infrared spectroscopic database: the computed spectra. Astrophysical J. Suppl. Ser.189 (2), 341–351. 10.1088/0067-0049/189/2/341
3
BauschlicherC. W.RiccaA.BoersmaC.AllamandolaL. J. (2018). The NASA Ames PAH IR spectroscopic database: computational version 3.00 with updated content and the introduction of multiple scaling factors. Astrophysical J. Suppl. Ser.234 (2), 32–41. 10.3847/1538-4365/aaa019
4
Buss JrR. H.CohenM.TielensA. G. G. M.WernerM. W.BregmanJ. D.WittebornF. C.et al (1990). Hydrocarbon emission features in the IR spectra of warm supergiants. Astrophysical J.365, L23–L26. 10.1086/185879
5
CamiJ.Bernard-SalasJ.PeetersE.MalekS. E. (2010). Detection of C 60 and C 70 in a young planetary nebula. Science329 (5996), 1180–1182. 10.1126/science.1192035
6
CataldoF. (2003). Fullerane, the hydrogenated C60Fullerene: properties and astrochemical considerations. Fullerenes, Nanotub. Carbon Nanostructures11 (4), 295–316. 10.1081/FST-120025852
7
CataldoF.Iglesias-GrothS. (2009). On the action of UV photons on hydrogenated fulleranes C60H36and C60D36. Mon. Notices R. Astronomical Soc.400 (1), 291–298. 10.1111/j.1365-2966.2009.15457.x
8
CerrigoneL.HoraJ. L.UmanaG.TrigilioC.HartA.FazioG. (2011). Identification of three new protoplanetary nebulae exhibiting the unidentified feature at 21 m. Astrophysical J.738 (2), 121–129. 10.1088/0004-637X/738/2/121
9
Díaz-LuisJ. J.García-HernándezD. A.ManchadoA.CataldoF. (2016). A search for hydrogenated fullerenes in fullerene-containing planetary nebulae. Astronomy Astrophysics589, 5–11. 10.1051/0004-6361/201527727
10
DuleyW. W.WilliamsD. A. (1981). The infrared spectrum of interstellar dust: surface functional groups on carbon. Mon. Notices R. Astronomical Soc.196 (2), 269–274. 10.1093/mnras/196.2.269
11
EvansA.van LoonJ. T.WoodwardC. E.GehrzR. D.ClaytonG. C.HeltonL. A.et al (2012). Solid-phase in the peculiar binary XX Oph?Mon. Notices R. Astronomical Soc. Lett.421 (1), L92–L96. 10.1111/j.1745-3933.2012.01213.x
12
FrischM. J. E. A. (2009). Gaussian 09, Revision d. 01, Gaussian. Wallingford CT: Inc, 201.
13
García-HernándezD. A.Iglesias-GrothS.Acosta-PulidoJ. A.ManchadoA.García-LarioP.StanghelliniL.et al (2011a). The formation of fullerenes: clues from new c 60, c 70, and (possible) planar c 24 detections in magellanic cloud planetary nebulae. Astrophysical J. Lett.737 (2), L30–L37. 10.1088/2041-8205/737/2/L30
14
García-HernándezD. A.ManchadoA.García-LarioP.StanghelliniL.VillaverE.ShawR. A.et al (2010). Formation of fullerenes in H-containing planetary nebulae. Astrophysical J. Lett.724 (1), L39–L43. 10.1088/2041-8205/724/1/L39
15
García-HernándezD. A.RaoN. K.LambertD. L. (2011b). Are C molecules detectable in circumstellar shells of R Coronae Borealis stars?Astrophysical J.729 (2), 126–131. 10.1088/0004-637X/729/2/126
16
García-LarioP.ManchadoA.UllaA.ManteigaM. (1999). Infrared space observatory observations of IRAS 16594–4656: a new Proto-planetary nebula with a strong 21 micron dust feature. Astrophysical J.513 (2), 941–946. 10.1086/306895
17
GehrzR. D.RoelligT. L.WernerM. W.FazioG. G.HouckJ. R.LowF. J.et al (2007). The NASA Spitzer space telescope. Rev. Sci. Instrum.78 (1), 011302. 10.1063/1.2431313
18
GielenC.CamiJ.BouwmanJ.PeetersE.MinM. (2011). Carbonaceous molecules in the oxygen-rich circumstellar environment of binary post-AGB stars: C60fullerenes and polycyclic aromatic hydrocarbons. Astronomy and Astrophysics536, A54–A62. 10.1051/0004-6361/201117961
19
GładkowskiM.SzczerbaR.SloanG. C.LagadecE.VolkK. (2019). 30-micron sources in galaxies with different metallicities. Astronomy and Astrophysics626, A92–A126. 10.1051/0004-6361/201833920
20
GlassgoldA. E.HugginsP. J. (1983). Atomic and molecular hydrogen in the circumstellar envelopes of late-type stars. Mon. Notices R. Astronomical Soc.203, 517–532. 10.1093/mnras/203.2.517
21
HrivnakB. J.VolkK.KwokS. (2000). 2-45 micron infrared spectroscopy of carbon-rich proto-planetary nebulae. Astrophysical J.535 (1), 275–292. 10.1086/308823
22
HrivnakB. J.VolkK.KwokS. (2009). A Spitzer study of 21 and 30m emission in several galactic carbon-rich protoplanetary nebulae. Astrophysical J.694 (2), 1147–1160. 10.1088/0004-637X/694/2/1147
23
Iglesias-GrothS.García-HernándezD. A.CataldoF.ManchadoA. (2012). Infrared spectroscopy of hydrogenated fullerenes (fulleranes) at extreme temperatures. Mon. Notices R. Astronomical Soc.423 (3), 2868–2878. 10.1111/j.1365-2966.2012.21097.x
24
JensenF. (2001). Polarization consistent basis sets: principles. J. Chem. Phys.115 (20), 9113–9125. 10.1063/1.1413524
25
JensenF. (2002). Polarization consistent basis sets. II. Estimating the Kohn–Sham basis set limit. J. Chem. Phys.116 (17), 7372–7379. 10.1063/1.1465405
26
KaboG. J.KarpushenkavaL. S.PaulechkaY. U. (2010). “Thermodynamic properties of fullerene hydrides C and equilibria,” in Fullerane: the hydrogenated fullerenes. Editors CataldoF.Iglesias-GrothS., 55–83.
27
KwokS., (1993). Proto-planetary nebulae. In Symp. - int. Astron. Union, symposium-international astronomical union. 155, 263–270. 10.1017/s0074180900170998Cambridge University Press
28
KwokS. (2004). The synthesis of organic and inorganic compounds in evolved stars. Nature430 (7003), 985–991. 10.1038/nature02862
29
KwokS. (2024). Planetary nebulae research: past, present, and future. Galaxies12 (4), 39–57. 10.3390/galaxies12040039
30
KwokS.VolkK.BernathP. (2001). On the origin of infrared plateau features in proto-planetary nebulae. Astrophysical J.554 (1), L87–L90. 10.1086/320913
31
KwokS.VolkK. M.HrivnakB. J. (1989). A 21 micron emission feature in four proto-planetary nebulae. Astrophysical J.345 (Oct. 1), L51–L54. 10.1086/185550
32
LiA.LiuJ. M.JiangB. W. (2013). On iron monoxide nanoparticles as a carrier of the mysterious 21 m emission feature in post-asymptotic giant branch stars. Astrophysical J.777 (2), 111–117. 10.1088/0004-637X/777/2/111
33
MatsuuraM.Bernard-SalasJ.Lloyd EvansT.VolkK. M.HrivnakB. J.SloanG. C.et al (2014). Spitzer Space Telescope spectra of post-AGB stars in the Large Magellanic Cloud–polycyclic aromatic hydrocarbons at low metallicities. Mon. Notices R. Astronomical Soc.439 (2), 1472–1493. 10.1093/mnras/stt2495
34
MattiodaA. L.HudginsD. M.BoersmaC.BauschlicherC. W.RiccaA.CamiJ.et al (2020). The NASA Ames PAH IR spectroscopic database: the laboratory spectra. Astrophysical J. Suppl. Ser.251 (2), 22–37. 10.3847/1538-4365/abc2c8
35
NeugebauerG.HabingH. J.Van DuinenR.AumannH. H.BaudB.BeichmanC. A.et al (1984). The infrared astronomical satellite (IRAS) mission. Astrophysical J.278 (March 1), L1–L6. 10.1086/184209
36
OtsukaM.KemperF.HyungS.SargentB. A.MeixnerM.TajitsuA.et al (2013). The detection of C in the well-characterized planetary nebula M1-11. Astrophysical J.764 (1), 77–96. 10.1088/0004-637X/764/1/77
37
PalotásJ.MartensJ.BerdenG.OomensJ. (2020). The infrared spectrum of protonated buckminsterfullerene C60H+. Nat. Astron.4 (3), 240–245. 10.1038/s41550-019-0941-6
38
PapoularR. (2011). Candidate carriers and synthetic spectra of the 21- and 30m proto-planetary nebular bands. Mon. Notices R. Astronomical Soc.415, 494–502. 10.1111/j.1365-2966.2011.18724.x
39
ParrR. G. (1985). “Density functional theory in Chemistry,” in Density functional methods in physics (Boston, MA: Springer US), 141–158.
40
RobertsK. R.SmithK. T.SarreP. J. (2012). Detection of C60 in embedded young stellar objects, a Herbig Ae/Be star and an unusual post-asymptotic giant branch star: C60 in young stellar and other objects. Mon. Notices R. Astronomical Soc.421 (4), 3277–3285. 10.1111/j.1365-2966.2012.20552.x
41
RosenbergM. J. F.BernéO.BoersmaC. (2014). Random mixtures of polycyclic aromatic hydrocarbon spectra match interstellar infrared emission. Astronomy Astrophysics566, L4–L10. 10.1051/0004-6361/201423953
42
RüchardtC.GerstM.EbenhochJ.CampbellE. E.TellgmannR.SchwarzH.et al (1993). Transfer hydrogenation and deuteration of buckminsterfullerene C60 by 9,10‐dihydroanthracene and 9,9′,10,10′[D4]Dihydroanthracene. Angew. Chemie-International Ed.32 (4), 584–586. 10.1002/anie.199305841
43
SadjadiS. A.ParkerQ. A. (2021). It remains a cage: ionization tolerance of C60 fullerene in planetary nebulae. FNCN29, 620–625. 10.1080/1536383x.2021.1876677
44
SellgrenK.WernerM. W.IngallsJ. G.SmithJ. D. T.CarletonT. M.JoblinC. (2010). C 60 IN REFLECTION NEBULAE. Astrophysical J. Lett.722 (1), L54–L57. 10.1088/2041-8205/722/1/L54
45
TielensA. G. G. M. (2008). Interstellar polycyclic aromatic hydrocarbon molecules. Annu. Rev. Astronomy Astrophysics46, 289–337. 10.1146/annurev.astro.46.060407.145211
46
Van WinckelH. (2003). Post-AGB stars. Annu. Rev. Astronomy Astrophysics41 (1), 391–427. 10.1146/annurev.astro.41.071601.170018
47
VolkK.HrivnakB. J.MatsuuraM.Bernard-SalasJ.SzczerbaR.SloanG. C.et al (2011). Discovery and analysis of 21 m feature sources in the magellanic clouds. Astrophysical J.735 (2), 127–154. 10.1088/0004-637X/735/2/127
48
VolkK.SloanG. C.KraemerK. E. (2020). The 21 m and 30 m emission features in carbon-rich objects. Astrophysics Space Sci.365, 88–25. 10.1007/s10509-020-03798-2
49
VolkK. M.KwokS. (1989). Evolution of protoplanetary nebulae. Astrophysical J.342 (July 1), 345–363. 10.1086/167597
50
WebsterA. (1995). The lowest of the strongly infrared-active vibrations of the fulleranes and astronomical emission band at a wavelength of 21m. Mon. Notices R. Astronomical Soc.277 (4), 1555–1566. 10.1093/mnras/277.4.1555
51
ZhangK.JiangB. W.LiA. (2009). On the carriers of the 21m emission feature in post-asymptotic giant branch stars. Mon. Notices R. Astronomical Soc.396, 1247–1256. 10.1111/j.1365-2966.2009.14808.x
52
ZhangY. (2020). Are fulleranes responsible for the 21 micron feature?Chin. J. Chem. Phys.33, 101–106. 10.1063/1674-0068/cjcp1911203
53
ZhangY.KwokS. (2011). DETECTION OF C60IN THE PROTOPLANETARY NEBULA IRAS 01005+7910. Astrophysical J.730 (2), 126–130. 10.1088/0004-637X/730/2/126
54
ZhangY.KwokS. (2015). On the viability of the PAH model as an explanation of the unidentified infrared emission features. Astrophysical J.798 (1), 37–42. 10.1088/0004-637X/798/1/37
55
ZhangY.KwokS.HrivnakB. J. (2010). A spitzer/infrared spectrograph spectral study of a sample of galactic carbon-rich proto-planetary nebulae. Astrophysical J.725 (1), 990–1001. 10.1088/0004-637X/725/1/990
56
ZhangY.SadjadiS.HsiaC. H. (2020). Hydrogenated fullerenes (fulleranes) in space. Astrophysics Space Sci.365 (4), 67–74. 10.1007/s10509-020-03779-5
57
ZhangY.SadjadiS.HsiaC. H.KwokS. (2017). Search for hydrogenated C60 (fulleranes) in circumstellar envelopes. Astrophysical J.845 (1), 76–91. 10.3847/1538-4357/aa71ac
Summary
Keywords
infrared, polycyclic aromatic hydrocarbons, fullerenes, fulleranes, protoplanetary nebulae, astrochemistry
Citation
Liao X-X, Zhang Y and Sadjadi S (2024) Exploring the possibility of polycyclic aromatic hydrocarbons and fulleranes as the carrier of the 21 micron emission feature. Front. Astron. Space Sci. 11:1489982. doi: 10.3389/fspas.2024.1489982
Received
02 September 2024
Accepted
25 September 2024
Published
11 November 2024
Volume
11 - 2024
Edited by
Ryszard Szczerba, Polish Academy of Sciences, Poland
Reviewed by
Franco Cataldo, Actinium Chemical Research Institute, Italy
Ransel Barzaga, Spanish National Research Council (CSIC), Spain
Updates
Copyright
© 2024 Liao, Zhang and Sadjadi.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yong Zhang, zhangyong5@mail.sysu.edu.cn
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.