The melting curve of Ni up to 100 GPa has been calculated using first principles methods based on density functional theory (DFT). We used two complementary approaches: i) coexistence simulations with a reference system and then free energy correctio
ns between DFT and the reference system, and ii) direct DFT coexistence using simulation cells including 1000 atoms. The calculated zero pressure melting temperature is slightly underestimated at $1637 pm 10$ K (experimental value is 1728 K), and at high pressure is significantly higher than recent measurements in diamond anvil cell experiments [Phys. Rev. B {bf 87}, 054108 (2013)]. The zero pressure DFT melting slope is calculated to be $30 pm 2$ K, in good agreement with the experimental value of 28 K.
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 pres
sure-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 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 b
y 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 kinetics of hydrogen absorption by magnesium bulk is affected by two main activated processes: the dissociation of the H$_2$ molecule and the diffusion of atomic H into the bulk. In order to have fast absorption kinetics both activated processed
need to have a low barrier. Here we report a systematic ab-initio density functional theory investigation of H$_2$ dissociation and subsequent atomic H diffusion on TM(=Ti,V,Zr,Fe,Ru,Co,Rh,Ni,Pd,Cu,Ag)-doped Mg(0001) surfaces. The calculations show that doping the surface with TMs on the left of the periodic table eliminates the barrier for the dissociation of the molecule, but the H atoms bind very strongly to the TM, therefore hindering diffusion. Conversely, TMs on the right of the periodic table dont bind H, however, they do not reduce the barrier to dissociate H$_2$ significantly. Our results show that Fe, Ni and Rh, and to some extent Co and Pd, are all exceptions, combining low activation barriers for both processes, with Ni being the best possible choice.
