First-principles energetics of water: a many-body analysis

Standard forms of density-functional theory (DFT) have good predictive power for many materials, but are not yet fully satisfactory for solid, liquid and cluster forms of water. We use a many-body separation of the total energy into its 1-body, 2-body (2B) and beyond-2-body (B2B) components to analyze the deficiencies of two popular DFT approximations. We show how machine-learning methods make this analysis possible for ice structures as well as for water clusters. We find that the crucial energy balance between compact and extended geometries can be distorted by 2B and B2B errors, and that both types of first-principles error are important.

Assessment of density functional theory for iron(II) molecules across the spin-crossover transition

Octahedral Fe$^{2+}$ molecules are particularly interesting as they often exhibit a spin-crossover transition. In spite of the many efforts aimed at assessing the performances of density functional theory for such systems, an exchange-correlation functional able to account accurately for the energetic of the various possible spin-states has not been identified yet. Here we critically discuss the issues related to the theoretical description of this class of molecules from first principles. In particular we present a comparison between different density functionals for four ions, namely [Fe(H$_2$O)$_6$]$^{2+}$, [Fe(NH$_3$)$_6$]$^{2+}$, [Fe(NCH)$_6$]$^{2+}$ and [Fe(CO)$_6$]$^{2+}$. These are characterized by different ligand-field splittings and ground state spin multiplicities. Since no experimental data are available for the gas phase, the density functional theory results are benchmarked against those obtained with diffusion Monte Carlo, one of the most accurate methods available to compute ground state total energies of quantum systems. On the one hand, we show that most of the functionals considered provide a good description of the geometry and of the shape of the potential energy surfaces. On the other hand, the same functionals fail badly in predicting the energy differences between the various spin states. In the case of [Fe(H$_2$O)$_6$]$^{2+}$, [Fe(NH$_3$)$_6$]$^{2+}$, [Fe(NCH)$_6$]$^{2+}$, this failure is related to the drastic underestimation of the exchange energy. Therefore quite accurate results can be achieved with hybrid functionals including about 50% of Hartree-Fock exchange. In contrast, in the case of [Fe(CO)$_6$]$^{2+}$, the failure is likely to be caused by the multiconfigurational character of the ground state wave-function and no suitable exchange and correlation functional has been identified.

The kinetics of homogeneous melting beyond the limit of superheating

Molecular dynamics simulation is used to study the time-scales involved in the homogeneous melting of a superheated crystal. The interaction model used is an embedded-atom model for Fe developed in previous work, and the melting process is simulated in the microcanonical $(N, V, E)$ ensemble. We study periodically repeated systems containing from 96 to 7776 atoms, and the initial system is always the perfect crystal without free surfaces or other defects. For each chosen total energy $E$ and number of atoms $N$, we perform several hundred statistically independent simulations, with each simulation lasting for between 500 ps and 10 ns, in order to gather statistics for the waiting time $tau_{rm w}$ before melting occurs. We find that the probability distribution of $tau_{rm w}$ is roughly exponential, and that the mean value $<tau_{rm w} >$ depends strongly on the excess of the initial steady temperature of the crystal above the superheating limit identified by other researchers. The mean $<tau_{rm w}>$ also depends strongly on system size in a way that we have quantified. For very small systems of $sim 100$ atoms, we observe a persistent alternation between the solid and liquid states, and we explain why this happens. Our results allow us to draw conclusions about the reliability of the recently proposed Z method for determining the melting properties of simulated materials, and to suggest ways of correcting for the errors of the method.

Melting properties of a simple tight-binding model of transition metals: I.The region of half-filled d-band

We present calculations of the free energy, and hence the melting properties, of a simple tight-binding model for transition metals in the region of d-band filling near the middle of a d-series, the parameters of the model being designed to mimic molybdenum. The melting properties are calculated for pressures ranging from ambient to several Mbar. The model is intended to be the simplest possible tight-binding representation of the two basic parts of the energy: first, the pairwise repulsion due to Fermi exclusion; and second, the d-band bonding energy described in terms of an electronic density of states that depends on structure. In addition to the number of d-electrons, the model contains four parameters, which are adjusted to fit the pressure dependent d-band width and the zero-temperature pressure-volume relation of Mo. We show that the resulting model reproduces well the phonon dispersion relations of Mo in the body-centred-cubic structure, as well as the radial distribution function of the high-temperature solid and liquid given by earlier first-principles simulations. Our free-energy calculations start from the free energy of the liquid and solid phases of the purely repulsive pair-potential model, without d-band bonding. The free energy of the full tight-binding model is obtained from this by thermodynamic integration. The resulting melting properties of the model are quite close to those given by earlier first-principles work on Mo. An interpretation of these melting properties is provided by showing how they are related to those of the purely repulsive model.

Comment on Molybdenum at High Pressure and Temperature: Melting from Another Solid Phase

There has been a major controversy over the past seven years about the high-pressure melting curves of transition metals. Static compression (diamond-anvil cell: DAC) experiments up to the Mbar region give very low melting slopes dT_m/dP, but shock-wave (SW) data reveal transitions indicating much larger dT_m/dP values. Ab initio calculations support the correctness of the shock data. In a very recent letter, Belonoshko et al. propose a simple and elegant resolution of this conflict for molybdenum. Using ab initio calculations based on density functional theory (DFT), they show that the high-P/high-T phase diagram of Mo must be more complex than was hitherto thought. Their calculations give convincing evidence that there is a transition boundary between the normal bcc structure of Mo and a high-T phase, which they suggest could be fcc. They propose that this transition was misinterpreted as melting in DAC experiments. In confirmation, they note that their boundary also explains a transition seen in the SW data. We regard Belonoshko et al.s Letter as extremely important, but we note that it raises some puzzling questions, and we believe that their proposed phase diagram cannot be completely correct. We have calculated the Helmholtz and Gibbs free energies of the bcc, fcc and hcp phases of Mo, using essentially the same quasiharmonic methods as used by Belonoshko et al.; we find that at high-P and T Mo in the hcp structure is more stable than in bcc or fcc.

Zero-temperature generalized phase diagram of the 4d transition metals under pressure

We use an accurate implementation of density functional theory (DFT) to calculate the zero-temperature generalized phase diagram of the 4$d$ series of transition metals from Y to Pd as a function of pressure $P$ and atomic number $Z$. The implementation used is full-potential linearized augmented plane waves (FP-LAPW), and we employ the exchange-correlation functional recently developed by Wu and Cohen. For each element, we obtain the ground-state energy for several crystal structures over a range of volumes, the energy being converged with respect to all technical parameters to within $sim 1$ meV/atom. The calculated transition pressures for all the elements and all transitions we have found are compared with experiment wherever possible, and we discuss the origin of the significant discrepancies. Agreement with experiment for the zero-temperature equation of state is generally excellent. The generalized phase diagram of the 4$d$ series shows that the major boundaries slope towards lower $Z$ with increasing $P$ for the early elements, as expected from the pressure induced transfer of electrons from $sp$ states to $d$ states, but are almost independent of $P$ for the later elements. Our results for Mo indicate a transition from bcc to fcc, rather than the bcc-hcp transition expected from $sp$-$d$ transfer.

Structural properties and enthalpy of formation of magnesium hydride from quantum Monte Carlo calculations

We have used diffusion Monte Carlo (DMC) calculations to study the structural properties of magnesium hydride (MgH$_2$), including the pressure-volume equation of state, the cohesive energy and the enthalpy of formation from magnesium bulk and hydrogen gas. The calculations employ pseudopotentials and B-spline basis sets to expand the single particle orbitals used to construct the trial wavefunctions. Extensive tests on system size, time step, and other sources of errors, performed on periodically repeated systems of up to 1050 atoms, show that all these errors together can be reduced to below 10 meV per formula unit. We find excellent agreement with the experiments for the equilibrium volume of both the Mg and the MgH$_2$ crystals. The cohesive energy of the Mg crystal is found to be 1.51(1) eV, and agrees perfectly with the experimental value of 1.51 eV. The enthalpy of formation of MgH$_2$ from Mg bulk and H$_2$ gas is found to be $0.85 pm 0.01$ eV/formula unit, or $82 pm 1$ kJ/mole, which is off the experimental one of $76.1 pm 1$ kJ/mole only by 6 kJ/mole. This shows that DMC can almost achieve chemical accuracy (1 kcal/mole) on this system. Density functional theory errors are shown to be much larger, and depend strongly on the functional employed.

Hydrogen Dissociation and Diffusion on Ni and Ti -doped Mg(0001) Surfaces

It is well known, both theoretically and experimentally, that alloying MgH$_2$ with transition elements can significantly improve the thermodynamic and kinetic properties for H$_2$ desorption, as well as the H$_2$ intake by Mg bulk. Here we present a density functional theory investigation of hydrogen dissociation and surface diffusion over Ni-doped surface, and compare the findings to previously investigated Ti-doped Mg(0001) and pure Mg(0001) surfaces. Our results show that the energy barrier for hydrogen dissociation on the pure Mg(0001) surface is high, while it is small/null when Ni/Ti are added to the surface as dopants. We find that the binding energy of the two H atoms near the dissociation site is high on Ti, effectively impeding diffusion away from the Ti site. By contrast, we find that on Ni the energy barrier for diffusion is much reduced. Therefore, although both Ti and Ni promote H$_2$ dissociation, only Ni appears to be a good catalyst for Mg hydrogenation, allowing diffusion away from the catalytic sites. Experimental results corroborate these theoretical findings, i.e. faster hydrogenation of the Ni doped Mg sample as opposed to the reference Mg or Ti doped Mg.