The first principles calculation of transport coefficients

We demonstrate the practical feasibility of calculating transport coefficients such as the viscosity of liquids completely from first principles using the Green-Kubo relations. Results presented for liquid aluminum are shown to have a statistical error of only ca. 5%. The importance of such calculations is illustrated by results for a liquid iron-sulfur alloy under Earth's core conditions, which indicate that the viscosity of the liquid outer core is not substantially higher than that of typical liquid metals under ambient conditions.Comment: four pages, including four figure

Ab initio statistical mechanics of surface adsorption and desorption: II. Nuclear quantum effects

We show how the path-integral formulation of quantum statistical mechanics can be used to construct practical {\em ab initio} techniques for computing the chemical potential of molecules adsorbed on surfaces, with full inclusion of quantum nuclear effects. The techniques we describe are based on the computation of the potential of mean force on a chosen molecule, and generalise the techniques developed recently for classical nuclei. We present practical calculations based on density functional theory with a generalised-gradient exchange-correlation functional for the case of H$_2$O on the MgO~(001) surface at low coverage. We note that the very high vibrational frequencies of the H$_2$O molecule would normally require very large numbers of time slices (beads) in path-integral calculations, but we show that this requirement can be dramatically reduced by employing the idea of thermodynamic integration with respect to the number of beads. The validity and correctness of our path-integral calculations on the H$_2$O/MgO~(001) system are demonstrated by supporting calculations on a set of simple model systems for which quantum contributions to the free energy are known exactly from analytic arguments.Comment: 11 pages, including 2 figure

The energetics of water on oxide surfaces by quantum Monte Carlo

Density functional theory (DFT) is widely used in surface science, but gives poor accuracy for oxide surface processes, while high-level quantum chemistry methods are hard to apply without losing basis-set quality. We argue that quantum Monte Carlo techniques allow these difficulties to be overcome, and we present diffusion Monte Carlo results for the formation energy of the MgO(001) surface and the adsorption energy of H$_2$O on this surface, using periodic slab geometry. The results agree well with experiment. We note other oxide surface problems where these techniques could yield immediate progress.Comment: 5 pages, 2 figure

Ivermectin - old drug, new tricks?

Ivermectin is one of the most important drugs in veterinary and human medicine for the control of parasitic infection and was the joint focus of the 2015 Nobel Prize in Physiology or Medicine, some 35 years after its remarkable discovery. Although best described for its activity on glutamate-gated chloride channels in parasitic nematodes, understanding of its mode of action remains incomplete. In the field of veterinary medicine, resistance to ivermectin is now widespread, but the mechanisms underlying resistance are unresolved. Here we discuss the history of this versatile drug and its use in global health. Based on recent studies in a variety of systems, we question whether ivermectin could have additional modes of action on parasitic nematodes

Thermodynamics of hexagonal-close-packed iron under Earth’s core conditions

The free energy and other thermodynamic properties of hexagonal-close-packed iron are calculated by direct ab initio methods over a wide range of pressures and temperatures relevant to the Earth’s core. The ab initio calculations are based on density-functional theory in the generalized-gradient approximation, and are performed using the projector augmented wave approach. Thermal excitation of electrons is fully included. The Helmholtz free energy consists of three parts, associated with the rigid perfect lattice, harmonic lattice vibrations, and anharmonic contributions, and the technical problems of calculating these parts to high precision are investigated. The harmonic part is obtained by computing the phonon frequencies over the entire Brillouin zone, and by summation of the free-energy contributions associated with the phonon modes. The anharmonic part is computed by the technique of thermodynamic integration using carefully designed reference systems. Detailed results are presented for the pressure, specific heat, bulk modulus, expansion coefficient and Grüneisen parameter, and comparisons are made with values obtained from diamond-anvil-cell and shock experiments

Ab-initio simulation of high-temperature liquid selenium

Ab initio molecular dynamics simulation is used to investigate the structure and dynamics of liquid Se at temperatures of 870 and 1370~K. The calculated static structure factor is in excellent agreement with experimental data. The calculated radial distribution function gives a mean coordination number close to 2, but we find a significant fraction of one-fold and three-fold atoms, particularly at 1370~K, so that the chain structure is considerably disrupted. The self-diffusion coefficient has values ($\sim 1 \times 10^{-8}$~m~s$^{-1}$) typical of liquid metals.Comment: 10 pages, 4 Poscript figures, uses REVTE

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.Comment: 1 page, 1 figure. submitted to Phys. Rev. Let