Issue |
A&A
Volume 682, February 2024
|
|
---|---|---|
Article Number | L18 | |
Number of page(s) | 7 | |
Section | Letters to the Editor | |
DOI | https://doi.org/10.1051/0004-6361/202349087 | |
Published online | 21 February 2024 |
Letter to the Editor
Hydrogenated amorphous carbon grains as an alternative carrier of the 9–13 μm plateau feature in the fullerene planetary nebula Tc 1
1
Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife, Spain
e-mail: magm@iac.es; agarcia@iac.es
2
Departamento de Astrofísica, Universidad de La Laguna, 38206 La Laguna, Tenerife, Spain
3
Consejo Superior de Investigaciones Científicas (CSIC), 28006 Madrid, Spain
Received:
22
December
2023
Accepted:
5
February
2024
Fullerenes have been observed in several astronomical objects since the discovery of C60 in the mid-infrared (mid-IR) spectrum of the planetary nebula (PN) Tc 1. It has been suggested that the carriers of the broad unidentified infrared (UIR) plateau features, such as the 9–13 μm emission feature (12 μm hereafter), may be related to the formation of fullerenes. In particular, their carriers have been suggested to be mixed aromatic or aliphatic hydrocarbons such as hydrogenated amorphous carbon (HAC-like hereafter) grains. For this study, we modeled the mid-IR emission of the C60-PN Tc 1 with a photoionization code, including for the first time the laboratory optical constants (n and k indices) of HAC-like dust at 300 K. Interestingly, we find that the broad 12 μm plateau feature in Tc 1 is well reproduced by using a distribution of canonical HAC grains, while at the same time they provide an important fraction of the IR dust continuum emission and are consistent with the other UIR features observed (e.g., the broad 6–9 μm plateau feature). This finding suggests that HAC-like grains may be possible carriers of the 12 μm plateau feature, being likely related to the fullerene formation mechanism in PNe. More laboratory experiments, to obtain the optical constants of HAC-like dust with several structures or a composition at different physical conditions, are strongly encouraged – that is, in order to extend this pilot study to more fullerene PNe, and to unveil the details of fullerene formation and of the potential carriers of the elusive UIR plateau features.
Key words: astrochemistry / stars: AGB and post-AGB / planetary nebulae: individual: PN Tc 1
© The Authors 2024
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
1. Introduction
Planetary nebulae (PNe) represent the late stages in the evolution of low- and intermediate-mass stars (LIM; ∼1–8 M⊙), the majority of stars in the Universe. On their way to the PN phase, LIM stars experience a strong mass loss during the preceding asymptotic giant branch (AGB; see e.g., Herwig 2005 for a review) phase and they chemically enrich the surrounding interstellar medium (ISM). The AGB stars are also major suppliers of dust grains and molecular species that are routinely seen in the local Universe (e.g., Ferrarotti & Gail 2006), which makes them (and their subsequent evolutionary stages) very important for our fundamental understanding of the enrichment and chemical composition of the ISM.
Fullerenes, among the most resistant and stable three-dimensional molecules that are only formed by C atoms, such as C60 (first discovered in the laboratory by Kroto et al. 1985), were supposed to be widely spread in space for decades. The presence of fullerenes in astrophysical environments was under debate until the infrared (IR) signatures of the C60 and C70 fullerenes were unambiguously detected in the Spitzer mid-IR spectrum of the young PN Tc 1 (Cami et al. 2010). After the first detection of fullerenes in a PN, the presence of C60 was found in several other astrophysical objects such as reflection nebulae (Sellgren et al. 2010), additional PNe (García-Hernández et al. 2010, 2011a, 2012; Otsuka et al. 2014), peculiar R Coronae Borealis (RCB) stars (García-Hernández et al. 2011b), post-AGB stars (Zhang & Kwok 2011; Gielen et al. 2011), and Herbig Ae/Be stars (Arun et al. 2023), among others, by its four strongest mid-IR emission features at ∼7.0, 8.5, 17.4, and 18.9 μm. However, the formation route of fullerenes in such H-rich objects, similar to the majority of PNe, is not well understood. Nowadays, the most suitable routes for the formation of fullerenes are1: i) the photochemical processing of hydrogenated amorphous carbon (HAC) grains or similar mixed aromatic or aliphatic hydrocarbons (i.e., same chemical composition but a different internal structure; hereafter HAC-like; García-Hernández et al. 2010) and ii) the photochemical processing of large polycyclic aromatic hydrocarbons (PAHs; Berné & Tielens 2012; Murga et al. 2022).
The C60 fullerenes are mainly detected toward PNe, and fullerene PNe are very young low-excitation (Teff ∼ 30 000–45 000 K) C-rich objects, which evolved from low-mass (∼1–3 M⊙) progenitors (García-Hernández et al. 2012; Otsuka et al. 2014). Remarkably, the IR Spitzer spectra (∼5–38 μm) of the fullerene-rich PNe are dominated by aliphatic hydrocarbon-rich dust emission features (superimposed on the underlying and featureless dust continuum emission), showing several broad unidentified IR (UIR) plateau emission features at 6–9 (hereafter 7 μm), 9–13 (hereafter 12 μm), 15–20, and 25–35 μm (e.g., García-Hernández et al. 2012). Characteristic aliphatic discrete features are located at 3.4 and 6.9 μm, arising from symmetric and asymmetric C-H stretching and bending modes of methyl and methylene groups attached to aromatic rings, respectively (see e.g., Kwok & Zhang 2011). The identification of the carriers of the broad (and discrete) UIR emission features widely observed in space (toward Solar System bodies, circumstellar envelopes of PNe, diffuse ISM, and remote galaxies, among others) is a long-standing problem in astrophysics (see e.g., Kwok 2016 for a review). Thus, the potential identification of the carrier(s) of any of the broad UIR plateau emission features in fullerene PNe would impact very different astronomical fields.
Nano- and micro-sized particles in the ISM and circumstellar environment play an important role in astrophysical processes that lead to the formation of such exotic UIR emissions seen in the spectra of C-rich PNe and other astrophysical objects. The optical constants, which are the complex refractive indices or the complex dielectric function, of different materials (e.g., graphite, silicates, SiC, MgS, and HAC) are commonly used, in conjunction with radiative transfer codes, to explain the origin of the IR emission in astronomical spectra (e.g., HACs, soots, SiC, and other hydrocarbon materials; Bernard-Salas et al. 2012; Gavilan et al. 2016; Gómez-Llanos et al. 2018; Dubosq et al. 2023 see respectively). Fortunately, the latter can be done thanks to the databases2 dedicated to compile optical properties of analog materials of cosmic dust in the wavelength range from the UV to the far-IR.
In particular, the broad UIR plateau emission feature at 12 μm was first observed in the low resolution spectrometer of the IR astronomical satellite (LRS IRAS) spectra by Cohen et al. (1985), and its carrier’s identification has been a long outstanding problem. Cohen et al. (1985) attributed the carrier to PAHs, while Blanco et al. (1988) and Buss et al. (1990) attributed both the 7 and 12 μm plateau features to a mixture of large PAHs and HAC grains. Later, Kwok et al. (2001) specifically suggested that the 7 and 12 μm plateaus are superpositions of in-plane and out-of-plane bending modes of aliphatic side groups. Despite these earlier suggestions, the broad UIR emission at 12 μm has been generally attributed, in the literature, to SiC (α-SiC; e.g., Speck et al. 2009; Jones et al. 2023). García-Hernández (2012) argued against the SiC identification of this feature in fullerene PNe because the spectral characteristics (position and shape) (see Table 4 in García-Hernández et al. 2012) of the 9–13 μm emission in C60-PNe significantly differ from those of the SiC 11.5 μm feature seen in the preceding AGB phase. This may suggest that the carrier could be different in AGBs and in fullerene PNe or even that the identification of SiC in AGB stars could be revisited.
The PN Tc 1 is the prototypical fullerene-rich PN (e.g., Cami et al. 2010; Aleman et al. 2019), being one of the best objects to study because of its simple round-elliptical morphology, the available high-quality and high-resolution IR Spitzer spectra, and the lack of PAH emission at 6.2, 7.7, 8.6, and 11.3 μm (Cami et al. 2010), which could strongly contribute to the fullerene and broad (and discrete) UIR emission features. In this Letter, we present an analysis of the high-quality Spitzer spectrum of PN Tc 1 in comparison with a photoionization model that includes, for the first time, the emission of HAC-like dust grains as measured in laboratory, showing that HAC-like grains are a possible explanation, different from the often assumed SiC, for the broad 12 μm plateau emission feature.
2. Infrared spectrum of the C60-PN Tc 1
We downloaded all the available mid-IR spectra from the Spitzer IR Spectrograph (IRS) (Werner et al. 2004; Houck et al. 2004) Heritage archive of the PN Tc 1 for which high-quality – a signal-to-noise ratio (S/N) > 50 – short-low (SL; 5.2–14.5 μm), short-high (SH; 9.9–19.6 μm), and long-high (LH; 18.7–37.2 μm) spectra were available (Program Identifier 3633) (see below and Perea-Calderón et al. 2009 for more data reduction details). We used the package SMART3 (v8.2.9; Higdon et al. 2004) to process the extracted wavelength and flux calibrated 1D spectra for each nod position, that is, to clean for bad data points, spurious jumps, and glitches and to combine and merge into one single spectrum per module (SL, SH, and LH). In order to get the final merged spectrum, we scaled the flux density of the SH spectra to match the overlapping wavelength region in the LH. Then, we scaled the SL spectra to match the overlapping wavelength region to the SH spectra. Finally, we scaled the flux density of the resulted combined spectra to that of the Wide-field IR Survey Explorer (WISE) band 4 (λc = 22 μm) photometric value (similarly to what was done in Otsuka et al. 2014). The resulting final Tc 1 Spitzer IR spectrum is shown in Fig. 1, where the emission features such as the C60 and C70 IR bands, and the broad 7 and 12 μm plateau emission features are indicated.
Fig. 1. Spitzer mid-IR spectrum of the PN Tc 1 (black line) compared with the best photoionization model spectrum including HAC-like dust grains (green line). The C60 (vertical dashed blue line) and C70 (vertical dashed red line) emission bands and UIR plateau emission features are indicated. We note that the narrow emission features seen in the best model spectrum (green) are just atomic nebular emission lines, whose detailed modeling is out of the scope of the present work. |
3. Cloudy simulations
To the best of our knowledge, the only previous attempt to model the 7 and 12 μm plateau emission features seen in the Spitzer spectrum of the C60-PN Tc 1 (Fig. 1) was made by Bernard-Salas et al. (2012). They used the Jones (2012) theoretical absorption coefficients (Qabs) of 3 nm HAC particles with a different H content (H/C = 0.23, 0.29, and 0.35) in conjunction with a 200 K blackbody (see Fig. 7 in Bernard-Salas et al. 2012). They could reproduce the 7 μm plateau emission feature – as seen in the Tc 1 continuum-subtracted Spitzer spectrum – with some degree of success, depending on the atom hydrogen fractions of the HAC particles. However, they could not reproduce the Tc 1 12 μm plateau emission feature. In particular, a too strong 11.3 μm feature was theoretically predicted for all H/C ratios, and they did not study or analyze the HAC contribution to the IR dust continuum emission. Bernard-Salas et al. (2012) concluded that their approximate reproduction of the 7 μm plateau emission feature may imply that fullerenes are formed by photochemical processing of HACs. Below we explain how we used, for the first time, a detailed photoionization model of Tc 1 together with the laboratory optical constants of HAC-like dust in order to model the 7 and 12 μm plateau features and the underlying dust continuum emission.
3.1. Modeling the 12 μm plateau emission
We used the stellar parameters and nebular abundances predicted in the detailed photoionization model developed by Aleman et al. (2019)4 (without graphite as the IR dust continuum source; see below) as input in CLOUDY (v22.02; Ferland et al. 2017) to model the mid-IR emission of the C60-PN Tc 1 by using, for the first time, the laboratory optical constants of HAC-like dust. We note that a change in the dust composition does not significantly change the emission line and Hβ fluxes predicted in the Aleman et al. model (a maximum flux difference of ∼5% is obtained), and so the stellar and nebular parameters and chemical abundances are unaffected. The laboratory optical constants, that is, the complex refractive indices n and k of HAC-like particles were obtained from Prof. W. Duley (priv. comm.; see Duley 1984 for a detailed description of the laboratory measurements), covering a spectral range from 0.1 to 17 μm with a gap between 0.65 to 2.5 μm due to a limitation of the laboratory instrumentation to measure over this wavelength range (Fig. 2). The former are the only available laboratory optical constants, prepared under simulated interstellar conditions, covering a sufficient range from the UV to the mid-IR. The laboratory n and k values were measured in HAC thin films, ∼50–200 nm thickness, deposited at 300 K and with slightly different hydrogen concentrations (H/C ∼ 0.3–0.4; Duley et al. 1998; Duley & Hu 2012). However, it should be noted that it is mostly the k index that is sensitive to the H content; this is because the concentration of CH and CHn in the HAC samples determines the IR k. For simplicity, we adopted the mean n and k values (presumably with H/C ∼ 0.35), which are shown in Fig. 2. Because CLOUDY needs a minimum wavelength mesh (∼0.001–96 μm), an extension of the n and k values to lower and upper wavelengths was carried out in a similar way as in Gómez-Llanos et al. (2018); that is, we extrapolated the n value to lower wavelengths on the basis of the BE1 amorphous carbon included in CLOUDY (Rouleau & Martin 1991), and for k we extrapolated using a power law of slope 2, whereas for higher wavelengths we extrapolated n as a constant value and k using a power law of slope −1. For the gap between 0.65 to 2.5 μm, a linear extrapolation was used.
Fig. 2. Laboratory average values of the complex refractive indices n and k (black and purple lines, respectively) of HAC-like particles at 300 K and with a hydrogen concentration of H/C ≃ 0.35 (Prof. W. Duley, priv. comm.; see text for more details). The errors of the complex refractive indices were measured to be ∼20% (see e.g., Duley 1984). |
For the model, we compiled the HAC optical constants with the grains command inside CLOUDY, using a canonical ISM grain distribution (amin = 0.001 μm to amax = 0.25 μm; 1–250 nm5, distributed in ten grain size logarithmic steps following a power of −3.5), and generated the opacities (with the integrated Mie code for spherical grains; see van Hoof et al. 2004) needed for the CLOUDY simulations. We note that CLOUDY takes into account the molar mass of the HAC (12.36 g mol−1), the band gap energy ( eV; Duley & Hu 2012), and the density of the material (1.5 g cm−3) (see Duley et al. 1998; Duley & Hu 2012). A reduced minimization was carried out for a grid of CLOUDY models – by varying the HAC grains’ relative content – to find a convergence between the model and the observed spectrum of Tc 1 between 10 and 14 μm (485 spectral points), yielding a ≃ 1.6 assuming a 10% error in the observations. This resulted in a C abundance (in units of log(C/H)) of −4.60 for the HAC-like grains, which is ∼3% of the total C abundance (gas + dust phase abundance; C ≃ −3.08), a mean dust temperature across the nebula of ∼215 K, and a dust-to-gas (D/G) ratio of 2.1 × 10−4.
Figure 1 also shows the best CLOUDY model (green line) obtained from the analysis applied to the observed spectrum of Tc 1. As it can be seen, the HAC-like grains at 300 K6 provide a very good fit to the broad 12 μm plateau feature in Tc 1, while at the same time they fully reproduce the IR dust continuum emission from 9 to 14 μm. For comparison, α-SiC dust grains were also used to similarly model the mid-IR spectrum of Tc 1, which resulted in a worst fit to the 12 μm feature (see Appendix A). Interestingly, it turns out that HAC-like grains may be an alternative explanation of the broad 12 μm feature because it is very well reproduced with a simple spherical grain size distribution using the laboratory HAC optical constants at a deposition temperature close to the fullerene-dust temperature in Tc 1.
In addition, there is a small contribution to the observed 7 μm plateau feature, which is also clearly seen in the model spectrum. However, we still would need to include other emitters and/or an additional dust component to reproduce the 7 μm plateau feature (see below) in addition to the shape and rise of the underlying dust continuum emission at 5–9 and 14–17 μm, respectively.
3.2. Modeling the 12 μm plateau emission and dust continuum
Thus, in conjunction with the HAC-like grains as above, we added graphite dust grains separately as a possible additional dust continuum component (as assumed by Aleman et al. (2019)) and we again applied a analysis to a grid of CLOUDY models in order to find the best fit by varying the HAC grains and graphite relative contents in the 9–16 μm range, yielding a value of ≃ 1.3. This resulted in C abundances of −3.74 and −4.90, which correspond to ∼18.7% and ∼1.3% of the total C abundance (gas + dust phase) as well as D/G of 1.5 × 10−3 and 1.1 × 10−4 for graphite and HAC grains, respectively.
We note that other carbon materials such as amorphous carbon (AC), generally used in the literature to model the IR dust continuum emission, provide worst fits because AC marginally contribute below 14 μm. Figure 3 displays the best photoionization model spectra for HAC-like grains and graphite, overplotted on the Tc 1 Spitzer spectrum; the HAC-like and graphite IR emission contributions are also shown separately for comparison. This way, we can reproduce the Tc 1 IR spectrum well from ∼5 to 17 μm (both the 12 μm plateau feature and the dust continuum) with the HAC-like grains providing an important fraction (∼20–40%) of the dust continuum emission, as evaluated from the integrated IR fluxes due to only graphite and HACs with respect to the total integrated IR continuum flux at three different wavelengths (5.6, 9.6, and 15.2 μm) representative of the dust continuum emission. Again, for comparison, α-SiC dust grains in conjunction with graphite dust grains were also used to similarly model the mid-IR spectrum of Tc 1, resulting in a worst fit of the red wing of the broad 12 μm feature (see the Appendix A). Finally, we stress that preliminary photoionization models with HAC-like grains deposited at 300 K also provide acceptable good fits to the 12 μm plateau feature (and dust continuum) in another fullerene PNe such as IC 418 and SMC 16, when playing with the model parameters (e.g., grain size and shape distribution; Gómez-Muñoz et al., in prep.).
Fig. 3. Spitzer mid-IR spectrum of the PN Tc 1 (black line) compared with the best photoionization model spectra for HAC-like (dashed gray line) and graphite (dashed magenta line) dust grains separately. The best photoionization model including HAC-like grains in conjunction with graphite dust grains (solid green line) is also shown for comparison. |
4. Astrophysical implications
To the best of our knowledge, this is the first time that HAC grains convincingly reproduce the 12 μm plateau emission feature in a fullerene PN (Tc 1), posing doubts on the generally accepted SiC identification in the literature as the main carrier of this feature. Remarkably, the HAC grains naturally provide an important fraction of the IR dust continuum emission and are consistent with the other UIR features observed such as the 7 μm plateau feature (see below). From our CLOUDY models, we determined that ∼1.3–3.0% (with and without graphite, respectively) of the C in Tc 1 is in the form of HAC-like dust. These values are comparable to the amount of C that is estimated to be in the form of fullerenes, such as C60, in Tc 1 (∼1.5%; see Cami et al. 2010).
The laboratory HAC-like grains (H/C ∼ 0.35) used in our models (Duley & Hu 2012) mainly contain aliphatic carbon chains (sp3 content), which represent 72% (i.e., aromatic/aliphatic ratio of ∼0.4) of the chemical composition, but also polyyne (-C≡C-)n and cumulenic (-C=C-)n chains. According to Duley & Hu (2012), these HAC-like grains can have more aromatic rings than aromatic chains. In fact, they suggest that above 10 μm, the features can be caused by a variety of vibrational modes associated with ring deformation and with CHn groups (i.e., out-of-plane CH vibrations). In particular, the spectral range of the broad 12 μm plateau emission may contain the features produced by small aromatic rings (pentagons and hexagons) with side chains such as dimethylnaphthalene or dimethylphenanthrene, but also those produced by solo, duo CH, and CHn groups as components of the hydrocarbon matrix (Duley & Hu 2012). Clearly, the broadening of the 12 μm plateau emission is the result of this coexistence of small aromatic rings and hydrocarbon chains; that is, the mixed aromatic and aliphatic hydrocarbon composition inside the HAC-like material is the key behind the broadening of the 12 μm feature. We remark that the only other laboratory optical constants available to us correspond to those measured in HAC films deposited at a lower temperature of 77 K (Prof. W. Duley, priv. comm.) and with slightly higher H concentrations (H/C ∼ 0.5–0.6), which have fewer aromatic rings in its composition (Duley et al. 1998). The 12 μm plateau feature in Tc 1 cannot be reproduced with the 77 K HAC-like sample (i.e., it is absent; Gómez-Muñoz et al., in prep.), suggesting that the aromatic ring/chain ratio, rather than the aromatic/aliphatic ratio, is a key parameter for the formation of the 12 μm plateau feature.
The presence of HAC-like grains in conjunction with fullerenes have been previously proposed in the literature as proof of the photochemical processing of HAC-like grains for fullerene formation (see e.g., García-Hernández et al. 2010, 2011a,b, 2012; Bernard-Salas et al. 2012). Naturally, our finding in Tc 1 suggests that the formation of fullerenes in PNe circumstellar environments may be related to the photochemical processing of HAC-like grains. However, the only experimental data on the decomposition of HAC-like grains are those from Scott et al. (1997), where they performed laser ablation (UV radiation-induced decomposition) of HAC films and found PAHs and fullerenes, among other species, but the chemical reactions are presently unknown. The lower aromatic content compared to more dehydrogenated hydrocarbons in our HAC sample suggest that the main process to generate the precursors to form C60, for example, the dehydrogenation of the HAC, is through radiative processing as aliphatic hydrocarbons have been proven to be efficiently dehydrogenated in the ISM by UV radiation (e.g., Muñoz Caro et al. 2001). We speculate that the presence of some small aromatic rings in our 300 K HAC-like sample might act as templates or catalysts for the formation of more complex Cn structures when the material is subjected to photochemical processing. More sophisticated laboratory experiments on the photochemical transformation of HAC-like grains at circumstellar conditions would be needed to unveil the details of fullerene formation via chemical reaction routes.
Another obvious implication of our work involves the carriers of the 7 μm plateau emission feature seen in the fullerene PN Tc 1. Our findings also suggest that HAC-like grains, with H/C = 0.35 and mainly composed by aliphatic carbon chains (72%; see above) with some small aromatic rings, are not a main contributor to such a 7 μm emission complex and that other still unidentified species should be the carriers. Fullerene-based species are natural candidate species as potential additional emitters at these wavelengths in fullerene PNe. Very recent quantum-chemistry simulations of metallofullerenes (neutral and charged endohedral and exohedral species) do indeed demonstrate that they strongly emit at 6–9 μm (Barzaga et al. 2023a,b; Hou et al. 2023) with no significant contribution in the 9–13 μm spectral region.
5. Concluding remarks
This pilot work leads us to conclude that HAC-like grains are a convincing alternative explanation (i.e., different from the often assumed SiC) for the broad 12 μm plateau emission feature seen in the fullerene PN Tc 1. In addition, our work suggests that HAC-like grains may be possible carriers of the 12 μm plateau feature, being likely related to the fullerene formation process (possibly as fullerene precursors) in the circumstellar environment of PNe.
Interestingly, HAC-like grains can also reproduce the far-ultraviolet (FUV) rise seen in the extinction curve (Gavilan et al. 2016, 2017) and very recently detected also in the International Ultraviolet Explorer (IUE) UV spectrum of Tc 1 (Gómez-Muñoz et al. 2024). In this particular case, Gómez-Muñoz et al. (2024) could reproduce the FUV rise toward Tc 1 very well by using HAC-like nanograins (1–15 nm). Although other carbonaceous materials such as hydrogenated carbon clusters and large PAHs, possibly among others, may also explain the FUV rise in extinction curves (see e.g., Dubosq et al. 2023; Lin et al. 2023); such an independent finding would be consistent with the results presented here.
Yet, we do not fully understand the exact chemical pathways from large HAC-like dust grains toward the most stable fullerene molecules as well as the origin of the other broad plateau emission features seen in fullerene PNe, such as 6–9, 15–20, and 25–35 μm. The James Webb Space Telescope offers the opportunity to observe these broad features in fullerene PNe at a much higher sensitivity and spectral resolution than Spitzer, permitting the variations in their positions and shapes to be studied precisely (even resolving a possible substructure) and providing new additional constraints for future modeling studies. Undoubtedly, more laboratory efforts are needed in order to definitively reveal the details of fullerene formation in PNe circumstellar environments as well as the carriers of the elusive UIR plateau features. Our work strongly encourages more laboratory experiments to obtain the refractive indices n and k (over a large wavelength range, from the UV to the far-IR) of HAC-like dust grains with several structures and a composition (e.g., H/C ratios) as well as at different physical conditions (e.g., temperatures or aromatic ring/chain ratios). Such laboratory measurements (not presently available) would be used to expand the pilot novel work presented here for Tc 1 to a larger sample of fullerene PNe.
We note that, more recently, Bernal et al. (2019) proposed the shock heating and ion bombardment induced processing of SiC grains as an alternative route toward fullerenes in the ISM, but such a strong energetic process seems to be unlikely to happen in the circumstellar environments of fullerene PNe where the atomic nebular emission lines can be well explained by photoionization and shocks seem to be unimportant for fullerene formation (see e.g., García-Hernández et al. 2012).
For example, the Jena database (https://www2.mpia-hd.mpg.de/HJPDOC/index.php) and the refractive index database (https://refractiveindex.info/).
Aleman et al. (2019) developed a detailed photoionization model to obtain accurate stellar and physical parameters of the star and the nebula, respectively, and in which they assumed graphite in order to model the IR dust continuum emission. It is available at https://github.com/Morisset/Tc1
We note that this value is close to the circumstellar dust grains’ equilibrium temperature and the excitation temperature (∼330 K) of the C60 emission in Tc 1 as derived by Cami et al. (2010), who assumed that C60 is attached to the dust grains.
Acknowledgments
We kindly acknowledge Prof. Walt Duley for supplying us with the laboratory optical constants of HAC-like particles and Dr. Christophe Morisset for earlier work and discussions on the IR spectral modeling, including HAC-like dust, with the photoionization code CLOUDY. We acknowledge the support from the State Research Agency (AEI) of the Spanish Ministry of Science and Innovation (MCIN) under grant PID2020-115758GB-I00/AEI/10.13039/501100011033. This article is based upon work from European Cooperation in Science and Technology (COST) Action NanoSpace, CA21126, supported by COST. This work is based on observations made with the Spitzer Space Telescope, which was operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA.
References
- Aleman, I., Leal-Ferreira, M. L., Cami, J., et al. 2019, MNRAS, 490, 2475 [NASA ADS] [CrossRef] [Google Scholar]
- Arun, R., Mathew, B., Manoj, P., et al. 2023, MNRAS, 523, 1601 [NASA ADS] [CrossRef] [Google Scholar]
- Barzaga, R., García-Hernández, D. A., Díaz-Tendero, S., et al. 2023a, ApJ, 942, 5 [NASA ADS] [CrossRef] [Google Scholar]
- Barzaga, R., García-Hernández, D. A., Díaz-Tendero, S., et al. 2023b, ApJS, 269, 26 [NASA ADS] [CrossRef] [Google Scholar]
- Bernal, J. J., Haenecour, P., Howe, J., et al. 2019, ApJ, 883, L43 [NASA ADS] [CrossRef] [Google Scholar]
- Bernard-Salas, J., Cami, J., Peeters, E., et al. 2012, ApJ, 757, 41 [NASA ADS] [CrossRef] [Google Scholar]
- Berné, O., & Tielens, A. G. G. M. 2012, Proc. Nat. Acad. Sci., 109, 401 [Google Scholar]
- Blanco, A., Bussoletti, E., & Colangeli, L. 1988, ApJ, 334, 875 [NASA ADS] [CrossRef] [Google Scholar]
- Buss, R. H., Jr, Cohen, M., Tielens, A. G. G. M., et al. 1990, ApJ, 365, L23 [NASA ADS] [CrossRef] [Google Scholar]
- Cami, J., Bernard-Salas, J., Peeters, E., & Malek, S. E. 2010, Science, 329, 1180 [NASA ADS] [CrossRef] [Google Scholar]
- Cohen, M., Tielens, A. G. G. M., & Allamandola, L. J. 1985, ApJ, 299, L93 [NASA ADS] [CrossRef] [Google Scholar]
- Dubosq, C., Pla, P., Dartois, E., & Simon, A. 2023, A&A, 670, A175 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Duley, W. W. 1984, ApJ, 287, 694 [NASA ADS] [CrossRef] [Google Scholar]
- Duley, W. W., & Hu, A. 2012, ApJ, 761, 115 [NASA ADS] [CrossRef] [Google Scholar]
- Duley, W. W., Scott, A. D., Seahra, S., & Dadswell, G. 1998, ApJ, 503, L183 [NASA ADS] [CrossRef] [Google Scholar]
- Ferland, G. J., Chatzikos, M., Guzmán, F., et al. 2017, Rev. Mex. A&A, 53, 385 [Google Scholar]
- Ferrarotti, A. S., & Gail, H. P. 2006, A&A, 447, 553 [CrossRef] [EDP Sciences] [Google Scholar]
- García-Hernández, D. A. 2012, Planetary Nebulae: An Eye to the Future, 283, 148 [Google Scholar]
- García-Hernández, D. A., Manchado, A., García-Lario, P., et al. 2010, ApJ, 724, L39 [CrossRef] [Google Scholar]
- García-Hernández, D. A., Iglesias-Groth, S., Acosta-Pulido, J. A., et al. 2011a, ApJ, 737, L30 [CrossRef] [Google Scholar]
- García-Hernández, D. A., Kameswara Rao, N., & Lambert, D. L. 2011b, ApJ, 729, 126 [CrossRef] [Google Scholar]
- García-Hernández, D. A., Villaver, E., García-Lario, P., et al. 2012, ApJ, 760, 107 [CrossRef] [Google Scholar]
- Gavilan, L., Alata, I., Le, K. C., et al. 2016, A&A, 586, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gavilan, L., Le, K. C., Pino, T., et al. 2017, A&A, 607, A73 [Google Scholar]
- Gielen, C., Cami, J., Bouwman, J., Peeters, E., & Min, M. 2011, A&A, 536, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gómez-Llanos, V., Morisset, C., Szczerba, R., García-Hernández, D. A., & García-Lario, P. 2018, A&A, 617, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gómez-Muñoz, M. A., García-Hernández, D. A., Manchado, A., Barzaga, R., & Huertas-Roldán, T. 2024, MNRAS, 528, 2871 [CrossRef] [Google Scholar]
- Herwig, F. 2005, ARA&A, 43, 435 [NASA ADS] [CrossRef] [Google Scholar]
- Higdon, S. J. U., Devost, D., Higdon, J. L., et al. 2004, PASP, 116, 975 [NASA ADS] [CrossRef] [Google Scholar]
- Hou, G.-L., Lushchikova, O. V., Bakker, J. M., et al. 2023, ApJ, 952, 13 [NASA ADS] [CrossRef] [Google Scholar]
- Houck, J. R., Roellig, T. L., van Cleve, J., et al. 2004, ApJS, 154, 18 [NASA ADS] [CrossRef] [Google Scholar]
- Jones, A. P. 2012, A&A, 542, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jones, O. C., Álvarez-Márquez, J., Sloan, G. C., et al. 2023, MNRAS, 523, 2519 [Google Scholar]
- Kroto, H. W., Heath, J. R., Obrien, S. C., Curl, R. F., & Smalley, R. E. 1985, Nature, 318, 162 [NASA ADS] [CrossRef] [Google Scholar]
- Kwok, S. 2016, A&ARv, 24, 8 [NASA ADS] [CrossRef] [Google Scholar]
- Kwok, S., & Zhang, Y. 2011, Nature, 479, 80 [NASA ADS] [CrossRef] [Google Scholar]
- Kwok, S., Volk, K., & Bernath, P. 2001, ApJ, 554, L87 [NASA ADS] [CrossRef] [Google Scholar]
- Lin, Q., Yang, X. J., & Li, A. 2023, MNRAS, 525, 2380 [NASA ADS] [CrossRef] [Google Scholar]
- Muñoz Caro, G. M., Ruiterkamp, R., Schutte, W. A., Greenberg, J. M., & Mennella, V. 2001, A&A, 367, 347 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Murga, M. S., Akimkin, V. V., & Wiebe, D. S. 2022, MNRAS, 517, 3732 [CrossRef] [Google Scholar]
- Otsuka, M., Kemper, F., Cami, J., Peeters, E., & Bernard-Salas, J. 2014, MNRAS, 437, 2577 [NASA ADS] [CrossRef] [Google Scholar]
- Perea-Calderón, J. V., García-Hernández, D. A., García-Lario, P., Szczerba, R., & Bobrowsky, M. 2009, A&A, 495, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Rouleau, F., & Martin, P. G. 1991, JRASC, 85, 201 [NASA ADS] [Google Scholar]
- Scott, A., Duley, W. W., & Pinho, G. P. 1997, ApJ, 489, L193 [NASA ADS] [CrossRef] [Google Scholar]
- Sellgren, K., Werner, M. W., Ingalls, J. G., et al. 2010, ApJ, 722, L54 [NASA ADS] [CrossRef] [Google Scholar]
- Speck, A. K., Corman, A. B., Wakeman, K., Wheeler, C. H., & Thompson, G. 2009, ApJ, 691, 1202 [NASA ADS] [CrossRef] [Google Scholar]
- van Hoof, P. A. M., Weingartner, J. C., Martin, P. G., Volk, K., & Ferland, G. J. 2004, MNRAS, 350, 1330 [NASA ADS] [CrossRef] [Google Scholar]
- Werner, M. W., Roellig, T. L., Low, F. J., et al. 2004, ApJS, 154, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Zhang, Y., & Kwok, S. 2011, ApJ, 730, 126 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Modeling the Tc 1 mid-IR spectrum using α-SiC grains
For comparison, Fig. A.1 shows the mid-IR spectrum of Tc 1 together with the best CLOUDY model spectrum obtained using α-SiC grains (as in Gómez-Llanos et al. 2018, for the case of the PN IC 418) with the same shape (spherical) and grain size distribution as when using HAC-like grains. As can be seen from the figure, when using only α-SiC (top panel), the best model spectrum only fits the weak 11.3 μm feature seen in Tc 1 and not the entire broad 12 μm feature. On the other hand, when using α-SiC+graphite (bottom panel), the dust continuum is very well fitted but the red wing of the broad 12 μm emission is poorly reproduced.
Fig. A.1. Spitzer mid-IR spectrum of the PN Tc 1 (black line) compared with the best photoionization model spectrum including α-SiC dust grains only (green line ) are shown (top panel). The C60 (blue line) and C70 (red line) emission bands and UIR plateau emission features are indicated. We note that the narrow emission features seen in the best model spectrum (green line) are just atomic nebular emission lines, whose detailed modeling is out of the scope of the present work. The best photoionization model including α-SiC grains in conjunction with graphite dust grains (green line) is also shown for comparison (bottom panel). |
We remark that it is already well known that much better fits with α-SiC can be obtained when playing with very different extreme shapes and grain size distributions (e.g., the case of PN IC 418; Gómez-Llanos et al. 2018) but here we are exploring if there is any other contribution to or explanation for the broad 12 μm plateau feature in PNe with fullerene-dominated IR spectra such as Tc 1, that is, different from the often assumed SiC.
It is important to mention that the 12 μm feature in IC 418 has different characteristics (position, width, and emission strength) than in Tc 1 (Otsuka et al. 2014). Gómez-Llanos et al. (2018) needed to use very different shapes (spheres and five extreme ellipsoid types) and grain size distributions along the nebula in order to fit the red wing of the 12 μm feature with α-SiC. Also, the IR spectrum of IC 418 is contaminated by other species (e.g., PAHs), displaying weaker C60 bands and a different fullerene and dust temperature than Tc 1.
All Figures
Fig. 1. Spitzer mid-IR spectrum of the PN Tc 1 (black line) compared with the best photoionization model spectrum including HAC-like dust grains (green line). The C60 (vertical dashed blue line) and C70 (vertical dashed red line) emission bands and UIR plateau emission features are indicated. We note that the narrow emission features seen in the best model spectrum (green) are just atomic nebular emission lines, whose detailed modeling is out of the scope of the present work. |
|
In the text |
Fig. 2. Laboratory average values of the complex refractive indices n and k (black and purple lines, respectively) of HAC-like particles at 300 K and with a hydrogen concentration of H/C ≃ 0.35 (Prof. W. Duley, priv. comm.; see text for more details). The errors of the complex refractive indices were measured to be ∼20% (see e.g., Duley 1984). |
|
In the text |
Fig. 3. Spitzer mid-IR spectrum of the PN Tc 1 (black line) compared with the best photoionization model spectra for HAC-like (dashed gray line) and graphite (dashed magenta line) dust grains separately. The best photoionization model including HAC-like grains in conjunction with graphite dust grains (solid green line) is also shown for comparison. |
|
In the text |
Fig. A.1. Spitzer mid-IR spectrum of the PN Tc 1 (black line) compared with the best photoionization model spectrum including α-SiC dust grains only (green line ) are shown (top panel). The C60 (blue line) and C70 (red line) emission bands and UIR plateau emission features are indicated. We note that the narrow emission features seen in the best model spectrum (green line) are just atomic nebular emission lines, whose detailed modeling is out of the scope of the present work. The best photoionization model including α-SiC grains in conjunction with graphite dust grains (green line) is also shown for comparison (bottom panel). |
|
In the text |
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.