# Melting properties by X-ray absorption spectroscopy: common signatures in binary fe-c, fe-o, fe-S and fe-Si systems

Silvia Boccato, Raffaella Torchio, Simone Anzellini, Eglantine Boulard, François Guyot, Tetsuo Irifune, Marion Harmand, Innokenty Kantor, Francesca Miozzi, Paraskevas Parisiades, et al.

► To cite this version:

Silvia Boccato, Raffaella Torchio, Simone Anzellini, Eglantine Boulard, François Guyot, et al. Melting properties by X-ray absorption spectroscopy: common signatures in binary fe-c, fe-o, fe-S and fe-Si systems. Scientific Reports, 2020, 10 (1), 10.1038/s41598-020-68244-3. hal-02934656

HAL Id: hal-02934656

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

Submitted on 9 Sep 2020

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
<figure><img src="image_2.png" /></figure>

OPEN

# Melting properties by X-ray absorption spectroscopy: common signatures in binary Fe–C, Fe–O, Fe–S and Fe–Si systems

Silvia Boccato<sup>1,2,✉</sup>, Raffaella Torchio<sup>2</sup>, Simone Anzullini<sup>3</sup>, Eglantine Boulard<sup>1</sup>, François Guyot<sup>1,4</sup>, Tetsuo Irifune<sup>5</sup>, Marion Harmand<sup>1</sup>, Innokenty Kantor<sup>2,6</sup>, Francesca Miozzi<sup>1</sup>, Paraskevas Parisiades<sup>1</sup>, Angelika D. Rosa<sup>2</sup>, Daniele Antonangeli<sup>1</sup> & Guillaume Morard<sup>7</sup>

X-ray absorption spectroscopy (XAS) is a widely used technique to probe the local environment around specific atomic species. Applied to samples under extreme pressure and temperature conditions, XAS is sensitive to phase transitions, including melting, and allows gathering insights on compositional variations and electronic changes occurring during such transitions. These characteristics can be exploited for studies of prime interest in geophysics and fundamental high-pressure physics. Here, we investigated the melting curve and the eutectic composition of four geophysically relevant iron binary systems: Fe–C, Fe–O, Fe–S and Fe–Si. Our results show that all these systems present the same spectroscopic signatures upon melting, common to those observed for other pure late 3d transition metals. The presented melting criterion seems to be general for late 3d metals bearing systems. Additionally, we demonstrate the suitability of XAS to extract melt compositional information in situ, such as the evolution of the concentration of light elements with increasing temperature. Diagnostics presented in this work can be applied to studies over an even larger pressure range exploiting the upgraded synchrotron machines, and directly transferred to time-resolved extreme condition studies using dynamic compression (ns) or fast laser heating (ms).

The coupling of X-ray absorption spectroscopy (XAS) with laser heated-diamond anvil cell (LH-DAC) is a challenging technical development only recently implemented in synchrotron facilities<sup>1-3</sup>. The opportunities that stem from this combination are of great interest in physics, chemistry and planetary science.

The element selectivity and the local structure sensitivity are two of the distinctive characteristics that make XAS a very informative and multipurpose technique. Combined together, they provide the possibility of microscopic investigation of the local atomic environment of a specific element in a compound or an alloy, irrespectively of the solid, liquid or amorphous state. The X-ray absorption near edge spectroscopy (XANES) region of the XAS spectrum, also probes the unoccupied electron density of states right above the Fermi level<sup>4-7</sup>. Recent high pressure/high temperature absorption experiments<sup>6,8-14</sup> and theoretical calculations<sup>8,9,11,15-19</sup> on pure metals directly correlate changes in the XANES spectra and/or in the pre-edge to phase transitions, but to which extent these findings could be generalized to multi-component systems remains to be assessed. XANES sensitivity to the electronic structure further grants the possibility to probe the chemical environment, through relatively simple approaches like the linear combination analysis (LCA), as demonstrated in many studies at ambient conditions (Benfatto et al.<sup>20</sup>, among others), but never attempted under extreme high-pressure, high-temperate conditions.

<footnote n="1"><sup>1</sup>Sorbonne Université, Muséum National d'Histoire Naturelle, UMR CNRS 7590, Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie (IMPMC), 75005 Paris, France. <sup>2</sup>ESRF - European Synchrotron Radiation Facility, Grenoble, France. <sup>3</sup>Diamond Light Source Ltd, Harwell Science and Innovation Campus, Didcot OX11 0DE, UK. <sup>4</sup>Institut Universitaire de France (IUF), Paris, France. <sup>5</sup>Geodynamics Research Center, Ehime University, 790-8577 Matsuyama, Japan. <sup>6</sup>Present address: Department of Physics, Technical University of Denmark, Kgs. Lyngby, Denmark. <sup>7</sup>Université Grenoble Alpes, Université Savoie Mont Blanc, CNRS, IRD, IFSTTAR, ISTerre, 38000 Grenoble, France. ✉email: silvia.boccato@upmc.fr</footnote>
In this paper, we show advanced application of XANES in combination with LH-DAC. Melting criteria are generalized and used to probe melting curves of Fe-X binary systems, where X = C, O, S and Si. Furthermore, the quality of the collected data allowed us to apply LCA-based methods for data analysis, obtaining a quantitative in situ determination of chemical composition of the temperature-quenched samples under pressure, during the heating runs. The here-presented experimental protocol and data analysis offer a valuable method alternative to more classic, and not always possible, time-demanding, ex-situ analysis of recovered samples.

Binary iron alloy systems have been here chosen as test case to prove the capabilities of XANES for the experi-
mental study of materials at high pressure and high temperature also in view of their very important implications
for the understanding of planetary cores. Position in temperature and in composition of the eutectic point in
a binary phase diagram anchors the reconstruction of the entire phase diagram. In the case of iron alloys, the
primary constituents of planetary cores, this allows to determine which phase crystallizes upon secular cool-
ing, yielding different crystallization paths and hence impacting the generation and sustainability of planetary
dynamo²¹. Moreover, the knowledge of the chemical composition is a primordial information required to know
the density contrast between solid and liquid, and thus crucial to determine whether the crystallizing solid is
gravitational stable, whether it sinks or floats and to which extent crystallization might power solutal convection.

In short, XANES is a powerful tool to determine the melting curve and to obtain direct information of the
electronic structure of probed materials. Some of the specific advantages of this technique are illustrated in this
work by providing the measurements of the eutectic temperatures and eutectic compositions over an extended
pressure and temperature range of Fe-X binary systems of primary geophysical interest, where X = C, O, S and Si.

**Analysis of XANES spectra**

**Melting criterion with XANES.**

**Melting criterion with XANES.** Classically, studies investigating metals at ambient pressure and high temperature use discontinuities in the temperature evolution of the absorption signal at a single photon energy to detect solid-solid phase transitions¹⁵, melting and undercooling of pure substances²²⁻²⁴. This method is also referred to as T-scan method. At simultaneous high pressures and temperatures, the energy and temperature resolutions are lowered by experimental constraints such as the need for fast acquisitions and/or the micrometric sample volumes in the central portion of the laser heated hot-spot area. The sample thickness, and consequently the absorption jump in the absorption-wavelength diagram, might vary during the experiment and thus the application of the T-scan method requires an accurate normalization. As already mentioned, the loss of features in the XANES upon increasing temperature has recently been proposed and validated as melting criterion for elemental 3d metals such as iron, nickel and cobalt in experiments very similar to the ones presented in this work in terms of the covered pressure-temperature range and heating method⁸⁻¹⁰,¹²,²⁵. In most cases results obtained by XAS experiments showed a very good agreement with those obtained by X-ray diffraction (XRD)¹⁰,¹²,¹³.

Pure Fe melts from a solid with a face-centered cubic (fcc) structure over a wide pressure range, from 15 to about 100 GPa¹³,²⁶. Melting from fcc structure over a comparable pressure range is common for all late 3d transition metals, i.e. metals with 3d orbitals more than half filled. The melting criterion in these systems is the flattening of the XANES shoulder and the disappearance of the first two oscillations after the edge, typical of an fcc structure. We observe that this melting criterion is also valid for the Fe–C and Fe–O systems (see left panel of Fig. 1 and “Methods” section for details on XAS data collection and normalization). We also note that in the investigated pressure range, Fe–C exhibits, before melting, a temperature-induced solid–solid phase transition, evidenced by a clear change in the XANES features at around 7130–7140 eV (violet spectrum in the top left panel of Fig. 1). The nature of this solid–solid transition will be discussed in the “Results” section. The high-temperature solid phase melts upon heating at 2550 K. Melt in the Fe–O system was observed at 50 GPa at 2920 K through a clear loss of features in the XANES (Fig. 1). XANES spectra collected on Fe–S and Fe–Si systems, when compared to Fe–C and Fe–O, show fewer features and less clear variations with increasing temperatures (right panel of Fig. 1). The solid–liquid transition is however still distinguishable from the flattening of the XANES features around the energies identified with B and D in Fig. 1. We therefore applied the T-scan method¹²,¹⁵,²²–²⁵ as a second criterion to double-check estimated melting temperature. In the inset of Fig. 1 we plot the variation of the absorption coefficient at a specific energy as a function of temperature, hereby called δα following the nomenclature in Di Cicco et al.²⁴. The transition from the solid to the liquid state leads to a discontinuous variation and an abrupt change of XANES features due to the loss of the long-range and medium-range order upon melting. In the insets of Fig. 1 the δα (A = 7134.0 eV for Fe–C, B = 7135.1 eV for Fe–S, C = 7130.5 eV for Fe–O and D = 7131.5 eV for Fe–Si) shows clear discontinuities that are associated to melting for the four binary systems. The energy for the T-scan method is chosen in order to maximize the contrast in the δα as a function of temperature. This value depends on the atomic structure of the last solid and the difference in the electronic density of states between the hot solid and the liquid²⁵.

The same changes observed in the XANES in this work were related to melting in previous experiments con-
ducted on pure late transition 3d metals such as iron⁸⁻¹⁰, nickel¹² and cobalt²⁵ and recently on Fe-Ni²⁷. In the case
of cobalt, ab-initio calculations showed that such changes are due to the appearance of multiple configurations
proper of the liquid phase and associated to the loss of crystalline order. The electron density of states also varies
at the energies where discontinuous changes are detected. The loss of XANES features, together with disconti-
nuities in the δα, can thus be used as a generalized melting criterion for late 3d transition metal based systems.

In the present study the adopted melting criterion was further verified through textural and chemical analysis of several recovered samples using the focused ion beam (FIB) and scanning electron microscope (SEM) approach detailed in Morard et al.²⁸. An example will be discussed later in the “Results” section and illustrated in Fig. 4.
<figure><img src="image_3.png" /><figcaption>**Figure 1.** XAS spectra at Fe K-edge collected for increasing temperature showing phase transitions and melting for Fe–C, Fe–O, Fe–S and Fe–Si. The spectra in black are at ambient temperature, the blue are solids, the violet (in the Fe–C system) is a solid with a different structure, the orange spectra are first appearance of liquid and the red are liquids. For clarity, sequences of spectra in the left panel are vertically shifted, while in the right panel they are superimposed. The insets show for each of the four binary systems the change of white line peak δα at the specific energies indicated with letters A, B, C and D respectively. Mixed phases are indicated with orange circles. A discontinuity in the single energy X-ray absorption is an indication of solid-solid phase transition (as for Fe–C) or melting. Horizontal dashed lines are a guide for the eyes.</figcaption></figure>

**Composition determined in situ by XANES.** As already mentioned, the XANES signal, due to its sensitivity to the electronic structure of the absorber, contains information about its chemical environment. The chemical composition of the sample can thus be determined by proper analysis of the spectra collected on the temperature-quenched samples, provided that XANES of reference compounds are available²⁰ (and preferably at close P-T conditions). In Fig. 2, the XANES spectra of pure Fe and stoichiometric Fe₃C and FeO are compared to the ones obtained in the binary Fe–C and Fe–O systems.

To quantify the relative amounts of end-members in Fe–C and Fe–O liquid binary systems, we used a linear combination analysis (LCA)²⁰ employing the program ATHENA³⁰. In this approach, the normalized spectra of temperature-quenched samples are fitted to a theoretical curve obtained as linear combination of the normalized spectra of stoichiometric samples.

A set of high-quality reference spectra at the relevant pressures is of crucial importance for obtaining reliable results. First of all because pressure can induce phase transitions and secondly, even when no transitions occur, the effect of pressure is to shift the EXAFS and the XANES features to higher energies. This shift was shown to start gradually few eV higher than the edge energy²⁵, therefore it is quite small in the edge region, while it becomes important in the EXAFS domain. Here a short range for the fitting around the edge has been chosen so to minimize any error due to small pressure offsets between the analyzed spectrum and the reference.

Reference spectra for iron and stoichiometric FeO at ambient temperature, quenched after heating at different pressures, were collected in this study and during previously published experiments (iron¹³ and FeO³¹). Due to the absence of high-pressure references, for the analysis of Fe₃C we used XANES measured at ambient conditions²⁹. We note however that for the examined compositional space of the Fe–C system and in the chosen energy range, the differences in the near edge features between the two end-members appear to be largely
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Absorption spectra at Fe K edge measured for the iron-bearing binary systems can be well reproduced by linear mixing of reference spectra of the corresponding end members, allowing the quantification of the amount of light element in each binary system. Example of linear combination analysis (LCA) of Fe–C (left) and Fe–O (right). Reference spectra of stoichiometric end-members are also represented (Fe₃C reference is from Yuan et al.²⁹, Fe is from Morard et al.¹³ and FeO is from Boulard et al.³¹), and their intensities are proportional to their weight resulted from the fit. Being E₀ = 7112 eV, the range for the fit was 7107–7130 eV for the Fe–C system and 7107–7147 eV for the Fe–O system, as indicated by dashed vertical lines. This range is chosen in order to have the best fit quality within the largest energy range. In the case of Fe and FeO the reference spectra are temperature-quenches collected after heating at the indicated temperature.</figcaption></figure>

dominant in respect of the pressure variations of Fe₃C XANES spectra, thus allowing to perform LCA analysis even though missing a set of stoichiometric Fe₃C spectra at relevant pressures.

An example of fit performed with LCA is shown in Fig. 2, together with the used reference spectra. Starting from the end-member proportions obtained from LCA, it is a straight exercise to quantify the eutectic composition, and specifically the amount of light elements in the eutectic liquid of the Fe–C and Fe–O iron binary systems, using the known content of carbon and oxygen in stoichiometric Fe₃C (6.7 wt% of carbon) and in stoichiometric FeO (22.3 wt% of oxygen).

## Results

**Eutectic melting curves.** All the investigated pressure–temperature points for the four different binary systems are shown in Fig. 3. Following the same color code as in Fig. 1, blue and violet spectra represent solid phases, orange spectra highlight the onset of melting and red spectra show liquid phases. In the case of Fe–S and Fe–Si, two slightly different starting compositions were used and are distinguished in Fig. 3 as circles and squares.

In all the four cases we observe an overall good agreement, within the pressure and temperature uncertainties, between our results and previous works, independently of the sample environment, the adopted detection technique or the starting composition. Over the pressure range of interest in this study, samples are most commonly pressurized and heated with LH-DAC, or multi anvil press for the lower pressures³²⁻³⁴. Melting detection methods include visual observation of sample surface movements³⁵, the plateau in the temperature versus laser power plot³⁶, appearance of liquid diffuse scattering in XRD measurements or disappearance of the solid peaks²⁸,³⁷⁻⁴⁴, and ex situ textural analysis (FIB coupled with SEM or microprobe) on recovered quenched samples²⁸,³²,³⁴,⁴⁰,⁴⁴. Starting compositions vary for Fe–C between 1.5 and 4.0 wt% C, for Fe–O between 1.2 and 11 wt% O, for Fe–S between 4 and 16.1 wt% S and for Fe–Si between 10 and 18 wt% Si. Same melting temperature for a system with different starting compositions is expected for incongruent melting, where the obtained liquid is at the eutectic composition.

Considering each individual system more in detail, in the top left panel of Fig. 3, blue and violet symbols indicate two different solid phases of Fe–C. Analysis of XANES suggests that the low-temperature phase (in blue) is dominantly hexagonal close packed (hcp) and the high-temperature phase (in violet) is dominantly fcc. The pressure–temperature conditions of this hcp to fcc phase transition are, within the error bars, compatible with the ones reported by Anzellini et al.²⁶ for pure Fe (thin black line in the figure), suggesting that the temperature induced hcp–fcc solid–solid transition is only slighted affected by C doping. In contrast, the melting temperature is significantly lowered by the presence of carbon. We compiled present XAS results and most recent literature
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Eutectic melting curves of four Fe–X binary systems are presented and compared to the literature: Fe–C<sup>28,36,39,44,45</sup>, Fe–O<sup>28,35,37</sup>, Fe–S<sup>28,32,34,40–42</sup> and Fe–Si<sup>28,33,38,43</sup>. Uncertainties on pressure and temperature determination are represented as black crosses for low and high P–T data points in the Fe–C system and are equivalent for all four binary systems at comparable pressure–temperature conditions. For reference, the melting curve of pure Fe as measured by XRD by Anzellini et al.<sup>26</sup> and recently confirmed by XAS<sup>13</sup> is shown as thin black line. The thick black curves represent a Simon–Glatzel fits to data points for which the first appearance of the liquid was observed in this study through XAS together with points from selected studies in literature (see main text for extended discussion). Fit parameters are reported in Table 1.</figcaption></figure>

results<sup>28,44</sup> and fit the whole data set to a Simon–Glatzel equation. The resulting melting curve is represented as a thick black curve in Fig. 3.

XAS data of Fe–O shown in the bottom left panel of Fig. 3 are in very good agreement with the most recent measurements<sup>28</sup> while a slight discrepancy can be noticed with the earlier data<sup>35,37</sup>. The joint Simon–Glatzel fit of our new data together with previous XRD data<sup>28</sup> results in the thick black curve that, in the range under analysis, is virtually indistinguishable from that from Morard et al.<sup>28</sup> (green thick curve).

In the top right panel of Fig. 3 the melting curves obtained in the Fe–S system, with and without the addition of 5 wt% of Ni, are compared (Fe-12 wt% S and Fe-5 wt% Ni-12 wt% S are represented with circles and squares respectively). The obtained melting temperatures are similar for both compositions, suggesting that the addition of 5 wt% of nickel has negligible effect on the melting temperature within current experimental uncertainties, which is in agreement with recent findings on the Fe–Ni melting curve<sup>27</sup>. The difference of about 125 K at 20 and 40 GPa previously reported by Stewart et al.<sup>34</sup> is indeed within present uncertainties and could therefore not be resolved. The melting curve for the Fe–S system, represented with a thick black curve in Fig. 3, was obtained by fitting a Simon–Glatzel equation to all our data points, irrespectively of the addition of nickel, together with data from Morard et al.<sup>28</sup> and Mori et al.<sup>40</sup>, the most recent literature results. The resulting melting curve is virtually indistinguishable from the one reported by Morard et al.<sup>28</sup>. Data from Fei et al.<sup>32</sup>, Stewart et al.<sup>34</sup> and Kamada et al.<sup>41,42</sup> are scattered around the fitted curve, thus in agreement within the uncertainties. Simon–Glatzel fit parameters for these three binary systems are reported in Table 1.

Previous studies reported a dependency of the melting temperature on the starting composition in the Fe–Si binary system, arguing for a eutectic behavior. However, in the present study we could not resolve a difference over the investigated pressure range (up to 40 GPa) for starting composition with 18 wt% and 10 wt% Si shown as circles and squares in Fig. 3 respectively. In view of the limited pressure range covered in this study, we did not attempt to fit the melting curve, but we stress that our melting temperatures are in very good agreement with melting curve proposed by Morard et al.<sup>28</sup> based on previous results<sup>38,43</sup>.
<table><caption>**Table 1.** Fit parameters of Simon–Glatzel equation: $T_m = T_0[(P_m - P_0)/a + 1]^{1/c}$.</caption><thead><tr><th></th><th>T<sub>0</sub> (K)</th><th>P<sub>0</sub> (GPa)</th><th>a</th><th>c</th></tr></thead><tbody><tr><td>Fe–C</td><td>1420</td><td>0</td><td>5.9</td><td>4.1</td></tr><tr><td>Fe–O</td><td>1800</td><td>0</td><td>21.6</td><td>3.3</td></tr><tr><td>Fe–S</td><td>1260</td><td>21</td><td>10.3</td><td>3.0</td></tr></tbody></table>

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Image on secondary electron secondary ion (SESI) detector (left panel) and elemental concentration maps from EDX-SEM (energy dispersive X-ray - scanning electron microscope) analysis (right panel) of a transverse section of quenched Fe–S sample of Fe-12 wt% S starting composition, recovered after laser heating at 34 ± 3 GPa and 1730 ± 120 K. Visible in the SESI image are the sample are the Fe–S sample (center, light grey area) and two KCl disks at the top and the bottom (dark grey). Clear textural differences can be distinguished in the sample. The homogeneous texture in the central area can be associated with a quenched liquid coexisting at high temperature with a solid formed at the rim of the laser heated spot (heterogeneous textures). The concentrations are shown for S (green), Fe (red), K (yellow) and Cl (violet) in the EDX-SEM maps. Chemical maps allow excluding melt contamination by K or Cl. Preferential liquid–solid partitioning of sulfur upon melting led to its depletion in the solid formed at the melt pool border.</figcaption></figure>

Out of the iron binary systems investigated in the present work, only Fe–Si exhibits a melting temperature comparable to that of pure iron. In contrast, the addition of O, C and S significantly lowers the melting temperatures, with differences reaching 500 K (O), 750 K (C) and 1000 K (S) at 90 GPa.

**Eutectic compositions.** The composition of the eutectic liquid was derived from XAS data obtained in situ at high pressures on samples quenched after melting and ex situ by microprobe measurements on recovered, FIB cut samples. In the case of ex situ chemical measurements we also determined the composition of the solid phase coexisting with the eutectic liquid. The compositions of the liquid and its temperature-quench are assumed to be the same. Laser heating in a DAC induces significant thermal gradients in the sample that also lead to chemical diffusion of elements, the so-called Soret effect$^{16}$. A liquid pocket forms at the center of the hotspot, surrounded by a coexisting solid phase (Fig. 4). These phases are in chemical equilibrium due to the fast diffusion coefficients at high temperatures for the elements present in the investigated systems. With increasing temperature above the eutectic point, the liquid portion will grow in size but will not significantly change in composition. This phenomenon was reported in previous LH-DAC studies on binary Fe–X systems$^{28,40,44,47}$.

For the Fe–S system, recovered samples were analyzed with EDX-SEM in order to check chemical contamination and for textural analysis (see Fig. 4). A global decreasing trend in S content in the eutectic liquid with increasing pressure was identified, with the actual S amount slightly higher compared to previous studies$^{32,34,40–42,50,53}$. The Fe–Si binary system is expected to melt with a 1:1 solid-liquid partitioning$^{54,55}$, so it was not considered in the following analysis.

Literature results for the Fe–O binary system show a large scatter of the liquid composition. The difference in the two most recent studies$^{28,47}$ can be rationalized by C contamination in Morard et al.$^{28}$. Our estimations, obtained by analysis of XANES data based on the available appropriate XANES references for both the end-members at high pressure, are in between, and show a trend compatible with a progressive increase in O content in the eutectic liquid with increasing pressure, as suggested by the non ideal model by Komabayashi et al.$^{48}$. For the Fe–O binary system, samples were not recovered.

In the case of Fe–C, the relative proportion between two end-members Fe and Fe$_3$C is evaluated in situ with LCA analysis and the total amount of carbon in the sample is obtained from the derived proportion of Fe$_3$C
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Eutectic liquid compositions from in situ XANES diagnostics on temperature-quenched samples ex situ and/or analysis techniques (microprobe analysis and EDX-SEM), compared to those of the literature for Fe–C<sup>28,36,44</sup>, Fe–O<sup>28,37,47–49</sup> and Fe–S<sup>32,34,40–42,50–52</sup>. XANES points are the average of the results from the analysis of several temperature-quenched data after melting at similar conditions. When the pressure of the spectrum to analyse is in between the pressure for which reference spectra were available, the chemical evaluation was performed with both lower and higher pressure references and the result averaged. The error bar is evaluated as the variance between the different results of analysis from individual spectra, as shown in the Supplementary information.</figcaption></figure>

phase plus the amount of carbon in solution in solid iron, which was estimated according to Mashino et al.<sup>44</sup>. Details on this calculation are presented in Supplementary information. For comparison, carbon content of the quenched samples was also evaluated ex situ by EDX-SEM. The eutectic liquid compositions obtained for the Fe–C binary system with XANES is in excellent agreement with those from EDX-SEM and with estimations from thermodynamics calculations by Fei and Brosh<sup>36</sup>. The eutectic liquid composition is observed to slightly evolve with pressure, in agreement with recent studies<sup>28,44</sup>, but in opposition to the strong negative trend proposed in early work<sup>36</sup>.

Textural analysis independently confirmed melting detection by XANES, but even more relevant is the clear agreement between ex situ and in situ eutectic composition measurements obtained in the present study (Fig. 5).

Discussion

The present work highlights the versatility and capability of XANES when applied in combination with laser heated diamond anvil cells. We demonstrate that solid-solid and solid-liquid phase transitions are unambiguously revealed by distinct changes in XANES features as a function of pressure and temperature (Fig. 1).

The criteria presented in this work for the solid-liquid phase transition are proved to be common not only to late 3d transition metals (Fe, Co, Ni, Cu)<sup>10,12,13,23–25</sup>, but also to several Fe-X systems under extreme conditions of pressure and temperature. Due to the structural and electronic similarities among the late 3d transition metals, we expect this criterion to be valid for a large number of late 3d transition metals bearing binary systems.

Focusing on in situ melting diagnostics, XAS can be used as a complementary technique to XRD. Both XAS and XRD are bulk techniques that allow to precisely track changes in material's properties, in contrast to the visual observation of the sample surface. In XRD, melting is related to the appearance of the diffuse scattering: the diffraction peaks disappear due to the loss of long range order, and are replaced by diffuse scattering signal, which contains information about the short range ordering. The changes detected in the XANES are also associated to intrinsic physical changes resulting from the phase transition from the solid to the liquid<sup>25</sup>: the smearing/loss of features at the absorption edge can be understood as the consequence of the coexistence of multiple configurations in the liquid phase, directly related to the modification of the electronic bands around the Fermi energy.
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Comparison of the XANES spectra at Fe K-edge of molten Fe-C, solid iron at similar pressure and temperatures and a solid iron at room temperature. The loss of features in the Fe-C XANES is a consequence of melting and cannot be ascribed to the damping of the oscillations due to the increase of temperature.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 7.** Fe₃C proportion derived using LCA on XANES measured on temperature-quenched samples, laser heated at the temperature indicated in the vertical axis. In the three cases, for three temperature runs at 38, 57 and 80 GPa, the starting Fe₃C proportion is ~ 22 wt%, corresponding to ~ 1.5 wt% C, well compatible with nominal composition and electron microprobe analysis of starting material. The decrease in the Fe₃C content at increasing temperature can be interpreted as the progressive dissolution of carbon in solid iron. The total amount of carbon has to be evaluated as the carbon in the Fe₃C phase plus the carbon in solid solution with Fe, as shown in Supplementary information. The yellow square in each run indicates the temperature at which the onset of melting was detected: this is the eutectic point (eutectic temperature and eutectic composition).</figcaption></figure>

These changes are in principle independent of the specifically considered system and thus give universality to the melting detection method. As the atomic vibrations become more important the higher is the temperature, the mean square relative displacement (Debye-Waller) increases with temperature, causing the damping of the EXAFS oscillation, in particular at higher energies. This oscillations damping is clearly visible when comparing the measured Fe K-edge spectra of solid iron at room temperature and at high temperature (Fig. 6). Differently, in a solid-liquid transition, sudden changes are visible at low energies, in the XANES region, while oscillations at higher energy still remain. The comparison between solid Fe and liquid Fe-C at very similar pressure and
temperature (also shown in Fig. 6) provides a further confirmation that the XANES changes attributed to melting are clear indication of solid-liquid phase transition and are not due to temperature increase. The oscillations still visible in the EXAFS contain information on the interatomic distance in the liquid phase, and can thus be further exploited¹⁴. However, because of the limited energy range and the data quality, the information extracted remains limited to the first neighbor distance. Thanks to the Extremely Brilliant Source (EBS) upgrade⁵⁷ with consequent renovation of the ID24 beamline⁵⁸, in the near future it will be possible to perform absorption measurements in LH-DAC over an even longer energy range in the EXAFS (~1000 eV) and using a much smaller X-ray beam (~1 µm full width half maximum, FWHM). Better data quality and more extreme conditions will be achieved allowing to reveal the structure of melts at planetary core conditions. In parallel, characterization of structural transformations using a supervised machine learning method, namely, artificial neural network, is ongoing, in order to extract information on the local structure starting from both XANES and EXAFS⁵⁹,⁶⁰. We believe that the combination of high quality data and the use of more sophisticated methods for the analysis will allow to fully explore all the structural and electronic information hidden in the XANES.

XANES measurements can be used to track changes in physical state as a function of temperature by analysis of the temperature quenches after each pulsed heating (here 1 s). LCA analysis of XANES spectra can be also used to obtain compositional information on the probed materials, including the critical case of quenched liquids (Fig. 2). The agreement between the eutectic liquid composition from in situ determination by XANES and from classical ex situ techniques such as EDX-SEM and electron microprobe, together with the overall agreement with previous literature, clearly shows the suitability and great potential of XANES as in situ probe of the eutectic melt composition at extreme conditions, complementary to the more time-consuming FIB ex situ analysis techniques (Fig. 4). Application of XANES is particularly relevant for light elements, difficult to be quantified in situ by other techniques such as XRD (whose sensitivity to the grain size and orientation hinders a quantitative analysis when approaching the melting temperature), and when samples cannot be recovered. Figure 7 shows the results of LCA analysis at three different pressures. The Fe₃C proportion was evaluated in situ on temperature-quenched samples after pulsed heating at the temperature reported in the y-axis. The decrease of Fe₃C proportion while increasing temperature is compatible with the dissolution of carbon in solid iron⁴⁴. Upon further temperature increase, the sample melts and the carbon content equilibrates to the eutectic composition. The carbon content can be evaluated as the one in the Fe₃C phase added to the one known to be in solid solution, as discussed in the Supplementary information. LCA method shown in this work allows to determine both the structural rearrangement and the composition in situ, thus tracking the pressure evolution of the eutectic point (yellow square in Fig. 7), which is a critical information to build phase diagram and to establish accurate thermodynamic models of the investigated systems.

In Aprilis et al.⁶¹ it was shown that chemical reactions of iron samples with diamonds to form iron carbide are independent of the continuous or pulsed nature of the laser heating, at least for a pulse duration down to 2 µs. Here we show that upon laser exposures of 1 or 2 s, chemical variations are detectable simultaneously with the phase transition. Reducing the laser pulse to shorter duration than the microseconds scale would allow identifying if there is a threshold in the pulse duration under which structural rearrangement or phase transitions happen faster than chemical changes. This is only possible thanks to the fast acquisitions proper of energy dispersive XAS.

Moreover, XAS is very suitable for measurements on iron-bearing samples at high pressure, where the thickness maximizing the signal to noise ratio ideally matches the LH-DAC requirements to minimize axial thermal gradients⁶². An initial sample thickness of 4 or 5 µm, allows to simultaneously determine the phase diagram and study the compression of the system in the liquid state. In the case of XRD, a good quality diffuse scattering from a liquid signal requires a thicker sample, and heating above melting temperature is often necessary to achieve large melting fraction versus remaining solid portion. Practically, two different sample thicknesses have to be chosen in XRD experiments dedicated to phase diagram determination (thinner samples) or for the study of the liquid properties (thicker samples).

## Conclusion

We report a systematic study on melting curves of iron-rich binary systems performed by X-ray absorption spectroscopy in combination with laser heated diamond anvil cell. Our experimental results allow confirming and generalizing the melting criterion, establishing the signature of melting in the XANES region common to late 3d metals and iron alloys. Specifically, the eutectic melting of iron alloyed with C, O, S and Si was determined and compared with the literature. The consistency between our in situ results by absorption spectroscopy, ex situ results by electron microprobe analysis, and the literature provided a further proof of the suitability of this technique to detect melting and allowed to refine melting curve and to update Simon-Glatzel fit parameters for Fe–C, Fe–O and Fe–S eutectic systems. We exploited the capability of XANES to probe both the atomic and the electronic structures, which are characteristic for each system, to track in situ changes in physical state together with changes in chemical composition of the investigated samples. In particular XANES have been used to quantify the variation of light elements content during the heating run and after complete melting. Accordingly, XANES can be used as complementary technique to XRD to measure the pressure evolution of the eutectic point (temperature and composition). We expect this approach to be used in future static or dynamic compression studies of iron binary systems, profiting of the recent development in the X-ray spectroscopy field⁹,¹¹. Linear combination analysis might be particularly important for pulsed laser heating, since it allows tracking possible changes in composition with temperature, and in all cases where samples cannot be recovered, such as in most of the dynamic compression experiments.
## Methods

**Starting materials.** Several binary systems have been investigated in this work, with nominal composition: Fe-1.5 wt% C, Fe-5 wt% O, Fe-12 wt% S, Fe-5 wt% Ni 12 wt% S, Fe-18 wt% Si and Fe-10 wt% Si. The Fe–O samples were synthesized by the DEPHIS Company, sputtering iron onto a glass slide under O₂ flow. The other samples were prepared by induction melting at the Institut de Chimie et des Matériaux de Paris-Est (ICMPE) in Thiais, France. The resulting samples, with a sub-micrometric texture and a thickness between 20 and 50 µm, were precompressed between two diamonds down to about 4 or 5 µm. The samples used for this work were from the same batches as the ones previously used in XRD experiments by Morard et al.²⁸,³⁸. For all samples, chemical analysis confirmed the nominal composition within the error bar.

**In situ laser heating XAS experiments in diamond anvil cell.** XAS experiments were performed at the energy dispersive beamline ID24 at the European Synchrotron Radiation Facility (ESRF). In the energy dispersive geometry the incident pink X-ray beam is dispersed by an elliptically bent Si (111) crystal and then focused to about 4x5 µm² (vertical × horizontal) FWHM on the sample position. XAS signal was measured for all the samples at the Fe K-edge 7112 eV. A range in energy between 7000 and 7500 eV allowed to record both X-ray absorption near edge spectroscopy (XANES) and extended X-ray absorption fine structure (EXAFS). The X-rays transmitted by the sample were recorded on a position sensitive detector, a FReLoN camera with an Hamamatsu chip⁶³, where the position was correlated to the energy using an iron standard reference.

High pressure conditions were obtained using membrane driven DAC equipped with either single crystal or nanopolycrystalline⁶⁴ diamonds, having culet diameters ranging from 150/300 µm (beveled diamonds) to 350 µm and a thickness of about 1.6 mm. 200 µm thick rhenium gaskets were pre-indented to a thickness of about one tenth of the culet size. An hole of one third of the culet diameter was laser drilled in the center of the pre-indentation to serve as sample chamber. The samples were then loaded in the DAC embedded in previously dried KCl disks. These disks were prepared at the Institut de Minéralogie, de Physique des Matériaux, et de Cosmochimie (IMPMC, Paris) by compressing KCl powder to the target thickness and then laser cut with a femtosecond-laser to the needed diameter as described in Ref.¹³. Samples were heated from both sides using the laser heating system installed at the ID24 beamline of the European Synchrotron Radiation Facility (ESRF)³, equipped with two CW Nd:YAG lasers (IPG photonics) of 1064 µm wavelength. Radial temperature gradients were reduced by slightly defocusing the lasers to a hot-spot size diameters of ~ 20 µm FWHM (this size was 3–4 times the size of the focused X-ray beam diameter). It should also be noted that the optimal initial sample thickness for iron binary systems is of few microns, which naturally matches DAC constraints for experimentation at Mbar pressure and minimizes the axial thermal gradients. As shown in Fig. 4, the sample thickness under pressure does not exceed 3–4 µm, granting temperature gradients below 200 K⁶², comparable to the error bars of spectroradiometric measurements.

At the target pressure, the sample was heated for about 1.1 up to 2.1 s at a selected laser's power. During this time the black body radiation emitted from the heated sample was collected using the beamline spectrometer and the temperature was measured by fitting the obtained black body radiation with a Planck function in the grey body approximation. Simultaneously, X-ray absorption spectra were collected on the FReLoN detector with an exposure time of about 100 ms per spectrum. The final XAS signal was obtained by averaging about 10–20 spectra to improve the signal to noise ratio. XAS spectra were collected at ambient temperature on the heated sample area after each heating pulse in order to monitor sample modifications and to identify potential chemical reactions. The reported error in the temperature corresponds to the standard deviation of the two-color fit as in Benedetti et al.⁶⁵. The pressure was measured before and after each heating run by ruby (Cr:Al₂O₃) fluorescence method, using the non-hydrostatic calibration by Mao et al.⁶⁶. The thermal pressure correction was applied to the data as in Boccato et al.¹². Overall, error bar on pressure is of ~ 10% of the pressure value.

More details about the measurement strategy and pressure and temperature metrology for these experiments are provided in Refs.¹²,¹³.

**Normalization of XANES spectra.** For the analysis of XAS data obtained in this work we first normalized the raw absorption spectra using two distinct energy points. The first point was located below the absorption edge (here referred to as Eb) and was set to zero. The second point was chosen after the absorption edge (here referred to as Ea) and was set to one. Eb is 7105.0 eV for the four systems, while Ea is 7142.8 eV for Fe–C, 7148.4 eV for Fe–O, 7146.1 eV for Fe–S and 7145.6 eV for Fe–Si, with a variability of about 0.1 eV due to the pixel to energy conversion in the different experiments. Eb and Ea are shown in Fig. 1 for the four systems. The choice of Ea, depends on the structure of the specific binary system. The two points have to be chosen such as their intensity is unaffected by the Debye Waller damping caused by the temperature increase.

**Analysis of recovered samples.** Several Fe–C and Fe–S samples were recovered for additional ex situ analysis after quenching from high pressure and temperature conditions. Fe–O and Fe–Si could not be recovered. A cross section was cut through the laser heated area using a focused ion beam (FIB) milling. Textural analysis was then performed using a Zeiss Neon40ESB SEM. The chemical composition determined by both electron microprobe and energy dispersive X-ray spectroscopy (EDX) elemental mapping in the SEM.

The carbon content of Fe–C recovered samples was determined using the Cameca SX100 electron microscope at CAMPARIS, Sorbonne Université, using pure Fe and stoichiometric Fe₃C as standard references. A cold catcher, cooled down with liquid nitrogen, and a micro-leak of oxygen were used to minimize the carbon contamination during the analysis. Operating conditions were 10 kV and 10 nA for a counting time of 45 s on peak and 20 s on background. To ensure electrical contact with the sample holder, the samples were coated with 3 nm of platinum (plasma coating for 10 s with a current of 60 mA). Carbon coating was avoided to prevent bias
in the carbon quantification. To ensure the stability of the calibration, pure iron was measured before and after each measurement on the sample.

The sulphur content in the Fe–S platinum coated samples was determined by EDX analysis on the SEM FEG Zeiss Ultra55 (IMPMC, Paris), using a 15 keV electron beam and measuring the emitted intensity on a silicon drift detector (SDD). The intensity calibration was performed on a copper reference sample. For quantitative analysis we used the Phi-Rho-Z method using pure iron and natural pyrite (FeS₂) as standard references. This method allowed to deconvolute the Fe and S peaks in natural pyrite and therefore provided an accurate quantification.

## Data availability

The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

Received: 3 March 2020; Accepted: 19 June 2020

Published online: 15 July 2020

## References

1. Boehler, R. et al. Portable laser-heating stand for synchrotron applications. *Rev. Sci. Instrum.* **80**, 045103 (2009).

2. Torchio, R. et al. Probing the local, electronic and magnetic structure of matter under extreme conditions of temperature and pressure. *High Press. Res.* **36**, 293–302 (2016).

3. Kantor, I., Marini, C., Mathon, O. & Pascarelli, S. A laser heating facility for energy-dispersive X-ray absorption spectroscopy. *Rev. Sci. Instrum.* **89**, 1–13 (2018).

4. Bunker, G. *Introduction to XAFS. A practical guide to X-ray absorption fine structure spectroscopy* (Cambridge University Press, Cambridge, 2010).

5. Fornasini, P. Introduction to X-ray absorption spectroscopy. in *XI School on Synchrotron Radiation*, 1–32 (2011).

6. Dorchies, F. et al. Unraveling the solid–liquid–vapor phase transition dynamics at the atomic level with ultrafast X-ray absorption near-edge spectroscopy. *Phys. Rev. Lett.* **107**, 245006 (2011).

7. Dorchies, F. & Recoules, V. Non-equilibrium solid-to-plasma transition dynamics using XANES diagnostic. *Phys. Rep.* **657**, 1–26 (2016).

8. Marini, C. et al. A microsecond time resolved x-ray absorption near edge structure synchrotron study of phase transitions in Fe undergoing ramp heating at high pressure. *J. Appl. Phys.* **115**, 1–5 (2014).

9. Harmand, M. et al. X-ray absorption spectroscopy of iron at multimegabar pressures in laser shock experiments. *Phys. Rev. B Condens. Matter Mater. Phys.* **92**, 1–7 (2015).

10. Aquilanti, G. et al. Melting of iron determined by X-ray absorption spectroscopy to 100 GPa. *Proc. Natl. Acad. Sci.* **112**, 12042–12045 (2015).

11. 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).

12. Boccato, S. et al. The melting curve of nickel up to 100 GPa explored by XAS. *J. Geophys. Res. Solid Earth* **122**, 1–10 (2017).

13. Morard, G. et al. Solving controversies on the iron phase diagram under high pressure. *Geophys. Res. Lett.* **45**, 11, 074–11, 082 (2018).

14. Boccato, S. et al. Compression of liquid Ni and Co under extreme conditions explored by x-ray absorption spectroscopy. *Phys. Rev. B* **100**, 180101. https://doi.org/10.1103/PhysRevB.100.180101 (2019).

15. Filipponi, A. et al. Single-energy x-ray absorption detection: a combined electronic and structural local probe for phase transitions in condensed matter. *J. Phys. Condens. Matter* **10**, 235–253 (1998).

16. Raji, A. T., Scandolo, S., Härting, M. & Britton, D. T. Probing the structure of iron at extreme conditions by X-ray absorption near-edge structure calculations. *High Press. Res.* **33**, 119–123 (2013).

17. Raji, A. & Scandolo, S. Theoretical x-ray absorption near-edge structure signatures of solid and liquid phases of iron at extreme conditions. *High Press. Res.* **34**, 250–258 (2014).

18. Mazevet, S. et al. Ab initio calculation of x-ray absorption of iron up to 3 Mbar and 8000 K. *Phys. Rev. B Condens. Matter Mater. Phys.* **89**, 1–5 (2014).

19. Yildirim, C. et al. Universal amorphous-amorphous transition in GexSe100-x glasses under pressure. *Sci. Rep.* **6**, 27317 (2016).

20. Benfatto, M. & Meneghini, C. *A close look into the low energy region of the XAS spectra: the XANES region* (Springer, Berlin, 2015).

21. Breuer, D., Rueckriemen, T. & Spohn, T. Iron snow, crystal floats, and inner-core growth: modes of core solidification and implications for dynamos in terrestrial planets and moons. *Prog. Earth Planet. Sci.* **2**, 39. https://doi.org/10.1186/s40645-015-0069-y (2015).

22. Filipponi, A. et al. An experimental station for advanced research on condensed matter under extreme conditions at the European Synchrotron Radiation Facility-BM29 beamline. *Rev. Sci. Instrum.* **71**, 2422–2432 (2000).

23. Di Cicco, A. & Trapananti, A. Study of local icosahedral ordering in liquid and undercooled liquid copper. *J. Non-Cryst. Solids* **353**, 3671–3678 (2007).

24. Di Cicco, A. et al. Local fivefold symmetry in liquid and undercooled ni probed by x-ray absorption spectroscopy and computer simulations. *Phys. Rev. B* **89**, 060102 (2014).

25. Boccato, S. Local Structure of Liquid 3d Metals Under Extreme Conditions of Pressure and Temperature. Ph.D. thesis, Université Grenoble Alpes (2017).

26. 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–6 (2013).

27. Torchio, R., Boccato, S., Miozzi, F., Rosa, A. D., Ishimatsu, N., Kantor, I. et al. Melting curve and phase relations of Fe-Ni alloys: implications for the Earth's core composition. *Geophys. Res. Lett.* **47**, e2020GL088169. https://doi.org/10.1029/2020GL088169 (2020).

28. Morard, G. et al. Fe-FeO and Fe-Fe3C melting relations at Earth's core-mantle boundary conditions: implications for a volatile-rich or oxygen-rich core. *Earth Planet. Sci. Lett.* **473**, 94–103 (2017).

29. Yuan, K. et al. Synergetic contribution of boron and F-Nx species in porous carbons toward efficient electrocatalysts for oxygen reduction reaction. *ACS Energy Lett.* **3**, 252–260 (2018).

30. Ravel, B. & Newville, M. ATHENA, ARTEMIS, HEPHAESTUS: data analysis for X-ray absorption spectroscopy using IFEFFIT. *J. Synchrotron Radiat.* **12**, 537–541. https://doi.org/10.1107/S0909049505012719 (2005).

31. Boulard, E. et al. Ferrous iron under oxygen-rich conditions in the deep mantle. *Geophys. Res. Lett.* **46**, 1348–1356 (2019).

32. Fei, Y., Li, J., Bertka, C. M. & Prewitt, C. T. Structure type and bulk modulus of Fe3S, a new iron-sulfur compound. *Am. Mineral.* **85**, 1830–1833 (2000).

33. Kuwayama, Y. & Hirose, K. Phase relations in the system FeFeSi at 21 GPa. *Am. Mineral.* **89**, 273–276 (2004).
34. Stewart, A. J., Schmidt, M. W., van Westrenen, W. & Liebske, C. Mars: a new core-crystallization regime. *Science* **316**, 1323–1325 (2007).

35. Boehler, R. Melting of the FeFeO and the FeFeS systems at high pressure: constraints on core temperatures. *Earth Planet. Sci. Lett.* **111**, 217–227 (1992).

36. Lord, O., Walter, M., Dasgupta, R., Walker, D. & Clark, S. Melting in the Fe–C system up to 70 GPa. *Earth Planet. Sci. Lett.* **284**, 157–167 (2009).

37. Seagle, C. T., Heinz, D. L., Campbell, A. J., Prakapenka, V. B. & Wanless, S. T. Melting and thermal expansion in the Fe–FeO system at high pressure. *Earth Planet. Sci. Lett.* **265**, 655–665 (2008).

38. Morard, G. et al. Melting of Fe Ni Si and Fe Ni S alloys at megabar pressures: implications for the core mantle boundary temperature. *Phys. Chem. Miner.* **38**, 767–776. https://doi.org/10.1007/s00269-011-0449-9 (2011).

39. Liu, J., Li, J., Hrubiak, R. & Smith, J. S. Origins of ultralow velocity zones through slab-derived metallic melt. *Proc. Natl. Acad. Sci.* **113**, 5547–5551 (2016).

40. Mori, Y. et al. Melting experiments on Fe–Fe3S system to 254 GPa. *Earth Planet. Sci. Lett.* **464**, 135–141 (2017).

41. Kamada, S. et al. Phase relationships of the Fe–FeS system in conditions up to the Earth’s outer core. *Earth Planet. Sci. Lett.* **294**, 94–100 (2010).

42. Kamada, S. et al. Melting relationships in the Fe–Fe3S system up to the outer core conditions. *Earth Planet. Sci. Lett.* **359–360**, 26–33 (2012).

43. Fischer, R. A. et al. Phase relations in the Fe–FeSi system at high pressures and temperatures. *Earth Planet. Sci. Lett.* **373**, 54–64 (2013).

44. Mashino, I., Miozzi, F., Hirose, K., Morard, G. & Sinmyo, R. Melting experiments on the Fe3C binary system up to 255 GPa: con-straints on the carbon content in the Earth's core. *Earth Planet. Sci. Lett.* **515**, 135–144. https://doi.org/10.1016/j.epsl.2019.03.020 (2019).

45. Wood, B. J. Carbon in the core. *Earth Planet. Sci. Lett.* **117**, 593–607. https://doi.org/10.1002/9781118992487.ch22 (1993).

46. Sinmyo, R. & Hirose, K. The Soret diffusion in laser-heated diamond-anvil cell. *Phys. Earth Planet. Inter.* **180**, 172–178 (2010).

47. Oka, K. et al. Melting in the Fe–FeO system to 204 GPa: implications for oxygen in Earth’s core. *Am. Miner.* **104**, 1603–1607. https://doi.org/10.2138/am-2019-7081 (2019).

48. Komabayashi, T. Thermodynamics of melting relations in the system Fe-FeO at high pressure: Implications for oxygen in the Earth's core. *J. Geophys. Res. Solid Earth* **119**, 4164–4177 (2014).

49. Tsuno, K., Ohtani, E. & Terasaki, H. Immiscible two-liquid regions in the Fe–O–S system at high pressure: implications for planetary cores. *Phys. Earth Planet. Inter.* **160**, 75–85 (2007).

50. Fei, Y., Bertka, C. M. & Finger, L. W. High-pressure iron–sulfur compound, Fe3S2, and melting relations in the Fe–FeS system. *Science* **275**, 1621–1623 (1997).

51. Chudinovskikh, L. & Boehler, R. Eutectic melting in the system Fe–S to 44 Å GPa. *Earth Planet. Sci. Lett.* **257**, 97–103 (2007).

52. Morard, G. et al. In situ determination of Fe–Fe3S phase diagram and liquid structural properties up to 65 GPa. *Earth Planet. Sci. Lett.* **272**, 620–626 (2008).

53. Morard, G. et al. Structure of eutectic Fe–FeS melts to pressures up to 17 GPa: implications for planetary cores. *Earth Planet. Sci. Lett.* **263**, 128–139. https://doi.org/10.1016/j.epsl.2007.09.009 (2007).

54. Fischer, R. A. et al. High pressure metal-silicate partitioning of Ni, Co, V, Cr, Si, and O. *Geochim. Cosmochim. Acta* **167**, 177–194. https://doi.org/10.1016/j.gca.2015.06.026 (2015).

55. Kuwayama, Y., Sawai, T., Hirose, K., Sata, N. & Ohishi, Y. Phase relations of iron–silicon alloys at high pressure and high temperature. *Phys. Chem. Miner.* **9**, 511–518. https://doi.org/10.1007/s00269-009-0296-0 (2009).

56. Fei, Y. & Brosh, E. Experimental study and thermodynamic calculations of phase relations in the Fe–C system at high pressure. *Earth Planet. Sci. Lett.* **408**, 155–162 (2014).

57. Andrault, D. et al. ESRF Upgrade Programme Phase II (2015–2022) Technical Design Study (2014). [Online; Accessed 04 May 2020].

58. Rosa, A. D. et al. Nano-polycrystalline diamond anvils: key devices for xas at extreme conditions: their use, scientific impact, present status and future needs. *High Press. Res.* https://doi.org/10.1080/08957959.2019.1700978 (2019).

59. Timoshenko, J., Lu, D., Lin, Y. & Frenkel, A. I. Supervised machine-learning-based determination of three-dimensional structure of metallic nanoparticles. *J. Phys. Chem. Lett.* **8**, 5091–5098 (2017).

60. Timoshenko, J. et al. Neural network approach for characterizing structural transformations by X-ray absorption fine structure spectroscopy. *Phys. Rev. Lett.* **120**, 225502 (2018).

61. Aprilis, G. et al. Comparative study of the influence of pulsed and continuous wave laser heating on the mobilization of carbon and its chemical reaction with iron in a diamond anvil cell. *J. Appl. Phys.* **125**, 095901. https://doi.org/10.1063/1.5067268 (2019).

62. Sinmyo, R., Hirose, K. & Ohishi, Y. Melting curve of iron to 290 GPa determined in a resistance-heated diamond-anvil cell. *Earth Planet. Sci. Lett.* **510**, 45–52 (2019).

63. Kantor, I. et al. A new detector for sub-millisecond EXAFS spectroscopy at the European Synchrotron Radiation Facility. *J. Synchrotron Radiat.* **21**, 1240–1246 (2014).

64. Ishimatsu, N. et al. Glitch-free X-ray absorption spectrum under high pressure obtained using nano-polycrystalline diamond anvils. *J. Synchrotron Radiat.* **19**, 768–772 (2012).

65. Benedetti, L. & Loubeyre, P. Temperature gradients, wavelength-dependent emissivity, and accuracy of high and very-high temperatures measured in the laser-heated diamond cell. *High Press. Res.* **24**, 423–445 (2004).

66. Mao, H. K. Static compression of simple molecular systems in the megabar range. In *Simple molecular systems at very high density* (eds Polian, P. L. A. & Boccara, N.) 221–236 (Plenum Press, New York, 1989).

## Acknowledgements

The authors are thankful to Jeroen Jacobs for technical support in the diamond anvil cell preparation and to Sakura Pascarelli, Nicolas Sévelin-Radiguet, Olivier Mathon, Florian Perrin et Sebastien Pasternak for support on the beamline. The authors wish to thank also Imène Estève for her help with sample analysis by SEM, and Michel Fialin and Nicolas Rividi for their help during microprobe analysis. The authors acknowledge the European Synchrotron Radiation Facility (ESRF) for provision of synchrotron radiation within Projects HC1307, ES358, HC2578, ES778. Femtosecond laser micromachining at the Institut de Minéralogie de Physique des Matériaux et de Cosmochimie (IMPMC), Paris, has been developed and realized by the “Cellule Project” with the financial support of ANR 2010-JCJC-604-01. The Focused Ion Beam (FIB) and Scanning Electron Micro-scope (SEM) facility at IMPMC is supported by Region Ile de France Grant SESAME 2006 N I-07-593/R, INSU-CNRS, Institute de Physique (INP)–CNRS, University Pierre et Marie Curie–Paris 6, and by the French National Research Agency (ANR) Grant ANR-07-BLAN-0124-01. This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (Grant
Agreement No. 724690). F. M. has received funding from the European Research Council (ERC) under the
European Union's Horizon 2020 research and innovation programme (Grant Agreement No. 670787).

**Author contributions**

G.M. conceived the experiments, G.M., S.B., R.T., S.A., E.B., F.G., M.H., I.K., F.M., P.P. and A.D.R. participated to the experiment(s), T.I. provided nanopolycrystalline diamonds, S.B. and G.M. performed ex situ compositional analysis, S.B. analyzed the XAS data, S.B., G.M., R.T. and D.A. worked on data interpretation and discussion of the results. 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-68244-3.

**Correspondence** and requests for materials should be addressed to S.B.

**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_10.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