Alexie
Boyer
,
Marius
Hervé
,
Audrey
Scognamiglio
,
Vincent
Loriot
and
Franck
Lépine
*
Univ Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622, Villeurbanne, France. E-mail: franck.lepine@univ-lyon1.fr
First published on 29th November 2021
Unraveling ultrafast processes induced by energetic radiation is compulsory to understand the evolution of molecules under extreme excitation conditions. To describe these photo-induced processes, one needs to perform time-resolved experiments to follow in real time the dynamics induced by the absorption of light. Recent experiments have demonstrated that ultrafast dynamics on few tens of femtoseconds are expected in such situations and a very challenging task is to identify the role played by electronic and nuclear degrees of freedom, charge, energy flows and structural rearrangements. Here, we performed time-resolved XUV-IR experiments on diamondoids carbon cages, in order to decipher the processes following XUV ionization. We show that the dynamics is driven by two timescales, the first one is associated to electronic relaxation and the second one is identified as the redistribution of vibrational energy along the accessible modes, prior to the cage opening that is involved in all fragmentation mechanisms in this family of molecules.
With the discovery of small nano-diamondoids in meteorites,9 diamondoids also gained a lot of attention from the astrophysical community. Due to isotope anomaly, these diamondoids are considered to have presolar origin9 and are even considered as the most abundant component of presolar grains.10 Regarding their high stability, they are also expected to be abundant in the interstellar medium. However, compared to polycyclic aromatic hydrocarbon (PAH) emission bands that dominate most of the mid-IR emission from astronomical objects, only few astronomical objects show spectroscopic evidences of the presence of diamondoids. For example, they are thought to be responsible for two broad emission bands (at 3.43 and 3.53 μm), observed in the spectrum of circumstellar disks around the stars Elias 1 and HD 97048 but their size is still a matter of debate.11–13 Regarding the strength of these bands, diamondoids are thought to represent 1–2% of the amount of cosmic carbon.14,15 The difficulty to detect emission spectra associated to diamondoids have been partially explained by the specific local conditions for their formation.16 However, the presolar origin of the diamondoids found in meteorites indicates sources of diamondoids outside the solar system. To understand the paucity of diamondoids in the interstellar medium, a large number of laboratory experiments and theoretical calculations have been developed, in particular to investigate the photostability of these molecules. Indeed, in space and for instance in Photo-Dissociation Regions, molecules are exposed to extreme ultraviolet (XUV) radiations above 10 eV. In that case, a low photostability could also explain the low detection of spectroscopic features. It was shown that upon ionization small diamondoids have low photostability, with efficient H-loss channels to form a stable closed-shell cation17 or fragmentation following the opening of the cage structure.18 The photostability of the doubly charged adamantane was investigated using XUV excitation, showing that the dissociation is driven by Coulomb repulsion.19 Although, these processes are expected to occur on ultrafast timescales, ranging from attosecond to picosecond, we note that no experiment has been able to investigate in a time-resolved fashion, the photoinduced evolution of diamondoids.
In this work, we investigate the electronic and vibrational relaxation dynamics of the two smallest diamondoids, adamantane (C10H16) and diamantane (C14H20), following XUV ionization. We use an ultrashort XUV pulse to ionize and excite the molecules and the induced dynamics is probe using an ultrashort IR pulse. The observed dynamics are associated to the electronic relaxation of the correlation bands (CBs), features created by electron correlation, and the following vibrational redistribution of the energy to lower electronic states. Similar temporal dynamics have been observed previously in polycyclic aromatic hydrocarbons (PAHs).20,21 Here, we investigate the effect of the 3D structure of these nanoscale diamonds on the energy redistribution and discuss its implication in understanding cage opening.
Mass spectra recorded from the interaction of neutral adamantane with either XUV or IR pulses are shown in Fig. 1a. These measurements were obtained with XUV photons energy centered around 17–23 eV or 30–33 eV and an IR intensity of about 5.5 TWcm−2. The spectra show a fragmentation pattern corresponding to the loss of CnHx groups. Similar patterns have been observed in collision induced dissociation (CID)24 or electron impact experiments.25,26 In CID experiments, the kinetic energy of the ions is converted into internal energy due to the collision, resulting in statistical fragmentation. Here, the main statistical fragments observed in CID correspond to the dominant peaks at groups of masses C6Hx+ (m/q ≈ 79) and C7Hx+ (m/q ≈ 93). The same behavior is observed in infra-red multi photon ionization (MPI) as shown in the mass spectrum measured with IR only (black line), where the number of photon is not sufficient to double ionize the molecule but the internal energy left in the monocation is sufficiently high to produce the same two groups of fragments. The same dominant peaks are also observed following the absorption of XUV photons of 17–23 eV (blue line). Since the absorption of 17–23 eV is not sufficient to double ionize the molecule, we conclude that most of the fragment ions observed in the corresponding mass spectrum are coming from the fragmentation of the excited monocation. In particular, some fragments such as the groups of masses C6Hx+ (m/q ≈ 79) and C7Hx+ (m/q ≈ 93) can be identify as statistical fragments. For the absorption of XUV photons of 30–33 eV (red line), the energy is sufficient to lead to the second ionization of the molecule. We note that no signature of the stable dication is measured at mass m/q = 68. This illustrates the low stability of the doubly charged adamantane. Indeed, one can see in Fig. 1a that some of the fragments detected in the mass spectrum at energy above IP2 (red line) are significantly different from the mass spectrum obtained at energy lower than the (blue line). For instance, the dominant peak for higher energy photons corresponds to the mass m/q = 41 (C3H5+). No significant amount of this fragment is observed in the statistical fragmentation of the monocation while it has been identified as part of the 2 and 3-body break-up channels of the dication.19 Some low masses can therefore be identified as coming mainly from the dissociation of the dication. Overall, the measured fragment ions in our experiment can come from two different mechanisms: statistical fragmentation of the monocation or dissociation of the unstable dication.
In the case of diamantane, only few experiments have been performed on the ionized molecule.27 To our knowledge, no study on the photostability of the doubly charged diamantane has been reported in literature. Its double ionization potential is also not known. However, regarding the similarity of adamantane and diamantane,2 it can be assumed to similar to IP2 of adamantane. The mass spectra recorded in diamantane are shown in Fig. 1b. The mass spectrum recorded from the interaction of the neutral diamantane with XUV photons energy centered around 20–23 eV (green line) exhibits the same fragmentation pattern than the one obtained in electron impact experiment.27 The dominant peaks are similar and correspond to the mass C7Hx+ (m/q ≈ 93). It means that the energy used in this experiment is probably not sufficient to significantly double ionize the molecule. Thus, the observed fragments mainly come from the statistical fragmentation of the monocation. When the XUV photon energy increases up to 33–39 eV, a significant number of small masses appears (purple line). Among these fragments, the masses m/q = 29 (C2H5+) and 41 (C3H5+) become dominant and correspond to fragments involved in the 2-body and 3-body break-up channels of the doubly charged adamantane. One can assume that, similarly to adamantane, these fragments come from the dissociation of the dication. Thus, the different fragments observed in the mass spectra arise either from statistical fragmentation of the monocation or from the dissociation of the dication. We also note that in the case of diamantane, the internal energy left in the monocation after the multi photon ionization by the IR is not sufficient to reach the fragmentation threshold of the monocation. Indeed, only the parent ion is observed in the IR only mass spectrum (black line).
Time-dependent signals are shown in Fig. 2. for adamantane (a–d) and diamantane (e–g). The signals were obtained using XUV pulses of 30–33 eV for adamantane and 20–23 eV for diamantane. The contribution from the XUV and IR pulses only are subtracted from the two color signal. Different temporal behaviors can be observed depending on the measured fragments. In the case of a fragment coming from the dissociation of the doubly charged molecules, C2H5+ (m/q = 29), the time-dependent signal exhibits a fast increase of the signal around the zero delay followed by an exponential decay of the yield (see Fig. 2b and f). In the case of major statistical fragments of the monocations, C7Hx+ (m/q ≈ 93), one can observe a population behavior with an increase of the signal after the zero delay following an exponential increase of the yield (see Fig. 2d and g). Thus, probing the dynamics by statistical fragmentation of the monocation gives access to a population timescale while probing the dynamics by the dissociation of the dication gives access to a decay dynamics. Depending on the observable, different ultrafast processes can be investigated.
In the case of diamantane, most of the studied fragments exhibit either a decay dynamics or a population dynamics as shown in Fig. 2f and g. In the case of adamantane, the number of fragments that exhibit both dynamics is high. The example of C5H7+ is given in Fig. 2c. However, despite the apparent complexity of these dynamics, our analysis led to the conclusion that only two time scales, τdecay and τpop, for the decay and population, respectively, are sufficient to describe the dynamics observed in all the fragments for adamantane and diamantane. This conclusion arises from a multistep procedure in our analysis starting from the use of free parameters (each transient is considered independently) to constrained parameters. Here we describe only the last step of our analysis. A multi-dimensional fitting procedure using the convolution between the cross-correlation of the XUV and IR pulses Icc and an exponential decay plus an exponential increase has been used:
(1) |
τ decay (fs) | τ pop (fs) | |
---|---|---|
Adamantane | 22 ± 5 | 28 ± 4 |
Diamantane | 29 ± 7 | 19 ± 6 |
The experiment was performed for different XUV photon energies and different IR intensities and no significant XUV and IR dependency of the timescales was observed. However, the amplitudes of the decay component Adecay and the population component Apop can change depending on the observed fragments, even going from positive to negative value in the case of adamantane. Moreover, the ratio of these two amplitudes also change with the XUV photon energy (without altering the measured timescales). This can be understood as a change in internal energy deposited into the molecule and the competition between statistical fragmentation of the monocation and direct dissociation of the dication. For instance, in the case of C5H7+ (m/q = 67) (Fig. 2c), the time-dependent signal is composed by a negative Adecay (orange dashed line) and a positive Apop (yellow dashed line). Since the two dynamics are measured from different observables, the negative Adecay means that the absorption of the IR pulse depopulated the states that should have led to the statistical fragmentation of the monocation towards C5H7+ to lead to the creation of dication. The formed dication later dissociates into other fragments.
A striking example of the two origins of the fragments is observed in the mass m/q = 41 (C3H5+) of adamantane. The 2D time-dependent mass spectrum is shown in Fig. 3a. One can see that signals surrounding mass 41 appear at the zero delay. It corresponds to the C3H5+ ions with high kinetic energy. Indeed, the presence of translational kinetic energy affects the arrival time of the ion on the detector, arriving earlier or later than the ion with low initial kinetic energy. This leads to the detection of signal around the m/q values. The corresponding time-dependent signals obtained by integrating the map at masses m/q = 40.6 ± 0.1, m/q = 41 ± 0.1 and m/q = 41.3 ± 0.1 are shown respectively in Fig. 3b–d. One can see that the signals of the fragments with high kinetic energy exhibit decay dynamics. Their high kinetic energy allows to identify them as purely coming from the dissociation of the dication following Coulomb repulsion. The signal at the central mass (ion with low kinetic energy) corresponds to mass 41 coming from statistical fragmentation only. It exhibits a negative amplitude of the decay, corresponding to a depletion of the statistical fragmentation channel towards the creation of dication. This shows that a given fragment can be reached following different pathways that can be separated as a function of the kinetic energy of the ions.
We note that the vibrational dynamics described above has been observed in PAHs in recent experiments. In PAHs, the population timescales were substantially longer than in the case of diamandoids. For instance, pyrene whose number of valence electron (N = 74) is similar to diamantane, has a population timescale of τpop = 101 ± 27 fs, while diamantane (N = 76) has a population timescale of τpop = 19 ± 6 fs. In the case of PAHs, ab initio calculations show that following CBs relaxation, mainly symmetric normal modes are excited. For instance, ab initio calculations performed on the smallest PAH, naphthalene, show that the frequency of the symmetric normal modes responsible for the electron-vibration energy redistribution are 514.30, 771.15 and 1050.83 cm−1 [ref. 20]. Since the normal modes are very sensitive to the geometrical structure of the system, one can expect that higher frequency modes (compared to PAHs) will be involved in the electronic relaxation of carbon cages such as diamondoids. For instance, the three lowest symmetric modes of adamantane have frequencies of 759.0, 1057.0 and 1496.0 cm−1 [ref. 29]. Consequently, we can assume that the initially populated vibrational modes have higher frequency in the case of diamondoids compared to PAHs. These higher frequency modes are considered as the starting point of the subsequent vibrational dynamics and therefore one can expect the redistribution of energy to occur faster. This could explain the shorter timescales measured in the case of diamondoids compared to PAHs.
We conclude that in the case of diamondoid cations, the vibrational energy is spread in few tens of fs over the vibrational degrees of freedom. With high energy in the nuclear degrees of freedom, diamondoids are unstable and fragment on longer timescales. Although there is not detailed theoretical or experimental investigations about these fragmentation mechanisms, any fragmentation of diamondoids starts with an essential first step of the molecular reorganization that is the opening of the carbon cage.18,19 Cage and ring openings are a general feature in organic molecules that recently attracted the attention of the ultrafast community.30 The energy barrier corresponding to the cage opening in adamantane was found to be 1.71 eV above the ionization potential of the molecule18 and therefore could be easily overcome in our experiment. Since the population dynamics is associated with the energy flow towards nuclear degrees of freedom and prior to dissociation, one can try to interrogate the link between τpop and the timescale of cage opening.
To answer this question, one should consider that the complete cage opening is expected to occur on femtosecond timescale and our measurements are rather connected to the first steps of this structural change. The comparison between the timescales measured in adamantane and diamantane provides a direct information on the size dependence of the process. One expects that the cage opening will take more time for a larger diamondoid (diamantane). However, the population timescales in diamantane is about 10 fs faster than in adamantane. This decrease of τpop with the size of the molecules is similar to what was observed in our previous work on PAHs.21 It is associated with the larger number of accessible degrees of freedom for a larger system, which leads to faster dynamics. We conclude that in our XUV-IR experiment the extracted timescale reveals the first instants of energy flow towards all accessible vibrational modes, which gives birth to the cage opening and subsequent fragmentation. Thus, our experiment is neither sensible to the cage opening of the molecule, neither to the fragmentation process itself which should occur on longer (picosecond) timescale. The vibrational dynamics we are measuring is an initial step of energetic relaxation, prior to any structural change and occurs in the intact parent ion. The larger the molecule, the faster this initial step even if the complete structural change itself is slower.
In contrast to the case of PAHs, the CBs dynamics was not observed by measuring the stable dication signal, but using fragments arising from the dissociation of unstable dications. We showed that the fast vibrational dynamics measured in diamondoids could be understood by the high frequency vibrational modes involved in the relaxation of the CBs in the case of 3D molecules. We show that in our experiment, the timescales are associated to the initial step of cage opening that is a general feature of diamondoid fragmentation. Our experimental results show that XUV induced ultrafast electronic and vibrational dynamics in diamondoids occurs on very short timescales of few tens of fs, reaching a threshold from where the molecule dissociates through complex structural rearrangements. Further time-resolved experiments, for instance using coulomb explosion imaging, could allow to measure XUV induced structural changes on longer timescales overall providing a multiscale “movie” of photoinduced processes in diamondoids.
This journal is © the Owner Societies 2021 |