# Probing local and electronic structure in Warm Dense Matter: single pulse synchrotron x-ray absorption spectroscopy on shocked Fe

Raffaella Torchio, Florent Occelli, Olivier Mathon, Arnaud Sollier, Emilien

Lescoute, Laurent Videau, Tommaso Vinci, Alessandra Benuzzi-Mounaix, Jon

Headspith, William Helsby, et al.

► To cite this version:

Raffaella Torchio, Florent Occelli, Olivier Mathon, Arnaud Sollier, Emilien Lescoute, et al. Probing local and electronic structure in Warm Dense Matter: single pulse synchrotron x-ray absorption spectroscopy on shocked Fe. Scientific Reports, 2016, 6, pp.26402. 10.1038/srep26402. hal-01339817

HAL Id: hal-01339817

https://hal.sorbonne-universite.fr/hal-01339817

Submitted on 30 Jun 2016

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
SCIENTIFIC REPORTS

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

OPEN

# Probing local and electronic structure in Warm Dense Matter: single pulse synchrotron x-ray absorption spectroscopy on shocked Fe

Received: 05 October 2015

Accepted: 28 April 2016

Published: 01 June 2016

Raffaella Torchio<fnref n="1" />,², Florent Occelli<fnref n="1" />, Olivier Mathon², Arnaud Sollier<fnref n="1" />, Emilien Lescoute<fnref n="1" />, Laurent Videau<fnref n="1" />, Tommaso Vinci³,⁴, Alessandra Benuzzi-Mounaix³,⁴, Jon Headspith⁵, William Helsby⁵, Simon Bland⁶, Daniel Eakins⁵,⁶, David Chapman⁶, Sakura Pascarelli² & Paul Loubeyre<fnref n="1" />

Understanding Warm Dense Matter (WDM), the state of planetary interiors, is a new frontier in scientific research. There exists very little experimental data probing WDM states at the atomic level to test current models and those performed up to now are limited in quality. Here, we report a proof-of-principle experiment that makes microscopic investigations of materials under dynamic compression easily accessible to users and with data quality close to that achievable at ambient. Using a single 100 ps synchrotron x-ray pulse, we have measured, by K-edge absorption spectroscopy, ns-lived equilibrium states of WDM Fe. Structural and electronic changes in Fe are clearly observed for the first time at such extreme conditions. The amplitude of the EXAFS oscillations persists up to 500 GPa and 17000 K, suggesting an enduring local order. Moreover, a discrepancy exists with respect to theoretical calculations in the value of the energy shift of the absorption onset and so this comparison should help to refine the approximations used in models.

Warm Dense Matter refers to a material state that occurs, under extreme thermodynamic conditions, at the intersection of condensed matter physics and plasma physics<fnref n="1" />. The limits of WDM are not sharply defined: it is the region where electron degeneracy, ion coupling and atomic physics all play a role in the material response, i.e. typically density and temperature are in units of solid density and 10⁴K (~eV) respectively. The most advanced calculations of WDM are made using first principle density functional theory²,³. It is a ground state electronic theory that has been very successfully extended to predict the properties of compressed matter at ambient temperature but the approximations for its extension in the eV temperature range, where the contributions of excited electronic states are important, still need to be validated.

Unfortunately, very little experimental data exist to test the predictions of these calculations at the atomic level. The equation of state of WDM is currently obtained at dynamic compression facilities (such as multi-stage gas guns, multi-mega ampere pulsed power generators and high energy lasers) by measuring the shock velocities. The link between the measured temperature–pressure–density relation and the microscopic state and atomic structure is the next step to progress in the validation of WDM models. Over the last decade, much effort has been devoted to the development of in-situ x-ray probes for atomic scale characterization of compressed matter in front of these various dynamic facilities. Most common techniques use x-ray back-lighters generated by lasers. Recently, significant results were achieved using x-ray absorption, on ramp compressed iron to around 500 GPa⁴ in the EXAFS (Extended X-ray Absorption Fine Structure) range, on aluminum shocked to 8 eV⁵ and SiO₂ to

<footnote n="1">¹CEA, DAM, DIF, F-91297, Arpajon, France. ²ESRF-The European Synchrotron 71, Avenue des Martyrs Grenoble, France. ³LULI - CNRS, Ecole Polytechnique, CEA: Université Paris-Saclay, France. ⁴UPMC Univ Paris 06: Sorbonne Universités - F-91128 Palaiseau cedex, France. ⁵Science and Technology Facilities Council, Daresbury Laboratory Sci-Tech Daresbury, Warrington WA4 4AD, UK. ⁶Institute of Shock Physics, Imperial College London, London SW7 2AZ, UK. Correspondence and requests for materials should be addressed to R.T. (email: torchio@esrf.fr)</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Schematic view of the experimental set-up. Panel (a) a curved Si-crystal focuses the polychromatic x-ray beam onto the sample. The beam transmitted by the sample diverges towards a position sensitive detector where energy is correlated to position; in this way, the full EXAFS spectrum is acquired at once by a position sensitive detector with Ge pixels (1024). A long laser pulse (10 ns, up to 35 J, at 1057 nm) focused at the sample position, drives the compression wave in Fe. The Fe target consists of a 3.5 µm iron layer sandwiched between two diamond windows. Panel (b) a series of typical single bunch XANES spectra obtained by changing the x-ray probe delay time with respect to the laser onset time, while keeping the same driving energy for each shot, shows that the compressed state is in the hcp phase and its thermodynamic conditions are stable over 2 ns at least. Panel (c) 1-D hydro-simulation of the shock in the Fe layer with a similar color scale for the pressure as for the XANES spectra.</figcaption></figure>

5 eV⁶ in the XANES (X-ray Absorption Near Edge Spectroscopy) range. However, these x-ray characterizations have not reached the energy resolution and level of detail achievable on synchrotron beamlines, even for samples under high static pressure-temperature conditions⁷.

Here, we propose an alternative line of progress by combining the best possible x-ray characterization methods available at synchrotrons to dynamic compression⁸,⁹. This is based on a similar strategy as the one used by the static high pressure community, namely minimizing the volume of the compressed sample to achieve extreme pressures and take full advantage of the high brightness (10¹⁴ photons/s/10 µm²), extreme stability (fluctuations of normalized intensity ≤0.002 at ~1 KHz frequency), high energy resolution (δE/E ~ 10⁻⁴), large spectral range (ΔE/E > 10%) and small spot size (typically less than 10 µm²) of the synchrotron x-ray beam. The shocked states were generated by a portable 35 J laser. The laser pulse was focused to a diameter of ~90 µm (high power density configuration) or ~350 µm (low power density configuration), resulting in similar compressions as those obtained by kJ lasers on mm size samples. Measurements were performed at the energy dispersive x-ray absorption beamline ID24¹⁰ of the European Synchrotron Radiation Facility, which has been recently upgraded. X-ray absorption fine structure, which includes XANES and EXAFS, is an element selective probe of the local electronic and atomic structure. Hence, it is well adapted for the characterization of WDM at the atomic level. Iron is a key constituent of terrestrial and exoplanets. Given its ubiquitous importance, Fe is often one of the first materials to be investigated with novel experimental techniques. With this in mind, the present data are thus compared to XAS measurements recently performed under dynamic compression at the OMEGA laser⁴ and LCLS free electron laser facilities¹¹.

## Results

The sketch of the experimental configuration is presented in Fig. 1 panel a. In the dispersive geometry, a fan of x-rays is dispersed and focused on the sample by a curved crystal, so that the energy of each ray varies as a function of convergence angle. The transmitted photons are then detected by a position sensitive detector allowing the simultaneous collection of the whole spectrum, up to a 300 eV range above the Fe K-edge (7.112 keV). Data were acquired using a single x-ray pulse of 100 ps duration, during a 4-bunch mode configuration of the synchrotron storage ring that assures the highest number of photons per pulse and a maximum separation between pulses (700 ns). The x-ray spot could be focused down to 5 µm × 7 µm FWHM (H × V) to ensure that the probed part of the compressed sample was not perturbed by edge effects. The position sensitive detector equipped with 1024 Ge pixels was operated with an integration time of 680 ns and phase locked to the ESRF radiofrequency system, ensuring integration of the signal generated by a unique bunch of electrons in the storage ring. The absorption spectrum so acquired for the undriven sample (in Fig. 1, panel b, the *ref* spectrum) shows fine structure features
<figure><img src="image_4.png" /><figcaption>**Figure 2.** EXAFS spectra, volume and temperature determination. Panel (a) Series of EXAFS spectra obtained by increasing the irradiation intensity at 7 ns of delay, driving the compression in the 350 µm phase plate configuration. The spectra under shock are compared to the reference spectra at ambient. The structural change between the bcc and hcp phase is clearly seen. Panel (b) Fit of the EXAFS for the reference spectrum and the most compressed hcp one. The compression factor and the temperature extracted for each intensity are indicated in panel (c) with comparison to the simulations outputs, literature Hugoniot curves<sup>27,28</sup> and shock data<sup>25,26</sup>. The pressure value for the EXAFS points is given by the average between the two hydrodynamic codes outputs.</figcaption></figure>

comparable to those of a spectrum acquired on an energy scanning XAS beamline. The <35 J energy of a 1057 nm laser delivered in a 10 ns square pulse was focused under vacuum on the target with an angle of 30°, launching an ablation driven compression in the target perpendicular to its surface. The x-ray beam probed the shocked Fe parallel to the compression axis with variable drive-probe delay.

The minimization of the compressed volume of Fe required a careful target design that was optimized using two 1D hydro-codes, ESTHER¹² and MULTI¹³. The target was made of a thin (3.5 µm) deposit of pure Fe sandwiched between two diamond windows, which act to confine the shocked state over a time interval much longer than the 100 ps x-ray pulse. A CH layer was deposited on the ablation surface to limit pre-heating effects. Two different kinds of targets were prepared for the two laser configurations used to cover a large thermodynamic range. In the first case, diamond windows of 40 (front) and 50 (rear) µm were used and the laser was focused with a phase plate of ~350 µm producing a maximum usable ablation laser intensity of 4 × 10¹² W/cm² and so a maximum pressure of ~80 GPa in Fe. The x-ray focal spot was adjusted to 5 × 100 µm (H × V) FWHM in this configuration. The shock launched in the front window generates reverberating compression waves in the iron foil due to the presence of the rear diamond window. The pressure reached in Fe is approximately the pressure of the diamond front window and because of the impedance matching between diamond and iron, the final state of iron is near its principal Hugoniot. A 2-D simulation with the DUED code¹⁴ was performed to evaluate rarefaction edge effects and the homogeneous compressed region in Fe was calculated to be 200 µm in diameter. The confinement time of the stable thermodynamic conditions in Fe was estimated by hydrodynamic simulations and measured directly from the modifications in the XANES spectra for various delays between the laser arrival on the target and the x-ray pulse. As seen in Fig. 1 panel b, a 4 ns confinement time is observed, in agreement with the hydro-simulations (panel c); within this interval, thermodynamic conditions are stable for at least 2 ns. The pressure was estimated by hydro-simulations using MULTI¹³ and ESTHER¹² codes and also double-checked using an empirical law giving the ablation pressure versus the laser intensity in the ablator¹⁵ (see the Supplementary material); the error in the pressure is estimated as the difference between the two hydrodynamic codes outputs.

The reproducibility of the measured spectra was good. We are thus confident that homogeneous well defined compressed states of Fe were measured.

The evolution of the K-edge EXAFS spectra at different pressures is shown in Fig. 2 (panel a); no deterioration of the spectra quality between the undriven and the shocked states is observed. The ambient spectrum was modeled with the bcc structure and those at higher pressure, from 40 GPa, with the hcp structure, as suggested by the changes in the shape of the spectral features. For example, the disappearance of the peak around 7.2 keV is the signature of the bcc-hcp transition¹⁶. A quantitative analysis of the EXAFS spectra enabled the determination of the volume and temperature of the compressed solid states⁴,¹⁷ (Fig. 2 panel c). The fits of normalized EXAFS signals, χ(k), of the unshocked and highest pressure shocked samples are shown in panel b of Fig. 2. The fit allows the lattice parameters $a_{bcc}$ and $a_{hcp}$ to be determined, and thus the volume ($V_{bcc} = a^3/2$ and $V_{hcp} = 0.433 a^2 c$). In the hcp phase, c/a was fixed to the value 1.61 since the k-range of the EXAFS was not sufficient for its precise adjustment; however, simulations and static measurements¹⁸⁻²² have shown that the c/a value remains almost
<figure><img src="image_5.png" /><figcaption>**Figure 3.** EXAFS spectra and structural changes. (Left) A series of EXAFS spectra obtained in the highly focused configuration of the laser. (Right). Zoom over the edge region.</figcaption></figure>

ideal (1.59–1.63) over the thermodynamic region covered here. The temperatures were also estimated from the EXAFS spectra using a correlated Debye model²³, where the Debye temperature as a function of compression was taken from ref. 24 (see the Methods session for further details on the fitting procedures). In panel c we compare V and T values obtained from the EXAFS analysis to the hydrodynamic codes outputs, literature EoS²⁷,²⁸ and shock data²⁵,²⁶. The compression data obtained by the EXAFS analysis are found quite close to previous shock data²⁵,²⁶, EoS²⁷,²⁸ and simulations values. The EXAFS temperature values are in agreement with the outputs from the ESTHER code, systematically higher than the Hugoniots. The difference found in the PT and PV outputs from the two codes origins from the different parameters and models used in the calculations and is a good estimation of the error made by this kind of calculations (see the Supplementary material for further details).

In order to achieve pressures in the 500 GPa range, the 350 µm phase plate was removed and the laser directly focused to a 90 µm spot with a Gaussian spatial profile. To maintain a homogeneous thermodynamic state and limit 2D effects, Fe targets with 25 µm thick diamond windows were used in that case. 2D-simulations give a 40 µm diameter homogeneous region, larger than the x-ray probe beam. An additional vertical focusing mirror was used to reduce the x-ray focal spot down to 5 × 7 µm (H × V) FWHM. Again, a homogeneous thermodynamic state could be probed, however, the reproducibility from shot to shot was poorer due to the request in alignment of a 10 µm precision coupled to the lower spatial homogeneity of the laser profile. As seen in Fig. 3, the EXAFS data quality obtained in this second configuration was not as good as for the spectra acquired with the phase plate and with a larger x-ray spot (the signal to noise ratio decreased by a factor ~2). The high frequency noise is related to lower photon statistics, because the additional mirror could not be put under vacuum. However, as seen in Fig. 3, clear trends can be followed in the spectra as a function of the laser intensity (hence pressure), both in the XANES and in the first EXAFS oscillations region. The compression and heating trend are revealed from the shift to higher energy and the broadening of the EXAFS oscillations respectively (black, purple and red spectra in Fig. 3). As pressure increases this trend is broken (orange spectrum) and the bump in the edge (at E = 7.12) starts to flatten, indicating the onset of melting, in agreement with the recent melting curve from Anzellini et al.¹⁸. The most extreme shocked state has been achieved with I = 5 × 10¹³ W/cm², with a corresponding pressure estimated around 500 GPa. Remarkably, weak EXAFS oscillations remain as liquid Fe is further compressed and heated to such an extreme state, indicating a persisting local order in WDM.

Figure 4 shows the states in the phase diagram of Fe that have been probed by XAS measurements under dynamic compression, namely our data (red full circles) and measurements on two other platforms⁴,¹¹. Although the thermodynamic domains of investigation are comparable in these three experiments, the use of the synchrotron facility enables to improve the XAS data. The experiment¹¹ at the x-ray free electron laser LCLS could not bring much insight on the local order in the liquid phase. The energy range (25 eV) was too narrow to cover even the first EXAFS oscillations and the data quality suffered from the intrinsic intensity fluctuations of the incident beam, leading to the necessity to average over few tens of shots to acquire satisfactory spectra. The second experiment⁴ at the OMEGA facility produced very good quality EXAFS data of Fe compressed to 560 GPa but the XANES region could not be exploited due to a lack in energy resolution²⁹.

The present PT data points are close to the Fe principal Hugoniot curve up to around 270 GPa and 6700 K (corresponding to a laser power of 3 × 10¹³ W/cm²) where the hydrodynamic simulations start to witness some preheating which makes the two last points higher in temperature (see the Supplementary material for further details on the simulations). Our data agree with the melting curve of Anzellini et al.¹⁸: our point at 170 GPa and 3800 K (red spectrum in Fig. 3), which is close to the new melting curve of Aquilanti et al.⁷ still shows a pronounced bump at the absorption edge, a signature of the solid phase. However, a precise determination of the melting curve is beyond the scope of this paper.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Iron phase diagram. Phase diagram of Fe showing our data (red full circles) in comparison to recent data from the literature. Full green triangles are from Harmand *et al.*¹¹; blue squares are from Ping *et al.*⁴; full is single shock, half horizontal is multiple shock with P₀ = 100 GPa and half vertical is multiple shock with P₀ = 150 GPa; black lines are melting lines from ref. 18 (full) and ref. 7 (dash). Grey curves indicate different estimations of Fe Hugoniot, from Brown²⁵, SESAME²⁷ and BLF²⁸.</figcaption></figure>

The evolution of the absorption at the K-edge of iron is shown in Fig. 5. The quality of the data is sufficient to detect the evolution of features (noted *a*, *b*, *c*, *d*) as predicted by ab-initio molecular dynamics calculations³,¹¹ and as observed in the laser heated DAC². In the compressed hcp solid (blue spectrum in Fig. 5 left bottom panel) feature *a* increases in intensity and feature *c* moves to higher energy due to the compression. At higher drive intensity up to I = 5 × 10¹³ W/cm², corresponding to around 500 GPa and 17000 K, feature *c* moves back to lower energies and feature *a* increases more in intensity as the edge is modified into a rounded shape: such drastic change is the sign of a molten phase. In fact, similar behavior has been recently observed in lower pressure laser heating experiments (<100 GPa) probing the melting of Fe⁷,³⁰. By comparing the available ab-initio molecular dynamics simulations of iron in similar conditions (Fig. 5 right panels) to our data, we observe that although the shape of the absorption onset is more structured than the experimental one, theory is capable of grasping the main changes in the data (indicated with arrows). However, whereas good agreement between the theoretical prediction and the data is observed in the energy shift of the absorption onset (7.112 keV, feature *d*) for the compressed solid (Fig. 5, bottom panels), an important discrepancy is found for the data recorded at the highest P, T values (Fig. 5, top panels): the theoretical shift is ~−2 eV, whereas the experimental shift is ~−0.5 eV.

Discussion

In the low power density configuration, the volume and the temperature could be extracted from the fit of the EXAFS oscillations. In the high power density configuration, the temperature induced broadening of the EXAFS oscillations does not enable volume extraction and local order structural refinement. However, a clear signature of the hcp structure could be followed up to the melting of Fe. This observation confirms the great stability of the hcp structure in Fe under very high pressure, as shown previously utilizing static laser heated DAC¹⁸ or the dynamic multi shock compression⁴.

At the most extreme states achieved (500 GPa and 17000 K), the electrons degeneracy factor is around 2 and the plasma coupling parameter around 6, and so thermodynamic conditions for Fe enter into the WDM regime¹. With further optimization of the technique, future measurements are expected to yield data up to 500 GPa and above of similar quality as in Fig. 2, thus allowing a finer comparison to theoretical models and a possibility to disclose the local order in WDM Fe.

The coupling between local order and the electronic state is also an important issue for models of WDM. Changes in the features of the absorption edge indeed reflect this coupling. The general trend of the modification of the K-edge spectra along the Hugoniot is well reproduced by DFT calculations with the exception of the value of the shift of the edge at 500 GPa (0.5 eV measured and 2 eV calculated). Since the K-edge shift is related to the energy difference, $E_F-E_{1s}$, between the Fermi energy and the 1s core energy level, it reflects the modification of the core energy level and of the electron density. In Fe, this K-edge shift along the Hugoniot is predicted to be much smaller than in Al⁵, in line with our observation. In the Warm Dense regime, this shift is particularly sensitive to changes in electronic shielding due to ionization, and to continuum lowering depression. The XANES data therefore provide important constraints on the electronic treatment of WDM within the DFT framework.
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Evolution of the absorption at the K-edge onset. Left: zoom over the edge region of a selection of spectra acquired in the highly focused configuration of the laser. a, b, c and d indicate regions where major changes are observed. Right: ab-initio molecular dynamic simulations (based upon private communication with Vanina Recoules from refs 3 and 11).</figcaption></figure>

In summary, we present data from a proof-of-principle experiment demonstrating that single pulse XAS synchrotron measurements will enable determination of structural and electronic properties of WDM with accuracy close to that on materials at ambient. Bringing the standard of synchrotron material studies towards the measurement of the atomic and microscopic properties of WDM will help validate models at the microscopic level and reveal new phenomena. We provide the first EXAFS data of Fe compressed to 500 GPa and heated to 17000 K, covering a region of the phase diagram not previously explored conclusively with x-ray diagnostics. Our data indicate a persistent local order in the dense Fe fluid. The experimentally measured value of the shift of the absorption edge at 500 GPa with respect to ambient is found to be smaller than the theoretical prediction, stimulating further work to improve the electronic treatment of WDM in the DFT framework. The present experimental approach could easily be extended to higher pressures and temperatures by using a more powerful (100–200 J) laser that would remain compatible with the beamline configuration.

## Methods

**Target preparation.** The targets were made of a thin (3.5 µm) deposit of pure Fe sandwiched between two diamond windows, whose role was to confine the shock for a time interval much larger than the synchrotron pulse duration (100 ps). The Fe deposit was realized by the DEPHIS company and the Fe initial density was measured to be that of the solid (7.87 ± 0.08 g/cc). This was confirmed by the obtained x-ray absorption jump, which is dependent on the density. The adherent Fe deposit was realized on the rear window, a good contact to the ablation window was assured by an interstitial thin layer of CH (<1 µm) and by the windows holder specifically designed to “close” the sandwich by exerting a slight mechanical pressure. Diamond windows of around 40 µm (front) and 50 µm (rear) were used when the laser was focused to 350 µm with the phase plate, whereas thinner windows (25 µm) were used when the laser was focused down to 90 µm to limit the 2D effects. The lateral dimension of the windows, and thus of the Fe deposit was 2 mm, and the holder opening was 1.2 mm. A 4 µm layer of CH was deposited on the ablation face to limit the pre-heating effect.

**Laser, synchronization and alignment.** The experiment was performed with the QUANTEL portable high power laser system provided by CEA. The maximum energy of the laser was 35 J, delivered in a square pulse of 10 ns. The diameter of the laser focal spot was around 350 µm and 90 µm with and without a phase plate respectively, as measured with a 30 degree incidence angle using a blade and sending the image to a Spiricon camera. The laser maximum repetition rate was 45 sec. All the laser clocks were locked to the ESRF radiofrequency system with the proper frequency divisions. The delay between the laser and the x-ray pulse was pre-aligned using a GaAs APD. A fine temporal adjustment was performed directly looking at the onset of the compressed hcp state
in the sample evident through the changes in the XANES. The spatial alignment between the laser and the x-ray beam was performed using the following procedure. First the x-ray beam position was determined by scanning a reference target containing a metallic foil. We then used a high magnification camera looking at the front face of the target to visualize the laser at low power and bring it to the x-ray beam position. Since the laser spot size was always much bigger than the x-ray spot, the error in the alignment is given by the accuracy in the determination of the x-ray beam position which we estimate to be around 10µm.

**Hydrodynamic simulation of the target.** 1D Hydrodynamic simulations using the Lagrangian ESTHER¹² code were performed using a mesh of around 900 units to describe the diamond-Fe-diamond sandwich. The EOS for plastic and diamond were taken from SESAME tables 7592 and 7834 respectively while a BLF²⁸ (Bushman, Lomonosov and Fortov) multiphase equation was used for iron. Opacities are given by the NOHEL code by A. Decoster and at low temperature by measurement from Henke³¹. The strength of diamond is taken into account by using a perfectly elastoplastic model, which includes a simple constitutive model where the elastic yield stress and the shear modulus are constant in the solid domain and vanish to zero when the diamond temperature reaches the melting temperature. The values of the elastic yield strength and shear modulus used to simulate the strength of diamond are 7.5 10⁹ (erg/cm³) and 4.77 10¹² (erg/cm³) respectively.

The 1D Hydrodynamic simulations using the MULTI¹³ code were performed using a mesh of around 550 units for the thin targets and 800 for the thick targets. The SESAME 7830 was used for the diamond, and the SESAME 2150 for iron. Opacities are taken from a SNOP model³². In both cases thermal conductivity and radiative transfer were included in the calculation.

**EXAFS analysis.** The EXAFS quantitative analysis was performed using the IFFEFIT package³³ for the spectra up to laser intensities of 4 × 10¹² W/cm². A five shells model was used for the fit to the ambient bcc spectrum; assuming an undistorted bcc structure, only one distance parameter was fitted. Temperature and Debye temperature were set to ambient values. The volume was then calculated as $V_{bcc} = a^3/2$. Other fitting parameters were the energy offset $\Delta E$ and the amplitude reduction factor $S_0^2$. The obtained value of $S_0^2 = 0.7 \pm 0.1$ was then fixed in the fits to the compressed hcp structures and its error used to evaluate the temperature error; in fact, temperature and $S_0^2$ are strongly correlated, therefore the main source of error in the temperature determination comes from the uncertainty over the $S_0^2$. A two shells model was used for the hcp structure, assuming a homogeneous compression between the two, so that only one distance parameter was fitted. The fit allows the lattice parameter $a_{hcp}$ and thus volume $V_{hcp} = 0.433 a^2 c$ to be obtained. The c/a parameter was fixed to the value 1.61 since the k-range of the EXAFS was not extended enough for its precise adjustment; however simulations and static measurements¹⁸⁻²² have shown that the c/a value remains almost ideal (1.59–1.63) over the thermodynamic region of interest. Other fitting parameters were the energy shift $\Delta E$ and the temperature $T$. The k-range of the collected data was not sufficient for the analysis to include anharmonic corrections. However, ab-initio molecular dynamics simulations suggest that these are small in hcp Fe up to inner core conditions²². A correlated Debye model was used for the temperature fit²³ which is appropriate in our temperature range⁴. The Debye temperature at different compressions in the hcp phase is taken from published measurements²⁴.

## References

1. Murillo, M. S. Strongly coupled plasma physics and high energy-density matter. *Physics of Plasmas* **11**, 2964–2971 (2004).

2. *Fontiers and challenges in warm dense matter.* (eds Graziani, F. et al.) (Springer 2014).

3. Mazevet, S. et al. Ab-initio calculations of x-ray absorption of iron up to 3 Mbar and 8000K. *Phys Rev B* **89**, 100103 (2014).

4. Ping, Y. et al. Solid iron compressed up to 560 GPa. *Phys. Rev. Lett.* **111**, 065501–5 (2013).

5. Benuzzi-Mounaix, A. et al. Electronic structure investigation of highly compressed aluminum with K edge absorption spectroscopy. *Phys. Rev. Lett.* **107**, 165006 (2011).

6. Denoeud, A. et al. Metallization of warm dense SiO₂ studied by XANES spectroscopy. *Phys. Rev. Lett.* **113**, 116404 (2014).

7. Aquilanti, G. et al. Melting of iron determined by X-ray absorption spectroscopy to 100 Gpa. *PNAS.* **112**, 12042–12045 (2015).

8. Luo, S. et al. Gas gun shock experiments with single-pulse x-ray phase contrast imaging and diffraction at the advanced photon source. *Rev. of Scient. Instrum.* **83**, 073903–10 (2012).

9. Milathianaki, D. et al. Femtosecond visualization of lattice dynamics in shock-compressed Matter. *Science* **342**, 220–223 (2013).

10. Pascarelli, S. et al. The Time resolved and Extreme conditions XAS (TEXAS) facility at the European Synchrotron Radiation Facility: the energy dispersive x-ray absorption spectroscopy beamline ID24. *Journal of Synchrotron Radiation* **23**, 353–368 (2016).

11. Harmand, M. et al. X-ray absorption spectroscopy of iron at multimegabar pressures in laser shock experiments. *Phys Rev B* **92**, 024108 (2015).

12. Colombier, J. P., Combis, P., Bonneau, F., Le Harzic, R. & Audouard, E. Hydrodynamic simulations of metal ablation by femtosecond laser irradiation, *Phys. Rev. B* **71**, 165406 1–6 (2005).

13. Ramis, R., Schmalz, R. & Meyer-Ter-Vehn, J. MULTI - A computer code for one-dimensional multigroup radiation hydrodynamics. *Comp. Phys. Comm.* **49**, 475–505 (1988).

14. Atzeni, S. et al. Fluid and kinetic simulation of inertial confinement fusion plasmas, *Comput. Phys. Commun.* **169**, 153–159 (2005).

15. Lindl, J. *Phys. Plasmas* **2**, 3933–4024 (1995).

16. Mathon, O. et al. Dynamics of the Magnetic and Structural $\alpha-\epsilon$ Phase Transition in Iron, *Phys. Rev. Lett.* **93**, 255503 (2004).

17. Yaakobi, B. EXAFS Measurement of Iron bcc-to-hcp Phase Transformation in Nanosecond-Laser Shocks. *Phys. Rev. Lett.* **95**, 075501 1–4 (2005).

18. Anzellini, S., Dewaele, A., Mezouar, M., Loubeyre, P. & Morard, G. Melting of iron at earth's inner core boundary based on fast X-ray diffraction. *Science* **340**, 464–466 (2013).

19. Tateno, S., Hirose, K., Ohishi, Y. & Tatsumi, Y. The structure of iron in Earth's inner core, *Science* **330**, 359 (2010).

20. Steinle-Neumann, G., Stixude, L., Cohen, R. E. & Gülseren, O. Elasticity of iron at the temperature of the Earth's inner core. *Nature* **413**, 57 (2001).

21. Belonoshko, A. B., Ahuja, R. & Johansson, B. Stability of the body-centred-cubic phase of iron in the Earth's inner core. *Nature* **424**, 1032 (2003).

22. Gannarelli, C. M. S., Alfe, D. & Gillan, M. J. The axial ratio of hcp iron at the conditions of the Earth's inner core. *Phys. Earth Planet. Inter.* **152**, 67 (2005).
23. Sevillano, E., Meuth, H. & Rehr, J. J. Extended x-ray absorption fine structure Debye-Waller factors. I. Monatomic crystals, *Phys. Rev. B* **20**, 4908 (1979).

24. Anderson, O. L., Dubrovinsky, L., Saxena, S. K. & LeBihan, T. Experimental vibrational Gruneisen ratio values for ε-iron up to 330 GPa at 300 K, *Geop. Res. Lett.* **28**, 399–402 (2001).

25. Brown, J. M. & McQueen, R. Phase transition, Gruneisen parameter and elasticity for shocked iron between 77 GPa and 400 GPa. *J. Geophys. Res.* **91**, 7485–7494 (1986).

26. McQueen, R. G., Marsh, S. P., Taylor, J. W., Fritz, J. N. & Carter, W. J. The equation of state of solids from shock wave studies in *High Velocity Impact Phenomena* (Ed. Kinslow, R.) 293–417 (New-York: Academic Press 1970).

27. Kerley, G. I. Multiphase equation of state for iron. *Sandia National Laboratories report N°SAND93-0027* (1993). Available at: http://prod.sandia.gov/techlib/access-control.cgi/1993/930027.pdf.

28. Bushman, A. V., Lomonosov, I. V. & Fortov, V. E. *Equations of State of Metals at High Energy Densities* (Institute of Chemical Physics, Russian Academy of Sciences, Chernogolovka, 1992).

29. Ping, Y. et al. A platform for x-ray absorption fine structure study of dynamically compressed materials above 1 Mbar. *Rev. of Sci. Instrum.* **84**, 123105 (2013).

30. Anzellini, S. Phase diagram of iron under extreme conditions measured with time resolved methods, PhD Thesis, Université Pierre et Marie Curie (2014).

31. Henke, B. L., Gullikson, E. M. & Davis, J. C. X-ray interactions: Photo-asorption, scattering, transmission, and reflection at E = 50–30000 eV, Z = 1–92, *Atomic Data and Nuclear Data Tables*, **54**, 181–342 (1993).

32. Tsakiris, G. D. & Eidmann, K. An approximate method for calculating Planck and Rosseland mean opacities in hot, dense plasmas. *Journal of Quantitative Spectroscopy and Radiative Transfer* **38**, 353–368 (1987).

33. Ravel, B. & Newville, M. ATHENA, ARTEMIS, HEPHAESTUS: data analysis for X-ray absorption spectroscopy using IFEFFIT, *J. Synchrotron Rad.* **12**, 537–541 (2005).

## Acknowledgements

We wish to thank Sebastien Paskernak, Florian Perrin, Gilles Berruyer, Marion Harmand and Katilin Mecsek for their assistance with the experiment. We thank Vanina Recoules for providing the DFT spectra used in Fig. 5 and for a critical reading of the manuscript. We acknowledge S. Atzeni, A. Marocchino, A. Schiavi for the use of 2D hydrodynamic code DUED. The ISP is supported by AWE Aldermaston plc and Imperial College London.

## Author Contributions

F.O., R.T., O.M., S.P. and P.L. designed the project. F.O. and O.M. performed all tests prior to the experiment. F.O., O.M., S.P., A.S., E.L. and R.T. conducted the experiment. R.T., P.L. and S.P. analyzed the data. A.S. and E.L. operated the laser. J.H., W.H. and O.M. designed, tested and operated the x-ray detector. R.T. and F.O. prepared the targets. R.T., L.V., T.V. and A.B.-M. performed hydro-simulations. All authors assisted with synchrotron experiments and discussed the results. R.T., S.P. and P.L. wrote the manuscript.

## Additional Information

**Competing financial interests:** The authors declare no competing financial interests.

**How to cite this article:** Torchio, R. et al. Probing local and electronic structure in Warm Dense Matter: single pulse synchrotron x-ray absorption spectroscopy on shocked Fe. *Sci. Rep.* **6**, 26402 (2016).

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

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/