# Ultrafast X-ray Diffraction Study of a Shock-Compressed Iron Meteorite above 100 GPa

Sabrina Tecklenburg, Roberto Colina-Ruiz, Sovanndara Hok, Cynthia Bolme, Eric Galtier, Eduardo Granados, Akel Hashim, Hae Ja Lee, Sébastien Merkel, Benjamin Morrow, et al.

► To cite this version:

Sabrina Tecklenburg, Roberto Colina-Ruiz, Sovanndara Hok, Cynthia Bolme, Eric Galtier, et al. Ultrafast X-ray Diffraction Study of a Shock-Compressed Iron Meteorite above 100 GPa. Minerals, 2021, 11 (6), pp.567. 10.3390/min11060567. hal-03256610

HAL Id: hal-03256610

https://hal.univ-lille.fr/hal-03256610

Submitted on 10 Jun 2021

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
Article

# Ultrafast X-ray Diffraction Study of a Shock-Compressed Iron Meteorite above 100 GPa

Sabrina Tecklenburg<fnref n="1" />,*, Roberto Colina-Ruiz <fnref n="2" />,†, Sovanndara Hok <fnref n="1" />, Cynthia Bolme <fnref n="3" />, Eric Galtier <fnref n="2" />,*,
Eduardo Granados <fnref n="2" />,†,*, Akel Hashim <fnref n="2" />,†,*, Hae Ja Lee <fnref n="2" />, Sébastien Merkel <fnref n="4" />,*, Benjamin Morrow <fnref n="3" />,*,
Bob Nagler <fnref n="2" />, Kyle Ramos <fnref n="3" />, Dylan Rittman <fnref n="1" />, Richard Walroth <fnref n="2" />, Wendy L. Mao <fnref n="1" /> and Arianna E. Gleason <fnref n="1" />,<fnref n="2" />,*

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

check for
updates

**Citation:** Tecklenburg, S.;
Colina-Ruiz, R.; Hok, S.; Bolme, C.;
Galtier, E.; Granados, E.; Hashim, A.;
Lee, H.J.; Merkel, S.; Morrow, B.; et al.
Ultrafast X-ray Diffraction Study of a
Shock-Compressed Iron Meteorite
above 100 GPa. *Minerals* **2021**, *11*, 567.
https://doi.org/10.3390/min11060567

Academic Editor: Martin Kunz

Received: 13 April 2021

Accepted: 24 May 2021

Published: 26 May 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Copyright:** © 2021 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** Natural kamacite samples (Fe₉₂.₅Ni₇.₅) from a fragment of the Gibeon meteorite were studied as a proxy material for terrestrial cores to examine phase transition kinetics under shock compression for a range of different pressures up to 140 GPa. In situ time-resolved X-ray diffraction (XRD) data were collected of a body-centered cubic (bcc) kamacite section that transforms to the high-pressure hexagonal close-packed (hcp) phase with sub-nanosecond temporal resolution. The coarse-grained crystal of kamacite rapidly transformed to highly oriented crystallites of the hcp phase at maximum compression. The hcp phase persisted for as long as 9.5 ns following shock release. Comparing the *c/a* ratio with previous static and dynamic work on Fe and Fe-rich Fe-Ni alloys, it was found that some shots exhibit a larger than ideal *c/a* ratio, up to nearly 1.65. This work represents the first time-resolved laser shock compression structural study of a natural iron meteorite, relevant for understanding the dynamic material properties of metallic planetary bodies during impact events and Earth's core elasticity.

**Keywords:** ultrafast X-ray diffraction; laser shock compression; iron meteorite

## 1. Introduction

Iron meteorites provide unique samples of early solar system material and offer
insight into planetary processes such as impact events and core formation originating
from differentiated planetesimals bodies [1,2]. The Gibeon meteorite, an iron meteorite
that was discovered in 1836 after falling in Namibia, Africa [3], is studied here due to
its compositional similarity to the Earth's core and commercially availability. It can be
used for studying how the physical properties of planetary core materials are affected by
dynamic processes that may have occurred during the formation and evolution—much
like the metallic asteroid 16-Psyeche [4]. During cosmic hit-and-run collisions, planetesimals
are stripped of their mantle material and leave behind a residual metal core that later
could break up further through collisions with other planetary bodies [5]. In addition,
exploring the effects of shock metamorphism on iron meteorites, as characterized by mi-
crostructure and microchemistry studies, provides insight into impact events and thermal
history [6]. Understanding the structural properties of meteorites under shock compression

<footnote n="1">¹ Department of Geological Sciences, Stanford University, Stanford, CA 94305, USA; stecklen@stanford.edu (S.T.); darahok@stanford.edu (S.H.); drittman13@gmail.com (D.R.); wmao@stanford.edu (W.L.M.)</footnote>

<footnote n="2">² SLAC National Accelerator Laboratory, Menlo Park, CA 94025, USA; colina@slac.stanford.edu (R.C.-R.); egaltier@slac.stanford.edu (E.G.); eduardo.granados@cern.ch (E.G.); ahashim@berkeley.edu (A.H.); haelee@slac.stanford.edu (H.J.L.); bnagler@slac.stanford.edu (B.N.); walroth@slac.stanford.edu (R.W.)</footnote>

<footnote n="3">³ Los Alamos National Laboratory, Los Alamos, NM 87545, USA; cbolme@lanl.gov (C.B.); morrow@lanl.gov (B.M.); kramos@lanl.gov (K.R.)</footnote>

<footnote n="4">⁴ Unité Matériaux et Transformations Cité scientifique, CNRS, INRAE, Centrale Lille, UMR 8207—UMET, F-59000 Lille, France; sebastien.merkel@univ-lille.fr</footnote>

<footnote>* Correspondence: ariannag@stanford.edu</footnote>

<footnote>† Now at CERN, 1211 Geneva 23, Switzerland.</footnote>

<footnote>‡ Now at Physics Department, University of California, Berkeley, CA 94704, USA.</footnote>
and during the timescale of dynamic compression may provide insight into the evolution
of these asteroids.

The Gibeon meteorite is made of an Fe-rich Fe-Ni alloy and belongs to the most
common structural class for iron meteorites called octahedrites (group IVA). This meteorite
contains three major phases, depending on the Fe/Ni ratio in its composition: an Fe-
rich/Ni-poor ferrite phase (<7.5% Ni) known as kamacite, which has a body-centered cubic
(bcc) structure; a second major phase taenite, also called austenite, which can incorporate
up to 25% Ni, crystallizing in a face-centered cubic (fcc) structure [7,8]; and plessite,
the third major meteorite texture, which is a mixture of kamacite and taenite, exhibiting
fine grains between larger bands of pure taenite and kamacite phases. The kamacite phase
is located in the Widmanstätten lamellae [9], with plessite and taenite making up the
material between the lamellae, producing complex structural orientations and boundaries
between the bcc and fcc crystals [10]. This is the dominant microstructure seen in Gibeon
meteorite samples [11].

The phase diagram of Fe has been studied extensively as a function of pressure and temperature and is of particular importance for Earth and planetary science as well as for industry. It shows transitions between the bcc or fcc structures transforming to a hexagonal close-packed (hcp) structure at high pressure, (e.g., [12–15]). Some reports include investigations on the Fe-Ni binary system, demonstrating that the fcc–hcp boundary is shifted toward higher pressure and lower temperature with respect to the pure Fe phase diagram. Here, the Ni stabilizes the fcc structure, showing that the Ni content could be crucial in stabilizing the two-phase region (fcc together with hcp) at extreme conditions [16].

The ratio of hcp crystal structure lattice parameters, *c/a*, is often tracked as a function of pressure (*P*) or volume (*V*) and temperature (*T*) to provide insight into atomic packing and examine deformation mechanisms in metals. At Earth's inner core pressure-temperature conditions, hcp iron is likely the stable phase. However, despite the importance of understanding the inner core structure and properties, the elastic constants of hcp iron at relevant core conditions are still poorly constrained. The ideal *c/a* ratio of 1.633 is derived from the assumption of non-interacting, identical hard spheres in an hcp structure. One way to track the elasticity of hcp iron at extreme conditions, and in particular, elastic anisotropy, is to measure departure from the ideal *c/a* ratio. For instance, alignment of hcp crystallites could be responsible for the measured seismic anisotropy in Earth's inner core. Substantially larger *c/a* ratios of 1.7, for instance, have been predicted for hcp Fe at inner core conditions [17], although static experiments do not agree with those predictions [18]. The axial *c/a* ratio can be used to characterize hcp Fe elasticity under static [17] or dynamic [19] compression. Similar to static compression measurements, dynamic compression experiments and simulations on polycrystalline bcc Fe show shock-induced transitions to hcp Fe and also reveal a *c/a* ratio below the ideal value [19]. This could be associated with defects, such as dislocations, or grain boundaries that act as sources of stress on the time scale of the measurements (i.e., nanoseconds) [19]. Furthermore, it has been demonstrated that hcp Fe-Ni and Fe-Ni-Si alloys only have slightly higher *c/a* ratios (less than 1% higher) than pure hcp Fe for the same pressure range [20] but that are all still lower than the ideal value. Tracking this parameter could provide insight into which important underlying physical or chemical properties (e.g., starting structure, composition, and microstructure like crystal orientation or texture) and/or processes (e.g., static vs. dynamic compression, including strain rate) dictate the deformation response in meteoritic material—including the typical timescales of phase transformation of meteorites exposed to extreme conditions.

High-pressure dynamic compression experiments represent a unique way to recreate
dynamic processes such as asteroid/planetary collisions in the lab and make it possible to
track the lattice-level response of materials using in situ time-resolved X-ray diffraction [21,22].
Dynamic compression techniques also provide access to combined high-pressure and high-
temperature states of a material in concert with characterization via ultrafast X-ray probes
that more traditional static compression techniques cannot access, e.g., [23]. Coupling the
laser-induced shock compression experiments with reliable structural studies to elucidate
phase transitions that occur at short timescales of hundreds of picoseconds to nanoseconds requires an ultrafast and tunable X-ray source such as X-ray free-electron lasers (XFELs) that can operate on femtosecond timescales with a high repetition rate, e.g., [24–26].

Herein, the first time-resolved structural study of a laser-shock-compressed natural meteorite is presented. Phase transformation kinetics of Fe<sub>92.5</sub>Ni<sub>7.5</sub> natural kamacite was investigated from the Gibeon meteorite using in situ time-resolved X-ray diffraction and compared to previous static compression work on related systems.

**2. Materials and Methods**

**2.1. *Sample Preparation and Characterization***

The starting material was a hand sample fragment (0.4 cm × 2 cm × 3 cm) of the Gibeon meteorite (Figure 1a) that was characterized by an electron microprobe to determine the precise metallic composition before shock compression (Figure 1b). Electron backscatter diffraction (EBSD, Philips XL30 FEG SEM using TSL Data Collection & Analysis software, MST-8, LANL, Los Alamos, NM, USA) was also performed to determine crystal orientations (Supplementary Materials, Figure S1). Individual, large-grained crystals of kamacite (bcc) are clear in the Widmanstätten lamellae, with Ni-rich taenite (fcc) and plessite concentrated in the grain boundaries. The hand sample fragment was sectioned and double-side parallel-polished into 34-µm-thick foils. Here, kamacite was isolated for primary investigation of this study due to its similarity to Earth's core composition. Individual kamacite lamellae were cut along the perimeter and diced into 1 mm × 2 mm foils. The composition within the selected lamellae was Fe<sub>92.5</sub>Ni<sub>7.5</sub>, with taenite and plessite (Ni content up to 25%) concentrated in the grain boundaries. The samples were aligned such that the length of the lamellae was perpendicular to the shock propagation direction. Dupont Kapton CB plastic film, 54 µm thick, was glued (glue layer thickness: 2–5 µm) to the meteorite foil, serving as the laser ablation layer. This plastic smooths out any inhomogeneities (e.g., hot spots) in the drive laser to generate a planar, uniform shock front with respect to the X-ray probed volume.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** (a) Photomicrograph of Gibeon meteorite starting material. (b) Electron microprobe map showing the Ni distribution (in wt %) in a typical portion of the sample. On average, the Ni content within the grains is 7.5%, while the grain boundaries are Ni rich (up to 25%).</figcaption></figure>

**2.2. *Experimental Setup***

A schematic of the setup for the dynamic compression experiments carried out at the Matter at Extreme Conditions (MEC) end station is shown in Figure 2. The XFEL probe was generated by the Linac Coherent Light Source (LCLS) at the SLAC National Accelerator Laboratory [27] through a quasi-monochromatic beam with high-energy reso-lution (ΔE/E~0.2%), fully transverse coherent at 6.9 keV, with a 60 fs temporal duration and an average of 10¹² photons per pulse. The XFEL probe beam was 20 µm in diameter on the sample. Phase plates on the optical drive laser, a frequency-doubled Nd: Glass laser
system (527 nm), set a 176-µm-diameter laser spot on target. This was used to achieve a focal spot intensity of ~10¹² W/cm² over a 12 ns quasi-flat top temporal pulse profile. The energy of the drive laser was varied from 6 to 40 J to achieve different peak pressures in the sample. The angle between the target normal and the XFEL probe was 65°. The optical laser and X-ray beam were spatially overlapped and operated in single-shot mode with absolute time zero corresponding to the overlap of their leading edges of the pulse. Prior to each laser shot, an X-ray-only pattern was first collected to confirm the starting sample area corresponding to the expected bcc structure. For each shot, a time delay (Δt) was selected for the XFEL pulse relative to the optical laser pulse with a jitter of 0.5 ns. This Δt was verified by oscilloscope traces captured for each shot. The pump-probe delay scans at several nanosecond intervals enabled collection of X-ray diffraction (XRD) time series, where the 2-dimensional (2D) patterns were captured by Cornell-SLAC Pixel Array Detectors (CSPADs) constructed of individual application-specific integrated circuits (ASICs) [28]. The combined azimuthal angle coverage was 180°. One target was shot per selected Δt.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Experimental setup at the Matter in Extreme Conditions (MEC) end station with the Linac Coherent Light Source (LCLS), an X-ray free-electron laser (XFEL) at the SLAC National Accelerator Laboratory. An MEC long-pulse laser was used to generate and launch a shock wave, via laser ablation, in the plastic (labeled 'CH' in the schematic), which then traversed the meteorite foil. At discrete time delays (Δt) during the shockwave propagation, the XFEL was used to probe the lattice structure by measuring XRD in transmission. Inset: cross section of the sample package. Meteorite sample thickness was 34 ± 0.5 µm.</figcaption></figure>

# 3. Results

## 3.1. X-ray Diffraction

Calibration for the sample-to-detector distance and orientation of the CSPADs were performed using internally developed software [29] and Dioptas [30]. Peak fitting was performed using Fityk [31]. An image of raw data from the CSPADs projected onto a single plane is shown in Figure 3a. The XRD patterns of the shock-compressed material showed new peaks identified as (100), (002), and (101) for high-pressure hcp Fe [18,32] (Figure 3b). This hcp Fe manifested as large spots with uneven intensity distribution in the XRD about the azimuthal coverage. A possible explanation for this diffraction character is highly oriented crystallites of hcp Fe, rather than a few large-grained crystals, since previous shock compression work on Fe demonstrated nano-crystalline hcp Fe formation at ul-
trafast timescales [13]. Due to the large angle between the shock propagation direction and the incident X-ray beam, some small volume (less than a few µm³) of unshocked sample was recorded in the data, identified as ambient-pressure bcc Fe. For each shot, the hcp (100), (002), and (101) peak positions were used to determine the lattice parameters and density (ρ), adjusted to include 7.5% Ni. A compositionally close Fe₉₀Ni₁₀ Hugoniot [33] was used to determine pressure for the measured density for in-shock data. Data collected after the cessation of the drive laser, and/or when the compressive wave reached the edge of the sample, were on release and no longer represented an on-Hugoniot measurement. Here, the quasi-isentropic release path from our peak pressure was estimated by extrapolating pressure and temperature from previously reported post-shock studies on iron [34,35]; see Supplemental Discussion, Section 1. Using the measured density from our XRD, we used these extrapolated release paths to constrain the pressure. Pressure uncertainty was derived from standard error propagation from XRD peak fitting. Velocimetry records were also collected corroborating pressure estimates within 8% using the velocity interferometry system for any reflector (VISAR) diagnostic. However, only transit times were available due to the low reflectivities of the natural samples. LiF as a VISAR window material was not used. Representative VISAR data are shown in Supplementary Materials, Figure S2.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Representative raw X-ray diffraction data for a sample shocked to 42 GPa. (a) 2D diffraction image from CSPADs projected onto a single plane. (b) Caked pattern of (a) displayed as an s function of the 2θ and azimuthal angle. The blue diamonds and dashed lines show the hcp peaks for the shock-compressed high-pressure phase. The red line and triangle indicate uncompressed starting bcc Fe. The large spots and uneven intensity distribution indicate the sample is coarse-grained and highly textured.</figcaption></figure>

Two types of XRD data series were collected, a pressure series and a time series on shock release from a peak pressure. A pressure series consisted of varied drive laser parameters to achieve a suite of peak pressure states, ~40 GPa to above 140 GPa, on the Hugoniot where the sample was probed at the state of maximum compressed sample volume, ~6.5 ns after onset of compression, for our target design. The 2D XRD data were recorded for each pressure. After a steady shock was no longer supported (i.e., the drive laser turned off or the shock front reached the end of the sample at the vacuum interface), the sample under-
went quasi-isentropic release, during which time the pressure and temperature decreased from some peak state. XRD patterns were collected at time delays of >6.5 ns to examine the metastability kinetics of hcp Fe. In this regime, though not totally released back to the starting bcc structure, the sample was still within the hcp phase stability region.

Spotty, high-intensity features in the 2D XRD made the hcp peaks challenging to identify. The hcp Fe transformation agreed with previous static compression experiments on pure Fe and Fe-rich Fe-Ni alloys up to 30% [36]. As expected, when Ni was added, the unit cell contracted to accommodate the slightly smaller Ni atoms (~1.24 Å) compared to the Fe atoms (~1.26 Å) [37]. As the drive laser energy increased, a suite of high-pressure states on the Hugoniot could be achieved. On compression, the hcp Fe peak positions shifted to a higher 2θ (or smaller d-spacing), as seen in azimuthally integrated XRD patterns at pressures ~40 GPa to above 140 GPa (Figure 4a).

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Integrated XRD patterns at 6.9 keV. The blue diamonds show hcp peak locations: light blue, mid-tone blue, and dark blue indicate hcp (100), hcp (002), and hcp (101) peaks, respectively. The red diamonds show the unshocked bcc (110) peak from the starting material at ambient pressure. (a) Pressure series for separate runs collected at full compression. Inset: zoomed-in 2θ region to illustrate a weak S/N of the hcp (100) and (200) peaks at 85 and 113 GPa, relative to the more intense (101) peak. (b) Representative time series showing XRD data collected for the same incident laser intensity, generating a 113 GPa peak pressure, but different XFEL time delays: at nearly full compression (bottom trace) and on release (top trace) at ~25 GPa.</figcaption></figure>

XRD data were also collected during shock release. After cessation of the drive laser and/or reaching the edge of the sample, peak compression was no longer supported and lower pressures and temperature could be accessed. Successively longer XFEL time delays probed the lattice structure as it reverted from the high-pressure phase back to the low-pressure state. The material eventually reached ambient conditions. An example of this shock release XRD data (Figure 4b) still showed retention of the hcp Fe phase at a pressure lower than the peak compression state.

## 3.2. Axial c/a Ratio

After identifying the hcp Fe peaks in the XRD for each run, least-squares fitting was performed to determine the hexagonal lattice parameters, *a* and *c*. Studying the *c/a* ratio as a function of *P*, *ρ*, and *T* is a common approach to track changes in elasticity of hexagonal metals at extreme conditions. As such, the axial *c/a* ratio was calculated and plotted
against volume to compare to previous static work on Fe-Ni alloys and previous dynamic compression of single- and poly-crystalline Fe (Figure 5); see Supplemental Discussion, Section 2. Previous static compression experiments on pure hcp Fe and Fe₉₀Ni₁₀ studied the effects of pressure, temperature, and composition on the axial ratio, in particular revealing an increase in *c/a* with increasing volume [38]. Most notably, these previous studies showed that the *c/a* ratio stays below the ideal 1.63 and is, on average, 1.61. Similarly, a shock compression study on polycrystalline Fe found that the axial ratio stays at 1.61 ± 0.01 regardless of pressure (up to 80 GPa) [19]. However, dynamic compression of single-crystal bcc Fe shocked along the [100] axis could produce a *c/a* that exceeds 1.7. Molecular dynamics simulations [39] show that in a compression regime imposing uniaxial strain, spacing between the atoms in the direction orthogonal to the shock front cannot change, and there can be a lack of plasticity from generation and propagation of dislocations. Specifically, the high strain rate and short timescale of these experiments inhibit defect and dislocation motion, such that the interplanar spacing of the (011) planes, which are perpendicular to the shock propagation direction, are completely unchanged during compression. In this case, the [100] of a bcc lattice remains less compressed and larger than expected as it becomes the *c*-axis in an hcp lattice and makes it possible to get an apparently large *c/a* ratio [40]. Our experimental data show a distribution in the *c/a* ratio—some instances with larger than, similar to, or less than the ideal ratio. The XRD data giving larger than ideal *c/a* ratios are cases where the *c*-axis is 2–3% larger and the *a*-axis is comparable to values of the same pressure under quasi-hydrostatic compression in a diamond-anvil cell [18]. This large *c*-axis value is consistent with model predictions, e.g., [39] and is the first experimental demonstration of a *c/a* ratio > 1.63.

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

**Figure 5.** Comparison of the axial *c/a* ratios determined in this study (black circles) to the ideal *c/a* ratio of 1.633 (dashed black line) and to previous static and dynamic compression work on Fe and Fe-rich Fe-Ni alloys [17]. Solid, colored lines are plotted to show the trend line for different groups from individual datasets from [17]. The shaded polygons show the spread of the data: shock-compressed polycrystalline Fe (gray) from [19], static-compressed polycrystalline Fe (pink) and Fe₉₀Ni₁₀ (blue) from datasets compiled in [38], pure Fe from [14,41–48], and Fe-Ni alloys from [36,47,49,50]. We plotted the static-compressed Fe and Fe₉₀Ni₁₀ data that fell in a similar temperature region as our study, 1000–3000 K. Error bars were calculated by propagating the uncertainty of the least-squares fit for hcp lattice parameters. The error bars on the data at 5.12 and 5.35 cm³ are smaller than the symbol size on the plot.
4. Discussion

Results from this work were compared with previous data within the same tem-
perature range as this study (1000–3000 K). Our c/a data show a large distribution as a
function of volume, with some values near previous work, e.g., [17] or near the ideal 1.63
and even above or below previous trends by as much as 3% (Figure 5). Though the data
distribution from all previous static experiments falls below the ideal c/a ratio, we present
the first experimental observation of a high c/a, at 1.65, which is even higher than pre-
vious ratios determined during shock compression on polycrystalline Fe. This suggests
that the requirements for achieving a higher than ideal c/a ratio include starting with an
oriented single crystal or highly textured polycrystalline bcc sample and using dynamic
compression to transform the bcc to the hcp phase. Though the starting orientation of
the exsolution lamellae remained constant shot to shot, meaning the long edge of the
lamellae was always placed parallel to the top of the sample holder, specific control of
intra-grain orientation for these natural large-grained samples was not possible. This could
explain the wide distribution in the c/a values for the Gibeon meteorite—some crystals
may have been oriented with the bcc [100] parallel to the shock propagation direction, as in
Run, reaching c/a = 1.65, whereas other samples may have had grains oriented along the
body diagonal or at high angles from the [100] direction. Recent work [51] systematically
tracked the starting orientation in single-crystal Fe during laser-driven shock compres-
sion, notably documenting higher than ideal c/a ratios as well. Continued effort in this
area is needed to potentially reveal how crystalline orientation trends support proposed
deformation mechanisms [13,39], as previously reported for Fe single crystals [52].

The pure Fe Hugoniot plotted on an equilibrium P–T phase diagram was used to compare our data with previous work, including the phase boundaries for pure Fe and the two-phase (fcc + hcp) region from a previous study on Fe₉₀.₃Ni₉.₇ [36] (Figure 6). The hcp–fcc boundary on the pure Fe phase diagram shifts to higher pressures and lower temperatures with increasing Ni content. We examined the composition dependence, ranging from pure Fe to Fe₉₀Ni₁₀, for the Hugoniot distribution across P–ρ and P–V, finding little change in our measured pressure regime, i.e., scatter in the datasets for each Ni content composition was within the scatter for pure Fe Hugoniot data [32,33]. We did not measure temperature directly and therefore used the Fe Hugoniot to estimate temperature. For XRD data collected after peak compression, on release ~7–16 ns, we estimated isentropic release paths using previously reported post-shock temperatures [34,35]. This helped to identify where, on the phase diagram, we were probing the lattice structure.

The hcp structure persisted to peak pressures on the Hugoniot, as expected. Interestingly, the hcp structure also remained in existence well past peak compression states to time delays as long as 16 ns. Here, samples follow an isentropic release path, similar to a natural impact or large-body collision process, with relevance to Psyche. However, it is important to consider that the starting conditions of our sample may be admittedly different than a natural iron meteorite in space or a metallic asteroid in our solar system at different points in the solar system evolution. For instance, starting temperatures can differ—our samples started at 300 K in this experiment, whereas initial temperatures of asteroids (e.g., Psyche) could have been substantially higher in the early history of the solar system. At a later time, the temperatures of iron asteroids in the main belt can be significantly colder than room temperature on Earth.

It is remarkable to study the physical properties of natural Fe-Ni samples at extreme
conditions to obtain insight into the rheological behavior of metal-rich and remanent core
planetary bodies for modeling impacts, collisions, and planet-forming processes. The next
steps include a systematic examination of oriented single-crystal Fe to further reveal
deformation mechanisms from c/a studies and the examination of further microstructure
complexities, such as variable and increased porosity and multiphase response during a
dynamic compression event.
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Phase diagram of Fe and results of shock-compressed Gibeon meteorite, kamacite lamellae, Fe₉₂.₅Ni₇.₅ on the Fe Hugoniot from this study. Shock compression data (blue circles) plotted on a pure Fe phase diagram: bcc boundaries from [53], melt curve from [14], along with the Fe₉₀.₃Ni₉.₇ phase boundary from [36] (pink section). XRD data collected at full compression are plotted along the Hugoniot for pure Fe (green curve) [54]. Guides for the eye of release isentropes (dashed lines) are derived from post-shock temperatures [34,35]. Our temperature uncertainties are 20%.</figcaption></figure>

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/min11060567/s1, Figure S1: Electron backscatter diffraction image for Gibeon meteorite sample, Figure S2: Representative VISAR data, Supplementary Discussion: Error Analysis and Shock Release Temperature Estimate.

**Author Contributions:** Project lead, A.E.G.; experiment participants, C.B., E.G. (Eric Galtier), E.G. (Eduardo Granados), A.H., H.J.L., S.M., B.N., K.R., D.R., W.L.M. and A.E.G.; sample characterization, B.M., A.E.G. and S.T.; data visualization software, C.B.; data analysis, S.T. and S.H.; writing—original draft preparation, S.T., R.C.-R. and A.E.G.; writing—review and editing, R.C.-R., R.W., W.L.M. and A.E.G.; All authors have read and agreed to the published version of the manuscript.

**Funding:** A.E.G. and W.L.M. were supported by NSF Geophysics (EAR0738873). A.E.G. acknowledges support from LANL Reines LDRD. A.E.G. and R.C.-R. acknowledge support from the DOE Office of Science, Fusion Energy Science ECA. This work was performed at the Matter at Extreme Conditions (MEC) instrument of the LCLS, supported by the US DOE Office of Science, Fusion Energy Science, under contract no. SF00515, and was supported by the LCLS, a National User Facility operated by Stanford University on behalf of DOE-BES. The Los Alamos National Laboratory is operated for the US DOE National Nuclear Security Administration under contract no. DE-AC52-06NA25396.

**Data Availability Statement:** All XRD and sample characterization data are available upon request.

**Acknowledgments:** We thank the LCLS and SLAC staff for assistance during/after the experiment.

**Conflicts of Interest:** The authors declare no conflict of interest.
## References

1. Asphaug, E. Interiors of small bodies and moons. *Nat. Commun.* **2020**, *11*, 10–12. [CrossRef]

2. Goldstein, J.I.; Scott, E.R.D.; Chabot, N.L. Iron meteorites: Crystallization, thermal history, parent bodies, and origin. *Geochemistry* **2009**, *69*, 293–325. [CrossRef]

3. Buchwald, V.F. *Handbook of Iron Meteorites. Their History, Distribution, Composition and Structure*; Danmarks Tekniske Hojskole: Lyngby, Denmark, 1975.

4. Oh, D.Y.; Goebel, D.; Polanskey, C.; Snyder, S.; Carr, G.; Collins, S.M.; Lantoine, G.; Landau, D.; Elkins-Tanton, L.; Lord, P.; et al. Psyche: Journey to a metal world. In *Proceedings of the 52nd AIAA/SAE/ASEE Joint Propulsion Conference, Salt Lake City, UT, USA, 25–27 July 2016*; p. 4541. [CrossRef]

5. Burbine, T.H.; Meibom, A.; Binzel, R.P. Mantle material in the main belt: Battered to bits? *Meteorit. Planet. Sci.* **1996**, *31*, 607–620. [CrossRef]

6. Yang, J.; Goldstein, J.I.; Scott, E.R.D.; Michael, J.R.; Kotula, P.G.; Pham, T.; Mccoy, T.J. Thermal and impact histories of reheated group IVA, IVB, and ungrouped iron meteorites and their parent asteroids. *Meteorit. Planet. Sci.* **2011**, *46*, 1227–1252. [CrossRef]

7. McKeehan, L.W. The crystal structure of iron-nickel alloys. *Phys. Rev.* **1923**, *21*, 402–407. [CrossRef]

8. Owen, E.A.; Yates, E.L.; Sully, A.H. An X-ray investigation of pure iron-nickel alloys. Part 4: The variation of lattice-parameter with composition. *Proc. Phys. Soc.* **1937**, *49*, 315–322. [CrossRef]

9. Purdy, G.R. Widmanstätten Structures. In *Encyclopedia of Materials: Science and Technology*; Elsevier: Amsterdam, The Netherlands, 2001; pp. 9575–9578.

10. Johnson, A.A.; Remo, J.L. A new interpretation of the mechanical properties of the Gibeon Meteorite. *J. Geophys. Res.* **1974**, *79*, 1142–1146. [CrossRef]

11. He, Y.; Godet, S.; Jacques, P.J.; Jonas, J.J. Crystallographic relations between face- and body-centred cubic crystals formed under near-equilibrium conditions: Observations from the Gibeon meteorite. *Acta Mater.* **2006**, *54*, 1323–1334. [CrossRef]

12. Jeanloz, R. The nature of the Earth's core. *Annu. Rev. Earth Planet. Sci.* **1990**, *18*, 357–386. [CrossRef]

13. Hawreliak, J.; Colvin, J.D.; Eggert, J.H.; Kalantar, D.H.; Lorenzana, H.E.; Stölken, J.S.; Davies, H.M.; Germann, T.C.; Holian, B.L.; Kadau, K.; et al. Analysis of the x-ray diffraction signal for the α- transition in shock-compressed iron: Simulation and experiment. *Phys. Rev. B Condens. Matter Mater. Phys.* **2006**, *74*, 184107. [CrossRef]

14. 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* **2013**, *340*, 464–466. [CrossRef] [PubMed]

15. Hwang, H.; Galtier, E.; Cynn, H.; Eom, I.; Chun, S.H.; Bang, Y.; Hwang, G.C.; Choi, J.; Kim, T.; Kong, M.; et al. Subnanosecond phase transition dynamics in laser-shocked iron. *Sci. Adv.* **2020**, *6*, 5132–5137. [CrossRef]

16. Mao, W.L.; Campbell, A.J.; Heinz, D.L.; Shen, G. Phase relations of Fe-Ni alloys at high pressure and temperature. *Phys. Earth Planet. Inter.* **2006**, *155*, 146–151. [CrossRef]

17. Steinle-Neumann, G.; Stixrude, L.; Cohen, R.E.; Gülseren, O. Elasticity of iron at the temperature of the Earth's inner core. *Nature* **2001**, *413*, 57–60. [CrossRef]

18. Mao, H.K.; Wu, Y.; Chen, L.C.; Shu, J.F.; Jephcoat, A.P. Static compression of iron to 300 GPa and Fe<sub>0.8</sub>Ni<sub>0.2</sub> alloy to 260 GPa: Implications for composition of the core. *J. Geophys. Res.* **1990**, *95*, 21737–21742. [CrossRef]

19. Hawreliak, J.A.; El-Dasher, B.; Lorenzana, H.; Kimminau, G.; Higginbotham, A.; Nagler, B.; Vinko, S.M.; Murphy, W.J.; Whitcher, T.; Wark, J.S.; et al. In situ x-ray diffraction measurements of the c/a ratio in the high-pressure ε phase of shock-compressed polycrystalline iron. *Phys. Rev. B Condens. Matter Mater. Phys.* **2011**, *83*, 1–6. [CrossRef]

20. Morrison, R.A.; Jackson, J.M.; Sturhahn, W.; Zhang, D.; Greenberg, E. Equations of State and Anisotropy of Fe-Ni-Si Alloys. *J. Geophys. Res. Solid Earth* **2018**, *123*, 4647–4675. [CrossRef]

21. Reagan, M.M.; Gleason, A.E.; Liu, J.; Krawczynski, M.J.; Van Orman, J.A.; Mao, W.L. The effect of nickel on the strength of iron nickel alloys: Implications for the Earth's inner core. *Phys. Earth Planet. Inter.* **2018**, *283*, 43–47. [CrossRef]

22. Schoelmerich, M.O.; Tschentscher, T.; Bhat, S.; Bolme, C.A.; Cunningham, E.; Farla, R.; Galtier, E.; Gleason, A.E.; Harmand, M.; Inubushi, Y.; et al. Evidence of shock-compressed stishovite above 300 GPa. *Sci. Rep.* **2020**, *10*, 1–8. [CrossRef] [PubMed]

23. Duffy, T.S.; Smith, R.F. Ultra-high pressure dynamic compression of geological materials. *Front. Earth Sci.* **2019**, *7*, 23. [CrossRef]

24. Coleman, A.L.; Gorman, M.G.; Briggs, R.; McWilliams, R.S.; McGonegle, D.; Bolme, C.A.; Gleason, A.E.; Fratanduono, D.E.; Smith, R.F.; Galtier, E.; et al. Identification of Phase Transitions and Metastability in Dynamically Compressed Antimony Using Ultrafast X-Ray Diffraction. *Phys. Rev. Lett.* **2019**, *122*, 255704. [CrossRef] [PubMed]

25. Tracy, S.J.; Smith, R.F.; Wicks, J.K.; Fratanduono, D.E.; Gleason, A.E.; Bolme, C.A.; Prakapenka, V.B.; Speziale, S.; Appel, K.; Fernandez-Pañella, A.; et al. In situ observation of a phase transition in silicon carbide under shock compression using pulsed x-ray diffraction. *Phys. Rev. B* **2019**, *99*, 214106. [CrossRef]

26. Gleason, A.E.; Bolme, C.A.; Lee, H.J.; Nagler, B.; Galtier, E.; Kraus, R.G.; Sandberg, R.; Yang, W.; Langenhorst, F.; Mao, W.L. Time-resolved diffraction of shock-released SiO<sub>2</sub> and diaplectic glass formation. *Nat. Commun.* **2017**, *8*, 1–6. [CrossRef] [PubMed]

27. Nagler, B.; Arnold, B.; Bouchard, G.; Boyce, R.F.; Boyce, R.M.; Callen, A.; Campell, M.; Curiel, R.; Galtier, E.; Garofoli, J.; et al. The Matter in Extreme Conditions instrument at the Linac Coherent Light Source. *J. Synchrotron Radiat.* **2015**, *22*, 520–525. [CrossRef]

28. Blaj, G.; Caragiulo, P.; Carini, G.; Carron, S.; Dragone, A.; Freytag, D.; Haller, G.; Hart, P.; Hasi, J.; Herbst, R.; et al. X-ray detectors at the Linac Coherent Light Source. *J. Synchrotron Radiat.* **2015**, *22*, 577–583. [CrossRef]
29. Bolme, C. *Bolmian XRD Projection Software*; LA-UR-15-25840, US DOE Copyright; Los Alamos National Laboratory: Los Alamos, NM, USA, 2015.

30. Prescher, C.; Prakapenka, V.B. DIOPTAS: A program for reduction of two-dimensional X-ray diffraction data and data exploration. *High Press. Res.* **2015**, *35*, 223–230. [CrossRef]

31. Wojdyr, M. Fityk: A general-purpose peak fitting program. *J. Appl. Crystallogr.* **2010**, *43*, 1126–1128. [CrossRef]

32. Brown, J.M.; McQueen, G. Phase transitions, Grüneisen parameter, and elasticity. *J. Geophys. Res.* **1986**, *91*, 7485–7494. [CrossRef]

33. Marsh, S.P. LASL Shock Hugoniot Data; University of California Press: Berkeley, CA, USA, 1980; ISBN 0-520-04008-2.

34. Huser, G.; Koenig, M.; Benuzzi-Mounaix, A.; Henry, E.; Vinci, T.; Faral, B.; Tomasini, M.; Telaro, B.; Batani, D. Temperature and melting of laser-shocked iron releasing into an LiF window. *Phys. Plasmas* **2005**, *12*, 1–4. [CrossRef]

35. Dorogokupets, P.I.; Dymshits, A.M.; Litasov, K.D.; Sokolova, T.S. Thermodynamics and Equations of State of Iron to 350 GPa and 6000 K. *Sci. Rep.* **2017**, *7*, 1–11. [CrossRef] [PubMed]

36. Komabayashi, T.; Hirose, K.; Ohishi, Y. In situ X-ray diffraction measurements of the fcc-hcp phase transition boundary of an Fe-Ni alloy in an internally heated diamond anvil cell. *Phys. Chem. Miner.* **2012**, *39*, 329–338. [CrossRef]

37. Wells, A.F. *Structural Inorganic Chemistry*; Oxford University Press: Oxford, UK, 1984; Volume 168.

38. Fischer, R.A.; Campbell, A.J. The axial ratio of hcp Fe and Fe-Ni-Si alloys to the conditions of Earth's inner core. *Am. Mineral.* **2015**, *100*, 2718–2724. [CrossRef]

39. Kalantar, D.H.; Belak, J.F.; Collins, G.W.; Colvin, J.D.; Davies, H.M.; Eggert, J.H.; Germann, T.C.; Hawreliak, J.; Holian, B.L.; Kadau, K.; et al. Direct observation of the α-ε Transition in shock-compressed iron via nanosecond X-ray diffraction. *Phys. Rev. Lett.* **2005**, *95*, 075502. [CrossRef] [PubMed]

40. Merkel, S.; Lincot, A.; Petitgirard, S. Microstructural effects and mechanism of bcc-hcp-bcc transformations in polycrystalline iron. *Phys. Rev. B* **2020**, *102*, 104103. [CrossRef]

41. Boehler, R.; Santamaría-Pérez, D.; Errandonea, D.; Mezouar, M. Melting, density, and anisotropy of iron at core conditions: New x-ray measurements to 150 GPa. *J. Phys. Conf. Ser.* **2008**, *121*. [CrossRef]

42. Fischer, R.A.; Campbell, A.J.; Shofner, G.A.; Lord, O.T.; Dera, P.; Prakapenka, V.B. Equation of state and phase diagram of FeO. *Earth Planet. Sci. Lett.* **2011**, *304*, 496–502. [CrossRef]

43. Yamazaki, D.; Ito, E.; Yoshino, T.; Yoneda, A.; Guo, X.; Zhang, B.; Sun, W.; Shimojuku, A.; Tsujino, N.; Kunimoto, T.; et al. P-V-T equation of state for ε-iron up to 80 GPa and 1900 K using the Kawai-type high pressure apparatus equipped with sintered diamond anvils. *Geophys. Res. Lett.* **2012**, *39*. [CrossRef]

44. Ono, S.; Kikegawa, T.; Hirao, N.; Mibe, K. High-pressure magnetic transition in hcp-Fe. *Am. Mineral.* **2010**, *95*, 880–883. [CrossRef]

45. Tateno, S.; Hirose, K.; Ohishi, Y.; Tatsumi, Y. The structure of iron in earth's inner core. *Science* **2010**, *330*, 359–361. [CrossRef] [PubMed]

46. Uchida, T.; Wang, Y.; Rivers, M.L.; Sutton, S.R. Stability field and thermal equation of state of ε-iron determined by synchrotron X-ray diffraction in a multianvil apparatus. *J. Geophys. Res. Solid Earth* **2001**, *106*, 21799–21810. [CrossRef]

47. Sakai, T.; Ohtani, E.; Hirao, N.; Ohishi, Y. Stability field of the hcp-structure for Fe, Fe-Ni, and Fe-Ni-Si alloys up to 3 Mbar. *Geophys. Res. Lett.* **2011**, *38*. [CrossRef]

48. Dewaele, A.; Loubeyre, P.; Occelli, F.; Mezouar, M.; Dorogokupets, P.I.; Torrent, M. Quasihydrostatic equation of state of Iron above 2 Mbar. *Phys. Rev. Lett.* **2006**, *97*, 215504. [CrossRef] [PubMed]

49. Tateno, S.; Hirose, K.; Komabayashi, T.; Ozawa, H.; Ohishi, Y. The structure of Fe-Ni alloy in Earth's inner core. *Geophys. Res. Lett.* **2012**, *39*. [CrossRef]

50. Lin, J.F.; Heinz, D.L.; Campbell, A.J.; Devine, J.M.; Mao, W.L.; Shen, G. Iron-nickel alloy in the Earth's core. *Geophys. Res. Lett.* **2002**, *29*, 1091–1093. [CrossRef]

51. Hok, S.; Gleason, A.E.; Bolme, C.; Rittman, D.; Merkel, S.; Ramos, K.J.; Morrow, B.; Lee, H.J.; Nagler, B.; Galtier, E.; et al. In situ strength measurement of shock-compressed iron via time-resolved X-ray diffraction. *AGUFM* **2019**, *2019*, DI23B-0043.

52. Dewaele, A.; Denoual, C.; Anzellini, S.; Occelli, F.; Mezouar, M.; Cordier, P.; Merkel, S.; Véron, M.; Rausch, E. Mechanism of the α-ε Phase transformation in iron. *Phys. Rev. B Condens. Matter Mater. Phys.* **2015**, *91*, 174105. [CrossRef]

53. Komabayashi, T.; Fei, Y. Internally consistent thermodynamic database for iron to the Earth's core conditions. *J. Geophys. Res. Solid Earth* **2010**, *115*. [CrossRef]

54. Wang, J.; Smith, R.F.; Eggert, J.H.; Braun, D.G.; Boehly, T.R.; Reed Patterson, J.; Celliers, P.M.; Jeanloz, R.; Collins, G.W.; Duffy, T.S. Ramp compression of iron to 273 GPa. *J. Appl. Phys.* **2013**, *114*. [CrossRef]