Revealing Exotic Nanophase Iron in Lunar Samples Through Impact-Driven Spatial Fingerprints

Nanophase iron (npFe) plays a crucial role in controlling the optical, chemical, and physical evolution of lunar regolith grains. While in-situ formation of npFe via reduction of native Fe-bearing minerals has long been considered a dominant pathway, recent mineralogical evidence from X.Zeng et al. (2025) reveals that the source of a significant fraction of npFe may be delivered directly by exotic micrometeoroid impacts (exotic npFe). Yet the atomic-scale processes governing how exotic np-Fe forms and survives during hypervelocity impacts remain largely unknown. To quantitatively compare in-situ and exotic delivery and formation of npFe, we perform a series of innovative atomistic modeling of micrometeoroid impacts with distinct projectile target compositions: (1) SiO$_2$ projectiles on Fe$_2$SiO$_4$ targets (in-situ formation), (2) Fe$_2$SiO$_4$ projectiles on SiO$_2$ targets (exotic delivery). Our results reveal distinct mechanistic fingerprints: in-situ np-Fe forms diffusely and radially around the impact site, whereas exotic np-Fe is efficiently retained and concentrated in asymmetric, momentum-aligned clusters. These contrasting spatial signatures provide a potential diagnostic criterion for distinguishing exotic versus in-situ np-Fe in returned lunar soils. In agreement with Chang’e-5 observations, our simulations demonstrate that exotic np-Fe production can be substantial, particularly in Fe-poor terrains such as highland regions. These findings highlight the need to account for exotic np-Fe when interpreting space weathering processes and remote-sensing data for the Moon and other airless bodies.

💡 Research Summary

This paper presents the first atomistic‐scale comparison of two competing pathways for nanophase iron (npFe) formation on airless bodies: (1) in‑situ generation of Fe from the lunar target material during a micrometeoroid impact, and (2) direct delivery of Fe by an Fe‑rich impactor (exotic npFe). Using reactive force‑field (ReaxFF) molecular dynamics, the authors simulate hypervelocity collisions at 12 km s⁻¹ and 45° incidence for two complementary cases: SiO₂ projectiles striking Fe₂SiO₄ (olivine) targets (in‑situ case) and Fe₂SiO₄ projectiles striking SiO₂ targets (exotic case). The simulations employ a hemispherical target geometry, a 0.1 fs timestep, and run for up to 30 ps, allowing the authors to capture the ultrafast heating, shock compression, bond breaking, and subsequent chemical reactions that occur during impact.

Key findings are:

• Retention efficiency – In the exotic case, about 91 % of the impactor’s Fe atoms remain in the near‑surface region (68 % embedded, 23 % redeposited). This is comparable to the retention of Si and O, indicating that mechanical entrainment dominates the early stages of delivery.
• Spatial fingerprints – In‑situ npFe nucleates diffusely around the impact point, forming a roughly spherical halo of nanometer‑scale metallic Fe particles within the melt‑quenched silicate matrix. By contrast, exotic Fe aggregates into asymmetric, momentum‑aligned clusters that are 2–5 nm in size and oriented along the impact direction.
• Velocity distribution – The velocity distribution function (VDF) of ejected Fe shows a dominant low‑velocity peak (0–10 m s⁻¹) corresponding to the clustered Fe that is well below lunar escape velocity (≈2.38 km s⁻¹). A high‑velocity tail (20–60 m s⁻¹) represents individual Fe‑bearing molecules (FeO₂, Fe₂SiO₂, etc.) generated by vaporization. The low‑velocity component implies that most exotic Fe will re‑deposit within tens of meters of the impact site, preserving a localized signature.
• Implications for lunar terrains – Because highland regions are Fe‑poor, the exotic delivery pathway can contribute a substantial fraction (30–40 %) of the total npFe budget, consistent with observations from Chang’e‑5 samples. This challenges the long‑standing assumption that npFe is produced solely by reduction of native lunar minerals.

Methodologically, the study demonstrates the power of ReaxFF MD to resolve bond‑order changes, charge transfer, and short‑lived intermediates (e.g., Fe²⁺ → Fe⁰ reduction, disproportionation reactions) that are inaccessible to classical non‑reactive potentials or laser‑ablation experiments, which neglect the projectile’s chemistry. The authors also develop a DBSCAN‑inspired clustering algorithm that incorporates bond‑length criteria to identify physically meaningful Fe clusters, and they compute VDFs to quantify escape versus retention.

Limitations include the relatively short simulation time (30 ps), which does not capture long‑term cooling, grain growth, or diffusion that may alter cluster size distribution over microseconds to seconds. Additionally, the ReaxFF parameter set, while validated for silicates, may have uncertainties under the extreme pressures and temperatures of hypervelocity impacts.

Future work suggested by the authors involves extending the temporal window with coarse‑grained heat‑diffusion models, exploring a broader range of impact velocities, angles, and projectile compositions, and coupling the atomistic results with laboratory impact experiments that combine laser heating and actual projectile delivery.

In summary, the paper provides compelling atomistic evidence that exotic npFe, delivered by Fe‑rich micrometeoroids, forms distinct, momentum‑aligned clusters that remain locally concentrated and low‑velocity. This spatial fingerprint offers a diagnostic tool for distinguishing exotic from in‑situ npFe in lunar samples and remote‑sensing data, and it necessitates the inclusion of exotic Fe delivery in models of space weathering for the Moon and other airless planetary bodies.