Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userFirst principles study of density, viscosity, and diffusion coefficients of liquid MgSiO3 at conditions of the Earths deep mantle
92
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Constant-pressure constant-temperature {it ab initio} molecular dynamics simulations at high temperatures have been used to study MgSiO$_3$ liquid, the major constituent of the Earths lower mantle to conditions of the Earths core-mantle boundary (CMB). We have performed variable-cell {it ab initio} molecular dynamic simulations at relevant thermodynamic conditions across one of the measured melting curves. The calculated equilibrium volumes and densities are compared with the simulations using an orthorhombic perovskite configuration under the same conditions. For molten MgSiO$_3$, we have determined the diffusion coefficients and shear viscosities at different thermodynamic conditions. Our results provide new constraints on the properties of molten MgSiO$_3$ at conditions near the core-mantle boundary. The volume change on fusion is positive throughout the pressure-temperature conditions examined and ranges from 5% at 88 GPa and 3500 K to 2.9% at 120 GPa and 5000 K. Nevertheless, neutral or negatively buoyant melts from (Mg,Fe)SiO$_3$ perovskite compositions at deep lower mantle conditions are consistent with existing experimental constraints on solid-liquid partition coefficients for Fe. Our simulations indicate that MgSiO$_3$ is liquid at 120 GPa and 4500 K, consistent with the lower range of experimental melting curves for this material. Linear extrapolation of our results indicates that the densities of liquid and solid perovskite MgSiO$_3$ will become equal near 180 GPa.
rate research
Read More
The temperature anomalies in the Earths mantle associated with thermal convection1 can be inferred from seismic tomography, provided that the elastic properties of mantle minerals are known as a function of temperature at mantle pressures. At present, however, such information is difficult to obtain directly through laboratory experiments. We have therefore taken advantage of recent advances in computer technology, and have performed finite-temperature ab initio molecular dynamics simulations of the elastic properties of MgSiO3 perovskite, the major mineral of the lower mantle, at relevant thermodynamic conditions. When combined with the results from tomographic images of the mantle, our results indicate that the lower mantle is either significantly anelastic or compositionally heterogeneous on large scales. We found the temperature contrast between the coldest and hottest regions of the mantle, at a given depth, to be about 800K at 1000 km, 1500K at 2000 km, and possibly over 2000K at the core-mantle boundary.
Calcium silicate perovskite (CaSiO$_3$) is one of the major mineral components of the lower mantle, but has been the subject of relatively little work compared to the more abundant Mg-based materials. One of the major problems related to CaSiO$_3$ that is still the subject of research is its crystal structure under lower mantle conditions - a cubic Pm$bar{3}$m structure is accepted in general, but some have suggested that lower-symmetry structures may be relevant. In this paper, we use a fully first-principles vibrational self-consistent field method to perform high accuracy anharmonic vibrational calculations on several candidate structures at a variety of points along the geotherm near the base of the lower mantle to investigate the stability of the cubic structure and related distorted structures. Our results show that the cubic structure is the most stable throughout the lower mantle, and that this result is robust against the effects of thermal expansion.
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.
The electronic state and transport properties of hot dense iron are of the utmost importance to geophysics. Combining the density functional and dynamical mean field theories we study the impact of electron correlations on electrical and thermal resistivity of hexagonal close-packed $epsilon$-Fe at Earths core conditions. $epsilon$-Fe is found to behave as a nearly perfect Fermi liquid. The quadratic dependence of the scattering rate in Fermi liquids leads to a modification of the Wiedemann-Franz law with suppression of the thermal conductivity as compared to the electrical one. This significantly increases the electron-electron thermal resistivity which is found to be of comparable magnitude to the electron-phonon one. The implications of this effect on the dynamics of Earths core is discussed.
Liquid metals at extreme pressures and temperatures are widely interested in the high-pressure community. Based on density functional theory molecular dynamics, we conduct first-principles investigations on the equation of state (EOS) and structures of four metals (Cu, Fe, Pb, and Sn) at 1.5--5 megabar conditions and 5$times10^3$--4$times10^4$ K. Our first-principles EOS data enable evaluating the performance of four EOS models in predicting Hugoniot densities and temperatures of the four systems. We find the melting temperature of Cu is 1000--2000 K higher and shows a similar Clapeyron slope, in comparison to those of Fe. Our structure, coordination number, and diffusivity analysis indicates all the four liquid metals form similar simple close-packed structures. Our results set theoretical benchmarks for EOS development and structures of metals in their liquid states and under dynamic compression.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
