Novel phases in the Fe-Si-O system at terapascal pressures

Using crystal-structure prediction and ab initio calculations, this study identifies three new metallic ternary Fe-Si-O compounds stable at terapascal pressures that adopt pseudo-binary structures and suggest a fundamentally different mechanism for iron incorporation in super-Earth mantles, potentially triggering silicate dissociation at pressures below ~3 TPa.

The Gist

Imagine the Earth's core as a giant, cosmic pressure cooker. Deep inside our planet, rocks are crushed so hard that they behave nothing like the sand or gravel we know on the surface. Now, imagine planets much bigger than Earth—called "Super-Earths." Inside these giants, the pressure is so extreme (trillions of times stronger than the air pressure at sea level) that the rules of chemistry and physics change completely.

This paper is like a cosmic treasure hunt to discover what new "mineral recipes" exist in these extreme environments, specifically looking at a mixture of Iron (Fe), Silicon (Si), and Oxygen (O).

Here is the story of what the scientists found, explained simply:

1. The Problem: We Didn't Know the Recipe

For a long time, scientists knew how Silicon and Oxygen behave under high pressure (making things like quartz or glass). They also knew how Iron and Oxygen behave. But they didn't know what happens when you mix them all together in the crushing depths of a Super-Earth. It's like knowing how to bake a chocolate cake and how to bake a vanilla cake, but having no idea what happens if you try to bake a "chocolate-vanilla hybrid" in a furnace that's 10,000 degrees hot.

2. The Method: The Digital Crystal Builder

Instead of trying to crush rocks in a lab (which is incredibly hard to do at these pressures), the scientists used a supercomputer. They used a digital tool called an "Adaptive Genetic Algorithm."

• The Analogy: Think of this like a digital "Lego Master." The computer tries millions of different ways to snap Iron, Silicon, and Oxygen atoms together. It builds a structure, checks if it falls apart, and if it does, it tries a new arrangement. It keeps evolving the best designs until it finds the ones that are strong enough to survive the pressure.

3. The Discovery: Three New "Super-Rocks"

The computer found three stable new compounds that act like a bridge between Iron Oxide and Silicon Oxide.

• The Characters: They found three new "characters" in the story: FeSiO₄, Fe₄Si₅O₁₈, and FeSi₂O₆.
• The Temperature Twist: Two of these rocks are stable when it's cold (like deep space), but the third one, FeSi₂O₆, only becomes the "king" when it gets very hot (around 2,000°C to 4,000°C). This is crucial because the insides of Super-Earths are incredibly hot.

4. The Structure: A New Way to Stack Bricks

In our Earth's mantle, Iron usually fits into a specific "cage" made of Oxygen atoms, like a guest sitting in a square chair.

• The New Twist: In these new high-pressure rocks, the Iron atoms are doing something weird. Some are sitting in six-sided cages, and others are squeezed into nine-sided cages.
• The Metaphor: Imagine a hotel where guests usually get a standard room. But in these Super-Earths, the hotel has to build brand new, weirdly shaped suites (some with 9 walls!) to fit the Iron guests in. This is a completely different way of organizing the atoms than we see on Earth.

5. The Big Implication: The "Dissociation" Domino Effect

This is the most exciting part. The scientists realized that these new rocks might change how we think about the entire interior of a Super-Earth.

• The Old Idea: We thought that in Super-Earths, Iron would just quietly hide inside the main silicate rocks (like a secret ingredient in a soup).
• The New Reality: These new findings suggest that at extreme pressures, the main silicate rocks might actually break apart (dissociate).
• The Analogy: Imagine a chocolate bar (the silicate rock). On Earth, if you squeeze it, it just gets smaller. But in a Super-Earth, the pressure is so high that the chocolate bar snaps into separate pieces: pure cocoa butter (Iron Oxide) and pure sugar (Silicon Oxide).
• Why it matters: If the rocks break apart, the planet's interior might become layered, like a cake with distinct frosting and sponge layers, rather than a smooth mix. This could stop the "convection" (the churning motion that drives plate tectonics and volcanoes), meaning these Super-Earths might be geologically dead or very different from Earth.

Summary

In short, this paper tells us that under the crushing pressure of Super-Earths, Iron and Silicon don't just mix quietly. They form brand new, weirdly shaped metallic rocks that only exist at high temperatures. These rocks might cause the planet's mantle to split apart into layers, fundamentally changing how these giant planets evolve and whether they could ever be habitable.

It's a reminder that in the extreme universe, the rules we learned in school are just the beginning, and the real magic happens when you turn the pressure dial all the way up.

Technical Summary

1. Problem Statement

The interior structure of terrestrial planets, particularly "Super-Earths" (planets with masses up to ~20 Earth masses), involves extreme pressures reaching the Terapascal (TPa) range (1 TPa = 10 Mbar). While the behavior of the dominant Mg-Si-O system in these environments has been extensively modeled, the role of Iron (Fe) remains poorly constrained.

• The Gap: In Earth's lower mantle (GPa range), iron primarily substitutes for Magnesium ([Fe]Mg) or Silicon ([Fe]Si) in silicate structures. However, it is unknown how Fe incorporates into silicates at TPa pressures characteristic of Super-Earth mantles.
• The Challenge: Existing models often assume Fe behaves similarly to Earth's mantle conditions. There is a lack of knowledge regarding stable ternary Fe-Si-O compounds at ~1 TPa, which are critical for understanding the thermodynamics, density, and seismic properties of deep planetary interiors.

2. Methodology

The authors employed a combination of advanced computational techniques to predict and validate new crystal structures:

• Crystal Structure Prediction: Utilized the Adaptive Genetic Algorithm (AGA) assisted by classical Embedded-Atom Method (EAM) potentials to explore the compositional space of the Fe-Si-O system. They searched various stoichiometries ($Fe_xSi_yO_z$) and pseudo-binary combinations.
• First-Principles Calculations: Performed ab initio calculations using Density Functional Theory (DFT) within the Quantum ESPRESSO package.
  • Functional: Non-spin-polarized Local Density Approximation (LDA) was used for structural optimization, with spin-polarized checks confirming quenched magnetic moments.
  • Pseudopotentials: Vanderbilt-type ultrasoft pseudopotentials were used for Fe, Si, and O.
• Thermodynamic Stability:
  • Phonon Calculations: Used the finite displacement method and PHONOPY code to ensure dynamical stability (absence of imaginary frequencies).
  • Quasi-Harmonic Approximation (QHA): Calculated Gibbs free energies at finite temperatures (up to ~10,000 K) by including vibrational and electronic entropy contributions.
• Convex Hull Analysis: Constructed the convex hull of formation enthalpies at 1 TPa and 0 K to identify thermodynamically stable ground states.

3. Key Contributions

The study identifies three novel ternary compounds in the Fe-Si-O system stable near 1 TPa:

1. $P\bar{3}$ FeSiO$_4$
2. $P\bar{3}$ Fe$_4$Si$_5$O$_{18}$
3. $P\bar{3}$ FeSi$_2$O$_6$

Key findings include:

• Structural Continuity: These phases are "pseudo-binaries" derived from the stable end-member oxides at 1 TPa: Fe$_2$P-type SiO$_2$ and Pnma-type FeO$_2$.
  • The structures arise from substituting Fe for Si in Fe$_2$P-type SiO$_2$ or Si for Fe in Pnma-type FeO$_2$.
  • $P\bar{3}$ FeSi$_2$O$_6$ acts as the fundamental structural unit, with the other two phases derivable via further atomic substitutions.
• Coordination Chemistry: Unlike FeO at similar pressures (which features 8-fold cubic coordination), these new phases accommodate Fe in 6-fold (octahedral) and 9-fold (tri-capped trigonal prism) coordination.
• Electronic Properties: All three identified phases are metallic and paramagnetic, with electronic density of states dominated by Oxygen and Iron states.

4. Results and Thermodynamic Behavior

• Temperature Dependence:
  • Low Temperature (< ~2000 K): $P\bar{3}$ FeSiO$_4$ and $P\bar{3}$ Fe$_4$Si$_5$O$_{18}$ are thermodynamically stable.
  • High Temperature (> ~2000 K): $P\bar{3}$ FeSi$_2$O$_6$ becomes the favored phase. Its stability at high temperatures is driven significantly by electronic entropy due to a high density of states at the Fermi level.
• Phase Diagrams:
  • The stability fields depend heavily on the Fe/Si ratio and planetary mass (which correlates with pressure/temperature at the Core-Mantle Boundary, CMB).
  • Fe/Si = 1.0: FeSiO$_4$ is stable in smaller Super-Earths (<4.5 $M_\oplus$) but dissociates into FeSi$_2$O$_6$ + FeO$_2$ in larger planets.
  • Fe/Si < 1: Complex recombination transitions occur. For example, at Fe/Si = 0.8, FeSiO$_4$ and SiO$_2$ coexist in smaller planets, transitioning to FeSi$_2$O$_6$ + FeO$_2$ in larger ones.
• Clapeyron Slopes (CS): The phase transitions involving these new compounds exhibit unusually large negative Clapeyron slopes (e.g., -109.9 mPa/K).
  • Implication: Negative CS inhibits convective flow across these boundaries, potentially inducing chemical or thermal layering in Super-Earth mantles, unlike the more continuous convection inferred for Earth.

5. Significance and Implications

• Redefining Iron Incorporation: The discovery suggests that at TPa pressures, Iron does not simply substitute into canonical Mg-silicates (like bridgmanite or olivine) as it does in Earth's mantle. Instead, Fe may preferentially form distinct FeO$_2$-SiO$_2$ pseudo-binary phases.
• Silicate Dissociation: The stability of these phases may trigger the dissociation of silicates into oxides at pressures below ~3 TPa:
  $$(Mg,Fe)_2(Si,Fe)O_4 \rightarrow 2(Mg,Fe)O + (Si,Fe)O_2$$
  The extent of this dissociation is governed by the iron content.
• Planetary Modeling: These findings necessitate a revision of models for Super-Earth interiors. The presence of metallic, paramagnetic Fe-Si-O phases could alter the electrical conductivity, seismic velocity profiles, and thermal evolution of these planets.
• Experimental Guidance: The identification of low-pressure analogs (e.g., NaLnF$_4$ structures) provides a pathway for laboratory synthesis and experimental validation of these high-pressure phases.

In conclusion, this work fundamentally shifts the understanding of iron chemistry in deep planetary interiors, moving from a substitutional model to a model where distinct ternary oxides dominate the mineralogy at Terapascal pressures.