No Arabic abstract
Fe, Mg, and O are among the most abundant elements in terrestrial planets. While the behavior of the Fe-O, Mg-O, and Fe-Mg binary systems under pressure have been investigated, there are still very few studies of the Fe-Mg-O ternary system at relevant Earths core and super-Earths mantle pressures. Here, we use the adaptive genetic algorithm (AGA) to study ternary Fe$_x$Mg$_y$O$_z$ phases in a wide range of stoichiometries at 200 GPa and 350 GPa. We discovered three dynamically stable phases with stoichiometries FeMg$_2$O$_4$, Fe$_2$MgO$_4$, and FeMg$_3$O$_4$ with lower enthalpy than any known combination of Fe-Mg-O high-pressure compounds at 350 GPa. With the discovery of these phases, we construct the Fe-Mg-O ternary convex hull. We further clarify the composition- and pressure-dependence of structural motifs with the analysis of the AGA-found stable and metastable structures. Analysis of binary and ternary stable phases suggest that O, Mg, or both could stabilize a BCC iron alloy at inner core pressures.
The solid inner core of the Earth is predominantly composed of iron alloyed with several percent Ni and some lighter elements, Si, S, O, H, and C being the prime candidates. There have been a growing number of papers investigating C and H as possible light elements in the core, but the results are contradictory. Here, using ab initio simulations, we study the Fe-C and Fe-H systems at inner core pressures (330-364 GPa). Using the evolutionary structure prediction algorithm USPEX, we have determined the lowest-enthalpy structures of possible carbides (FeC, Fe2C, Fe3C, Fe4C, FeC2, FeC3, FeC4 and Fe7C3) and hydrides (Fe4H, Fe3H, Fe2H, FeH, FeH2, FeH3, FeH4) and have found that Fe2C (Pnma) is the most stable iron carbide at pressures of the inner core, while FeH, FeH3 and FeH4 are stable iron hydrides at these conditions. For Fe3C, the cementite structure (Pnma) and the Cmcm structure recently found by random sampling are less stable than the I-4 and C2/m structures found here. We found that FeH3 and FeH4 adopt chemically interesting thermodynamically stable structures, in both compounds containing trivalent iron. The density of the inner core can be matched with a reasonable concentration of carbon, 11-15 mol.percent (2.6-3.7 wt.percent) at relevant pressures and temperatures. This concentration matches that in CI carbonaceous chondrites and corresponds to the average atomic mass in the range 49.3-51.0, in close agreement with inferences from the Birchs law for the inner core. Similarly made estimates for the maximum hydrogen content are unrealistically high, 17-22 mol.percent (0.4-0.5 wt.percent), which corresponds to the average atomic mass in the range 43.8-46.5. We conclude that carbon is a better candidate light alloying element than hydrogen.
We report on the thermal and electrical conductivities of two liquid silicon-oxygen-iron mixtures (Fe$_{0.82}$Si$_{0.10}$O$_{0.08}$ and Fe$_{0.79}$Si$_{0.08}$O$_{0.13}$), representative of the composition of the Earths outer core at the relevant pressure-temperature conditions, obtained from density functional theory calculations with the Kubo-Greenwood formulation. We find thermal conductivities $k$ =100 (160) W m$^{-1}$ K$^{-1}$, and electrical conductivities $sigma = 1.1 (1.3) times 10^6 Omega^{-1}$ m$^{-1}$ at the top (bottom) of the outer core. These new values are between 2 and 3 times higher than previous estimates, and have profound implications for our understanding of the Earths thermal history and the functioning of the Earths magnetic field, including rapid cooling rate for the whole core or high level of radiogenic elements in the core. We also show results for a number of structural and dynamic properties of the mixtures, including the partial radial distribution functions, mean square displacements, viscosities and speeds of sound.
The transport properties of iron under Earths inner core conditions are essential input for the geophysical modelling but are poorly constrained experimentally. Here we show that the thermal and electrical conductivity of iron at those conditions remains high even if the electron-electron-scattering (EES) is properly taken into account. This result is obtained by ab initio simulations taking into account consistently both thermal disorder and electronic correlations. Thermal disorder suppresses the non-Fermi-liquid behavior of the body-centered cubic iron phase, hence, reducing the EES; the total calculated thermal conductivity of this phase is 220 Wm$^{-1}$K$^{-1}$ with the EES reduction not exceeding 20%. The EES and electron-lattice scattering are intertwined resulting in breaking of the Matthiessens rule with increasing EES. In the hexagonal close-packed iron the EES is also not increased by thermal disorder and remains weak. Our main finding thus holds for the both likely iron phases in the inner core.
Being a lithophile element at ambient pressure, magnesium is long believed to be immiscible with iron. A recent study by Gao et al. [1] showed that pressure turns magnesium into a siderophile element and can produce unconventional Fe-Mg compounds. Here, we extend the investigation to exoplanetary pressure conditions using an adaptive genetic algorithm-based variable-composition structural prediction approach. We identify several Fe-Mg phases up to 3 TPa. Our cluster alignment analysis reveals that most of the predicted Fe-Mg compounds prefer a BCC packing motif at terapascal pressures. This study provides a more comprehensive structure database to support future investigations of the high-pressure structural behavior of Fe-Mg and ternary, quaternary, etc. compounds involving these elements.
The stability, structure and properties of carbonate minerals at lower mantle conditions has significant impact on our understanding of the global carbon cycle and the composition of the interior of the Earth. In recent years, there has been significant interest in the behavior of carbonates at lower mantle conditions, specifically in their carbon hybridization, which has relevance for the storage of carbon within the deep mantle. Using high-pressure synchrotron X-ray diffraction in a diamond anvil cell coupled with direct laser heating of CaCO$_{3}$ using a CO$_{2}$ laser, we identify a crystalline phase of the material above 40 GPa $-$ corresponding to a lower mantle depth of around 1,000 km $-$ which has first been predicted by textit{ab initio} structure predictions. The observed $sp^{2}$ carbon hybridized species at 40 GPa is monoclinic with $P2_{1}/c$ symmetry and is stable up to 50 GPa, above which it transforms into a structure which cannot be indexed by existing known phases. A combination of textit{ab initio} random structure search (AIRSS) and quasi-harmonic approximation (QHA) calculations are used to re-explore the relative phase stabilities of the rich phase diagram of CaCO$_{3}$. Nudged elastic band (NEB) calculations are used to investigate the reaction mechanisms between relevant crystal phases of CaCO$_{3}$ and we postulate that the mineral is capable of undergoing $sp^{2}$-$sp^{3}$ hybridization change purely in the $P2_{1}/c$ structure $-$ forgoing the accepted post-aragonite $Pmmn$ structure.