# Non-thermal structural transformation of diamond driven by x-rays

Heimann Phil, Nicholas J Hartley, Inoue Ichiro, Tkachenko Victor, Antoine Andre, Dorchies Fabien, Falcone Roger, Gaudin Jérôme, Höppner Hauke, Inubushi Yuichi, et al.

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

## To cite this version:

Heimann Phil, Nicholas J Hartley, Inoue Ichiro, Tkachenko Victor, Antoine Andre, et al. Non-thermal structural transformation of diamond driven by x-rays. Structural Dynamics, 2023, 10 (5), pp.054502. 10.1063/4.0000193. hal-04262209

HAL Id: hal-04262209

https://hal.science/hal-04262209v1

Submitted on 5 Nov 2023

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Non-thermal structural transformation of diamond driven by x-rays

Philip Heimann ⚙️ ID ; Nicholas J. Hartley ID ; Ichiro Inoue; Victor Tkachenko ⚙️ ID ; Andre Antoine ID ;

Fabien Dorchies ID ; Roger Falcone ID ; Jérôme Gaudin ID ; Hauke Höppner ID ; Yuichi Inubushi ID ;

Konrad J. Kapcia ID ; Hae Ja Lee ID ; Vladimir Lipp ID ; Paloma Martinez ID ; Nikita Medvedev ID ;

Franz Tavella ID ; Sven Toleikis ID ; Makina Yabashi ID ; Toshinori Yabuuchi ID ; Jumpei Yamada ID ;

Beata Ziaja ID

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

Check for updates

Struct. Dyn. 10, 054502 (2023)

https://doi.org/10.1063/4.0000193

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

View
Online

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

Export
Citation

CrossMark

<figure><img src="image_6.png" /></figure>
# Non-thermal structural transformation of diamond driven by x-rays

Cite as: Struct. Dyn. **10**, 054502 (2023); doi: 10.1063/4.0000193

Submitted: 18 May 2023 · Accepted: 10 October 2023 ·

Published Online: 27 October 2023

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

View Online

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

Export Citation

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

CrossMark

Philip Heimann,<fnref n="1" />,ª ID Nicholas J. Hartley,<fnref n="2" /> ID Ichiro Inoue,<fnref n="3" /> Victor Tkachenko,<fnref n="4" />,<fnref n="5" />,ª ID Andre Antoine,<fnref n="6" /> ID Fabien Dorchies,<fnref n="7" /> ID Roger Falcone,<fnref n="8" /> ID Jérôme Gaudin,<fnref n="7" /> ID Hauke Höppner,<fnref n="9" /> ID Yuichi Inubushi,<fnref n="3" />,<fnref n="10" /> ID Konrad J. Kapcia,<fnref n="5" />,<fnref n="11" /> ID Hae Ja Lee,<fnref n="1" /> ID Vladimir Lipp,<fnref n="5" />,<fnref n="12" /> ID Paloma Martinez,<fnref n="7" /> ID Nikita Medvedev,<fnref n="13" />,<fnref n="14" /> ID Franz Tavella,<fnref n="1" /> ID Sven Toleikis,<fnref n="15" /> ID Makina Yabashi,<fnref n="3" />,<fnref n="10" /> ID Toshinori Yabuuchi,<fnref n="3" />,<fnref n="10" /> ID Jumpei Yamada,<fnref n="3" /> ID and Beata Ziaja<fnref n="5" />,<fnref n="12" /> ID

## AFFILIATIONS

## ABSTRACT

Intense x-ray pulses can cause the non-thermal structural transformation of diamond. At the SACLAC XFEL facility, pump x-ray pulses triggered this phase transition, and probe x-ray pulses produced diffraction patterns. Time delays were observed from 0 to 250 fs, and the x-ray dose varied from 0.9 to 8.0 eV/atom. The intensity of the (111), (220), and (311) diffraction peaks decreased with time, indicating a disordering of the crystal lattice. From a Debye-Waller analysis, the rms atomic displacements perpendicular to the (111) planes were observed to be significantly larger than those perpendicular to the (220) or (311) planes. At a long time delay of 33 ms, graphite (002) diffraction indicates that graphitization did occur above a threshold dose of 1.2 eV/atom. These experimental results are in qualitative agreement with XTANT+ simulations using a hybrid model based on density-functional tight-binding molecular dynamics.

© 2023 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/). https://doi.org/10.1063/4.0000193

## I. INTRODUCTION

In metals, thermal melting induced by short pulses of light occurs on picosecond time scales, triggered by electron-phonon coupling.<fnref n="1" /> In covalently bonded materials, exposure to a short burst of either x-ray

or optical radiation<fnref n="2" /> may induce a non-thermal phase transition. For a non-thermal phase transition, calculations predict that the excitation of a few percent of the valence band electrons leads to a modification of the potential energy surface, which results in a change in the atomic

<footnote n="1">¹Linac Coherent Light Source, SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA</footnote>

<footnote n="2">²SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA</footnote>

<footnote n="3">³RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo, Hyogo 679-5148, Japan</footnote>

<footnote n="4">⁴European XFEL GmbH, Holzkoppel 4, 22869 Schenefeld, Germany</footnote>

<footnote n="5">⁵Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestr. 85, 22607 Hamburg, Germany</footnote>

<footnote n="6">⁶University of Michigan, 500 S State St, Ann Arbor, Michigan 48109, USA</footnote>

<footnote n="7">⁷University Bordeaux, CNRS, CEA, CELIA, UMR 5107, F-33500 Talence, France</footnote>

<footnote n="8">⁸Department of Physics, University of California, Berkeley, California 94720, USA</footnote>

<footnote n="9">⁹Helmholtz-Zentrum Dresden-Rossendorf, 01328 Dresden, Germany</footnote>

<footnote n="10">¹⁰Japan Synchrotron Radiation Research Institute, 1-1-1 Kouto, Sayo, Hyogo 679-5198, Japan</footnote>

<footnote n="11">¹¹Institute of Spintronics and Quantum Information, Faculty of Physics, Adam Mickiewicz University in Poznań, ul. Uniwersytetu Poznańskiego 2, 61-614 Poznań, Poland</footnote>

<footnote n="12">¹²Institute of Nuclear Physics, Polish Academy of Sciences, Radzikowskiego 152, 31-342 Krakow, Poland</footnote>

<footnote n="13">¹³Institute of Physics, Czech Academy of Sciences, Na Slovance 1999/2, 182 21 Prague 8, Czech Republic</footnote>

<footnote n="14">¹⁴Institute of Plasma Physics, Czech Academy of Sciences, Za Slovankou 3, 182 00 Prague 8, Czech Republic</footnote>

<footnote n="15">¹⁵Deutsches Elektronen-Synchrotron DESY, Notkestr. 85, 22607 Hamburg, Germany</footnote>

<footnote>ª) Authors to whom correspondence should be addressed: paheim@slac.stanford.edu and victor.tkachenko@xfel.eu</footnote>
structure on femtosecond timescales.³,⁴ With x-ray pump pulses focused on diamond, previous work has observed a non-thermal structural disordering⁷ and a transition to another crystalline phase, graphite, at lower pump intensities.⁶,⁷ For the solid-to-solid transformation, there is indirect evidence that the phase transition is non-thermal from the measurement of the optical transmission.⁷ Calculations have also predicted the x-ray-induced phase transition from diamond to graphite⁴ although more recent simulations indicate partial graphitization.⁸

Reitze et al. conducted optical pump-optical probe measurements of graphite and diamond and derived optical properties for liquid carbon.⁹ In comparison with the optical pump-optical probe technique, the x-ray pump-probe method involves a 'simpler' interaction of the x-rays with materials (with photoabsorption being the only channel for the absorption of the x-ray energy) and an uniform absorption with depth (in contrast to the absorption of optical radiation, which is strongly non-uniform with depth).¹⁰ High energies of the probing x-ray photons can provide very high spatial resolution for the structural changes observed.¹¹ On the other hand, with the x-ray pump-probe measurements, there is a spatial gradient of pulse intensity and a spatial transport of the energetic photo and Auger electrons. The latter can significantly decrease the energy deposition in the beam focus. Both these effects imply that the observed diffraction signal is volume integrated, i.e., it collects the contributions from differently illuminated regions of the material. This makes a comparison between experimental and theory predictions not straightforward, as the latter should also include volume integration of the predicted diffraction signal. Gaudin et al.⁶ examined the threshold fluence of graphitization in irradiated diamond for a variety of extreme ultraviolet wavelengths and demonstrated from *ex situ* Raman measurements that a phase transition does occur. Using soft x-ray irradiation, Tavella et al.⁷ observed the graphitization phase transition indirectly through changes in the optical transmission.

The interaction of diamond with high-intensity x-rays has applications in many areas of high energy density physics. One such area of interest is found in inertial confinement fusion (ICF); recent experiments performed at the National Ignition Facility have achieved a burning plasma state using nanocrystalline diamond as the ablator material.¹² Also, many devices used in x-ray free-electron laser (FEL) beamlines involve diamond and silicon.¹³⁻¹⁵ Thus, a good understanding of the damage mechanisms within these materials is of great practical utility for the operation of x-ray FELs. These studies can give further prospects for the ultrafast processing of materials with intense x-ray pulses.¹⁶ In general, the elucidation of diamond's response to x-ray radiation will provide us with fundamental knowledge related to how changes in a material's electronic structure can affect its atomic structure and, hence, provide new insight into several aspects of high energy density physics and materials science.

Under femtosecond x-ray irradiation, covalently bonded materials, such as diamond, undergo a sequence of processes.¹⁷,¹⁸ First, photoabsorption promotes electrons from the bound states of the deep atomic shells (K-shell for carbon) or valence band to the conduction band. The deep-shell holes can then decay through Auger processes, which are the predominant relaxation channels for low-Z (light) elements.⁴ The Auger decay results in more electronic excitations from the valence band to the conduction band, following the relaxation of the core holes into the valence band. The released photo- and Auger electrons scatter further via inelastic (impact ionization of valence

band or deep-shell electrons) and elastic channels. In the case of 10 keV photoelectrons, the impact ionization cascading is completed on a 100-femtosecond timescale. It finishes when the electrons do not have sufficient energy to perform more impact ionizations.¹⁹ Meanwhile, initiated by 10 keV photoelectrons, the spatial distribution of the secondary electron cloud expand to 600 nm.²⁰

Ziaja and co-workers have developed a theoretical framework for modeling x-ray-induced phase transitions using hybrid models. The resulting simulation tools are (i) the XTANT code based on tight-binding molecular dynamics⁴,²¹,²² and (ii) the XTANT+ code based on density-functional tight-binding molecular dynamics.⁸,¹¹ These models include photoabsorption, Auger decay, and collisional scattering processes occurring in the electronic system.⁴,⁸ Because of these processes, the occupation of electronic bands transiently changes, causing the evolution of the potential energy surface. In diamond, populating a sufficient number of antibonding states could lead to an ultrafast rearrangement of atoms from sp³ to sp² bonding, causing its transition to graphite-like structures. The present experiment provides a test of the accuracy of these theoretical methods. This can validate their applicability to obtain predictions on the x-ray interaction with materials. For the present experiment, we use the XTANT+ tool, first introduced in 2022.⁸

For the observation of graphitization, this experiment was performed at the SACLA x-ray free electron laser facility²³ operating in the two-color mode.²⁴ The first x-ray pulse triggered a phase transition. The second x-ray pulse probed the structural change through x-ray diffraction. The intensity of the diamond reflections was monitored with respect to the time separation between the two x-ray pulses. These experimental results are then compared with XTANT+ calculations.

## II. EXPERIMENTAL

The SACLA x-ray Free Electron Laser was operated in the split undulator mode generating x-ray pump and probe pulses.²⁴ The two pulses consisted of pump pulses at a photon energy of 7 keV and probe pulses at 10.5 keV, both ~6 fs in duration.²⁵,²⁶ The jitter between the x-ray pump and probe pulses is less than 1 fs.²⁴ The pump-probe time delay was varied between 0 and 250 fs. A diamond spectrometer measured the x-ray spectrum and, in particular, determined the pump and probe pulse energies.²⁷ At the sample, the average pulse energies were 23 µJ for the pump pulses and 14 µJ for the probe pulses. The pump and probe pulse energies varied following the self-amplified spontaneous emission process. Events were sorted according to the pump pulse energy into four bins each with a width of ±3 µJ. At BL3 experimental hutch 5,²⁸ the Kirkpatrick-Baez mirrors focused the x-rays to spot size of 140 × 160 nm (full-width at half-maximum). To vary the x-ray dose, the sample was translated along the x-ray propagation direction, which resulted in increased x-ray focal areas up to 1.06 × 1.05 µm². The average profile of both pump and probe beams was measured at a set of sample translations and shown to be two-dimensional Gaussians. The pump pulse was selected by kicking the electron beam at the magnetic chicane preventing emission of the probe pulse from the downstream undulators. The probe pulse was isolated with a 1.6-mm-thick silicon filter, attenuating the pump and probe by 9 × 10⁻¹⁵ and 5 × 10⁻⁵, respectively.²⁹ The spatial separation between the pump and probe beams was observed to be ~20 nm. By a continuous motion of the sample along a serpentine path, an undamaged location was provided for each x-ray pump-probe exposure of the sample.
<figure><img src="image_10.png" /><figcaption>FIG. 1. (a) The experimental setup at the SACLA experimental hutch 5, (b) a selected diffraction image, and (c) azimuthally integrated diffraction patterns at 0 and 250 fs time delay.</figcaption></figure>

There were between 6000 and 25 000 exposures for each x-ray spot size and time delay, where the largest number of shots was performed with the smallest spot size and the smallest number with the largest spot size.

Figure 1 shows (a) the experimental setup, (b) a probe diffraction pattern on the Multi-Port Charge-Coupled Device (MPCCD) detector,³⁰ and (c) azimuthally averaged diffraction patterns at 0 and 250 fs time delay and 3.1 eV/atom x-ray dose. The left side of the detector was covered by a 50 µm thick Cu filter to select the pump diffraction and the right side by a 600 µm thick Al filter to select the probe diffraction. The samples were 20 µm thick nanocrystalline free-standing diamond (Diamond Materials). The thickness is more than an order of magnitude less than the x-ray attenuation length at the pump photon energy of 7 keV (428 µm),²⁹ implying that the x-rays are uniformly absorbed through the depth of the sample. From the Scherrer equation, the average nanocrystal size is estimated as 20 nm. From scanning electron microscope measurements, diamond films fabricated by the same method were analyzed to have an average grain size of 30 nm.³¹ The samples were inside a small helium enclosure. The nanocrystalline diamond produced a powder diffraction pattern as seen in Fig. 1(b). The octal MPCCD collected the (111), (220), and (311) diffraction peaks. From a comparison of the diffraction intensity distribution along the azimuthal angle, $\chi$, with the polarization factor, it is concluded that the nanocrystal orientation is nearly isotropic. It is not

expected that the dimensions of the nanocrystals would contribute dynamics on the 250 fs time range of this experiment; in nanocrystal-line gold films, the propagation of melting from grain boundaries has been observed but with time scales of ~50 ps.³²

### III. RESULTS

Three diffraction peaks are observed, and the diamond reflections are (111), (220), and (311). To extract the peak intensities, the peaks were fit using a Lorentzian function, $\frac{I}{\pi w(1 + (\frac{\lambda - \lambda_c}{w})^2)} + B$, where $I$ is the intensity, $w$ is the half-width at half-maximum (HWHM), $\lambda_c$ is the scattering angle ($2\theta$), and $B$ is the background. The background is a linear function including a constant and slope. A first normalization was made to the probe pulse intensities observed by a spectrometer in an upstream hutch and then a second normalization to the diffraction intensity at the pump-probe time delay of 0 fs. Figure 2 shows the diffraction intensities between 0 and 250 fs at a low, medium, and high x-ray dose. The effective x-ray dose $D_{eff}$ is calculated as follows:

$$D_{eff} = \frac{E}{\sqrt{w_x^2 + \lambda_e^2} \sqrt{w_y^2 + \lambda_e^2 \lambda_x \rho_a}}, \quad (1)$$

where $E$ is the pump pulse energy, $w_x$ is the pump horizontal x-ray focus (width at the 1/e level), $w_y$ is the pump vertical x-ray focus

<figure><img src="image_11.png" /><figcaption>FIG. 2. The time-dependent intensities of the (111), (220), and (311) diffraction peaks at the effective x-ray doses of (a) 0.9, (b) 3.1, and (c) 6.1 eV/atom. The curves are fits to an exponential function.</figcaption></figure>
(width at the 1/e level), $\lambda_e$ is the electron cascade size at the pump photon energy of 7 keV (0.41 µm according to Refs. 20 and 33), $\lambda_x$ is the x-ray penetration depth at the pump photon energy, and $\rho_a$ is the atomic density. The pump pulse energy is evaluated from the pump pulse intensity observed by the spectrometer, which is then normalized by the signal of a beamline intensity monitor and corrected by the reflectivity of the Kirkpatrick-Baez focusing mirrors. The inclusion of $\lambda_e$ in Eq. (1) represents an approximation because $\lambda_e$ is a maximum distance at which the electrons lose their energy below a cutoff.²⁰ The electron cascade size increases the volume of energy deposition especially for the smaller x-ray spot sizes, i.e., higher effective doses. In addition, the pump and probe beams have Gaussian intensity profiles, which results in different regions of the sample being excited by a range of x-ray doses. Therefore, one cannot expect a one-to-one correspondence between the theoretical (uniform) dose and the experimental effective dose.

The curves in Fig. 2 represent fits with an exponential function. Figure 3 shows the time constants and the diffraction intensities at 250 fs from the exponential fits. In Fig. 3, selected error bars are shown. For each reflection and at each x-ray dose, the diffraction intensities decrease with time, while no new peaks were observed, as shown in Fig. 1(c). It is concluded that there is an ultrafast disordering of the crystal lattice, but the diamond symmetry is preserved. As expected, the diffraction intensity at the maximum time delay, 250 fs, decreases with increasing x-ray dose as the sample becomes more disordered. A general trend is observed where the time constant for the changing diffraction intensity becomes faster with the increasing x-ray dose. In addition, the (111) time constants are observed to be longer than those of the (220) and (311) reflections.

The decreasing diffraction intensities can be understood further from a Debye-Waller analysis. In this model, the diffraction intensity is reduced by $\exp(-q^2\langle u^2 \rangle)$, where q is the scattering vector, $q = 4\pi\sin\theta/\lambda$, $2\theta$ is the scattering angle, $\lambda$ is the x-ray wavelength, and $\sqrt{\langle u^2 \rangle}$ is the root mean square (rms) atomic displacement perpendicular to a particular lattice plane. At zero pump-probe time delay, $\sqrt{\langle u_0^2 \rangle} = 0.043$ Å corresponds to the rms displacement of carbon atoms in room temperature diamond.³⁴ The Debye-Waller theory is valid for small atomic displacements.³⁵ In the case of the ultrafast melting of silicon, it has been concluded that the Debye-Waller theory was correct for displacements $\sqrt{\langle u^2 \rangle}$ up to ~1 Å,³⁶ compensating for the

relative dimensions of the diamond and silicon unit cells and for the projection of $u$ onto a direction hkl gives $\sqrt{\langle u_{hkl}^2 \rangle} \sim 0.4$ Å as an upper limit. An additional confirmation of the Debye-Waller analysis can be made by calculating the rms displacements directly from the XTANT+ simulation's atomic positions and by evaluating the rms displacements using the Debye-Waller model with the simulation's diffraction intensities.³⁷ For all three x-ray doses 0.75, 2.50, and 3.88 eV/atom, the atomic displacements from the atomic positions and from the Debye-Waller model agree within 10%.

Figure 4 shows the experimental and simulated rms atomic displacements calculated for the three diffraction peaks. At early time delays, the experimental atomic displacements derived from the three reflections are the same within the uncertainties, which indicates that the motion is isotropic. At all the x-ray fluences and later time delays, the experimental atomic displacements perpendicular to the (111) planes are substantially higher than that perpendicular to the (220) and (311) planes. It is noted that the (111) direction is along the carbon-carbon bond, which is weakened by the excitation of electrons into the conduction band.³⁸ The atomic displacements from the simulations, Figs. 4(a)-4(c), are discussed in Sec. IV.

At SACL A, x-ray pulses are delivered to beamline 3 at a repetition rate of 30 Hz. By accepting two pump pulses without probe pulses at individual sample locations, one measures the x-ray diffraction at a long time delay, 33 ms. Figure 5(a) shows the azimuthally integrated diffraction pattern and Fig. 5(b) the intensity of the diamond (111) and graphite (002) diffraction peaks with the maximum intensities normalized to 1. The sample was translated after each double x-ray pump exposure, with a sequential motion. There were between 401 and 546 exposures for each x-ray spot size. The graphite (002) peak appears at a $2\theta$ of 29.1°, which is somewhat shifted from the nominal angle of 30.7°. The peak is also broad with a width of ~2.3° (fwhm). Interestingly, shock compressed graphite-diamond mixtures display a broad graphite (002) peak at a $2\theta$ angle equivalent to 29.0° at the pump photon energy, 7 keV.³⁹ These observations suggest that the graphite crystals are affected by stress.

The observed threshold for graphitization was found to be at $D_{eff}$ of 1.2 ± 0.1 eV/atom. The threshold was evaluated by fitting the graphite (002) and diamond (111) intensity data to constants and lines and then finding the intersection points between them, as shown in Fig. 5. The error bar is derived from the standard deviation of the fits. This threshold is similar to the non-thermal ones obtained from the

<figure><img src="image_12.png" /><figcaption>FIG. 3. The (a) time constants and (b) intensities at 250 fs from the exponential fits of the (111), (220), and (311) diffraction intensity curves.</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>FIG. 4. The time-dependent root mean square displacements calculated from the experimental (111), (220), and (311) diffraction peak intensities at the effective x-ray doses of (a) 0.9, (b) 3.1, and (c) 6.1 eV/atom and the root mean square displacements from the simulations at the x-ray doses of (a) 0.75, (b) 2.50, and (c) 3.88 eV/atom.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>FIG. 5. (a) The azimuthally integrated diffraction pattern and (b) the intensities of the diamond (111) and graphite (002) diffraction peaks measured at a delay of 33 ms.</figcaption></figure>

tight-binding molecular dynamics (~1.0 eV/atom with XTANT) and the density-functional tight-binding molecular dynamics (1.3 eV/atom with XTANT+) calculations.⁸ It exceeds the graphitization damage threshold observed experimentally by Gaudin *et al.*⁶ (~0.7 eV/atom). The threshold from Gaudin *et al.* has the advantage that at XUV photon energies, the electron cascade size in diamond is in the nm range²⁰ and does not have a significant influence on the x-ray dose. The long timescale required to observe the graphite diffraction peak implies that graphitization may occur by a thermal process, rather than non-thermal one.⁶,⁷

At a high x-ray dose, about 4 eV/atom and above, the intensity of the graphite (002) and diamond (111) peaks both decrease, indicating a regime of ablation or permanent structural disordering, similar to that predicted theoretically.⁸ In a limited postmortem analysis of the samples, scanning probe microscopy showed indentations 1 µm wide and 50–100 nm deep, indicating that some ablation does occur.

## IV. SIMULATIONS AND DISCUSSION

Figure 4 displays the root mean square atomic displacements from XTANT+ simulations and Fig. 6 the corresponding atomic

snapshots. The XTANT+ calculations were performed with 512 atoms in a simulation box with periodic boundaries, at 0.05 fs molecular dynamics time step, and with 100 000 Monte Carlo iterations. Five molecular dynamics realizations were performed at each x-ray dose. The application of averaged electron distributions removes strong fluctuations from the molecular dynamics. Figure 4 shows the rms atomic displacements from the simulations at three x-ray doses: 0.75, 2.5, and 3.88 eV/atom with homogeneous x-ray absorption. The simulations were performed at x-ray doses similar to the effective x-ray doses used in the experiment. In the calculations, all the atoms experience a constant x-ray dose. In contrast, in the experiment, because of the Gaussian profiles of the pump and probe beams, the sample is excited by a range of x-ray doses. The difficulty of comparing simulations with x-ray pump-probe experiments has been discussed in depth in Ref. 40. In Fig. 4, a measure of simulation error is given as bands from the minimum to maximum rms atomic displacement chosen from all the MD realizations performed. They are plotted together with curves showing the rms displacements averaged over the number of MD realizations.

In the XTANT+ atomic displacements shown in Fig. 4, dynamics are observed on several time scales. The atomic displacements remain
<figure><img src="image_15.png" /><figcaption>FIG. 6. XTANT+ atomic snapshots at x-ray doses of (a) 0.75, (b) 2.5, and (c) 3.88 eV/atom and time delay of 250 fs.</figcaption></figure>

nearly constant until ~10 fs after the pump pulse similar to what was previously observed by Inoue et al.⁵,¹¹ The delayed onset for atomic disordering is not observable in the experiment because of the lack of sufficient data points at short time delays. Then, the atomic displacements increase sharply until ~40 fs. Both the delayed onset and the fast increase in the atomic displacements occur during the timescale of the electron cascading, during which the collisions of photoelectrons and Auger electrons with atoms result in excitations of the valence electrons to the conduction band. Then, a slower increase in the atomic displacements occurs. This behavior at later times takes place after the electron cascading is finished, and the atoms move in the new potential energy surface. In Fig. 4, good agreement is found between the experimental and simulation results for the atomic displacements in the (220) and (311) directions. The experimentally observed difference between the atomic displacements in the (111) direction compared with those in the (220) or (311) directions is not found in the XTANT+ simulation results. While performing a volume integration of simulations at different x-ray doses according to the x-ray pump and probe profiles would improve the quantitative agreement between the experiment and simulations, it was outside the scope of the present study because of the required computer resources. Overall, the simulations in Fig. 4 provide qualitative agreement with the experimental results in terms of time scales and trends.

Figure 6 shows XTANT+ atomic snapshots from individual molecular dynamics realizations. This figure has been generated using VESTA 3.⁴¹ The disorder is seen to increase with the x-ray dose. At 2.5 eV/atom, evidence of graphite planes is observed in the upper left of the simulation box, although in the other molecular dynamics realizations, the graphite planes were less clearly seen. For the short-range order, carbon pair distribution functions predicted by XTANT+⁸ show a well-defined peak for the second nearest-neighbor at an x-ray dose of 1.0 eV/atom. At 3.0 eV/atom, there is only a well-defined peak for the first nearest-neighbor.

The experiment does not observe extra diffraction peaks from graphite on the femtosecond timescale although the graphite (002) reflection is detected at the pump-probe time delay of 33 ms. From the XTANT+ atomic structures using the XSINC code,²⁷ one may derive diffraction intensities for particular reflections including graphite (002). The graphite peak is orders of magnitude lower than the diamond peaks. In addition, in Fig. 6, the graphite structure is seen over a fraction of the simulated atoms indicating small graphite domains, which would result in broad diffraction peaks. These weak broad diffraction peaks would not be observable in the experiment. The experiment provides evidence for

disordering on an ultrafast timescale and graphitization on a long timescale, which may be of thermal origin. The simulations show partial graphitization but predominantly disordering on the ultrafast timescale.

## ACKNOWLEDGMENTS

The experiments were performed with the approval of the Japan Synchrotron Radiation Research Institute (JASRI, Proposal Nos. 2019B8070 and 2021A8039). P.H. and H.J.L. were supported by the DOE Office of Science, Fusion Energy Science, under FWP No. 100106. N.J.H. was funded by the DOE Office of Science, Fusion Energy Science, under FWP No. 100182. R.W.F. acknowledges support from the U.S. Department of Energy, Office of Science, Office of Fusion Energy Sciences, Award No. DE-AC02-05CH11231, and National Nuclear Security Administration, Award No. DE-NA0003842. I.I. gratefully acknowledges the financial support from the Japan Society for the Promotion of Science KAKENHI (Grant Nos. 19K20604 and 22H03877). K.J.K. acknowledges the Polish National Agency for Academic Exchange for funding in the frame of the Bekker programme (No. PPN/BEK/2020/1/00184). N.M. gratefully acknowledges the financial support from the Czech Ministry of Education, Youth and Sports (Grant Nos. LTT17015, LM2018114, and EF16_013/0001552).

## AUTHOR DECLARATIONS

### Conflict of Interest

The authors have no conflicts to disclose.

### Author Contributions

**Philip Heimann:** Writing – original draft (equal). **Yuichi Inubushi:** Investigation (equal). **Konrad J. Kapcia:** Formal analysis (equal). **Hae Ja Lee:** Investigation (equal). **Vladimir Lipp:** Formal analysis (equal). **Paloma Martinez:** Investigation (equal). **Nikita Medvedev:** Formal analysis (equal). **Franz Tavella:** Investigation (equal). **Sven Toleikis:** Investigation (equal). **Makina Yabashi:** Investigation (equal). **Toshinori Yabuchi:** Investigation (equal). **Nicholas Hartley:** Formal analysis (equal). **Jumpei Yamada:** Investigation (equal). **Beata Ziaja:** Formal analysis (equal). **Ichiro Inoue:** Methodology (equal). **Victor Tkachenko:** Formal analysis (equal). **Andre Antoine:** Investigation (equal). **Fabien Dorchies:** Investigation (equal). **Roger Falcone:** Investigation (equal). **Jérôme Gaudin:** Investigation (equal). **Hauke Höppner:** Investigation (equal).
DATA AVAILABILITY

The data that support the findings of this study are available from
the corresponding authors upon reasonable request.

REFERENCES

<footnote n="1">¹M. Z. Mo, Z. Chen, R. K. Li, M. Dunning, B. B. L. Witte *et al.*, "Heterogeneous to homogeneous melting transition visualized with ultrafast electron diffraction," *Science* **360**, 1451 (2018).</footnote>

<footnote n="2">²A. M. Lindenberg, J. Larsson, K. Sokolowski-Tinten, K. J. Gaffney, C. Blome *et al.*, "Atomic-scale visualization of inertial dynamics," *Science* **308**, 392 (2005).</footnote>

<footnote n="3">³E. S. Zijlstra, A. Kalitsov, T. Zier, and M. E. Garcia, "Fractional diffusion in sili-con," *Adv. Mater.* **25**, 5605 (2013).</footnote>

<footnote n="4">⁴N. Medvedev, H. O. Jeschke, and B. Ziaja, "Nonthermal phase transitions in semiconductors induced by a femtosecond extreme ultraviolet laser pulse," *New J. Phys.* **15**, 015016 (2013).</footnote>

<footnote n="5">⁵I. Inoue, Y. Inubushi, T. Sato, K. Tono, T. Katayama *et al.*, "Observation of femtosecond x-ray interactions with matter using an x-ray-x-ray pump-probe scheme," *Proc. Natl. Acad. Sci. U. S. A.* **113**, 1492 (2016).</footnote>

<footnote n="6">⁶J. Gaudin, N. Medvedev, J. Chalupsky, T. Burian, S. Dastjani-Farahani *et al.*, "Photon energy dependence of graphitization threshold for diamond irradiated with an intense XUV FEL pulse," *Phys. Rev. B* **88**, 060101 (2013).</footnote>

<footnote n="7">⁷F. Tavella, H. Höppner, V. Tkachenko, N. Medvedev, F. Capotondi *et al.*, "Soft x-ray induced femtosecond solid-to-solid phase transition," *High Energy Density Phys.* **24**, 22 (2017).</footnote>

<footnote n="8">⁸V. Lipp, V. Tkachenko, M. Stransky, B. Aradi, T. Frauenheim *et al.*, "Density functional tight binding approach utilized to study x-ray-induced transitions in solid materials," *Sci. Rep.* **12**, 1551 (2022).</footnote>

<footnote n="9">⁹D. H. Reitze, H. Ahn, and M. C. Downer, "Optical-properties of liquid carbon measured by femtosecond spectroscopy," *Phys. Rev. B* **45**, 2677 (1992).</footnote>

<footnote n="10">¹⁰B. C. Stuart, M. D. Feit, A. M. Rubenchik, B. W. Shore, and M. D. Perry, "Laser-induced damage in dielectrics with nanosecond to subpicosecond pulses," *Phys. Rev. Lett.* **74**, 2248 (1995).</footnote>

<footnote n="11">¹¹I. Inoue, V. Tkachenko, K. J. Kapcia, V. Lipp, B. Ziaja *et al.*, "Delayed onset and directionality of x-ray-induced atomic displacements observed on subatomic length scales," *Phys. Rev. Lett.* **128**, 223203 (2022).</footnote>

<footnote n="12">¹²A. L. Kritcher, C. V. Young, H. F. Robey, C. R. Weber, A. B. Zylstra *et al.*, "Design of inertial fusion implosions reaching the burning plasma regime," *Nat. Phys.* **18**, 251 (2022).</footnote>

<footnote n="13">¹³F. Uhlen, D. Nilsson, A. Holmberg, H. M. Hertz, C. G. Schroer *et al.*, "Damage investigation on tungsten and diamond diffractive optics at a hard x-ray free-electron laser," *Opt. Express* **21**, 8051 (2013).</footnote>

<footnote n="14">¹⁴C. David, S. Gorelick, S. Rutishauser, J. Krzywinski, J. Vila-Comamala *et al.*, "Nanofocusing of hard x-ray free electron laser pulses using diamond based Fresnel zone plates," *Sci. Rep.* **1**, 57 (2011).</footnote>

<footnote n="15">¹⁵I. Inoue, T. Osaka, T. Hara, T. Tanaka, T. Inagaki *et al.*, "Generation of narrow-band x-ray free-electron laser via reflection self-seeding," *Nat. Photonics* **13**, 319 (2019).</footnote>

<footnote n="16">¹⁶T. H. Dinh, N. Medvedev, M. Ishino, T. Kitamura, N. Hasegawa *et al.*, "Controlled strong excitation of silicon as a step towards processing materials at sub-nanometer precision," *Commun. Phys.* **2**, 150 (2019).</footnote>

<footnote n="17">¹⁷N. Medvedev, U. Zastrau, E. Forster, D. O. Gericke, and B. Rethfeld, "Short-time electron dynamics in aluminum excited by femtosecond extreme ultraviolet radiation," *Phys. Rev. Lett.* **107**, 165003 (2011).</footnote>

<footnote n="18">¹⁸P. Lorazo, L. J. Lewis, and M. Meunier, "Thermodynamic pathways to melting, ablation, and solidification in absorbing solids under pulsed laser irradiation," *Phys. Rev. B* **73**, 134108 (2006).</footnote>

<footnote n="19">¹⁹B. Ziaja, R. A. London, and J. Hajdu, "Unified model of secondary electron cascades in diamond," *J. Appl. Phys.* **97**, 064905 (2005).</footnote>

<footnote n="20">²⁰V. Lipp, I. Milov, and N. Medvedev, "Quantifying electron cascade size in various irradiated materials for free-electron laser applications," *J. Synchrotron Radiat.* **29**, 323 (2022).</footnote>

<footnote n="21">²¹H. O. Jeschke, M. E. Garcia, and K. H. Bennemann, "Microscopic analysis of the laser-induced femtosecond graphitization of diamond," *Phys. Rev. B* **60**, R3701 (1999).</footnote>

<footnote n="22">²²H. O. Jeschke, M. E. Garcia, and K. H. Bennemann, "Theory for the ultrafast ablation of graphite films," *Phys. Rev. Lett.* **87**, 015003 (2001).</footnote>

<footnote n="23">²³T. Ishikawa, H. Aoyagi, T. Asaka, Y. Asano, N. Azumi *et al.*, "A compact x-ray free-electron laser emitting in the sub-angstrom region," *Nat. Photonics* **6**, 540 (2012).</footnote>

<footnote n="24">²⁴T. Hara, Y. Inubushi, T. Katayama, T. Sato, H. Tanaka *et al.*, "Two-colour hard x-ray free-electron laser with wide tunability," *Nat. Commun.* **4**, 2919 (2013).</footnote>

<footnote n="25">²⁵I. Inoue, T. Hara, Y. Inubushi, K. Tono, T. Inagaki *et al.*, "X-ray Hanbury Brown-Twiss interferometry for determination of ultrashort electron-bunch duration," *Phys. Rev. Accel. Beams* **21**, 080704 (2018).</footnote>

<footnote n="26">²⁶I. Inoue, Y. Inubushi, T. Osaka, J. Yamada, K. Tamasaku *et al.*, "Shortening x-ray pulse duration via saturable absorption," *Phys. Rev. Lett.* **127**, 163903 (2021).</footnote>

<footnote n="27">²⁷K. Tamasaku, Y. Inubushi, I. Inoue, K. Tono, M. Yabashi *et al.*, "Inline spectrometer for shot-by-shot determination of pulse energies of a two-color x-ray free-electron laser," *J. Synchrotron Radiat.* **23**, 331 (2016).</footnote>

<footnote n="28">²⁸H. Yumoto, Y. Inubushi, T. Osaka, I. Inoue, T. Koyama *et al.*, "Nanofocusing optics for an x-ray free-electron laser generating an extreme intensity of 100 EW/cm² using total reflection mirrors," *Appl. Sci.* **10**, 2611 (2020).</footnote>

<footnote n="29">²⁹B. L. Henke, E. M. Gullikson, and J. C. Davis, "X-ray interactions: Photoabsorption, scattering, transmission, and reflection at $E = 50–30 000$ ev, $Z = 1–92$," *At. Data Nucl. Data Tables* **54**, 181 (1993).</footnote>

<footnote n="30">³⁰T. Kameshima, S. Ono, T. Kudo, K. Ozaki, Y. Kirihara *et al.*, "Development of an x-ray pixel detector with multi-port charge-coupled device for x-ray free-electron laser experiments," *Rev. Sci. Instrum.* **85**, 033110 (2014).</footnote>

<footnote n="31">³¹C. Dawedeit, S. O. Kucheyev, S. J. Shin, T. M. Willey, M. Bagge-Hansen *et al.*, "Grain size dependent physical and chemical properties of thick CVD diamond films for high energy density physics experiments," *Diamond Relat. Mater.* **40**, 75 (2013).</footnote>

<footnote n="32">³²T. A. Assefa, Y. Cao, S. Banerjee, S. Kim, D. Kim, H. Lee *et al.*, "Ultrafast x-ray diffraction study of melt-front dynamics in polycrystalline thin films," *Sci. Adv.* **6**, 2445 (2020).</footnote>

<footnote n="33">³³N. Medvedev, F. Akhmetov, R. A. Rymzhanov, R. Voronkov, and A. E. Volkov, "Modeling time-resolved kinetics solids induced by extreme electronic excitation," *Adv. Theory Simul.* **5**, 2200091 (2022).</footnote>

<footnote n="34">³⁴V. F. Sears and S. A. Shelley, "Debye-Waller factor for elemental crystals," *Acta Crystallogr., Sect. A* **47**, 441 (1991).</footnote>

<footnote n="35">³⁵B. E. Warren, *X-Ray Diffraction* (Dover, 1990).</footnote>

<footnote n="36">³⁶T. Zier, E. S. Zijlstra, A. Kalitsov, I. Theodonis, and M. E. Garcia, "Signatures of nonthermal melting," *Struct. Dyn.* **2**, 054101 (2015).</footnote>

<footnote n="37">³⁷M. M. Abdullah, Z. Jurek, S. K. Son, and R. Santra, "Calculation of x-ray scattering patterns from nanocrystals at high x-ray intensity," *Struct. Dyn.* **3**, 054101 (2016).</footnote>

<footnote n="38">³⁸I. Inoue, Y. Deguchi, B. Ziaja, T. Osaka, M. M. Abdullah *et al.*, "Atomic-scale visualization of ultrafast bond breaking in x-ray-excited diamond," *Phys. Rev. Lett.* **126**, 117403 (2021).</footnote>

<footnote n="39">³⁹D. K. Potter and T. J. Ahrens, "Shock consolidation of diamond and graphite mixtures to fused polycrystalline diamond," *J. Appl. Phys.* **63**, 910 (1988).</footnote>

<footnote n="40">⁴⁰V. Tkachenko, M. M. Abdullah, Z. Jurek, N. Medvedev, V. Lipp *et al.*, "Limitations of structural insight into ultrafast melting of solid materials with x-ray diffraction imaging," *Appl. Sci.* **11**, 5157 (2021).</footnote>

<footnote n="41">⁴¹K. Momma and F. Izumi, "VESTA 3 three-dimensional visualization crystal, volumetric morphology data," *J. Appl. Crystallogr.* **44**, 1272 (2011).</footnote>