Dario
Campisi
*a,
Frederik Doktor S.
Simonsen
b,
John D.
Thrower
b,
Rijutha
Jaganathan
b,
Liv
Hornekær
bc,
Rocco
Martinazzo
d and
Alexander G. G. M.
Tielens
*a
aLeiden Observatory, Leiden University, Niels Bohrweg 2, 2333 CA Leiden, The Netherlands. E-mail: campisi@strw.leidenuniv.nl
bDepartment of Physics and Astronomy, Aarhus University, Ny Munkegade 120, DK-8000 Aarhus C, Denmark
cInterdisciplinary Nano-Science Centre (iNano), Aarhus University, Denmark
dDepartment of Chemistry, Università degli Studi di Milano, Via Golgi 19, 20133, Milan, Italy
First published on 16th December 2019
Investigating the hydrogenation of carbonaceous materials is of interest in a wide range of research areas including electronic device development, hydrogen storage, and, in particular, astrocatalytic formation of molecular hydrogen in the universe. Polycyclic Aromatic Hydrocarbons (PAHs) are ubiquitous in space, locking up close to 15% of the elementary carbon. We have used thermal desorption measurements to study the hydrogenation sequence of pentacene from adding one additional H to the fully hydrogenated pentacene species. The experiments reveal that hydrogenated species with an even number of excess H atoms are highly preferred over hydrogenated species with an odd number of H atoms. In addition, the experiments show that specific hydrogenation states of pentacene with 2, 4, 6, 10, 16 and 22 extra H atoms are preferred over other even numbers. We have investigated the structural stability and activation energy barriers for the superhydrogenation of pentacene using Density Functional Theory. The results reveal a preferential hydrogenation pattern set by the activation energy barriers of the hydrogenation steps. Based on these studies, we formulate simple concepts governing the hydrogenation that apply equally well for different PAHs.
In this work, we focus, specifically, on one specific kind of linear catacondensed PAH, pentacene (C22H14), representative of a class called acenes. These compounds, consisting of a linear row of condensed six membered rings,17 are characterized by zigzag edges and they are representative of the smallest graphene nanoribbon. Zigzag edges are a structural repetition of carbon atoms already bound to a single hydrogen atom, well known in the literature as “solo sites”,18 and bridge carbon atoms that are not bound to any hydrogen atoms. The zigzag edges are important in astrophysics since they have strong infrared signatures that have been clearly detected in astronomical spectra.19 Specifically, the strong C–H out-of-plane bending mode at 11.2 μm is characteristic of C–H groups that belong to zigzag edges20–22 and this band dominates the observed IR spectra of HII regions, YSOs and reflection nebulae (RNe) as well as the diffuse ISM.18 The presence of zigzag edges is also of interest to the chemistry community since they allow charge localization at the edges of the PAHs, and that increases their reactivity.
In contrast to catacondensed PAHs such as the acenes, some pericondensed PAHs (e.g. coronene) have no zigzag edges and their structure is characterized by duo carbons (two adjacent hydrogenated carbon sites) with strong CH modes in the 12–13 μm range.20,21 Both catacondensed and pericondensed PAHs are susceptible to hydrogenation by hydrogen atoms under interstellar conditions, but the hydrogenation sequence differs due to their different structures.23,24 For instance, a recent study demonstrated that the hydrogenation sequence of coronene, a pericondensed PAH, starts from the duo carbons24 with delayed hydrogenation of the other carbon sites. For catacondensed species, the hydrogenation sequence begins from the central aromatic ring, saturating both inner and edge carbons.23,25
Superhydrogenated PAHs are important since they might be responsible for the formation of molecular hydrogen in Photodissociation Regions (PDRs).26 PDRs are regions where the ultraviolet radiation from nearby bright stars influences the chemical and physical properties of the present molecular structures.27 It is well appreciated that, because of the high UV flux, molecular hydrogen has to be rapidly reformed and that the normal routes considered for H2 formation in the ISM do not suffice to account for its observed abundance.28,29 As astronomical observations have revealed a spatial correlation between the H2 formation rate and interstellar PAHs emission in these regions, PAHs have been suggested as possible catalysts for H2 formation.29,30 Different mechanisms have been studied in order to understand the formation of molecular hydrogen on PAHs, including reactions between two chemisorbed atomic hydrogen atoms through Langmuir–Hinshelwood (LH) and Eley–Rideal (ER) mechanisms. These mechanisms have been studied for pyrene31 and for coronene,3,5 demonstrating that Eley–Rideal abstraction is a viable pathway for PAH catalyzed molecular hydrogen formation. Moreover, Rasmussen25 studied the correlation of the increase in hydrogen binding energy with increasing zigzag edge length, using several acenes with different length as models. Therefore, it is crucial to study the aromatic–aliphatic character of PAHs when they are in a hydrogen rich environment in order to shed light on the sequence that leads to the so called “Magic Numbers” in the hydrogenation sequence of these species.32
This study aims to determine and understand the hydrogenation sequence of pentacene that leads to complete superhydrogenation of the molecule, using a combination of experimental and theoretical techniques. Specifically, temperature programmed desorption (TPD) is employed to directly determine the sequence of products of hydrogen addition to pentacene while density functional theory (DFT) calculations are used to investigate the hydrogenation sequence. This also serves as a test of simple concepts that may be applicable to other aromatic species featuring a bipartite lattice (e.g. carbon sp2 structured as a graphene lattice).
The monolayer of pentacene on HOPG was then exposed to a beam of atomic hydrogen (H) using a hydrogen atom beam source (HABS).33 The HABS was operated at ≈2200 K, producing a hot atomic H beam which was subsequently cooled collisionally, through a quartz nozzle, to an estimated temperature of ≈1400 K.
TPD measurements were performed with a linear ramp of β = 1 K s−1 until reaching 1100 K. To mass-selectively detect the species that leave the surface a quadrupole mass spectrometer (QMS) (Extrel CMS LLC) was used, scanning the range 258–340 atomic mass units (amu) per charge. The TPD spectra within the given mass range are integrated and averaged over the relevant desorption temperature range to give the presented mass spectra.
The theoretical part of this work is divided into two sections: the first describes the sequence of H addition that leads to fully hydrogenated pentacene while the second describes the hydrogenation sequence of superhydrogenated species based on a binding energy and energy barrier analysis.
The binding energies are defined as:
Ebind = [(EH + Ereagent) − EnH-pentacene] |
The energy barriers are defined as:
Ebarr = [ETS − (EH + Ereagent)] |
Finally, the formation energy is defined as:
Ef = [Eproduct − (EH + Ereagent)] |
In order to gain broader insights into the hydrogenation process, the sequence of hydrogenation was rationalized with simple concepts derived from a tight-binding model where the most reactive site is related to the smallest π-coordination number (number of carbon atoms connected to a carbon center) and the highest hypercoordination number (number of second carbon atom neighbours that have the same coordination number) for closed-shell systems.42 Attaching a hydrogen on a closed-shell system (i.e. all odd extra hydrogens chemisorbed) generates an unpaired electron (mid-gap state) that localizes preferentially in the ortho- or para-positions respect to the already extra hydrogen attached. Therefore, the odd hydrogenations were predicted using the coordination and hyper-coordination rule, whereas the even hydrogenations were predicted using the ortho–para localization. Finally, these simple rules were confirmed through the binding energy and energy barrier analysis. The latter are correlated by the Bell–Evans–Polanyi principle according to which larger energy barriers are associated with lower binding energies.43,44
Fig. 1(b–e) show the mass distribution of superhydrogenated pentacene after hydrogenation times of (b) 180 s, (c) 900 s, (d) 1800 s and (e) 54000 s (17 hours). Fig. 1(b) shows a small degree of hydrogenation for both the pentacene and the oxidized species. The species with even masses such as 280, 282 and 284 are seen to be more prominent than the odd. This is consistent with the fact that superhydrogenated pentacene molecules with an odd mass will be open-shell, and hence exhibit a radical behaviour and therefore be more reactive towards the addition of a subsequent H atom. By contrast, the even mass species are less prone to further hydrogenation. The odd–even oscillation is clearly visible in Fig. 1(c) where the mass distribution covers the entire range from pristine pentacene, mass 278, to fully superhydrogenated pentacene, C22H36, with mass 300. While the odd–even oscillation is clearly visible, it is also evident that some specific species (with specific superhydrogenation degrees) appear more prominently than others. The most prominent species are: 280 (+2H), 282 (+4H), 284 (+6H), 288 (+10H), 294 (+16H) and 300 (+22H), also indicated in Fig. 1(c).
After extended H atom exposure, see Fig. 1(d–e), the mass distribution is pushed towards mass 300, corresponding to fully superhydrogenated pentacene C22H36. In Fig. 1(e) the peak at mass 300 completely dominates the spectrum. The smaller peaks seen at masses 301–302 and 298–299 are ascribed to the 13C isotope and fragmentation in the QMS, respectively. Hence, the fact that it is possible to drive the entire mass distribution to fully superhydrogenated pentacene for long H deposition times (see Fig. 1(e)) reveals that H atom addition reactions dominate over H atom abstraction for superhydrogenated pentacene.
Fig. 2 Optimized molecular structure of pentacene in balls and sticks. The numerical labels reported indicate the full hydrogenation sequence of pentacene. |
The hydrogenation sequence follows a simple rule: first hydrogenation of the carbon in the solo position (CH bonded H where adjacent C atoms do not have a H), followed by hydrogenation of the neighboring inner carbons (bridge carbons). The latter step reflects that an inner carbon is under-coordinated upon hydrogenation of a neighbour carbon site. This is in line with the study carried out by Klærke et al.48 and Bonfanti et al.42 confirming that it is easier to hydrogenate edge carbons than bridge ones.
The highest binding energy found is associated with the chemisorption of the extra hydrogen in position 2 (Fig. 2), which is para with respect to carbon 1. Once carbons 1 and 2 are hydrogenated, two topologically disconnected systems are generated. The molecule now consists of two naphthalene units, bound to the hydrogenated non-aromatic ring, outlined in bold in Fig. 4(a). These two naphthalene units are equivalent and therefore, we considered the chemisorption of the third hydrogen atom on C3 (Fig. 4(a)) which has the highest hyper-coordination. Fig. 4(a) also reports graphically the propensity of the carbon atoms to be hydrogenated (details reported in the ESI†) where it is clearly shown that C3 and C15 have higher affinity (high hyper-coordination) compared to the inner carbon site C11 as the latter has the highest coordination number, while, C16 and C4 have lower propensity to be hydrogenated since these carbon sites have lower hyper-coordination. The fourth hydrogen atom, then, chemisorbs on the least coordinated carbon site (labeled C4 in Fig. 4(a)) that is a bridge carbon. The sequence, schematically shown in Fig. 2, continues with the system tending to preserve the aromaticity of individual rings as long as possible. Moreover, during the hydrogenation, disconnected systems are created, as naphthalene and benzene rings (shown in bold in Fig. 4(b)), and the hydrogenation sequence might change accordingly. A complete overview of the sequence is reported in the ESI.†
Fig. 4 (a) Graphical representation of disconnected naphthalene units and their propensity to be hydrogenated. (b) Graphical representation of disconnected benzene units. |
Fig. 5 Side view and top view of the optimized molecular structure of syn-22H-pentacene in balls and sticks. |
As it may be relevant for the hydrogenation of gaseous pentacene, we also considered binding the hydrogen atoms to the bridging carbon atoms of the molecule on alternating sides of the molecular plane (trans-hydrogenation), referred to here as “anti position”, and shown in Fig. 6. The anti conformation in the case of pentacene is relevant only in carbon sites where no hydrogen atoms are already bound. In fact, in the case of edge carbons where two hydrogen atoms are attached, the hydrogen atoms are topologically equivalent.
The hydrogenation sequence for the anti path begins with the 4th hydrogenation since the hydrogen atoms already chemisorbed are topologically equivalent. Specifically, the hydrogen atom may chemisorb on carbon 4 (Fig. 2) from the bottom side with respect to the plane of the molecule. The addition of the next hydrogen atom to carbon 5, proceeds from the top of the molecule with respect to the molecular plane. Hence, if the hydrogen atom chemisorbs on the bottom side of the molecule, specifically on carbon atoms 4, 8, 11 and 20, we can observe a more linear structure with all rings in arm-chair conformation. A complete description of why those hydrogenated species are formed, based on the analysis of binding energies and energy barriers, is reported in the next section and a more complete description of the sequence is reported in the ESI.†
Fig. 6 Side view and top view of optimized molecular structure of anti-22H-pentacene in balls and sticks. |
Fig. 7 Front view and side view, respectively, of the optimized molecular structure in balls and sticks of quasi-nanotube formed by the chemisorption of 12th hydrogen on a less reactive carbon site of syn-11H-pentacene (see section 12th hydrogenation in the ESI†). |
Table 1 provides the exoergic formation energies, indicating that the formation of quasi-nanotubular structures is thermodynamically favourable, especially during the addition of the 16th hydrogen atom in a less reactive carbon site (more details are reported in ESI†). Moreover, the formation of these nanotubular structures occurs only in the gas-phase and involves an intramolecular homo-coupling (formation of a carbon bond within a molecule). In a similar study, the formation of intermolecular homo-coupling (formation of a carbon bond between two molecules) between two pentacene molecules on Au(111), resulting from superhydrogenation, has been investigated.49 In this case, intermolecular homo-coupling dominates on a surface. This process is not important in the gas phase as collisions are rare. Instead, in the gas phase, intramolecular nanotubular formation dominates.
Hn | ΔEf (eV) |
---|---|
10 | −2.39 |
12 | −1.20 |
14 | −2.50 |
16 | −2.70 |
syn-21H-pentacene can react through a barrierless reaction with another radical to produce the fully hydrogenated pentacene with a mass of 300 amu. The peak corresponding to mass 300 amu becomes prominent as exposure times to atomic hydrogen increases. The presence of magic numbers in the hydrogenation of pentacene reflects that the energy barriers are particularly high whenever a new aromatic ring has to be attacked.
The values of binding energies reported in Fig. 9 suggest that all even hydrogenated species are more favourable than the odd hydrogenated species and have binding energies between 3 and 4.5 eV. The species with an odd number of extra hydrogen atoms attached display smaller binding energies between 2.5 eV and 1 eV. This trend was studied and confirmed on several neutral PAHs,3,25,31 while for the cation the opposite is found as expected for an open shell species.32,50 Within the odd and even sequences there are small variations in binding energy which are due to high affinity of these carbon sites toward the hydrogenation. However, these differences are small and, at the temperature of the experiment, are unlikely to be relevant for the discussion of the magic numbers.
Fig. 9 Trend of binding energies for the hydrogenation of pentacene for two different paths: syn and anti. H number indicates the number of extra hydrogens chemisorbed on pentacene. |
Our study of the binding energies reveals that superhydrogenation is energetically favorable for all species. The odd–even pattern and the presence of magic numbers in the even sequence reflects the pattern of energy barriers in the production of odd species from even ones. The lowest energy barriers found in this work are related to the chemisorption of the hydrogen on aliphatic carbon and, in particular, for the case of the formation of syn-9H, syn-13H, syn-15H and syn-19H-pentacene.
Table 2 reports the difference in energy barrier between anti and syn. The trend revels significant differences between syn and anti that cannot explain the magic numbers found in our experiment. Specifically, based on this trend we should not expect the presence of anti-6H-pentacene due to the lower barrier with respect to the syn counterpart, whereas we might expect the anti-20H-pentacene due to the difficulties associated with further hydrogenation. For a complete overview, see Fig. 8.
Hn | ΔEantibar − ΔEsynbar (eV) |
---|---|
7 | −0.070 |
11 | 0.039 |
13 | 0.010 |
15 | 0.019 |
21 | −0.072 |
The anti path has the same trend in binding energies as the syn path. We show that anti-6H, anti-10H, anti-12H, anti-14H, anti-16H and anti-20H-pentacene have large binding energies (larger than 3.4 eV), while anti-11H and anti-17H have lower binding energies by 1 eV. Table 3 shows the difference of binding energies between syn and anti hydrogenation, indicating that the anti chemisorption of hydrogen atoms might be more favourable for some hydrogenations.
Hn | ΔEantibind − ΔEsynbind (eV) |
---|---|
5 | 0.12 |
10 | 0.24 |
14 | 0.25 |
16 | 0.17 |
20 | 0.16 |
22 | 0.30 |
Some hydrogenation steps are not included in Table 3 as the energy differences are small (∼0.007–0.09 eV). The observed pattern of magic numbers in the hydrogenation of pentacene adsorbed on graphite is fully consistent with our expectation that, because of steric reasons, hydrogenation will occur following the syn sequence. Specifically, we observed a dominant peak at the mass of 6H-pentacene and only a minor peak corresponding to 20H-pentacene in accordance with the syn-sequence, but not the anti-sequence.
Analysis of this data has revealed that H2 must be very efficiently made in these surface layers.29 This is not understood as the dust in PDRs is too warm to allow efficient H2 formation.28 This observed efficient H2 formation in PDR surfaces has been ascribed to the catalytic activity of PAHs and detailed studies have appeared in the astronomical literature.8,52–54 The present study is geared towards elucidating the interaction of atomic H with PAHs and quantitative measuring of the molecular parameters involved that are needed for these astronomical models. Our results show that addition of the first two hydrogens will be very efficient as the energy barriers involved are small and the kinetic temperature of the gas is quite high. However, adding the third hydrogen has a prohibitive energy barrier (Fig. 8) even at gas temperatures relevant for PDRs (200–500 K). To estimate the importance of quantum mechanical tunneling, we have evaluated the cross-over temperatures for the reactions involved and these are reported in ESI.† Typically these vary from 136 to 213 K and hence, quantum mechanical tunneling is not expected to play a role in the warm gas of PDRs but it may play more of a role in cold molecular clouds. We note also that the binding energies of these extra H's are relatively low and these are readily lost in strong radiation field of a PDR.8,53 Detailed astronomical modeling will be required to assess the balance of reactions of atomic H with PAHs and the UV photolysis in a PDR environment as well as the importance of hydrogenation of PAHs in PDRs for the formation of H2. The experimental and theoretical results reported here will form the basis for such a modeling effort. We defer this to a follow up study that will be reported in an appropriate journal.
Specifically, we found that even hydrogenated species of pentacene are more preferred than the odd hydrogenated species since chemisorbing odd hydrogens on pentacene leads to the formation of reactive open-shell radicals. These extremely reactive species interact with further hydrogen atoms forming preferred species (species with magic numbers seen to dominate the mass spectra). The analysis of binding energies and energy barriers carried out in this study have confirmed the validity of the hypercoordination and coordination rule (discussed in detail in the ESI†), demonstrating that these simple concepts are applicable to small condensed and linear systems. These concepts are validated for systems formed by hexagonal aromatic rings and are applicable during the hydrogenation sequence for all benzene-like structures that still present bipartite lattices. However, we emphasize that these simple concepts are useful tools and can guide the choice of hydrogenation sites, but they cannot replace a deeper analysis of binding energies and energy barriers. Specifically, energy barriers calculation provide insights into the hydrogenation sequence, showing which species are more likely to form. We found that the energy barrier required for hydrogenating an aliphatic carbon is lower (<0.13 eV) than the energy required for hydrogenating aromatic carbons (up to 0.15 eV). Hence, the presence of magic numbers in the mass spectra are explained by the energy barriers rather than the binding energies. In fact, the absence in the mass spectra of the peaks corresponding to the masses of 12H, 14H, 18H and 20H-pentacene was confirmed by our theoretical calculations since a further hydrogenation (11H, 13H, 17H and 19H-pentacene) requires the highest energy barrier, up to 0.16 eV. Our theoretical results also suggest that both anti and syn confirmations are possible in the gas phase. The sequence found for the syn-path in our calculation matches the magic numbers found in the mass spectra for the hydrogenation of pentacene, whereas the anti-path might be more relevant for the hydrogenation sequence that will occur in the gas phase.
Furthermore, we noticed differences in the hydrogenation sequence that occurs in pentacene with respect to pericondensed PAHs such as coronene.24 The hydrogenation scheme tends at first, in these aromatic molecules, to break the aromaticity of the central aromatic ring of pentacene through addition to solo sites. In contrast in coronene, where such sites are not present, this occurs for the aromatic rings located at the edges. Moreover, pentacene is more reactive at the first hydrogenation with a barrier of only 109 meV with respect to coronene.24 As a last difference, coronene has a barrierless reaction for the addition of the last odd hydrogen. Therefore, zigzag edges are significantly more reactive for the first hydrogen addition reaction with respect to pericondensed PAHs.
We also discovered a pathway that leads to the formation of quasi-nanotubular structures during the hydrogenation of pentacene in the gas phase. This is due to the formation of two topologically separated structures in the molecule. The two unpaired electrons located in two different systems cannot be shared any longer. Therefore, a spontaneous curling up of the structure occurs with the formation of a quasi-nanotube.
We discovered that the flexibility of the molecule provides extra structural stability with respect to more rigid molecules as coronene or graphene. Zigzag edges present on linear and flexible molecules may provide extra reactivity that may correlate the presence of H2 with the formation of superhydrogenated species.
Footnote |
† Electronic supplementary information (ESI) available: Sequence analysis of the hydrogenation of pentacene. Crossover temperatures and Cartesian coordinates of the optimized geometries of minima and transition states of the hydrogenation of syn and anti-nH-pentacene. See DOI: 10.1039/c9cp05440e |
This journal is © the Owner Societies 2020 |