# Dissociation dynamics of the diamondoid adamantane upon photoionization by XUV femtosecond pulses

Sylvain Maclot, Jan Lahl, Jasper Peschel, Hampus Wikmark, Piotr Rudawski, Fabian Brunner, Hélène Coudert-Alteirac, Suvasthika Indrajith, Bernd Huber, Sergio Díaz-Tendero, et al.

► To cite this version:

Sylvain Maclot, Jan Lahl, Jasper Peschel, Hampus Wikmark, Piotr Rudawski, et al. Dissociation dynamics of the diamondoid adamantane upon photoionization by XUV femtosecond pulses. Scientific Reports, 2020, 10 (1), pp.2884. 10.1038/s41598-020-59649-1. hal-03668782

HAL Id: hal-03668782

https://hal.science/hal-03668782

Submitted on 5 Dec 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
OPEN

# Dissociation dynamics of the diamondoid adamantane upon photoionization by XUV femtosecond pulses

Sylvain Maclot<sup>1,2*</sup>, Jan Lahl<sup>1</sup>, Jasper Peschel<sup>1</sup>, Hampus Wikmark<sup>1</sup>, Piotr Rudawski<sup>1</sup>, Fabian Brunner<sup>1</sup>, Hélène Coudert-Alteirac<sup>1</sup>, Suvasthika Indrajith<sup>3</sup>, Bernd A. Huber<sup>3</sup>, Sergio Díaz-Tendero<sup>4,5,6</sup>, Néstor F. Aguirre<sup>7</sup>, Patrick Rousseau<sup>3</sup> & Per Johnsson<sup>1*</sup>

This work presents a photodissociation study of the diamondoid adamantane using extreme ultraviolet femtosecond pulses. The fragmentation dynamics of the dication is unraveled by the use of advanced ion and electron spectroscopy giving access to the dissociation channels as well as their energetics. To get insight into the fragmentation dynamics, we use a theoretical approach combining potential energy surface determination, statistical fragmentation methods and molecular dynamics simulations. We demonstrate that the dissociation dynamics of adamantane dications takes place in a two-step process: barrierless cage opening followed by Coulomb repulsion-driven fragmentation.

Diamondoids are a class of carbon nanomaterials based on carbon cages with well-defined structures formed by C(sp³)-C(sp³)-hybridized bonds and fully terminated by hydrogen atoms. All diamondoids are variants of the adamantane molecule, the most stable among all of the isomers with the formula C₁₀H₁₆, shown in Fig. 1. On Earth, diamondoids are naturally found in petroleum deposits and natural gas reservoirs, and their most common applications are for the characterization of petroleum and gas fields, offering possibilities to e.g. trace the source of oil spills<fnref n="1" />. Today, diamondoids are attracting increasing interest for use as an applied nanomaterial in e.g. nano- and optoelectronics as well as in biotechnology and medicine due to their high thermal stability and well-defined structure in combination with no known toxicity². In space, diamondoids have been found to be the most abundant component of presolar grains³, and due to their high stability they are thus also expected to be abundant in the interstellar medium⁴.

However, when compared to laboratory measurements based on infrared spectroscopy⁵, astronomical observations show a deficiency of diamondoids in the interstellar medium which to date is not completely understood⁶. The first ionization limit in diamondoids lies around 8–9 eV with a maximum in the ionization yield between 10 and 11 eV⁷, close to the hydrogen Lyman-α line, and the efficient production of cations followed by dissociation has been suggested as a possible explanation for the apparent lack of diamondoids in the interstellar medium. Steglich et al. investigated the stability of diamondoid cations using ultraviolet irradiation, finding that rapid loss of a neutral hydrogen followed ionization⁸. Since then, further studies have suggested that small hydrocarbons are also created as dissociation products. In a recent work at the Swiss Light Source, vacuum ultraviolet radiation (9–12 eV) was used in combination with threshold photoelectron and photoion coincidence detection to determine the appearance energies and branching ratios of the resulting photofragments of the singly charged adamantane cation⁹. The study reveals, in addition to the expected hydrogen loss, dissociation via a number of parallel channels which all start with an opening of the carbon cage and hydrogen migration indicating that the

<footnote n="1">¹Department of Physics, Lund University, P.O. Box 118, 22100, Lund, Sweden. ²Biomedical and X-Ray Physics, Department of Applied Physics, AlbaNova University Center, KTH Royal Institute of Technology, SE-10691, Stockholm, Sweden. ³Normandie Université, ENSICAEN, UNICAEN, CEA, CNRS, CIMAP, 14000, Caen, France. ⁴Departamento de Química, Módulo 13, Universidad Autónoma de Madrid, 28049, Madrid, Spain. ⁵Condensed Matter Physics Center (IFIMAC), Universidad Autónoma de Madrid, 28049, Madrid, Spain. ⁶Institute for Advanced Research in Chemical Sciences (IAdChem), Universidad Autónoma de Madrid, 28049, Madrid, Spain. ⁷Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico, 87545, USA. *email: smaclot@gmail.com; per.johnsson@fysik.lth.se</footnote>
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Structure of adamantane.</figcaption></figure>

<figure><img src="image_3.png" /><figcaption>**Figure 2.** XUV spectrum and the first three ionization thresholds of adamantane.</figcaption></figure>

low photostability of adamantane could explain its deficiency in astronomical observations. While this study was recently complemented by a first time-resolved study¹⁰, to date no results have been published on the dissociation dynamics of multiply charged adamantane molecules.

In this work, we study the fragmentation dynamics of the adamantane dication after ionization by extreme ultraviolet (XUV) femtosecond pulses, the use of which ensures prompt and well-defined ionization. The experimental technique used in this study is based on correlated ion and electron spectroscopy, enabling the characterization of the charged products of interaction (identification and energetics). The support of various theoretical methods such as molecular dynamics simulations, potential energy surface determination and statistical fragmentation models, helps to unravel the fragmentation dynamics of such a complex molecular system.

## Methods

**Experiments.** The high-intensity XUV beamline at the Lund Attosecond Science Centre provides trains of attosecond pulses in the XUV spectral region using the high-order harmonic generation (HHG) technique¹¹,¹². This is achieved by focusing an intense infrared (IR) pulse (high-power Ti:Sapphire chirped pulse amplification laser with a pulse energy of 50 mJ, a central wavelength of 810 nm, a pulse duration of 45 fs and a repetition rate of 10 Hz) into an argon gas medium (6 cm long cell) in a loose focusing geometry (~8 m focal length)¹³. The train of XUV attosecond pulses, with a total duration of 20 fs, contains around 15 attosecond pulses with an estimated individual pulse duration of approximately 300 as, spaced by 1.35 fs. The photon energy spectrum of the produced XUV light is a characteristic harmonic comb spanning from ~20 to 45 eV (Fig. 2). Then, the XUV light is micro-focused (~5 × 5 µm²) on target via a double toroidal mirror¹⁴ and with a pulse train energy on target of around 10 nJ this leads to intensities of the order of 10¹² W·cm⁻².

Adamantane molecules, C₁₀H₁₆ (powder from Aldrich with >99% purity), are produced in the gas phase by a pulsed Even-Lavie valve¹⁵,¹⁶, heated to 100 °C and using He as a carrier gas, in the form of a cold and collimated supersonic jet.

The photon-molecule interaction leads to the formation of highly excited singly and doubly charged adamantane molecules (single and double ionization thresholds IT₁~9.2 eV and IT₂~23.9 eV - see Fig. 2). The trication is only produced in negligible amounts since the triple ionization threshold is IT₃~43.6 eV (Fig. 2) and thus is not discussed in the following. The charged products of interaction are detected by a double-sided velocity map imaging (VMI) spectrometer¹⁷ giving access to the kinetic energy distributions of ions and electrons on a shot-to-shot basis. In addition, the ion side of the spectrometer can measure the time-of-flight (TOF) of the ions, providing the mass spectrum. Despite the high count rates (several tens of counts per pulse), the use of the partial covariance technique¹⁸ enables us to disentangle the contributions of the different fragmentation channels. For instance, applying this technique on single-shot ion TOF spectra gives the possibility to produce ion-ion correlation maps that reveal the dissociation dynamics of the doubly charged adamantane molecules. In addition, the use of this
technique on single-shot ion TOF spectra and the single-shot ion VMI data gives access to the kinetic energy distribution of specific ionic fragments.

**Theory.** Three different theoretical methods have been used: (i) molecular dynamics (MD) simulations in the framework of the density functional theory (DFT) and the density functional tight binding (DFTB) method, (ii) exploration of the potential energy surface (PES) employing the DFT, and (iii) statistical fragmentation using the Microcanonical Metropolis Monte-Carlo (M3C) method.

*Molecular dynamics simulations.* The molecular dynamics approach, using the DFTB method<sup>19,20</sup>, has been used to compute the lifetime of doubly-ionized and excited adamantane. To this end, we have considered double ionization of adamantane in a Franck-Condon way; that is, our molecular dynamics simulations start from the optimized geometry of the neutral adamantane molecule after removal of two electrons. We assume that the electronic excitation energy is rapidly redistributed into the nuclear degrees of freedom and thus we run these simulations in the electronic ground state of the dication Ada<sup>2+</sup>, taking the equilibrium geometry of the neutral ground state as initial point. We have taken four values of excitation energy corresponding to the relative energy of the highest order of harmonics in the XUV spectrum (Fig. 2) with respect to IT<sub>2</sub> (8.46, 11.50, 14.63 and 17.80 eV, which correspond to temperature values 2520, 3450, 4360 and 5300 K, respectively). The used E<sub>exc</sub> values can be considered as the upper limits of the remaining electronic excitation energy after ionization with the four highest energy harmonics, assuming that the kinetic energy of the ejected electrons is zero, thus E<sub>ph</sub> = E<sub>exc</sub>. This excitation energy is randomly distributed into the nuclear degrees of freedom in each trajectory (i.e., in the velocity components v<sub>x</sub><sup>i</sup>, v<sub>y</sub><sup>i</sup>, v<sub>z</sub><sup>i</sup> of each atom i). Molecular dynamics trajectories are propagated up to a maximum time of t<sub>max</sub> = 1, 5, 10, 20 and 100 ps. For each value of excitation energy and propagation time a set of 1000 independent trajectories are considered (that is, the initial conditions are separately established in each set of trajectories). Statistics are then carried out over these trajectories to obtain information on the survival time of the doubly-ionized adamantane with different excitation energies. To ensure adiabaticity in the simulations a time step of Δt = 0.1 fs is used. These simulations have been carried out with the deMonNano code<sup>21</sup> and the results show that even with these relatively high excitation energies the dication does not fully fragment until after tens to hundreds of picoseconds making it computationally too heavy to perform full MD simulations at the DFT level (see SI).

*Potential energy surface calculations.* We have explored the PES using DFT, in particular we have employed the B3LYP functional<sup>22-24</sup> in combination with the 6-31G(d) basis set. This part of the study provides useful energetic and structural information of the experimentally measured exit channels. In order to identify the most relevant stationary points of the PES we have adopted the following strategy:

(1) We have first performed molecular dynamics simulations at the same DFT-B3LYP/6-31G(d) level, to mimic the evolution of the system during the first femtoseconds after the ionization and excitation. This part of the simulations also starts by computing the energy required to doubly ionize adamantane in a Franck-Condon transition from the optimized neutral structure. This is our starting point for the dynamics. 160 trajectories were carried out using the ADMP method with a maximum propagation time of 500 fs, and considering a time step of Δt = 0.1 fs, and a fictitious mass of μ = 0.1 au. Thus, after propagation of the doubly charged excited adamantane, we have obtained the evolution of the system in the first femtoseconds.

(2) Then, using the last step in the dynamics as an initial guess, we have optimized the geometry of the produced species.

(3) Finally, geometry optimization of fragments observed in the experiments has also been computed. To this end, we have considered several structures for each C<sub>n</sub>H<sub>x</sub>H<sub>x</sub><sup>q+</sup> fragment, thus obtaining the relative energy of the exit channels observed in the experiment. Harmonic frequencies have been computed after geometry optimization to confirm that the obtained structures are actual minima in the PES (no imaginary frequencies) and to evaluate the Zero-Point-Energy correction. These calculations were carried out with the Gaussian09 package<sup>25</sup>. The proposed strategy was used in the past with success to study the fragmentation dynamics of ionized biomolecules (see e.g.<sup>26-28</sup>).

Notice that we are not performing a complete exploration of the full PES. The very large number of degrees of freedom makes such study unfeasible. However, with the followed strategy, we are confident that the most populated fragmentation channels are taken into account, and the minima of the relevant energy paths are located.

*Statistical fragmentation simulations.* We study statistical fragmentation of doubly ionized adamantane with the recently developed M3C method<sup>29,30</sup>, using the constrained approach presented in<sup>31</sup>. The key aspect of this methodology is that it provides a random way to move in the phase space (the so-called Markov chain) until a region of maximum entropy is reached, where the physical observables are computed. This description should be equivalent to an MD simulation in the infinite integration time limit. In this work, we focus on two observables: (1) the probability of each fragmentation channel as a function of internal energy (the so-called breakdown curves), and (2) the distribution of the internal energy of the system in its components. This method was successfully used in the past to describe the fragmentation of carbonaceous species<sup>29,31-34</sup>.

The main ingredients of the M3C simulations, i.e. structures, energies, and vibrational frequencies of the fragments, are those obtained in the PES exploration at the B3LYP/6-31G(d) level. In total, 148 molecules are included in the fragmentation model (see SI for details). Geometries are available as an additional file in the SI and in the M3C-store project database<sup>35</sup>.
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Mass spectrum of cationic products resulting from the interaction between XUV pulses and neutral adamantane molecules (purple). The dashed gray lines correspond to the mass spectrum obtained by electron impact at 70 eV³⁹.</figcaption></figure>

The statistical simulations have been carried out such that the sum of angular momenta from all fragments exactly compensates the orbital momentum resulting in a total angular momentum equal to zero. We set the radius of the system to 30.0 Å. Implementation of larger radii implies similar Coulomb interaction among fragments, but requires increased sampling to achieve convergence; on the other hand, a smaller radius results in an artificial overestimation of the fragments' angular momentum. We have performed a scan of the internal energy from 0 to 10 eV. 10000 numerical experiments for each value of internal energy have been carried out in order to estimate the error in the computed observables (and thus to use them as convergence criteria). The numerical experiments each differ from one another in their initial values for vibrational energy, angular momentum, and molecular orientation, which were randomly chosen. All numerical experiments start from the most stable structure of the doubly charged $C_{10}H_{16}^{2+}$. The sequence 5*V, T, R, S: 0, 5*V, T, R, S: -1:1 has been used as a Markov chain, including a total of 2000 events; among them 10% have been used as a burn-in period (see ref. ²⁹ for details).

In summary, a complete picture of the fragmentation of excited doubly-charged adamantane is obtained with the theoretical simulations: dynamic, energetic and entropic approaches allow us to infer the main factors governing the experimentally observed processes, also providing complementary information.

## Results and Discussion

The experimental total mass spectrum of the charged products of interaction (Fig. 3) is obtained by calibration of the TOF spectrum recorded after 275000 laser shots. The most intense peak (excluding helium) corresponds to the singly charged parent ion at *m/z* = 136. The loss of one hydrogen atom is observed at a mass-to-charge ratio of *m/z* = 135 with an intensity of 9% of the parent ion. The losses of two and three hydrogen atoms are also observed, however they have an intensity two orders of magnitude lower than that of the single hydrogen loss. As a general feature, the production of a wide distribution of $C_nH_x^+$ fragments resulting from dissociation of singly and doubly charged adamantane molecules is observed. The most intense peaks of each *C_n* group are attributed to $CH_3^+$, $C_2H_5^+$, $C_3H_5^+$, $C_4H_7^+$, $C_5H_7^+$, $C_6H_7^+$, $C_7H_9^+$, $C_8H_{11}^+$ and $C_9H_{13}^+$ (respectively *m/z* = 15, 29, 41, 55, 67, 79, 93, 107 and 121) and are rather similar to the ones observed in the case of ionization by electron impact at 70 eV (dashed line in Fig. 3). On the other hand, the main fragments of the *C_n* groups *n* = 3, 4 and 8 are different from the ones found by Candian et al.⁹ (photodissociation around first ionization threshold), i.e. $C_3H_7^+$, $C_4H_8^+$ and $C_8H_{12}^+$ respectively, demonstrating that the dynamics of fragmentation is sensitive to the ionization/excitation energy.

Most of the fragments indicate strong intramolecular rearrangements with multiple hydrogen migrations and/or hydrogen losses. Some of these rearrangements may occur before fragmentation and lead to the cage opening of adamantane cations. The corresponding cage opening of the singly charged cation was already studied by Candian et al. using DFT and RRKM simulations⁹. In the case of the doubly charged adamantane, the results of our *ab initio* molecular dynamics calculations (ADMP using DFT-B3LYP up to 500 fs) followed by PES exploration are summarized in Fig. 4. The three lowest energy configurations (Fig. 4(a)), appearing ≈ 4 eV below the double ionization threshold, have an open-cage geometry and have at least one hydrogen migration (CH₃ termination).

**Dication fragmentation pattern.** In order to help to understand the complex energetic picture of the dication dynamics, Fig. 5 shows a schematic dissociative potential energy curve for the dication, including the different energetic quantities that are useful for the discussion. We consider the vertical ionization (Franck-Condon region) from the neutral ground state with a photon of energy *hν* (purple arrow), such as *hν* > *IT*₂. The excess energy after photoionization is defined as *E*<sub>ph</sub> = *hν* − *IT*₂. We already know that the relaxation from the double ionization threshold *IT*₂ to the ground state of the dication has a fixed energy of *E*<sub>relax</sub> ≈ 4.1 eV (see Fig. 4(a)). The kinetic energy of the electron pair involved in the double ionization is called *E*<sub>2e<sup>-</sup></sub> and the internal energy of the dication ground state is denoted *E*<sub>int</sub>. Under our assumptions, the internal energy is the sum of the rapidly redistributed electronic excitation energy (*E*<sub>exc</sub>) and the relaxation energy (*E*<sub>relax</sub>), and thus *E*<sub>int</sub> ≥ *E*<sub>relax</sub> ≈ 4.1 eV. The final ionic products of the dissociation have a total energy equal to the sum of their internal energy (*E*<sub>intfrag</sub>) and the kinetic energy release (KER). This total energy also corresponds to the sum of the
<figure><img src="image_5.png" /><figcaption>**Figure 4.** Key energy levels of adamantane dication processes. (a) Double ionization threshold and lowest energy configurations for doubly charged adamantane found in the PES exploration. (b,c) Final energy levels of the fragmentation channels of the adamantane dication (2- and 3-body breakups) corresponding to the ones in Table 1. Energy levels in gray are not explicitly labelled but can be found in Table 1. The energy values are relative to the neutral ground state in units of eV.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 5.** Schematic of the dissociation of adamantane dication displaying the different energy quantities involved.</figcaption></figure>

initial internal energy of the dication ($E_{int}$) and the energy difference between the dication ground state and the energy levels of the exit channel.

Experimentally, the intact doubly charged parent ion is not observed in the mass spectrum (Fig. 3) at the timescale of the detection (a few microseconds). Moreover, no doubly charged fragments are detected. While the total mass spectrum is dominated by the fragments of singly charged adamantane molecules, the use of the partial covariance technique$^{18,36}$ on the ion TOF spectrum enables to case correlate singly charged ions coming from the dissociation of the dication of adamantane using an ion-ion correlation map representation (Fig. 6). Correlation islands in this map give the ion pairs that are summarized in Table 1 as well as the branching ratios (BR) of these fragmentation channels. In addition, the PES exploration provides the final energy levels of the dication fragmentation channels referred to the neutral ground state (ΔE), which are represented in Fig. 4 and given in Table 1. The energy levels of the 2-body breakup channels appear at lower energies than most of the 3-body breakup ones. It is interesting to notice that all the energy levels are below the double ionization threshold, meaning that the dication of adamantane is metastable and will spontaneously dissociate (if no barrier, *i.e.* transition state, at higher energy than the double ionization threshold is involved).

**Statistical fragmentation simulations.** Assuming a low internal energy after the ionization, the fragmentation time of doubly charged adamantane is expected to be very long (≥ 100 ps, according to our DFTB simulation, see SI) and therefore we cannot afford to carry out simulations with ab initio molecular dynamics. Thus, we use the M3C statistical method to obtain complementary information. This method was developed to study the fragmentation of molecular systems based on entropic criteria (see$^{29–31}$ for details).
<figure><img src="image_7.png" /><figcaption>Figure 6. Ion-ion correlation map resulting from the fragmentation of adamantane dications.</figcaption></figure>

<table><caption>**Table 1.** List of correlated singly charged fragments coming from the dissociation of adamantane dication observed experimentally. In the case of n-body breakups with n > 2, the neutral losses are given in mass losses such that the chemical formulae have to be seen as chemical element indicators and not necessarily as fragments. BR stands for branching ratio and is given in percent. ΔE is the calculated final energy level of the dication fragmentation channels (2- and 3-body breakups) referred to the neutral ground state (in eV).</caption><thead><tr><th colspan="2">Fragment 1</th><th colspan="2">Fragment 2</th><th colspan="2">Neutral loss</th><th>BR</th><th>ΔE</th></tr><tr><th>Formula</th><th>m (a.u.)</th><th>Formula</th><th>m (a.u.)</th><th>Formula</th><th>m (a.u.)</th><th>%</th><th>(eV)</th></tr></thead><tbody><tr><td rowspan="4">CH<sub>3</sub><sup>+</sup></td><td rowspan="4">15</td><td>C<sub>8</sub>H<sub>9</sub><sup>+</sup></td><td>105</td><td>CH<sub>4</sub></td><td>16</td><td>2.8 ± 0.2</td><td>19.27</td></tr><tr><td>C<sub>7</sub>H<sub>9</sub><sup>+</sup></td><td>93</td><td>C<sub>2</sub>H<sub>4</sub></td><td>28</td><td>2.7 ± 0.2</td><td>20.32</td></tr><tr><td>C<sub>7</sub>H<sub>7</sub><sup>+</sup></td><td>91</td><td>C<sub>2</sub>H<sub>6</sub></td><td>30</td><td>4.3 ± 0.2</td><td>19.92</td></tr><tr><td>C<sub>6</sub>H<sub>6</sub><sup>+</sup></td><td>78</td><td>C<sub>3</sub>H<sub>7</sub></td><td>43</td><td>2.3 ± 0.2</td><td>22.24</td></tr><tr><td rowspan="5">C<sub>2</sub>H<sub>5</sub><sup>+</sup></td><td rowspan="5">29</td><td>C<sub>8</sub>H<sub>11</sub><sup>+</sup></td><td>107</td><td>—</td><td>—</td><td>23.3 ± 0.3</td><td>17.21</td></tr><tr><td>C<sub>8</sub>H<sub>9</sub><sup>+</sup></td><td>105</td><td>H<sub>2</sub></td><td>2</td><td>6.3 ± 0.2</td><td>18.11</td></tr><tr><td>C<sub>7</sub>H<sub>8</sub><sup>+</sup></td><td>92</td><td>CH<sub>3</sub></td><td>15</td><td>5.4 ± 0.2</td><td>20.18</td></tr><tr><td>C<sub>7</sub>H<sub>7</sub><sup>+</sup></td><td>91</td><td>CH<sub>4</sub></td><td>16</td><td>6.1 ± 0.2</td><td>18.06</td></tr><tr><td>C<sub>6</sub>H<sub>7</sub><sup>+</sup></td><td>79</td><td>C<sub>2</sub>H<sub>4</sub></td><td>28</td><td>7.1 ± 0.2</td><td>19.05</td></tr><tr><td rowspan="6">C<sub>3</sub>H<sub>5</sub><sup>+</sup></td><td rowspan="6">41</td><td>C<sub>7</sub>H<sub>11</sub><sup>+</sup></td><td>95</td><td>—</td><td>—</td><td>2.5 ± 0.2</td><td>18.26</td></tr><tr><td>C<sub>6</sub>H<sub>8</sub><sup>+</sup></td><td>80</td><td>CH<sub>3</sub></td><td>15</td><td>5.3 ± 0.2</td><td>20.31</td></tr><tr><td>C<sub>5</sub>H<sub>7</sub><sup>+</sup></td><td>67</td><td>C<sub>2</sub>H<sub>4</sub></td><td>28</td><td>5.3 ± 0.2</td><td>19.78</td></tr><tr><td>C<sub>4</sub>H<sub>7</sub><sup>+</sup></td><td>55</td><td>C<sub>3</sub>H<sub>4</sub></td><td>40</td><td>4.6 ± 0.2</td><td>21.06</td></tr><tr><td>C<sub>4</sub>H<sub>6</sub><sup>+</sup></td><td>54</td><td>C<sub>3</sub>H<sub>5</sub></td><td>41</td><td>6.6 ± 0.2</td><td>21.42</td></tr><tr><td>C<sub>7</sub>H<sub>10</sub><sup>+</sup></td><td>94</td><td>—</td><td>—</td><td>2.8 ± 0.2</td><td>18.90</td></tr><tr><td rowspan="3">C<sub>3</sub>H<sub>6</sub><sup>+</sup></td><td rowspan="3">42</td><td>C<sub>6</sub>H<sub>7</sub><sup>+</sup></td><td>79</td><td>CH<sub>3</sub></td><td>15</td><td>6.4 ± 0.2</td><td>20.27</td></tr><tr><td>C<sub>5</sub>H<sub>6</sub><sup>+</sup></td><td>66</td><td>C<sub>2</sub>H<sub>4</sub></td><td>28</td><td>3.7 ± 0.2</td><td>20.80</td></tr><tr><td>C<sub>6</sub>H<sub>9</sub><sup>+</sup></td><td>81</td><td>—</td><td>—</td><td>2.3 ± 0.2</td><td>18.39</td></tr></tbody></table>

**Breakdown curves.** In Fig. 7, showing the breakdown curves, the gray areas mark the inaccessible regions of internal energy ($E_{int} < E_{relax}$). The probabilities for the 2-body channels (Fig. 7(a)) all peak below $E_{relax}$, with tails reaching into the accessible internal energy region. This is consistent with the fact that they appear in the lowest energy region in the PES (Fig. 4). Above $E_{relax}$, only the $C_2H_5^+/C_8H_{11}^+$ has a significant probability, while the other 2-body breakup channels $C_3H_5^+/C_7H_{11}^+$, $C_3H_6^+/C_7H_{10}^+$ and $C_4H_7^+/C_6H_9^+$ are almost not populated. The 3-body breakup channels (Fig. 7(b)) all peak around 6 eV, also consistent with their higher energy levels according to the PES in Fig. 4.
<figure><img src="image_8.png" /><figcaption>**Figure 7.** Breakdown curves of the channels observed in the experiments for 2-body breakups (a) and 3-body breakups (b). The errors (shaded areas around curves) correspond to the standard deviation. The black dashed lines indicate the minimum of internal energy that we can reach in our case, meaning $E_{relax}$, and the gray areas mark the regions of internal energy that are inaccessible in the experiment.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 8.** The average of the energy components $E_{int}$, $E_{vib}$, and $E_{rot}$ is shown as a function of the internal energy. The errors (shaded areas around curves) correspond to the standard deviation.</figcaption></figure>

Comparing the calculated breakdown curves with the experimentally measured branching ratios in Table 1, suggests that the internal energy of the adamantane dication is close to $E_{relax}$ under the current conditions, since in this region the $C_2H_2/C_8H_{11}$ channel dominates the other 2-body breakup channels, in good agreement with the experiment (BR > 23% vs. BR < 3%). If the internal energy was higher, the 3-body breakup channels would start to dominate, which is not observed in the experiment.

Under our assumption that the electronic excitation energy is considered to be rapidly redistributed into the nuclear degrees of freedom, this implies that most of the dication remain in the ground state or low excited states after ionization, and that the internal energy of the system mainly corresponds to the relaxation energy.

*Energy storage in the fragments.* We have further obtained additional valuable information by analyzing how the internal energy is distributed after ionization using the statistical simulations with the M3C code. Note that this analysis is based on the ergodic assumption, i.e. at infinite time when the system has reached equilibrium and the maximum entropy region in the phase space is populated.

We can decompose the internal energy as

$$E_{int} = E_{kin} + E_{pot} + E_{vib} + E_{rot}$$

where $E_{kin}$, $E_{pot}$, $E_{vib}$, and $E_{rot}$ represent the kinetic or translational, the potential, the vibrational, and the rotational energy components, respectively. The potential energy $E_{pot}$ is the energy difference between different geometrical configurations, in this case the exit channel and the most stable dication structure (see Fig. 4). While the potential energy is important for the total available energy, we now focus on how the latter is distributed between the remaining degrees of freedom, i.e. between $E_{kin}$, $E_{vib}$ and $E_{rot}$. Fig. 8 shows the average and the standard deviation of these energy components as a function of the internal energy $E_{int}$. It is clear that in the considered internal energy range the vibrational contribution is most prominent, with smaller contributions of the rotational and
<figure><img src="image_10.png" /><figcaption>**Figure 9.** VMI images obtained after filtering using the partial covariance method on the TOF peaks correspond to the fragments C₂H₅⁺ (a) and C₈H₁₁⁺ (b) (left part: raw data and right part: inverted data). Artefacts are coming from intense signal (helium and parent ion) that were not filtered out by covariance analysis (see SI for more details). (c,d) Ion kinetic energy distributions of the respective fragments obtained by angular integration of the inverted data avoiding the artefacts signal and energy calibrated using ion trajectory simulations (SIMION⁴⁰). (e) Kinetic energy release distribution (KERd) for the channel C₂H₅⁺ / C₈H₁₁⁺ obtained by convolution of the kinetic energy distributions of the two fragments.</figcaption></figure>

kinetic energy components. As already shown, the fragmentation is the dominant process in the relaxation of the adamantane dication; thus, the available energy of the system is primarily absorbed by the vibrational component, i.e. the produced fragments can store a large amount of energy in nuclear degrees of freedom (mainly vibrations). Assuming an internal energy of ~4.1 eV (corresponding to the relaxation energy), ~70% of the available energy is stored in vibration and ~15% in rotation while the remaining 15% are shared among the other components.

**The main fragmentation channel C₂H₅⁺/C₈H₁₁⁺ - experimental evidence of a two-step process.**

As we have seen, the channel C₂H₅⁺/C₈H₁₁⁺ is strongly dominating the fragmentation dynamics of the dication with a branching ratio ~23.3% whereas the other channels are at least three times less intense (Fig. 6 and Table 1). This can be roughly interpreted by looking at the energy levels of the different fragmentation channels obtained from the exploration of the PES of the adamantane dication (Fig. 4 and Table 1). This main channel is energetically favorable since it has the lowest energy level at around ~7 eV below the double ionization threshold. In addition, we have seen from the M3C calculations that the breakdown curve (Fig. 7(a)) of this channel was dominant at low accessible internal energy (~4 – 5 eV), thus being also entropically favorable in this energy region.
<figure><img src="image_11.png" /><figcaption>**Figure 10.** (a) Total photoelectron spectrum (purple). The black dashed lines indicate the position of helium photoelectrons. The gray line is the calculated sum of helium photoelectrons (Gaussian) and of the photoelectrons coming from the singly charged adamantane (taken from³⁸) considering our XUV spectrum and the photoabsorption cross section of adamantane⁸. (b) Estimated photoelectron spectrum coming from double ionization of adamantane obtained as the difference between the purple and the gray curves of panel (a).</figcaption></figure>

In order to have a better insight into the energetics of the main fragmentation channel, we can regard the ion kinematics of the dissociation process, particularly the kinetic energy release distribution (KERd). Performing partial covariance analysis between the ion TOF and the ion VMI data gives the "mass-selected" velocity map images displayed on the left parts of panels (a,b) in Fig. 9. The right parts in Fig. 9(a,b) show the result of Abel inversion using an iterative method³⁷. The angular integration of the inverted images gives, after energy calibration, the kinetic energy distributions of the fragments C₂H₅⁺ and C₈H₁₁⁺ (Fig. 9(c,d)). The intense signal close to zero corresponds to the contribution from the dissociation of the singly charged adamantane. At higher energies, clear peaks show contributions at ~2.5 eV and ~0.7 eV respectively due to the Coulomb repulsion of the 2-body breakup of the dication and verify the momentum conservation principle. It is possible to obtain the KERd of the channel C₂H₅⁺/C₈H₁₁⁺ by convolution of the two individual kinetic energy distributions (Fig. 9(e)). The peak at ~3 eV indicates the main energy contribution of this channel.

The small value of the KER reflects complex fragmentation dynamics with strong molecular rearrangement before the charge separation takes place: In a first step, a cage opening leading, most probably, to one of the structures in Fig. 4; then a second molecular reorganization producing both charged fragments; and finally charge repulsion between them. The energy difference between the open structures in Fig. 4 and the C₂H₅⁺/C₈H₁₁⁺ exit channel is between ~2.6 and ~2.9 eV, which is consistent within a multiple step fragmentation as the one presented here. Considering the Coulomb repulsion between the two positive charges, the energy released is C/R, with C = 14.4 eV · Å and R the initial inter-charge distance in [Å]. An energy release of ~3.0 eV, as measured in the experiment, corresponds to an inter-charge distance of 4.8 Å. However, the maximum distance between two C atoms in the closed-cage adamantane structure is ~3.5 Å, discarding the assumption of an instantaneous double ionization followed by prompt fragmentation. Thus, structural rearrangements before the charge separation would produce a considerable extension of the structure, thus increasing the inter-charge distance up to 4.8 Å (as inferred from the experiment). This further confirms the multi-step processes with cage opening preceding Coulomb repulsion.

**Photoelectron spectrum.** Figure 10(a) (purple line) shows the total photoelectron energy spectrum obtained by angular integration of the electron VMI data after inversion using an iterative method³⁷. Electrons coming from the ionization of the helium buffer gas can be seen around the dashed lines and were used for energy calibration. The contribution to the photoelectron spectrum associated with single ionization by the harmonics is found in the higher energy part of the total spectrum. The expected contribution, shown by the gray line Fig. 10(a), is calculated using the measured photoelectron spectrum for singly charged adamantane³⁸, the harmonic spectrum (see Fig. 2), and the photoabsorption cross section of adamantane⁸. In this estimated contribution, the calculated photoelectron spectrum from ionization of helium has also been included, and by subtracting this from the total photoelectron spectrum, the expected spectrum of the photoelectrons associated with double ionization is obtained and shown in Fig. 10(b).

Although the absolute signal in the resulting spectrum is sensitive to the scaling of the calculated spectrum, it is clear that the photoelectrons primarily occupy the low-energy part of the spectrum (<10 eV) compared to the photoelectrons from single ionization. From the harmonic spectrum (Fig. 2) the maximum total electron pair energy is E₂e⁻ ≈ 21 eV, assuming that the dications remain in the ground state or in low excited states after ionization. Thus, the observed cut-off at ~10 eV, suggests that the available energy is shared rather evenly between the two electrons. All-in-all, the shape of the photoelectron energy distribution supports the assumption of low excitation, apart from the unexpectedly strong contribution at energies below 2 eV. The latter feature is a possible indication of excitation of higher-lying electronic states, resulting in low energy electrons. Possible candidates
for such states are found through a calculation of the excited states of the adamantane dication (see SI for more details) exhibiting a dense band of excited states between 35 and 40 eV (relative to the ground state of the neutral) that would result in large internal energy values $E_{\text{int}}$ between 15 and 20 eV, and values of $E_{2e-}$ in the 5–10 eV range following ionization by the cut-off harmonics. Such excitation could not be inferred from the experimental fragmentation pattern, nor is included in the current level of theory, which calls for further investigation.

## Conclusion

We have performed a detailed study of the photodissociation of adamantane, focused on the fragmentation dynamics of the dication. By combining the experimental analysis with multiple theoretical methods we unraveled key processes governing charge and energy distribution after the ultrafast photoionization and the subsequent fragmentation dynamics.

We found that the most stable structures of the dication of adamantane present an open-cage geometry, appearing at ~4 eV below the double ionization threshold, that can be reached in a few tens of femtoseconds after the ionization. However, these structures are metastable and evolve producing several fragments in a Coulomb repulsion process. Much like other carbonaceous species, adamantane dications are quite efficient energy reservoirs, being able to store a large amount of energy in particular in vibrational modes. This occurs when the internal energy generated in the photoionization, together with the energy produced in the Coulomb explosion, is redistributed among the nuclear degrees of freedom of the produced fragments.

Among the different fragmentation pathways, the most populated channel $C_2H_5^+/C_8H_{11}^+$ is the lowest in energy in the PES (~2.6 eV below the most stable structure of the adamantane dication) leading us to conclude that the fragment distribution is mainly governed by energetic criteria. We have shown that the ion KERd of this channel peaks at ~3 eV with a width of 1 eV allowing us to experimentally confirm that the cage opening takes place prior to fragmentation. A qualitative comparison between the experimental branching ratios and the results of statistical fragmentation simulations suggests that the internal energy of the dication largely consists of the relaxation energy from the cage opening, with only minor contribution from redistribution of electronic excitation energy. Finally, measurements of the photoelectron kinetic energy spectrum largely confirms these observations, but also indicate the possible existence of electronic excitation which could not be further investigated in the current experiment.

The presented results highlight the complexity of, and provides pieces of information on, the fragmentation of multiply charged diamondoids. While this study was able to identify the dominant fragmentation pathways and shine light on the redistribution of energy and charge, the further elucidation of the ultrafast excitation dynamics, and in particular the timescales involved, calls for time-resolved experiments. Such experiments can be envisaged using ultrashort single wavelength XUV pulses, e.g. from free electron lasers, in combination with multicoincidence ion-electron spectroscopy techniques.

Received: 11 December 2019; Accepted: 29 January 2020;

Published online: 19 February 2020

## References

1. Schwertfeger, H., Fokin, A. & Schreiner, P. Diamonds are a chemist's best friend: Diamondoid chemistry beyond adamantane. *Angewandte Chemie Int. Ed.* **47**, 1022–1036, https://doi.org/10.1002/anie.200701684 (2008).

2. Stauss, S. & Terashima, K. *Diamondoids: Synthesis, Properties, and Applications* (CRC Press, 2017).

3. Anders, E. & Zinner, E. Interstellar grains in primitive meteorites: Diamond, silicon carbide, and graphite. *Meteorit.* **28**, 490–514, https://doi.org/10.1111/j.1945-5100.1993.tb00274.x (1993).

4. Henning, T. & Salama, F. Carbon in the universe. *Sci.* **282**, 2204–2210, https://doi.org/10.1126/science.282.5397.2204 (1998).

5. Oomens, J. *et al*. Infrared spectroscopic investigation of higher diamondoids. *J. Mol. Spectrosc.* **238**, 158–167, https://doi.org/10.1016/j.jms.2006.05.001 (2006).

6. Acke, B. & van den Ancker, M. E. A survey for nanodiamond features in the 3 micron spectra of herbig ae/be stars. *Astron. & Astrophys.* **457**, 171–181, https://doi.org/10.1051/0004-6361:20065612 (2006).

7. Lenzke, K. *et al*. Experimental determination of the ionization potentials of the first five members of the nanodiamond series. *The J. Chem. Phys.* **127**, 084320, https://doi.org/10.1063/1.2773725 (2007).

8. Steglich, M., Huisken, F., Dahl, J. E., Carlson, R. M. K. & Henning, T. Electronic spectroscopy of fuv-irradiated diamondoids: a combined experimental and theoretical study. *The Astrophys. J.* **729**, 91, https://doi.org/10.1088/0004-637x/729/2/91 (2011).

9. Candian, A., Bouwman, J., Hemberger, P., Bodi, A. & Tielens, A. G. G. M. Dissociative ionisation of adamantane: a combined theoretical and experimental study. *Phys. Chem. Chem. Phys.* **20**, 5399–5406, https://doi.org/10.1039/C7CP05957D (2018).

10. Marciniak, A. *et al*. Electron correlation driven non-adiabatic relaxation in molecules excited by an ultrashort extreme ultraviolet pulse. *Nat. communications* **10**, 1–8, https://doi.org/10.1038/s41467-018-08131-8 (2019).

11. Ferray, M. *et al*. Multiple-harmonic conversion of 1064 nm radiation in rare gases. *J. Phys. B: At. Mol. Opt. Phys.* **21**, L31–L35, https://doi.org/10.1088/0953-4075/21/3/001 (1988).

12. McPherson, A. *et al*. Studies of multiphoton production of vacuum-ultraviolet radiation in the rare gases. *J. Opt. Soc. Am. B* **4**, 595–601, https://doi.org/10.1364/JOSAB.4.000595 (1987).

13. Manschwetus, B. *et al*. Two-photon double ionization of neon using an intense attosecond pulse train. *Phys. Rev. A* **93**, 061402, https://doi.org/10.1103/PhysRevA.93.061402 (2016).

14. Coudert-Alteirac, H. *et al*. Micro-focusing of broadband high-order harmonic radiation by a double toroidal mirror. *Appl. Sci.* **7**, 1159, https://doi.org/10.3390/app7111159 (2017).

15. Even, U. Pulsed supersonic beams from high pressure source: Simulation results and experimental measurements. *Adv. Chem.* **2014**, 11, https://doi.org/10.1155/2014/636042 (2014).

16. Even, U. The even-lavie valve as a source for high intensity supersonic beam. *EPJ Tech. Instrumentation* **2**, 17, https://doi.org/10.1140/epjti/s40485-015-0027-5 (2015).

17. Rading, L. *et al*. A versatile velocity map ion-electron covariance imaging spectrometer for high-intensity xuv experiments. *Appl. Sci.* **8**, 998, https://doi.org/10.3390/app8060998 (2018).

18. Frasinski, L. J. Covariance mapping techniques. *J. Phys. B: At. Mol. Opt. Phys.* **49**, 152004, https://doi.org/10.1088/0953-4075/49/15/152004 (2016).
19. Porezag, D., Frauenheim, T., Köhler, T., Seifert, G. & Kaschner, R. Construction of tight-binding-like potentials on the basis of density-functional theory: Application to carbon. *Phys. Rev. B* **51**, 12947–12957, https://doi.org/10.1103/PhysRevB.51.12947 (1995).

20. Seifert, G., Porezag, D. & Frauenheim, T. Calculations of molecules, clusters, and solids with a simplified Icao-dft-Ida scheme. *International Journal of Quantum Chemistry*, **58**, 185–192, 10.1002/(SICI)1097-461X(1996)58:2<185::AID-QUA7>3.0.CO;2-U (1996).

21. Heine, T. et al. demon-nano Jacobs University Bremen: Bremen, Germany (2009).

22. Lee, C., Yang, W. & Parr, R. G. Development of the colle-salvetti correlation-energy formula into a functional of the electron density. *Phys. Rev. B* **37**, 785–789, https://doi.org/10.1103/PhysRevB.37.785 (1988).

23. Miehlich, B., Savin, A., Stoll, H. & Preuss, H. Results obtained with the correlation energy density functionals of becke and lee, yang and parr. *Chem. Phys. Lett.* **157**, 200–206, https://doi.org/10.1016/0009-2614(89)87234-3 (1989).

24. Becke, A. D. Density-functional thermochemistry. iii. the role of exact exchange. *The J. Chem. Phys.* **98**, 5648–5652, https://doi.org/10.1063/1.464913 (1993).

25. Frisch, M. J. et al. Gaussian 09 Revision E.01 (2013). Gaussian Inc. Wallingford CT 2013.

26. Maclot, S. et al. Dynamics of glycine dications in the gas phase: ultrafast intramolecular hydrogen migration versus coulomb repulsion. *The J. Phys. Chem. Lett.* **4**, 3903–3909, https://doi.org/10.1021/jz4020234 (2013).

27. Piekarski, D. G. et al. Unusual hydroxyl migration in the fragmentation of β-alanine dication in the gas phase. *Phys. Chem. Chem. Phys.* **17**, 16767–16778, https://doi.org/10.1039/C5CP01628B (2015).

28. Maclot, S. et al. Determination of energy-transfer distributions in ionizing ion-molecule collisions. *Phys. Rev. Lett.* **117**, 073201, https://doi.org/10.1103/PhysRevLett.117.073201 (2016).

29. Aguirre, N. F., Díaz-Tendero, S., Hervieux, P.-A., Alcamí, M. & Martín, F. M3C: A Computational Approach To Describe Statistical Fragmentation of Excited Molecules and Clusters. *J. Chem. Theory Comput.* **13**, 992–1009, https://doi.org/10.1021/acs.jctc.6b00984 (2017).

30. Aguirre, N. F., Díaz-Tendero, S., Hervieux, P.-A., Alcamí, M. & Martín, F. M3C Software Package, https://github.com/nfaguirrec/M3C, https://doi.org/10.5281/zenodo.3533987 (2019).

31. Aguirre, N. F. et al. Fully versus constrained statistical fragmentation of carbon clusters and their heteronuclear derivatives. *The J. Chem. Phys.* **150**, 144301, https://doi.org/10.1063/1.5083864 (2019).

32. Martinet, G. et al. Fragmentation of highly excited small neutral carbon clusters. *Phys. Rev. Lett.* **93**, 063401, https://doi.org/10.1103/PhysRevLett.93.063401 (2004).

33. Díaz-Tendero, S., Hervieux, P.-A., Alcamí, M. & Martín, F. Statistical fragmentation of small neutral carbon clusters. *Phys. Rev. A* **71**, 033202, https://doi.org/10.1103/PhysRevA.71.033202 (2005).

34. Erdmann, E., Labuda, M., Aguirre, N. F., Díaz-Tendero, S. & Alcamí, M. Furan fragmentation in the gas phase: New insights from statistical and molecular dynamics calculations. *The J. Phys. Chem. A* **122**, 4153–4166, https://doi.org/10.1021/acs.jpca.8b00881 (2018).

35. Aguirre, N. F. et al. M3C-store dataset, https://github.com/nfaguirrec/M3C-store, https://doi.org/10.5281/zenodo.2632188 (2019).

36. Frasinski, L. J., Codling, K. & Hatherly, P. A. Covariance mapping: A correlation method applied to multiphoton multiple ionization. *Sci. 246*, 1029–1031, https://doi.org/10.1126/science.246.4933.1029 (1989).

37. Vrakking, M. J. J. An iterative procedure for the inversion of two-dimensional ion/photoelectron imaging experiments. *Rev. Sci. Instruments* **72**, 4084–4089, https://doi.org/10.1063/1.1406923 (2001).

38. Schmidt, W. Photoelectron spectra of diamondoid molecules: Adamantane, silamantane and urotropine. *Tetrahedron* **29**, 2129–2134, https://doi.org/10.1016/0040-4020(73)80153-X (1973).

39. Wallace, W. E. "Mass Spectra" by NIST Mass Spectrometry Data Center in NIST Chemistry WebBook, NIST Standard Reference Database Number 69, Eds. P.J. Linstrom and W.G. Mallard (National Institute of Standards and Technology, Gaithersburg MD, 20899, 2019).

40. Dahl, D. A. Simion for the personal computer in reflection. *Int. J. Mass Spectrom.* **200**, 3–25 (2000).

## Acknowledgements

A special thanks is given to A. Persson for his support for the high-power laser and to M. Gisselbrecht for critical reading of the manuscript. This research was supported by the Swedish Research Council, the Swedish Foundation for Strategic Research and the Crafoord Foundation. This project has received funding from the European Union's Horizon 2020 research and innovation program under the Marie Skłodowska-Curie Grant Agreement No. 641789 MEDEA and under grant agreement no 654148 Laserlab-Europe. The research was conducted in the framework of the International Associated Laboratory (LIA) Fragmentation DYNAmics of complex MOlecular systems - DYNAMO, funded by the Centre National de la Recherche Scientifique (CNRS). The authors acknowledge the generous allocation of computer time at the Centro de Computación Científica at the Universidad Autónoma de Madrid (CCC-UAM). This work was partially supported by the project CTQ2016-76061-P of the Spanish Ministerio de Economía y Competitividad (MINECO). Financial support from the MINECO through the "María de Maeztu" Program for Units of Excellence in R&D (MDM-2014-0377) is also acknowledged.

## Author contributions

P.R., P.J. and S.M. designed the experiments, S.M., J.L., J.P., H.W., P.R., P.R., S.I., F.B., H.C.-A. and B.A.H. performed the experiment, S.M. and J.L. analysed the results. S.D.-T. performed the potential energy surface and the molecular dynamics calculations, and N.F.A. the M3C simulations. S.M., S.D.-T. and P.J. wrote the original draft, S.D.-T. and N.F.A. wrote the theoretical sections and part of the discussion and all authors reviewed the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41598-020-59649-1.

Correspondence and requests for materials should be addressed to S.M. or P.J.

Reprints and permissions information is available at www.nature.com/reprints.
**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_12.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020