No Arabic abstract
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.
The Earth acts as a gigantic heat engine driven by decay of radiogenic isotopes and slow cooling, which gives rise to plate tectonics, volcanoes, and mountain building. Another key product is the geomagnetic field, generated in the liquid iron core by a dynamo running on heat released by cooling and freezing to grow the solid inner core, and on chemical convection due to light elements expelled from the liquid on freezing. The power supplied to the geodynamo, measured by the heat-flux across the core-mantle boundary (CMB), places constraints on Earths evolution. Estimates of CMB heat-flux depend on properties of iron mixtures under the extreme pressure and temperature conditions in the core, most critically on the thermal and electrical conductivities. These quantities remain poorly known because of inherent difficulties in experimentation and theory. Here we use density functional theory to compute these conductivities in liquid iron mixtures at core conditions from first principles- the first directly computed values that do not rely on estimates based on extrapolations. The mixtures of Fe, O, S, and Si are taken from earlier work and fit the seismologically-determined core density and inner-core boundary density jump. We find both conductivities to be 2-3 times higher than estimates in current use. The changes are so large that core thermal histories and power requirements must be reassessed. New estimates of adiabatic heat-flux give 15-16 TW at the CMB, higher than present estimates of CMB heat-flux based on mantle convection; the top of the core must be thermally stratified and any convection in the upper core driven by chemical convection against the adverse thermal buoyancy or lateral variations in CMB heat flow. Power for the geodynamo is greatly restricted and future models of mantle evolution must incorporate a high CMB heat-flux and explain recent formation of the inner core.
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.
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.
Molecular dynamics simulations were performed to understand the role of twin boundaries on deformation behaviour of body-centred cubic (BCC) iron (Fe) nanopillars. The twin boundaries varying from one to five providing twin boundary spacing in the range 8.5 - 2.8 nm were introduced perpendicular to the loading direction. The simulation results indicated that the twin boundaries in BCC Fe play a contrasting role during deformation under tensile and compressive loadings. During tensile deformation, a large reduction in yield stress was observed in twinned nanopillars compared to perfect nanopillar. However, the yield stress exhibited only marginal variation with respect to twin boundary spacing. On the contrary, a decrease in yield stress with increase in twin boundary spacing was obtained during compressive deformation. This contrasting behaviour originates from difference in operating mechanisms during yielding and subsequent plastic deformation. It has been observed that the deformation under tensile loading was dominated mainly by twin growth mechanism, due to which the twin boundaries offers a negligible resistance to slip of twinning partials. This is reflected in the negligible variation of yield stress as a function of twin boundary spacing. On the other hand, the deformation was dominated by nucleation and slip of full dislocations under compressive loading. The twin boundaries offer a strong repulsive force on full dislocations resulting in the yield stress dependence on twin boundary spacing. Further, it has been observed that the curved twin boundary can acts as a source for full dislocation. The occurrence of twin-twin interaction during tensile deformation and dislocation-twin interaction during compressive deformation were presented and discussed.
Radiation damage in body-centered cubic (BCC) Fe has been extensively studied by computer simulations to quantify effects of temperature, impinging particle energy, and the presence of extrinsic particles. However, limited investigation has been conducted into the effects of mechanical stresses and strain. In a reactor environment, structural materials are often mechanically strained, and an expanded understanding of how this strain affects the generation of defects may be important for predicting microstructural evolution and damage accumulation under such conditions. In this study, we have performed molecular dynamics simulations in which various types of homogeneous strains are applied to BCC Fe and the effect on defect generation is examined. It is found that volume-conserving shear strains yield no statistically significant variations in the stable number of defects created via cascades in BCC Fe. However, strains that result in volume changes are found to produce significant effects on defect generation.