domain: phys chim

text:
Clément Dubosq

Paula Pla

Emmanuel Dartois

Aude Simon
email: aude.simon@irsamc.ups-tlse.fr

Spectroscopic investigation of interstellar hydrogenated carbon clusters Infrared and optical spectra of C 24 H n populations

Keywords:

Context. The assignment of the mid-infrared (mid-IR) emission features and plateaus observed in C-rich and H-rich regions of the interstellar medium (ISM) is still debated. Such mid-IR assignments must also be related to their contribution to the extinction curve in our galaxy and to the ultraviolet (UV) bump. Aims. The aim of this work is to investigate the influence of hydrogenation rate on the mid-IR spectra of populations of carbon clusters in order to constraint the nH /nC ratios in regions of the ISM where carbon is an important component. Their potential contribution to the extinction curve and in particular to the UV bump is also investigated. Methods. The absorption IR and optical spectra of tens of thousands of C24Hn (n=0,6,12,18,24) isomers classified into structural families -namely flakes, branched, pretzels, and cages-were computed using the density functional based tight binding electronic structure method and its time-dependent version, respectively. Final spectra were obtained by averaging the spectra of many individual isomers. Results. The shapes and the relative intensities of the bands centered at ∼ 3.25 and 3.40 µm and assigned to the C-H stretch of sp 2 and sp 3 carbon atoms, respectively, present a clear dependence on the nH /nC ratio. From a comparison with the astronomical spectrum from the Orion bar H2S1, the most interesting emitting candidates would pertain to the flakes population; this is the most energetically favorable family of clusters, possessing a high content of five and six carbon rings and being mostly planar, with no sp 1 carbon atoms and with a nH /nC ratio of lower than 0.5. The same conclusion is drawn when comparing the computed IR features in the [4-20 µm] region with the observed plateaus from some C-rich and H-rich planetary nebulae objects of the Small Magellanic Cloud. The contribution of the same family could be considered for the UV bump. When nH /nC increases, only a contribution to the high-energy part of the continuum due to σ → π excitations can reasonably be considered. Conclusions. These results bring some constraints on the structural features and on the nH /nC ratio of the hydrogenated carbon populations emitting in the mid-IR domain in interstellar objects such as protoplanetary and reflection nebulae. The flakes population, with a low nH /nC ratio, is an interesting candidate for the carbon population emitting in these regions, but not for that absorbing in the diffuse ISM. None of the populations studied in the present work can account for the UV bump, but they would contribute to a broad extinction rise in this domain. The computed features reported in this article could be used to interpret future astronomical data provided by the James Webb Space Telescope.

Introduction

Due to its allotropy, a large fraction of interstellar carbon is expected to be included in large molecules or in dust grains with various degrees of organization and aromaticity, including fullerenes, polycyclic aromatic hydrocarbons (PAHs), and amorphous possibly hydrogenated carbon grains (Dartois 2019). PAHs have attracted considerable and sustained interest since the proposal in the mid-1980s, based on their photophysical emission properties, that they are the carriers of the aromatic infrared bands (AIBs), a set of mid-IR emission bands observed at about 3.3, 6.2, 7.7, 8.6, 11.3, and 12.7 µm in many regions of the interstellar medium (ISM) (Allamandola et al. 1985;Léger & Puget 1984). However, it is only very recently that specific PAH molecules (two cyano-naphthalene isomers and indene) have been successfully detected (McGuire et al. 2021;Cernicharo et al. 2021). Some years earlier, the detection of C 60 and C 70 fullerenes in various nebulae was confirmed (Cami et al. 2010;Sellgren et al. 2010). In general, the identification of the specific molecular carriers of the AIBs is challenging, as they are presumably contributed by many unknown species that emit in the same regions of the spectra, therefore leaving many questions still open (Hansen et al. 2022).

Starting in the mid-1980s, the identification of the carriers of the AIBs motivated experimental and theoretical studies of individual molecules. In particular, quantum

Article number, page 1 of 16 chemical calculations, mainly based on density functional theory (DFT) calculations, proved to be very powerful in computing the infrared (IR) features of many individual PAHs of different sizes and shapes (Bauschlicher et al. 2009;Ricca et al. 2012), neutral or charged (Banhatti et al. 2021), with varying hydrogenation degree (Mackie et al. 2015;Pla et al. 2020), and incorporating heteroatoms (Hudgins et al. 2005;Oliveira et al. 2020) or complexed to a metal atom (Simon & Joblin 2010). Such calculations were also used to build spectroscopic databases of PAHs (Malloci et al. 2007a;Bauschlicher et al. 2018). However, the AIBs not only include discrete bands but also broad plateaus, which are a collection of broad blended features whose carriers remain elusive (Shannon et al. 2018). For instance, fullerene-rich planetary nebulae were observed to exhibit a broad plateau with substructure in the 6-9 µm range (García-Hernández et al. 2010, 2011b,a, 2012;Bernard-Salas et al. 2012). Several types of carriers were proposed, including hydrogenated amorphous carbon (HAC) compounds (Bernard-Salas et al. 2012) and very small grains (Tielens 2008;Buss et al. 1993). Very recently, intensive quantum calculations showed that a so-called cages population of C 60 -containing hundreds of thousands of isomers characterized by a low asphericity and a high fraction of sp 2 carbon atoms-was a good candidate to carry the 6-9 µm plateau (Dubosq et al. 2019), while the planar polyaromatic population (PPAs, or flakes) was shown to possess features in the 8-10 µm range. However, the so-far unexplained astronomical features of planetary nebulae in the 10-15 µm range could not be reproduced by any of the studied populations of C 60 , that is, neither the more organized (cages, PPAs or flakes) nor the more amorphous ones (pretzels and branched), addressing the possible contribution of hydrogenated carbon populations.

The 3 µm region of the IR spectra is especially relevant for hydrogenated carbon compounds as it accounts for the C-H stretching features. In many regions of the ISM, the 3.3 µm emission is dominant and associated to aromatic C-H stretching in polyaromatic systems. This feature is often accompanied by weaker emission bands at around 3.4 µm on top of a fainter emission plateau, which was assigned to methylene (-CH 2 ) and methyl (-CH 3 ) aliphatic resonances, but the exact identification of the carriers remains under debate. It was indeed proposed to be due to the emission from highly excited superhydrogenated PAHs (Bernstein et al. 1996;Wagner et al. 2000), which are expected be stable and to form efficiently in the ISM while also having a catalytic role in the formation of H 2 (Thrower et al. 2012). Methylated PAHs were also proposed (Wagner et al. 2000). The analysis of the 3.4 µm band with respect to experiments and computations on pure aliphatic carbon molecules (alkanes) revealed that the -CH 2 /-CH 3 fraction depended on the object and that, in dark clouds, -CH 3 groups were suppressed while the -CH 2 groups remained abundant (Duley et al. 2005). In parallel, a 3.4 µm absorption band with a different profile has been widely observed in the diffuse ISM (Sandford 1991;Godard et al. 2012;Li & Draine 2012). This band was assigned to a hydrogenated amorphous carbon dust component and presents substructures at 3. 38, 3.48, 3.42, and 3.50 µm that were assigned to asymmetric and symmetric C-H stretching modes of -CH 3 and -CH 2 groups, respectively (see Godard et al. (2012) and references therein). A contribution of molecular PAHs with excess hydrogens has also been proposed (Steglich et al. 2013). In the ISM, the stronger the UV radiation field, the lower the 3.4/3.3 µm intensity ratio is, which could be interpreted as a consequence of the photodestruction of the more fragile aliphatic bonds (Pilleri et al. 2015). Interestingly, the effect of cosmic rays on this band was studied experimentally (Godard et al. 2011;Dartois et al. 2017), showing the evolution of the band when the aliphatic C-H bands are destroyed and, through an appropriate modeling, that the destruction of just these grains by cosmic rays cannot account for the disappearance of the 3.4 µm IR feature in denser regions.

Determining the exact nature of the carbon compounds contributing to the ultraviolet (UV) bump (Stecher 1965) has also been of longstanding interest. The UV bump is a broad ultraviolet absorption bump observed on the ISM extinction curves in systems ranging from the Milky Way to high-redshift galaxies, with a position centered at 217.5 nm (5.77 eV or 4.65 µm -1 ) stable to ∼2 % with relatively large variations in width (0.75 to 1.3 µm -1 ) in different lines of sight (Fitzpatrick & Massa 1986, 2007;Fitzpatrick et al. 2019). Interestingly, the variations in the width and shape of the UV bump have been shown to be related to the presence of defects within the sp 2 conjugated carbon network and several theoretical and experimental studies were subsequently carried out to unravel the structural origin of these variations. The contribution of disordered species under the form of HAC (Jones et al. 1990;Gadallah et al. 2011) and soot particles (see Gavilan et al. (2017) and references therein) was proposed. In an attempt to establish correlations between carbon-grain morphology and spectral data, Rotundi et al. (1998) proposed that carbon nanostructures ordered on the micrometer scale could be better candidates to explain the UV bump than graphitic particles. Recent quantum modeling also showed that a population of carbon spherical isomers with a high fraction of conjugated sp 2 atoms, the so-called cages population, could contribute to the astronomical UV bump (Dubosq et al. 2020).

In this work, we aim to investigate the contribution of populations of C 24 H n (n=0,6,12,18,24) encompassing tens of thousands of isomers to emission IR features in the [3-3.5 µm] domain and to the AIB plateaus. Their UV-visible spectral signatures are also investigated to probe their possible contribution to the UV bump. To this end, intensive quantum chemical calculations were performed using an approximate DFT scheme, namely self-consistent charge density functional based tight binding (SCC-DFTB) (Elstner et al. 1998), which is a similar approach to that taken by Dubosq et al. (2019Dubosq et al. ( , 2020) ) to study pure carbon clusters C n=24,42,60 . Indeed, as mentioned above, some astronomical features could not be accounted for by large pure carbon clusters, revealing the need to consider the role of hydrogenation. Therefore, in this article, we present the computed IR and UV-visible spectra of populations of C 24 H n clusters (n=0, 6, 12, 18, 24) that were determined and characterized in our previous work (Pla et al. 2021), as described in section 2.1. The methodology used in this work to determine the individual IR and UV-visible (optical) spectra of hundreds of thousands of isomers is described in section 2.2. The results are presented in section 3 and the astrophysical implications are discussed in section 4.

Methodology

Isomer landscapes of C 24 H n populations

In our previous work, the potential energy surfaces of C 24 H n=0,6,12,18,24 clusters (tens of thousands of isomers per stoichiometry; see Table 1) were investigated using a genetic algorithm (Larsen et al. 2017) in combination with the density functional based tight binding (DFTB) potential (Porezag et al. 1995;Seifert et al. 1996;Elstner et al. 1998). Global exploration was achieved using the zerothorder DFTB (DFTB0) Hamiltonian (Porezag et al. 1995;Seifert et al. 1996) and subsequent local optimization was performed using the second-order (self-consistent charge or SCC-) DFTB with the mio set of parameters (Elstner et al. 1998), and an empirical dispersion correction contribution (Rapacioli et al. 2009). The structural diversity was analyzed using order parameters, in particular the sum of the numbers of five-and six-carbon rings R 5/6 and the asphericity parameter β, which we used to classify the populations of isomers into families (see e.g., the 2D map for C 24 H 12 on top of Fig. 1, along with an example of one isomer of each family for C 24 H n (n=6,12,18, 24) at the bottom of Fig. 1).

The most abundant and lowest-energy population, which we designated the flakes population, is made up of isomers of various shapes and large R 5/6 values (see Table 1). Flakes are characterized by an increasing number of spherical isomers when n H /n C increases (Pla et al. 2021). The second-most-abundant population, designated the branched population, contains the isomers with the smallest R 5/6 values and remains the highest-energy population for all n H /n C ratios. The fraction of the pretzels population, which is made up of spherical isomers with smaller R 5/6 values, is smaller but increases with n H /n C . Finally, the fraction of the cages population (largest R 5/6 values), which was the major fraction and contained the greatest number of stable isomers in the case of C 60 , remains extremely minor (see Table 1). For all C 24 H n=0,6,12,18,24 clusters, the evolution of the carbon-ring size distribution with energy clearly shows a correlation between stability and the number of six-carbon rings. In addition, a result of astrophysical interest, especially regarding the role of carbon species in the gas ionization, is that the computed average values of the ionization potentials of all populations were found to decrease when n H /n C increases, ranging from 7.9 down to 6.4 eV. This trend was correlated to geometric and electronic factors, in particular to carbon hybridization (Pla et al. 2021).

Infrared and optical spectra of C 24 H n populations

In this work, the IR absorption spectra of all independent isomers of C 24 H n=0,6,12,18,24 clusters were computed at the SCC-DFTB level of theory in the double harmonic approximation. The raw vibrational frequencies were scaled to match the best reference data in a two-step procedure: first, the SCC-DFTB frequencies were scaled in order to match DFT results (see Appendix A). Contrary to pure carbon clusters (Dubosq et al. 2019), determining scaling factors for hydrogenated carbon clusters was a challenge because of the resonance of C-H and C-C modes, for which scaling factors are different in the same regions of the spectrum. This leads to some uncertainties in the shape of the presented et al. (2018) for the same functional and a similar basis set as the ones used in the present work. In the absence of clear information about the underlying thermodynamics, uniform populations were chosen to reconstruct the typical IR spectra of the four structural families from the individual contributions of each member. We are aware that such uniform distributions are an oversimplification of the energetics and this is a limitation of our approach. By summing over thousands of individual topologically distinct isomers, the spectra are necessarily broad and cannot account for sharp bands usually assigned to individual abundant molecules using more traditional methods.

Among these structures, a fraction was selected (see numbers in italics in Table 1) for electronic spectroscopy in order to save computational burden. These configurations were selected based on their values of β, and R 5/6 values lying near the average values of the entire sample for this family and within two standard deviations of the corresponding distribution. When the number of configurations meeting this criterion was considered too large with respect to computational time limitations, some were randomly removed. The time-dependent (TD) DFTB method in its linear response formulation derived by Niehaus et al. (2001) was used to calculate vertical excitation spectra for all minima in our SCC-DFTB sample. This method was carefully benchmarked in our previous article (Dubosq et al. 2020) and no additional benchmark was performed for the present study. The electronic spectrum was computed up to about 15 eV for all stoichiometries, which corresponds to the computation of 600 up to 900 roots for C 24 to C 24 H 24 . However, as explained in our previous study, high excitations are likely to involve orbitals such as Rydberg orbitals (Malloci et al. 2007b), which are not present in our theoretical method due to the use of a confined minimal basis set. Such orbitals should not be involved in transition below 8 eV, that is to say that the transitions possibly contributing to the UV bump are not expected to be impacted. From the former benchmark by Dubosq et al. (2020), all computed transition energies are scaled by 0.988.

Results

Infrared spectra

In this section, we discuss the evolution of the IR spectra as a function of the hydrogenation rate, focusing on the most abundant and stable population, which encompasses the flakes isomers. The IR spectra for all hydrogenation rates are reported in Fig. 2 (c), with a zoom in the C-H stretching region in Fig. 2 (a). For the next-highest energy and less abundant family, that is the branched population, the spectra are reported in Fig. 2 (b) and (d). Interestingly, the features for these latter two populations are similar. The IR spectra of the pretzels population are reported in Fig. B.1 (a) and (c). The fraction of the pretzels population is small, except in the case of C 24 H 24 where the energy distribution of the isomers is similar to that of the flakes population. Finally, the IR spectra of the cages isomers, which are much less abundant, are reported in Fig. B.1 (b) and (d). In order to interpret the IR features and their evolution in terms of structural and electronic characteristics, the ratio of C(sp x )H n (x=1-3, n=1-3) patterns with respect to all hydrogenated carbon atoms was determined for all stoichiometries and all families. The results are reported in Fig. 3 for the flakes and branched families, and in Fig. B.2 for the pretzels and cages families, which are expected to be less relevant based on energetic criteria and number of isomers, respectively. The precise assignment of the modes is based on an analysis of a reduced sample of isomers.

C-H stretching modes: The [2.9-3.5 µm] region

The C-H stretching features of the flakes population of C 24 H n (n=6,12,18,24) are reported in Fig. 2 (a). For all hydrogenation states, IR features were found to resonate in three main distinct wavelength (energy) domains: [2.99-3.00 µm], [3.16-3.30 µm] , and [3.37-3.56 µm]. These globally correspond to the resonance of C-H stretches involving mainly sp 1 , sp 2 , and sp 3 carbon atoms, respectively. As can be seen in Fig. 3 largest among all sp 2 carbon atoms. Regarding sp 3 carbon atoms, the C(sp 3 )H 2 fraction, which is much smaller than that of C(sp 2 )H, is the dominant one. The fraction of C(sp 1 )H moieties remains small (below 20%).

A first weak band assigned to C(sp 1 )-H stretch is located at 3.00 µm (3331 cm -1 ) for n=6 and becomes slightly blueshifted up to 2.99 µm (3384 cm -1 ) for n=24 when increasing n H /n C . A weak band is also observed at 3.17 µm (3155 cm -1 ) except for n=6 and its intensity increases with n H /n C . This band was assigned to antisymmetric (as) C-H stretching modes of terminal C(sp 2 )H 2 groups. Two broad bands are observed at higher wavelengths. A first one is observed in the [3.22-3.30 µm] range ([3029-3105 cm -1 ]) for all values of n, with the intensity of this band increasing with n H /n C . The maximum is located at 3.26 µm for n=6,12,18, a blue shoulder appearing for n=18, whereas two maxima appear at 3.24 (3086 cm -1 ) and 3.28 µm (3048 cm -1 ) for n=24. This band includes (C 2 )C(sp 2 )-H stretches but also C(sp 3 )-H 3 as stretches, the fraction of the latter remaining minor with respect to that of the former for all n H /n C ratios and for all families (see Fig. 3 and Fig. B.2). Finally, a broad and structured band is observed in the [3.37-3.56 µm] range ([2805-2965 cm -1 ]) with two maxima at ∼3.39 µm (2947 cm -1 ) and ∼3.45 µm (2896 cm -1 ). This band was assigned, by decreasing order of energy, to (C 2 )C(sp 2 )-H stretches, (C 2 )C(sp 3 )-H 2 as and symmetric (s) stretches, C(sp 3 )-H 3 s stretches, and (C 3 )C(sp 3 )-H stretches.

As can be seen in Fig. 2 (a), the positions of the bands fluctuate when n H /n C increases, but the most striking evolution concerns the relative intensities of the different energy domains, namely [2.99-3.00 µm], [3.16-3.30 µm], and [3.37-3.56 µm], which reflect the fractions of sp 1 , sp 2 , and sp 3 carbon atoms, respectively. This is quantified by the intensity ratios reported in Fig. 4. For all families, the I 2.99-3.00µm /I 2.99-3.56µm ratio was shown to decrease following an exponential law as a function of n H /n C (see Fig. 4 (a)). On the contrary, the I 3.37-3.56µm /I 2.99-3.56µm ratio was shown to increase with n H /n C following logarithmic (flakes and branched) or linear laws (pretzels). The evolution of I 3.16-3.30µm /I 2.99-3.56µm could not be fitted by a simple function. For the flakes and pretzels families, this ratio decreases as a function of n H /n C whereas its behavior is not as clear for the branched families: it increases up to n H /n C =3/4 and is smaller when n H /n C =1. The evolution of the I 3.37-3.56µm /I 3.16-3.30µm ratio was also reported in Fig. 4 (d). This ratio increases with n H /n C for the three families, the major evolution occurring for the branched family, and the minor evolution occurring for the pretzels family. In between, the data points for the flakes family could be fitted by an exponential function.

These trends can be interpreted based on the evolution of the ratio -with respect to the total number of hydrogenated carbons-of C(sp x )H n moieties (x=1-3,n=1-3) as a function of n H /n C (see Figs. 3 andB.2).The C(sp 2 )H fraction remains the largest one, but it decreases when n H /n C increases (84 % down to 44 % for the flakes family) while the C(sp 2 )H 2 fraction increases but remains small (2% to 7 % for the flakes family). This is in line with the trend observed for the evolution of the I 3.16-3.30µm /I 2.99-3.56µm ratio. atoms, the major fractionis that of C(sp 3 )H 2 , which increases when n H /n C increases (from 5 % (n H =6) to 25 % (n H =24). The C(sp 3 )H fraction also increases, although remaining smaller (from 3 % to 13 % for the flakes family), as does the C(sp 3 )H 3 fraction, which remains a relatively minor fraction (0.1 % to 5 % for the flakes family).

n H /n C C sp 1H 1 C sp 2H 1 C sp 2H 2 C sp 3H 1 C sp 3H 2 C sp 3H 3
These results are in line with the increase in the I 3.37-3.56µm /I 2.99-3.56µm ratio and also contribute to explaining the decreases in I 2.99-3.00µm /I 2.99-3.56µm and I 3.16-3.30µm /I 2.99-3.56µm when n H /n C increases. Although the trends are relatively similar for all families, we note some differences. The C(sp 1 )H fraction remains small for all stoichiometries (between 5 % and 7 % for the flakes families depending on the stoichiometry), indicating a small number of terminal C-H bonds. However, it is larger for the branched (15 % down to 5 %) and pretzels families (14 % down to 8 %) reflecting the larger I 2.99-3.00µm /I 2.99-3.56µm ratio (Fig. 4 (a)). Interestingly, the C(sp 3 )H fraction is more important in the cages family than in the other families (from 9 % to 29 %) because of the spherical shape of the clusters. Globally, the increase in the fraction of sp 3 carbon atoms to the detriment of sp 2 carbon atoms when n H /n C increases is illustrated by the evolution of the I 3.37-3.56µm /I 3.16-3.30µm ratio reported in Fig. 4 (d). Fig. 4: Intensity ratios I 2.99-3.00µm /I 2.99-3.56µm (a), I 3.37-3.56µm /I 2.99-3.56µm (b), I 3.16-3.30µm /I 2.99-3.56µm (c), and I 3.37-3.56µm /I 3.16-3.30µm (d) as a function of the n H /n C fraction for the flakes family (black crosses), branched family (red squares), and pretzels family (blue triangles). The dashed lines either result from the analytical data fitting (when starting below x=0.25; see Table B1 in Appendix B for the coefficients) or were drawn to guide the eye.

C-C and C-H modes: The [4-30 µm] region

The IR features of the families of pure C n carbon clusters (n=24,42,60) were described previously by Dubosq et al. (2019). Focusing on the flakes population of C 24 (Fig. 2), an intense narrow band can be found in the [4.47-5.14 µm] region. This is due to the resonance of the C-C stretch of sp 1 carbon atoms. The intensity of this band clearly decreases when n H /n C increases, reflecting the decrease in the fraction of sp 1 carbon atoms. For all families except that of the cages, this band remains intense up to n H /n C =0.5 and nearly disappears for n H /n C =1, while its intensity decreases for n H /n C =0.75. The maximum of this band is found at 4.80 µm for C 24 (flakes family), at 4.77 µm for C 24 H 6 and C 24 H 18 , and at 4.78 µm for C 24 H 12 and C 24 H 24 . A wide IR band is observed in the [5.4-12.8 µm] wavelength range for C 24 and its shape depends on the family. Globally, the [5.4-10.2 µm] wavelength range features are assigned to the C-C stretches of sp 2 carbon atoms and also of sp 3 carbon atoms at the longest wavelengths. When decreasing the energy up to ∼ 12.8 µm, softer CCC in-plane bending modes are found. Above 15 µm, more collective deformation modes of the carbon skeleton occur.

As mentioned in section 2, the difficulty in addressing the spectra of hydrogenated carbon clusters is caused by a difference of scaling between the C-H bending modes and some CC stretches and CCC bends that occur in the same energy domains; that is, above ∼ 10 µm. In this region, C-H bending modes are expected to be dominant, especially when n H /n C increases, and therefore we used the scaling law for the C-H bends determined in the Appendix A for the [10.2-22.7 µm] wavelength range (CC scaling), the band in the latter domain being finally shifted to the [10.0-15.3 µm] domain (see Fig. 2 and Fig. B.1). In the remainder of this paper, the results are presented using this final wavelength domain.

When n H /n C increases, the relative intensity in the [5.4-10.2 µm] range decreases, whereas that in the [10.0-15.3 µm] region increases for all families. As can be seen in Fig. 5 (a), the I 10.0-15.2µm /I 5.4-10.3µm ratio globally increases with n H /n C , reaching similar values (between 0.4 and 0.45) for the flakes, branched, and pretzels families when n H /n C =1. The features in the [5.4-10.2 µm] wavelength range are due to pure C-C stretches at the smaller wavelength but when the wavelength increases, the resonant vibrational modes include combinations of C-C stretches and C-H bends. The resonance of C-H bends appears above 10.0 µm. Regarding the flakes family, the maximum of the band in the [10.0 -15.2 µm] region is found at 12.4, 12.3, 12.1, and 11.8 µm from n H =6 to 24. This band appears globally slightly blueshifted, with a blue tail starting at 11.2 µm for n H =6 and at 10.8 µm for the higher hydrogen ratios; its shape is also modified when n H /n C increases: its lower-energy part becomes more intense to the detriment of the higher-energy part, which contains a higher fraction of features due to CC modes. Because of the difference in scaling between C-H and CC modes in this region, the profile of the band has to be considered with caution. However, it is clear that the presence of C-H bonds leads to new features in the [11.0 -13 µm] domain with respect to pure carbon clusters. As in the case of pure carbon clusters, less intense and delocalized soft modes occur at higher energies. Interestingly, we notice Article number, page 7 of 16 a slight bump at 14.8 µm for n H /n C =0.5 and 0.75. From our analysis, this appears to be due to the C-H bend of terminal C(sp 1 )-H bonds. As could have been expected, this bump is more intense in the case of the branched family for all stoichiometries (see Fig. 2 (d)). In the case of the branched families, hydrogenation also clearly leads to new features in the [10.0 -15 µm], range with more intensity in the lower energy domain (maximum at 12.4, 12.1, 11.9 and 11.8 µm for n H from 6 to 24, respectively). The red tail of the band finishes by a small bump at 14.8 µm. Contrary to the flakes family, a small bump is observed at 9.5, 9.3, and 9.3 µm for n H from 12 to 24. The evolution of the spectra for the pretzels families is similar to that of the branched families (see Fig. B.1 (c)) although the features due to the C-H bending modes start at a higher wavelength (above 11.0 µm and the broad band in the [11.0 -15 µm] region has a maximum at ∼ 12.5, 12.3, 12.1 and 12.1 µm for n H from 6 to 24). The features for the cages families are reported in Fig. B.1 (d). However, due to the lack of isomers, these cannot be considered to interpret the astronomical plateaus.

UV-visible spectra

The computed UV-visible spectra for C 24 H n (n=0, 6,12,18,24) are reported in Fig. 6 for the flakes and branched populations, and in Fig. C.1 for the pretzels and cages populations. Regarding the flakes population, all spectra are characterized by a wide plateau starting at ∼ 1 eV for C 24 and at higher energy, up to ∼ 2 eV, for C 24 H 24 when n H /n C increases. The occurrence of such low-energy transitions is expected for unsaturated carbon clusters such as C 24 . C n (n=24,42,60) clusters were indeed shown to have their lowest energy transitions below 1 eV and a high density of lowenergy electronic states, all the more as the cluster size increases (Lacinbala et al. 2022). These transitions correspond to the promotion of the higher energy π (delocalized) electrons to the lower energy π orbitals. When n H /n C increases, the fraction of electrons in covalent σ bonds increases to the detriment of π bonds, leading to relatively less intense features in the lower-energy part of the spectrum.

The shape of the plateau depends on the hydrogenation degree. In the case of C 24 , it presents two bumps centered at ∼ 4.75 eV and ∼ 8.40 eV. Although our method of calculation does not allow us to trust the spectrum above 10 eV, it appears that transitions are still expected at higher energies and that the spectrum resembles a continuum. Interestingly, good agreement is found between this spectrum and that of the C 24 population obtained from a set of structures determined using a replica exchange molecular dynamic (REMD) algorithm (Dubosq et al. 2020).

In the case of C 24 H 6 , this structured plateau is replaced by one single wide band starting at ∼ 1.7 eV with a maximum at ∼ 6.8 eV. As in the case of C 24 but to a lesser extent, transitions are expected at energies higher than 10 eV and the blue tail of the bump slopes gently downwards. For C 24 H 12 , the bump starts at the same energy as for C 24 H 6 with a redshift of the maximum down to ∼ 6 eV. This could be expected, as the maximum in the average spectrum of 40 PAHs (below 10 eV) is found at ∼ 6 eV Malloci et al. (2007b). Its shape is more symmetric than in the case of C 24 H 6 , with a steeper blue tail. Regarding the overhydrogenated molecules C 24 H 18 and C 24 H 24 , their UVvisible spectra are very similar and present a bump shifted towards higher energies with a maximum of ∼ 7.4 eV for both clusters.

Article number, page 8 of 16

The electronic transitions below 9 eV involve π→π transitions, mixed with σ→π transitions when the energy increases (see some examples in Fig. C.2). The blueshift of the maximum for C 24 H 18 and C 24 H 24 with respect to C 24 H 12 is consistent with the more important fraction of sp 3 C atoms leading to σ→π transitions. Transitions to σ orbitals are found at higher energies (above 9 eV). Interestingly, for the flakes population, when integrating the spectra between 2 and 9 eV, a clear decrease in the oscillator strength can be seen in this region when increasing n H /n C (see Fig. 7, black crosses).

Regarding the branched families, all spectra are also characterized by a wide plateau that, as in the case of the flakes family, appears at ∼ 1 eV for C 24 and at higher energies up to ∼ 2 eV for C 24 H 24 when n H /n C increases. The difference in the shapes of the plateaus with n H /n C is less obvious than in the case of the flakes. As in the case of the flakes family, the 2-9 eV oscillator strength decreases when n H /n C increases. The oscillator strength in the higher-energy part of the spectrum (6-9 eV) increases, whereas that in the 2-6 eV domain decreases when n H /n C increases, with a maximum at approximately 8.2 eV up to n H /n C =0.5, and 7.8 and 7.5 eV for n H /n C =0.75 and 1, respectively.

Regarding the pretzels family, the trends are similar (see Fig . C.1 (a) and Fig. 7). The maximum of the plateau is at approximately 8.3 eV for C 24 , 7.6 eV for C 24 H 6 and C 24 H 12 , and at a lower energy for C 24 H 18 and C 24 H 24 . The evolution of the UV-visible spectra of the cages (see Fig. C.1 (b)) is not discussed as they cannot be considered as a population due to the small number of isomers.

Astrophysical implications

The astronomical IR spectra reported in the present work are emission spectra and the computed spectra are absorption spectra. Simulating a full radiative emission cascade approach for comparison is out of the scope of this article. To first order, for the relatively large systems studied here, the emitted spectral positions should in any case remain within the error bars of the scaling of our absorption spectra.

Regarding the C-H stretching region, the error bars on the scaling factors may amount to 0.1 µm due to the high concentration of data on a small wavelength range. However, the relative intensities in the three different wavelength regions assigned to sp n (n=1-3) carbon atoms could be a good indicator of the structural characteristics and of the n H /n C ratio of the IR band carriers in the C-H region in the ISM (see Fig. 4). For instance, comparing the spectra computed in this work with the astronomical spectrum of the Orion bar (a so-called class A in the AIBs astronomical spectral classification; i.e., the more aromatic AIBs class), it appears that isomers with the characteristics of the flakes family (large number of five-and six-carbon cycles and mostly planar) with a n H /n C ratio smaller than 0.5 with no sp 1 carbon atom, are the most interesting candidates.

In the view of interpreting future astronomical observations, the I 3.37-3.56µm /I 3.16-3.30µm intensity ratio reported in Fig. 4 (d) could be used to bring constraints on the n H /n C fraction in regions where hydrogenated carbons are expected. From the dependence of the C(sp x )H n (x=1-3, n=1-3) fraction on n H /n C reported in Fig. 3, we derived the C(sp 3 )H 3 /C(sp 3 )H 2 fraction as a function of n H /n C for the flakes, branched, and pretzels families (see Fig. 8 (a)). As can be seen in Fig. 8 (a), this fraction remains small for all families (lower than 0.3) but increases with n H /n C . Interestingly, the shape of the [3.3-3.6 µm] band for the flakes and a n H /n C ratio of 0.5 is very similar to the astronomical features recorded for the protoplanetary nebula CRL 618 (Chiar et al. 1998). Chiar et al. (1998) deduced from their analysis that the 3.4 µm band was due to the presence of relatively short chains of aliphatic hydrocarbons, with a CH 3 /CH 2 fraction estimated at between 0.4 and 0.5, which is larger than our estimation for low n H /n C fractions.

Article number, page 9 of 16 In the diffuse ISM, the C(sp 3 )H 3 /C(sp 3 )H 2 fraction was estimated to be in the range of 2-2.5 (Matrajt et al. 2005;Sandford et al. 1991;Pendleton et al. 1994), and therefore the populations considered in this work are not good candidates for the HAC in these regions. On the contrary, the C(sp 3 )H 1 /C(sp 3 )H 2 ratio is larger and remains above 0.4 (see Fig. 8 (b)). For the lower values of n H /n C , the smallest C(sp 3 )H 1 /C(sp 3 )H 2 ratio is obtained for the flakes populations.

Regarding the IR features above 4 µm, the sharp band at ∼4.77 µm due to C≡C stretch is not observed in astronomical spectra. For instance, similarly to the 3.00 µm C(sp 1 -H) stretching band, the 4.77 µm band is not observed in the astronomical emission spectrum of Orion bar H2S1. Therefore, among all the isomers considered here, those possessing C≡C bonds are not observed in astronomical emission spectra. If they are formed in the ISM, then their lifetime is expected to be insufficient to be detected. The reason that can be invoked is that C≡C bonds are expected to be more reactive than C=C bonds. Consequently, the isomers possessing sp 1 carbon atoms are expected to isomerize under the action of either photons or collisions, leading to most stable isomers possessing a majority of sp 2 atoms. The isomers containing sp 1 carbon atoms could also be reactive intermediates in the growth of larger carbon clusters considering a bottom-up process. A shortcoming of the present work is that all isomers were taken into account without any energetic consideration, and we show that the higher the fraction of sp 1 atoms, the less stable the isomers are (see Fig. 10 in the article by Pla et al. (2021)). The intensity ratios reported in Fig. 5 could be used to estimate the n H /n C ratio in regions of the ISM where the presence of hydrogenated carbon clusters can be expected. It is clear that I 10.0-15.2µm /I 5.4-10.3µm (Fig. 5 (a)) and I 10.0-15.2µm /I 4.47-10.3µm (Fig. 5 (b)) increase with n H /n C for all families, leading to similar values of n H /n C =1 -namely 0.40-0.45 for the former and 0.3-0.4 for the latter-for all families. However, these evolutions could not be fitted by a simple law. The I 5.41-10.3µm /I 4.47-5.41µm and could also be used to constrain the hydrogen rate and to a lesser extent the structural features of astrophysical candidates. Interestingly, the I 5.41-10.3µm /I 4.47-5.41µm ratio increase for the pretzels families is less pronounced than the same increase for the branched and flakes families, which appear similar.

As the computed spectra are absorption spectra, and the profile of the band in the [10.0 -15 µm] range remains uncertain due to the difference of scaling factors for C-C and C-H modes, only trends can be compared to astronomical data. However, we compared the computed IR features to the astronomical emission spectrum recorded in a C-rich and H-rich planetary nebula of the Small Magellanic Cloud (SMC), where both fullerenes and PAHs have been detected (García-Hernández et al. 2010). Interestingly, it appears that the plateaus observed in the [6-10] and [10 -15 µm] ranges could be accounted for by the resonance of the C-C stretches and C-H bends, respectively. The best agreement in terms of the shape of the bands and their intensity ratio qualitatively corresponds to the flakes family populations, which possess a n H /n C rate below 0.5. This conclusion is similar to the one drawn after the analysis of the C-H stretching features in the Orion bar Article number, page 10 of 16 H2S1. However, we are aware that the relative intensities of these regions should be in favor of the lower energy range in astronomical emission spectra.

This study, based on the analysis of the relative intensities of the IR features in different domains, brings some constraints on the chemical nature of hydrogenated carbon clusters (structural organization and hybridization of carbon atoms) and on the n H /n C ratio of these clusters in the regions of the ISM where they are detected. Our results are clearly in favor of the flakes population cleared from sp 1 atoms with low n H /n C ratio. It is interesting to compare these trends with the experimental results obtained by Dartois, E. et al. (2020). In this study, interstellar-relevant dust particles were produced by high-energy mechanochemical synthesis from carbon-based solids (graphene flakes, fullerenes, nanotubes) under hydrogen atmosphere. From the experimental measurement of the n H /n C ratio and from the comparison of the IR-measured intensities of the produced grains and the AIB features of the so-called class C-D spectral family, which contains aromatic and aliphatic moieties, these latter authors concluded that the n H /n C ratio within polyaromatic-emitting structures should be very low, between 3 % and 7 %, which is one order of magnitude lower than the estimations provided in the present work. However, the size of the produced species in these experiments is much larger than that of the carbonaceous systems studied in this work. Regarding pure carbon clusters, quantum chemical calculations show that when size increases, the cages population becomes energetically more favorable than the flakes population (Dubosq et al. 2019). However, interestingly, a recent theoretical study by Perez-Mellor et al. (2022) shows that both entropy and addition of a moderate amount of hydrogen favors planar structures analogous to the flakes populations studied in this work. Typically, regarding the hydrogenated isomers of C 60 , the flakes population was determined as the most favorable for C 60 H n with 12 ≤ n ≤ 60 stoichiometries for T between 0 and 3500 K. Therefore, we may reasonably assume that the IR features for larger clusters such as C 60 H n with 12 ≤ n ≤ 60 are similar to the ones obtained for C 24 H 6 . However, for a hydrogen fraction lower than 0.20, the most stable population for C 60 H n is roughly that of the cage, which is also the most stable but encompasses fewer isomers for carbon clusters with less than 30 atoms. It is therefore hazardous to extrapolate our results to clusters of a larger size with a larger number of carbon atoms and a lower n H /n C ratio than the one considered in the present study.

When comparing the computed UV-visible spectra of the flakes families (this work) to the average interstellar medium extinction curve (Fitzpatrick & Massa 2007), it appears that the C 24 H 12 family is the one contributing the most to the UV bump, and C 24 H 6 is that contributing the least, its maximum being blueshifted with respect to the astronomical data. Furthermore, larger pure carbon clusters pertaining to the cages -and to a lesser extent the flakes-family were shown to be good candidates (Dubosq et al. 2020). When increasing the hydrogenation degree, only a contribution to the high-energy part of the absorption continuum can be considered. As for the spectrum of C 24 , it could also contribute to the continuum considering carbon-rich regions, but not to the UV bump, and this is also the case for the branched and pretzels populations.

Conclusions

In this work, we report the IR and optical absorption spectra of populations of C 24 H n (n=0,6,12,18,24) clusters containing up to several tens of thousands of isomers. These were obtained by summing individual DFTB (IR) and TD-DFTB (optical) spectra of individual isomers. We also present correlations between spectral features and the n H /n C ratio for different structural families.

Regarding the IR features, the computed bands in the [2.99-3.00 µm], [3.16-3.30 µm], and [3.37-3.56 µm] domains are assigned to C(sp x )-H stretches with n=1, 2, and 3, respectively. The intensity ratios within these three domains could be an indicator of the n H /n C fractions in the regions of the ISM where the presence of carbon species is expected. The comparison of the computed spectra with the Orion bar H2S1 spectrum suggests that the flakes population, which is deprived of sp 1 atoms with an n H /n C lower than 0.5, could be at the origin of these IR features. However, error bars due to the limited wavelength range where scaling factors were determined are expected on the absolute positions of the C-H stretching modes. Hydrogenation of carbon clusters also leads to a plateau in the [10-15 µm] region -which is more intense as n H /n C increases-, where no feature is observed for pure carbon clusters.We clearly show that the ratio between the intensities in the [10.0-15.2 µm] domain with respect to the [5.41-10.3 µm] increases with n H /n C . This could also be used to analyse future astronomical IR data to provide a constraint on the n H /n C fraction. The comparison between the computed IR features in the [10-15 µm] and [6-9 µm] regions and the observational data of a C-rich and H-rich planetary nebula indicates that hydrogenated flakes with a low n H /n C ratio (lower than 0.5) are good candidates to be the carriers of these plateaus.

However, an intense narrow band resonating at ∼ 4.8µm and due to C≡C stretch is observed in all calculated spectra, with its relative intensity decreasing when n H /n C increases. To our knowledge, such a feature has not been observed in the astronomical IR spectra of the ISM in C-rich and H-rich regions. However, in emission spectra, the intensity of this band is expected to be considerably decreased. As discussed in this article, sp 1 carbon atoms are expected to be more reactive than sp 2 carbon atoms and are expected not to be stable over time. A question that will be addressed in the future is the evolution of the spectra of populations with a low n H /n C ratio when only the most stable isomers -which also possess the fewest sp 1 carbon atoms-are taken into account in order to improve the agreement with astronomical features. In the present work, neither energetic constraints nor reaction processes were taken into account. The results reported here could be of interest for the interpretation of new JWST data, in particular in the framework of the PDRs4All program.

Regarding the optical features, a wide band in the case of the flakes and pretzels families, and a sort of continuum in the case of the branched family, are observed in the [1-9 eV] range. These were assigned to π, σ → π transitions. The total intensity decreases when n H /n C increases. Comparison with the astronomical spectrum Article number, page 11 of 16 suggests that the most stable C 24 H 12 flakes family (and to a lesser extent the C 24 H 6 flakes family) could contribute to the UV bump. When n H /n C increases, more higher-energy features are observed with an increasing contribution of σ → π transitions and these could contribute to the high-energy part of the continuum of the extinction curve.

Appendix A: IR spectra benchmark

Benchmarking was performed against well-known DFT harmonic spectra presented in the literature. This was done for a set of 15 isomers of C 24 H 1 , 13 isomers of C 24 H 2 , and 15 isomers of C 24 H 12 . All of these geometries were relaxed at the B3LYP/pc-1 level with the Gaussian 09 software suite. Corresponding IR spectra are computed at the corresponding level of theory. To this set of structures, we added the results of a previous benchmark performed over a set of 50 C 24 isomers (Dubosq et al. 2019). The benchmarking was performed as follows. For a given isomer, for each nondegenerate C-H vibrational mode and most remarkable C-C vibrational modes at the B3LYP/pc-1 level, we identify the same modes at the SCC-DFTB level. Doing this, we can correlate the vibrational energy computed at B3LYP/pc-1 level with the corresponding SCC-DFTB value (see Fig. For soft and collective C-C vibrational modes (ν < 1750 cm -1 ), we determined a scaling factor of 0.9046 ± 0.0021, while for harder C-C stretching modes (ν > 1750 cm -1 ), the scaling factor was found to be 0.9983±0.0020. A distinction between soft and hard C-C vibrational modes occurs around 1750 cm -1 with an overlap of ∼ 100 cm -1 around this value. To allow a smooth transition between these two trends, we proposed the following equation to express the scaling factor depending of the vibrational energy:

Scaling f actor = A erf x -1750 α + B ,
where constant A and B are determined by asymptotic limit. The α value is defined to have a 75% scaling factor of the soft C-C vibrational mode at 1700 cm -1 and a 75% scaling factor of the hard C-C vibrational mode at 1800 cm -1 .Values of these three parameters are given below:

A = 0.04685 ; B = 1.9029 2 × A ; α = 104.835808.
Concerning the C-H vibrational modes, three energy regions corresponding to three different types of vibrational modes were distinguished and fitted with affine functions a × x + b. The first range of C-H vibrational modes (500 -850 cm -1 ) corresponds to out-of-plane C-H vibration (affine law : 0.594 ± 0.063 × x + 376 ± 45 cm -1 ). The second range of C-H vibrational modes [2440 -2990 cm -1 ] is associated to C-H stretching modes of -CH 2 type. For such modes, the 0.581 ± 0.016 × x + 1346 ± 47 cm -1 affine law was determined. These modes were mostlyobserved in our set for the C 24 H 12 . The last vibration modes correspond to C-H stretching modes for C-H sp 2 and sp 1 atoms (ν > 2950 cm -1 ) and the 0.843 ± 0.016 × x + 627 ± 48 cm -1 was obtained. Due to the difficulty in distinguishing C-H out-of-plane vibrational modes from C-C soft modes of vibration from a pure energetic point of view, we do not take into account the affine law determined for the C-H out-of-plane vibrational mode at first to plot the spectra, but the law is applied to the [450-1000 cm -1 ] region afterwards, without proposing a smooth transition between the two regions, as this seemed beyond our capabilities within the scope of this work, where a tremendous number of isomers and modes are involved.

However, due to the clear distinction we have for the highest-energy C-C and C-H vibrational modes, we were able to propose a smooth transition of the scaling between the two types of C-H stretching modes, and we propose the following form: 6,12,18,24) for the pretzels (a) and cages (b) families. The average Galactic extinction curve adapted from Fitzpatrick & Massa (2007) is reported in dashed lines. Article number, page 16 of 16

Fig. 2 :

 2 Fig. 2: Infrared spectra for the flakes and branched families and comparison with astronomical spectra. Panels (a) and (b) show a zoom into the [2.8 -4 µm] region. Panels (c) and (d) show the features represented with dashed lines in the [10.2 -22.7 µm] region rescaled using the appropriate factors for C-H bend (see main text and Appendix A). The final shifted band in the [10.0 -15.3 µm] domain is represented in plain lines. The emission spectra from the Orion bar H2S1(Van Diedenhoven et al. 2004) and from the SMP SMC 16 planetary nebula (García-Hernández et al. 2010) are added for comparison. Error bars are reported in panels (a) and (c) for the most intense bands of the IR spectrum of C 24 H 6 flakes only, but can be applied to similar modes for all stoichiometries and families. Article number, page 5 of 16

Fig. 3 :

 3 Fig. 3: Evolution of the C(sp x )H n (x=1-3, n=1-3) fraction as a function of n H /n C for the flakes (a) and branched (b) families

  Fig. 5: Intensity ratios I 10.0-15.2µm /I 5.4-10.3µm (a), I 10.0-15.2µm /I 4.47-10.3µm (b), and I 5.41-10.3µm /I 4.47-5.4µm (c) as a function of the n H /n C fraction for the flakes family (black crosses), branched family (red squares), and pretzels family (blue triangles). The dashed lines in plots (a) and (b) were drawn to guide the eye, whereas in plot (c), these lines show the exponential fitting of the data (see Appendix B for the details of the fit).

  Fig. 6: TD-DFTB optical spectra of C 24 H n (n=0, 6,12,18,24) for the flakes (a) and branched (b) families. The average Galactic extinction curve adapted from Fitzpatrick & Massa (2007) is reported in dashed lines.

Fig. 7 :

 7 Fig.7: Oscillator strength ratio in the [2-9 eV] region between C 24 H n (n=6,12,18,24) and C 24 for the flakes (black crosses), branched (red squares), and pretzels (blue triangles) families. The data points were fitted by a linear law (see Appendix B).

Fig. 8 :

 8 Fig. 8: Evolution of the C(sp 3 )H 3 /C(sp 3 )H 2 (a) and C(sp 3 )H 1 /C(sp 3 )H 2 (b) fraction as a function of n H /n C for the flakes, branched, and pretzels families

  Fig. A.1: Correlation graph between SCC-DFTB and B3LYP/pc-1 wavenumbers. The black dotted line corresponds to a perfect correspondence between SCC-DFTB and B3LYP/pc-1 vibrational energy.

  Fig. B.1: Infrared spectra for the pretzels and cages families. (a), (b) Zoom into the [2.8 -4 µm] region. The intensities were multiplied by 4 for C 24 H 6 and by 2 for C 24 H 12 to make the structures of the bands visible. (c), (d) The features represented with dashed lines in the [10.2 -22.7 µm] region were rescaled using the appropriate factors for C-H bend (see main text and Appendix A). The final shifted band in the [10.0 -15.3 µm] domain is represented in plain lines. The emission spectra from the Orion bar H2S1(Van Diedenhoven et al. 2004) and from the SMP SMC 16 planetary nebula (García-Hernández et al. 2010) are added for comparison.

  (a) : f (x) = a0 * exp(-b0 * x), Fig. 4 (b) : f (x) = a1 * log(b1 * x) for flakes and branched, f (x) = a1 + b1 * x for pretzels, Fig. 4 (d): f (x) = a * exp(b * x) for flakes with a=0.39 ± 0.03 and b=2.21 ± 0.09 Fig. 5 (c) : f (x) = a2 * exp(b2 * x), Fig. 7: f (x) = a3 + b3 * x

Fig. B. 2 :Fig

 2 Fig. B.2: Evolution of the C(sp x )H n (x=1-3, n=1-3) ratio as a function of n H /n C for the pretzels (a) and cages (b) families.

Fig

  Fig. C.2: Representation of some relevant transitions for some isomers of the flakes families of C 24 H n (n=6,12,18,24).

Cage Flake Pretzel Branched C 24 H 6 C 24 H 12 C 24 H 18 C 24 H 24

		R 5/6	β	C 24	C 24 H 6		C 24 H 12		C 24 H 18	C 24 H 24
	flakes	3-8	1 0.3-1.0 9449 (62%) 49989 (73%) 41264 (79%) 17889 (69%) 4924 (45%)	10 3
				9049	8175			8446		4522	4556
	branched	0-2	0.8 0.3-1.0 5626 (37%) 18191 (26%) 9722 (19%)	5578 (21%) 2889 (27%)	10 2
				5242	4842			8833		5105	2701
						0.6
	pretzels	0-8	0.0-0.3	95(0.6%)	β 704 (1 %)		1224 (2.3%)	2293 (9%)	3009 (28%)	10 1
				94	667	0.4		1154		2160	2866
	cage	≥ 9 0.0-1.0	3 3	12 10	0.2		5 3		259 (1%) 249	72 72	10 0
	total			15173	68896		52213		26019	10894
						0	0	2	4	6	8	10	12	14	16	10 -1
											R 5/6
						cages			flakes		pretzels	branched
					Fig. 1: Adapted from Pla et al. (2021). Two-dimensional
					distribution of the isomer population for C 24 H 12 as a func-
					tion of R 5/6 and β (top) and examples of isomers of C 24 H n
					for the four families (bottom).
					spectra in the [10-16 µm] region. Second, the DFT frequen-
					cies were scaled using the 0.9731 scaling factor previously
											Article number, page 3 of 16

Table 1 :

 1 Total number of structures per family and in total for each C 24 H n system taken into account to compute the IR (roman) and optical spectra (italics). The percentage of structures of that family of the total number of structures for a particular C 24 H n system is reported in parentheses. R 5/6 and β give, for each family, the related to the sum of the numbers of five-and six-carbon rings and the asphericity parameter, respectively.

	determined in the [4-30 µm] region (Dubosq et al. 2019).
	Due to the larger anharmonicity of C-H stretching modes,
	these were scaled by 0.961 as determined by Bauschlicher

  for the flakes and branched populations, and in Fig. B.2 for the pretzels and cages populations, the fraction of C(sp 2 )H moiety is by far the

					Wavenumber (cm -1 )								Wavenumber (cm -1 )
			3500 3400 3300 3200 3100 3000 2900 2800 2700 2600 2500			3500 3400 3300 3200 3100 3000 2900 2800 2700 2600 2500
		20	(a) Flakes			C(sp 3 )-H 1-3				20	C(sp 1 )-H (b) Branched	C(sp 2 )-H	C(sp 3 )-H 1-3
					C(sp 2 )-H							C 24 H 24
			C(sp 1 )-H
		15	C 24 H 24									15
	Intensity (km.mol -1 /cm -1 )	5 10	C 24 H 12 C 24 H 18						x2		Intensity (km.mol -1 /cm -1 )	5 10	C 24 H 12 C 24 H 18		x6 x4
													C 24 H 6		x6
			C 24 H 6						x4
		0	Orion bar H2S1						5	Fv (Jy)	0	Orion bar H2S1	5	Fv (Jy)
										0					0
		2.8	3	3.2		3.4	3.6	3.8	4.0		2.8	3	3.2	3.4	3.6	3.8	4.0
					Wavelength (µm)								Wavelength (µm)
					Wavenumber (cm -1 )
			3000	2000 1500		1000		500
										0.02
	Fv (Jy) Orion bar	100 200	(c) Flakes Orion bar H2S1 SMC 16					0 0.01	Fv (Jy) SMC16
		0
			C 24 H 24
		8										8
	)		C 24 H 18										C 24 H 18
	-1 -1 /cm	6										6
	Intensity (km.mol	4	C 24 H 6 C 24 H 12									4	C 24 H 6 C 24 H 12		x0.7
		2										2
			C 24										C 24		x0.5
		0										0
			3	4	5 6	8	10 12 15	20 25 30					3	4	5 6	8	10 12 15	20 25 30
					Wavelength (µm)
	Article number, page 4 of 16

  Regarding the sp 3 carbon

		80 100	(a) flakes				C sp 1H 1 C sp 2H 1 C sp 2H 2 C sp 3H 1 C sp 3H 2 C sp 3H 3
		60
	%C sp xH n	40
		20
		0
			0.2	0.3	0.4	0.5	0.6	0.7	0.8	0.9	1
							n H /n C
		100	(b) branched
		80
		60
	%C sp xH n	40
		20
		0
			0.2	0.3	0.4	0.5	0.6	0.7	0.8	0.9	1

  Article number, page 14 of 16 Clément Dubosq et al.: Spectroscopic investigation of interstellar hydrogenated carbon clusters Table B.1: Coefficients of the analytical fitting equations for the intensity ratios as a function of n H /n C for Fig. 4

	family	a0	b0	a1	b1
	flakes	0.28 ± 0.02	3.49 ± 0.28	0.32 ± 0.03	10.0 ± 1.6
	branched	0.66 ± 0.06	3.15 ± 0.26	0.23 ± 0.02	19.0 ± 3.0
	pretzel	0.65 ± 0.08	3.98 ± 0.41	0.20 ± 0.01	0.62 ± 0.02
		a2	b2	a3	b3
	flakes	0.14 ± 0.02	3.74 ± 0.15	0.98 ± 0.02	-0.53 ± 0.04
	branched	0.056 ± 0.020	3.73 ± 0.39	0.999 ± 0.001	-0.489 ± 0.001
	pretzel	0.056 ± 0.010	4.90 ± 0.19	1.03 ± 0.03	-0.61 ± 0.04

Article number, page 2 of 16

Acknowledgements. The authors gratefully acknowledge financial support by the Agence Nationale de la Recherche (ANR) Grant No. ANR-16-CE29-0025 and the Spanish MICINN (PID2019-110091GB-I00). They also thank the computing mesocenter CALMIP ("CALcul en MIdi Pyrénées", UMS CNRS 3667) for generous allocation of computer resources (p0059 project).

Appendix B: IR features

Article number, page 13 of 16


Johannes Scheller
email: johannes.scheller@laplace.univ-tlse.fr

Karl-Joseph Rizzo
email: kjrizzo@laplace.univ-tlse.fr

Eric Duhayon
email: eric.duhayon@laplace.univ-tlse.fr

Jean-François Rouchon
email: rouchon@laplace.univ-tlse.fr

Marianna Braza
email: marianna.braza@imft.fr

Design and development of a hybrid NACA 4412 morphing airfoil

The present paper will describe the design and development of a hybrid NACA 4412 morphing airfoil equipped with both Shape Memory Alloy Actuators providing large deformations over a limited frequency range as well as Macro Fiber Composite Actuators providing small deformations with a large bandwidth. The design choices will be highlighted, the electrical and mechanical actuation performances characterized and finally the static and dynamic effects on the flow will be evaluated and compared to the unactuated airfoil.

Introduction

Today's wing geometries are usually a trade-off between the optimal shapes at different phases of the flight. Control surfaces, while modifying the aerodynamic characteristics of the airfoil, usually exhibit poor aerodynamic performance and efficiency [START_REF] Ursache | Morphing Winglets for Aircraft Multi-Phase Improvement[END_REF]. Adaptive or morphing structures hold the potential to solve this problem and studies on wing deformation are subject of much interest in the aerospace domain. Recent advances made in the field of smart-materials have renewed this interest [START_REF] Weisshaar | Morphing aircraft systems: Historical perspectives and future challenges[END_REF][START_REF] Rivas Mcgowan | Aeroservoelastic and structural dynamics research on smart structures conducted at NASA langley research center[END_REF]. The Electro-active morphing for micro-airvehicles (EMMAV) research program, which was created as part of the French foundation of «Sciences et Technologies pour l'Aéronautique et l'Espace»'s effort to develop micro-and nano-air-vehicles and is composed of three French laboratories (IMFT, LAPLACE, ISAE), aims at optimizing the performance of micro-air-vehicles in realistic environments via electro-active morphing [START_REF] Rouchon | Activation d'une voilure déformable par des câbles d'AMF répartis en surface[END_REF].

As part of this collaborative effort a prototype NACA 0012 wing was developed with embedded Shape memory alloys (SMA) and trailing-edge piezoelectric actuators. This allowed both large deformations (~10% of the chord) at limited frequency (≤ 1Hz) and small deformations (several µm) at high frequencies (≤100Hz) [START_REF] Chinaud | Hybrid Electroactive Wings Morphing for Aeronautic Applications[END_REF]. The characteristics of the SMA technology, which were activated using the well understood Joule effect [START_REF] Manzo | Analysis and design of a hyperelliptical cambered span morphing aircraft wing[END_REF], make it especially suitable to optimize the shape of the wing and to control the flight [START_REF] Barbarino | Morphing trailing edges with shape memory alloy rods[END_REF][START_REF] Song | Robust control of a shape memory alloy wire actuated flap[END_REF]. The high-frequent but low amplitude piezoelectric technology on the other hand is useful to produce trailing-edge vortex breakdown [START_REF] Prechtl | Design of a high efficiency, large stroke, electromechanical actuator[END_REF][START_REF] Hall | Design and testing of a double X-frame piezoelectric actuator[END_REF][START_REF] Straub | Development and whirl tower test of the SMART active flap rotor[END_REF].

Whereas the hybridization of different smart-materials is a first to the authors' knowledge various previous studies have already addressed the use of SMAs [START_REF] Barbarino | Morphing trailing edges with shape memory alloy rods[END_REF][START_REF] Manzo | Analysis and design of a hyperelliptical cambered span morphing aircraft wing[END_REF][START_REF] Song | Robust control of a shape memory alloy wire actuated flap[END_REF] and piezoelectric actuators [START_REF] Bent | Active fiber composites for structural actuation[END_REF][START_REF] Bilgen | Macro-fiber composite actuated simply supported thin airfoils[END_REF][START_REF] Magrab | Vibrations of Elastic Systems: With Applications to MEMS and NEMS[END_REF][START_REF] Park | Analytical development of single crystal macro fiber composite actuators for active twist rotor blades[END_REF] separately. Like in the previously mentioned hybrid prototype the materials used in these studies were fairly idealistic. To address the issues emerging when dealing with less compliant materials a second prototype was conceived addressing more industrial concerns. This paper will describe the design and development of this second prototype, which was designed based on a NACA 4412 airfoil, and highlight the integration issues encountered. This remainder of this abstract will briefly introduce the designed hybrid prototype, illustrate the design choices and used materials and provide an outlook on the full paper.

Prototype

As mentioned in the introduction the proposed hybrid NACA 4412 prototype was conceived in order to address industrial concerns such as the material choice, wing geometry, material distribution, etc. An illustration of the designed wing can be seen in Figure 1. As can be seen in the illustration the deformable area was limited to ~50% of the chord. Further constraints included the use of aluminium or carbon fiber instead of more compliant materials as well as the absence of a compliant skin. In order to conduct flow visualizations the chord size of the prototype was limited to 45 cm. Due to the limited available space the high frequency vortex breakdown is to be done using Macro Fiber Composite (MFC) actuators. These piezoelectric fiber actuators provide a high frequency of actuation as well as a suitable displacement range.

Conclusions and Outlook

The work to be presented is developed twofold: first the actuation performance of the hybrid prototype from both a mechanical and electrical point of view is highlighted, including the maximum achievable deflection as well as the energy consumption. This data is then compared to the achievable gains both macroscopically in terms of lift increase and/or drag decrease as well as mesoscopically in terms of vortex breakdown.

Fig. 1 :

 1 Fig. 1: NACA 4412 hybrid prototype