No Arabic abstract
It has long been assumed the Earths solid inner core started to grow when molten iron cooled to its melting point. However, the nucleation mechanism, which is a necessary step of crystallization, has not been well understood. Recent studies found it requires an unrealistic degree of undercooling to nucleate the stable hexagonal close-packed (hcp) phase of iron, which can never be reached under the actual Earths core conditions. This contradiction leads to the inner core nucleation paradox [1]. Here, using a persistent-embryo method and molecular dynamics simulations, we demonstrate that the metastable body-centered cubic (bcc) phase of iron has a much higher nucleation rate than the hcp phase under inner-core conditions. Thus, the bcc nucleation is likely to be the first step of inner core formation instead of direct nucleation of the hcp phase. This mechanism reduces the required undercooling of iron nucleation, which provides a key factor to solve the inner-core nucleation paradox. The two-step nucleation scenario of the inner core also opens a new avenue for understanding the structure and anisotropy of the present inner core.
In a first approximation the Earths interior has an isotropic structure with a spherical symmetry. Over the last decades the geophysical observations have revealed, at different spatial scales, the existence of several perturbations from this basic structure. Some of them are situated in the neighborhood of the inner core boundary (ICB). One of the best documented perturbations is the asymmetry at the top of the inner core (ATIC) characterized by faster seismic wave velocity in the eastern hemisphere than in the western hemisphere. All existing explanations are based on a hemispheric variation of the material properties near ICB inside the inner core. Using numerical simulations of the seismic ray propagation, we show that the ATIC can be explained as well by the displacement of the inner core towards east in the equatorial plane tens of kilometers from the Earths center, without modifying the spherical symmetry in the upper inner core. The hypothesis of a displaced inner core is also sustained by other observed hemispheric asymmetries at the top of the inner core and at the bottom of the outer core. A displaced inner core would have major implications for many mechanical, thermal, and magnetic phenomena in the Earths interior.
The crystal structure of iron in the Earths inner core remains debated. Most recent experiments suggest a hexagonal-close-packed (hcp) phase. In simulations, it has been generally agreed that the hcp-Fe is stable at inner core pressures and relatively low temperatures. At high temperatures, however, several studies suggest a body-centered-cubic (bcc) phase at the inner core condition. We have examined the crystal structure of iron at high pressures over 2 million atmospheres (>200GPa) and at high temperatures over 5000 kelvin in a laser-heated diamond cell using microstructure analysis combined with $textit{in-situ}$ x-ray diffraction. Experimental evidence shows a bcc-Fe appearing at core pressures and high temperatures, with an hcp-bcc transition line in pressure-temperature space from about 95$pm$2GPa and 2986$pm$79K to at least 222$pm$6GPa and 4192$pm$104K. The trend of the stability field implies a stable bcc-Fe at the Earths inner core condition, with implications including a strong candidate for explaining the seismic anisotropy of the Earths inner core.
We employ state-of-the-art ab initio simulations within the dynamical mean-field theory to study three likely phases of iron (hexogonal close-packed, hcp, face centered cubic, fcc, and body centered cubic, bcc) at the Earths core conditions. We demonstrate that the correction to the electronic free energy due to correlations can be significant for the relative stability of the phases. The strongest effect is observed in bcc Fe, which shows a non-Fermi liquid behaviour, and where a Curie-Weiss behaviour of the uniform susceptbility hints at a local magnetic moment still existing at 5800 K and 300 GPa. We predict that all three structures have sufficiently high magnetic susceptibility to stabilize the geodynamo.
We have studied the body-centered cubic (bcc), face-centered cubic (fcc) and hexagonal close-packed (hcp) phases of Fe alloyed with 25 at. % of Ni at Earths core conditions using an ab initio local density approximation + dynamical mean-field theory (LDA+DMFT) approach. The alloys have been modeled by ordered crystal structures based on the bcc, fcc, and hcp unit cells with minimum possible cell size allowing for the proper composition. Our calculations demonstrate that the strength of electronic correlations on the Fe 3d shell is highly sensitive to the phase and local environment. In the bcc phase the 3d electrons at the Fe site with Fe only nearest neighbors remain rather strongly correlated even at extreme pressure-temperature conditions, with the local and uniform magnetic susceptibility exhibiting a Curie-Weiss-like temperature evolution and the quasi-particle lifetime {Gamma} featuring a non-Fermi-liquid temperature dependence. In contrast, for the corresponding Fe site in the hcp phase we predict a weakly-correlated Fermi-liquid state with a temperature-independent local susceptibility and a quadratic temperature dependence of {Gamma}. The iron sites with nickel atoms in the local environment exhibit behavior in the range between those two extreme cases, with the strength of correlations gradually increasing along the hcp-fcc-bcc sequence. Further, the inter-site magnetic interactions in the bcc and hcp phases are also strongly affected by the presence of Ni nearest neighbors. The sensitivity to the local environment is related to modifications of the Fe partial density of states due to mixing with Ni 3d-states.
In a first approximation, the Earths interior has an isotropic structure with a spherical symmetry. Over the last decades the geophysical observations have revealed, at different spatial scales, the existence of several perturbations from this basic structure. In this paper we discuss the hemispheric perturbations induced to this basic structure if the inner core is displaced from the center of mass of the Earth. Using numerical simulations of the observed hemispheric asymmetry of the seismic waves traveling through the upper inner core, with faster arrival times and higher attenuation in the Eastern Hemisphere, we estimate that the present position of the inner core is shifted by tens of kilometers from the Earths center eastward in the equatorial plane. If the only forces acting on the inner core were the gravitational forces, then its equilibrium position would be at the Earths center and the estimated displacement would not be possible. We conjecture that, due to interactions with the flow and the magnetic field inside the outer core, the inner core is in a permanent chaotic motion. To support this hypothesis we analyze more than ten different geophysical phenomena consistent with an inner core motion dominated by time scales from hundreds to thousands of years.