First principles modelling of magnesium titanium hydrides

Mixing Mg with Ti leads to a hydride Mg(x)Ti(1-x)H2 with markedly improved (de)hydrogenation properties for x < 0.8, as compared to MgH2. Optically, thin films of Mg(x)Ti(1-x)H2 have a black appearance, which is remarkable for a hydride material. In this paper we study the structure and stability of Mg(x)Ti(1-x)H2, x= 0-1 by first-principles calculations at the level of density functional theory. We give evidence for a fluorite to rutile phase transition at a critical composition x(c)= 0.8-0.9, which correlates with the experimentally observed sharp decrease in (de)hydrogenation rates at this composition. The densities of states of Mg(x)Ti(1-x)H2 have a peak at the Fermi level, composed of Ti d states. Disorder in the positions of the Ti atoms easily destroys the metallic plasma, however, which suppresses the optical reflection. Interband transitions result in a featureless optical absorption over a large energy range, causing the black appearance of Mg(x)Ti(1-x)H2.Comment: 22 pages, 9 figures, 4 table

First principles simulations of liquid Fe-S under Earth's core conditions

First principles electronic structure calculations, based upon density functional theory within the generalized gradient approximation and ultra-soft Vanderbilt pseudopotentials, have been used to simulate a liquid alloy of iron and sulfur at Earth's core conditions. We have used a sulfur concentration of $\approx 12$wt, in line with the maximum recent estimates of the sulfur abundance in the Earth's outer core. The analysis of the structural, dynamical and electronic structure properties has been used to report on the effect of the sulfur impurities on the behavior of the liquid. Although pure sulfur is known to form chains in the liquid phase, we have not found any tendency towards polymerization in our liquid simulation. Rather, a net S-S repulsion is evident, and we propose an explanation for this effect in terms of the electronic structure. The inspection of the dynamical properties of the system suggests that the sulfur impurities have a negligible effect on the viscosity of Earth's liquid core.Comment: 24 pages (including 8 figures

The Structure, Dynamics and Electronic Structure of Liquid Ag-Se Alloys Investigated by Ab Initio Simulation

Ab initio molecular-dynamics simulations have been used to investigate the structure, dynamics and electronic properties of the liquid alloy Ag(1-x)Se(x) at 1350 K and at the three compositions x=0.33, 0.42 and 0.65. The calculations are based on density-functional theory in the local density approximation and on the pseudopotential plane-wave method. The reliability of the simulations is confirmed by detailed comparisons with very recent neutron diffraction results for the partial structure factors and radial distribution functions (RDF) of the stoichiometric liquid Ag2Se. The simulations show a dramatic change of the Se-Se RDF with increasing Se content. This change is due to the formation of Se clusters bound by covalent bonds, the Se-Se bond length being almost the same as in pure c-Se and l-Se. The clusters are predominantly chain-like, but for higher x a large fraction of 3-fold coordinated Se atoms is also found. It is shown that the equilibrium fractions of Se present as isolated atoms and in clusters can be understood on a simple charge-balance model based on an ionic interpretation. The Ag and Se diffusion coefficients both increase with Se content, in spite of the Se clustering. An analysis of the Se-Se bond dynamics reveals surprisingly short bond lifetimes of less than 1 ps. The changes in the density of states with composition arise directly from the formation of Se-Se covalent bonds. Results for the electronic conductivity obtained using the Kubo-Greenwood approximation are in adequate agreement with experiment for l-Ag2Se, but not for the high Se contents. Possible reasons for this are discussed.Comment: 14 pages, Revtex, 14 Postscript figures embedded in the tex

Surface Structure of Liquid Metals and the Effect of Capillary Waves: X-ray Studies on Liquid Indium

We report x-ray reflectivity (XR) and small angle off-specular diffuse scattering (DS) measurements from the surface of liquid Indium close to its melting point of $156^\circ$C. From the XR measurements we extract the surface structure factor convolved with fluctuations in the height of the liquid surface. We present a model to describe DS that takes into account the surface structure factor, thermally excited capillary waves and the experimental resolution. The experimentally determined DS follows this model with no adjustable parameters, allowing the surface structure factor to be deconvolved from the thermally excited height fluctuations. The resulting local electron density profile displays exponentially decaying surface induced layering similar to that previously reported for Ga and Hg. We compare the details of the local electron density profiles of liquid In, which is a nearly free electron metal, and liquid Ga, which is considerably more covalent and shows directional bonding in the melt. The oscillatory density profiles have comparable amplitudes in both metals, but surface layering decays over a length scale of $3.5\pm 0.6$ \AA for In and $5.5\pm 0.4$ \AA for Ga. Upon controlled exposure to oxygen, no oxide monolayer is formed on the liquid In surface, unlike the passivating film formed on liquid Gallium.Comment: 9 pages, 5 figures; submitted to Phys. Rev.

First principles simulations of direct coexistence of solid and liquid aluminium

First principles calculations based on density functional theory, with generalised gradient corrections and ultrasoft pseudopotentials, have been used to simulate solid and liquid aluminium in direct coexistence at zero pressure. Simulations have been carried out on systems containing up to 1000 atoms for 15 ps. The points on the melting curve extracted from these simulations are in very good agreement with previous calculations, which employed the same electronic structure method but used an approach based on the explicit calculation of free energies [L. Vo\v{c}adlo and D. Alf\`e, Phys. Rev. B, {\bf 65}, 214105 (2002).]Comment: To appear in Phys. Rev.

Composition Dependence of the Structure and Electronic Properties of Liquid Ga-Se Alloys Studied by Ab Initio Molecular Dynamics Simulation

Ab initio molecular dynamics simulation is used to study the structure and electronic properties of the liquid Ga-Se system at the three compositions Ga$_2$Se, GaSe and Ga$_2$Se$_3$, and of the GaSe and Ga$_2$Se$_3$ crystals. The calculated equilibrium structure of GaSe crystal agrees well with available experimental data. The neutron-weighted liquid structure factors calculated from the simulations are in reasonable agreement with recent neutron diffraction measurements. Simulation results for the partial radial distribution functions show that the liquid structure is closely related to that of the crystals. A close similarity between solid and liquid is also found for the electronic density of states and charge density. The calculated electronic conductivity decreases strongly with increasing Se content, in accord with experimental measurements.Comment: REVTeX, 8 pages and 12 uuencoded PostScript figures, submitted to Phys. Rev. B. corresponding author: [email protected]

Iron under Earth's core conditions: Liquid-state thermodynamics and high-pressure melting curve

{\em Ab initio} techniques based on density functional theory in the projector-augmented-wave implementation are used to calculate the free energy and a range of other thermodynamic properties of liquid iron at high pressures and temperatures relevant to the Earth's core. The {\em ab initio} free energy is obtained by using thermodynamic integration to calculate the change of free energy on going from a simple reference system to the {\em ab initio} system, with thermal averages computed by {\em ab initio} molecular dynamics simulation. The reference system consists of the inverse-power pair-potential model used in previous work. The liquid-state free energy is combined with the free energy of hexagonal close packed Fe calculated earlier using identical {\em ab initio} techniques to obtain the melting curve and volume and entropy of melting. Comparisons of the calculated melting properties with experimental measurement and with other recent {\em ab initio} predictions are presented. Experiment-theory comparisons are also presented for the pressures at which the solid and liquid Hugoniot curves cross the melting line, and the sound speed and Gr\"{u}neisen parameter along the Hugoniot. Additional comparisons are made with a commonly used equation of state for high-pressure/high-temperature Fe based on experimental data.Comment: 16 pages including 6 figures and 5 table

Structure and Dynamics of Liquid Iron under Earth's Core Conditions

First-principles molecular dynamics simulations based on density-functional theory and the projector augmented wave (PAW) technique have been used to study the structural and dynamical properties of liquid iron under Earth's core conditions. As evidence for the accuracy of the techniques, we present PAW results for a range of solid-state properties of low- and high-pressure iron, and compare them with experimental values and the results of other first-principles calculations. In the liquid-state simulations, we address particular effort to the study of finite-size effects, Brillouin-zone sampling and other sources of technical error. Results for the radial distribution function, the diffusion coefficient and the shear viscosity are presented for a wide range of thermodynamic states relevant to the Earth's core. Throughout this range, liquid iron is a close-packed simple liquid with a diffusion coefficient and viscosity similar to those of typical simple liquids under ambient conditions.Comment: 13 pages, 8 figure

Half-metallicity and Slater-Pauling behavior in the ferromagnetic Heusler alloys

Introductory chapter for the book "Halfmetallic Alloys - Fundamentals and Applications" to be published in the series Springer Lecture Notes on Physics, P. H. Dederichs and I. Galanakis (eds). It contains a review of the theoretical work on the half-metallic Heusler alloys.Comment: Introductory chapter for the book "Halfmetallic Alloys - Fundamentals and Applications" to be published in the series Springer Lecture Notes on Physics, P. H. Dederichs and I. Galanakis (eds