Monthly Archives: December 2013

Anharmonic vibrational properties of solids

Bartomeu Monserrat
, Neil Drummond, Richard Needs

The vibrational properties of solids are usually studied within the harmonic approximation, which is valid when the motion of the atomic nuclei is restricted to the neighbourhood of their equilibrium positions. This approximation is usually very good when investigating standard solids under normal conditions. However, for the lighter nuclei and at high temperatures, atoms explore regions far away from their equilibrium positions, and the harmonic approximation is no longer valid.

We have worked on a methodology to study the vibrational properties of solids beyond the harmonic approximation. We map the Born-Oppenheimer energy surface on which the atoms move along phonon modes, and then solve the resulting equations self-consistently within a mean-field formalism. We then calculate phonon expectation values from the anharmonic phonon wave function. This allows us to investigate physical properties such as the electronic band gap renormalisation due to electron-phonon coupling, or thermal expansion using the stress tensor [PRB 87, 144302 (2013)].

The first applications of this methodology have included the study of systems comprising of the lightest elements: hydrogen and helium. For hydrogen we have investigated the dissociation of molecular solid hydrogen at high pressure with quantum Monte Carlo calculations for the electronic part, and anharmonic vibrational calculations for the nuclear part. We have predicted a dissociation pressure of 370 GPa, within reach of static diamond anvil cell experiments in the near future. This work has been done in collaboration with Sam Azadi and Matthew Foulkes. For helium, we have calculated the effects of electron-phonon coupling on the metallization transition, predicting a strong zero-point correction and temperature dependence, both increasing the static lattice metallization pressure. This result has implications in the area of white dwarf cooling. This work has been done in collaboration with Chris Pickard.

Discovering new crystalline forms of metal oxides


Matthew J. Lyle, Chris J. Pickard, and Richard J. Needs

High performance computing is an essential part of many areas of scientific research. Materials discovery, in particular, is especially suited to computational simulation. Here we use an ab-initio random structure searching (AIRSS) algorithm to generate thousands of possible atomic arrangements for our system of choice. Constraints such as symmetry and stoichiometry can be incorporated to help our search towards particular structures we may be after. These candidate structures are energetically ranked, usually using density functional theory (DFT), and low-energy or other interesting structures are further investigated.

My projects have focussed on finding new crystalline forms of metal oxides, including alumina (Al2O3) and titania (TiO2). Alumina exists in a number of distinct crystalline forms that are widely used in industry; annual global production is 100 Mtonne, valued at approximately U.S.$32 billion. Conventional crystalline aluminas exhibit negligible porosity, though many of their industrial applications require high porosities and surface areas. This is usually achieved through the fashioning of mesoporous channels, though often at the expense of crystallinity which is known to increase alumina stability. Moreover, crystalline microporous silicates and aluminosilicates, e.g. zeolites, have long been used as adsorbents, ion exchangers, and catalysts.


In this project we performed an extensive first-principles computational search for new crystalline forms of alumina that are stable at ambient pressure. We found 147 unique structures with energies intermediate between those of the conventional alumina phases and identified a new class of crystalline aluminas which are up to 43% less dense. We attribute these low densities to significant amounts (up to 100%) of fivefold Al coordination and the formation of highly microporous zeolite-like channels. These are the first crystalline aluminas exhibiting such extensive fivefold Al coordination and surface areas, suggesting a new paradigm for the processing of alumina for industrial use.


MAGE: Materials Automated Generation and Exploration

Gareth Conduit

Through the stone, bronze, and iron ages the discovery of new materials has chronicled human history. The coming of each age was sparked by a chance discovery of a new material. Despite the central importance of materials in enabling new technologies, even today the only way to develop new materials is through experiment driven trial and error. We have developed the first tool [Patent GB1302743.8], MAGE, Materials Automated Generation and Exploration: computationally designs materials with specified physical properties. We have proven its accuracy by predicting five new alloys [Patent GB1307533.8, Patent GB1307535.3, Rolls-Royce Group plc invention submissions NC12261, NC13006 & Acta Materialia, 61, 3378 (2013)] that were subsequently experimentally verified. MAGE can help not only materials scientists, but moreover design engineers. At present engineers must design new objects and products around the shortcomings of pre-existing but non-ideal materials. The capability to develop materials computationally would allow engineers to instantly optimize bespoke materials for their application, bringing materials into the heart of the design process.

The alloys developed include a new nickel-based alloy for turbine discs in jet engines, a novel alloy for lining the combustion chamber of a jet engine, and a molybdenum based forging alloy. Each alloy has eight individual physical properties that are predicted to match or exceed commercially available alternatives including fracture toughness, oxidation resistance, yield strength, creep resistance and processibility. These properties are calculated using a variety of techniques including ab initio calculations, physically based models, and interpolation of existing experimental data. Several properties for each alloy have subsequently been experimentally verified so the alloys are now undergoing compliance testing by Rolls Royce plc.

In the future we plan to further incorporate first principles methods into MAGE. Working with Samsung Electronics Co Ltd through a Global Research Outreach Grant we will first focus on designing new GaInN-based LEDs materials, exploiting Density Functional Theory to predict properties from first principles. This investigation should allow us to not only design potentially important semiconductor materials, but moreover act as a launch-pad to guide the future design of new materials from first principles.

A Stochastic Approach to Battery Design

Martin Mayo, Hugh Glass, Andrew J. Morris.

LIBs are the rechargeable battery of choice for electric vehicles and portable electronic devices. There is substantial interest in enhancing the capacity of LIBs, driven by the economic and environmental advantages of increasing the range of electric vehicles, and enabling longer-life portable electronic devices. The traditional anode is composed of graphite but recently silicon has been proposed as an alternative which has a theoretical capacity some ten-times larger. Silicon anodes are not fully understood however and suffer from large volume expansion on charge and hysteresis on charge cycling. Phosphorous is another such candidate for high-capacity anodes due to its multiple allotropes (white, red and black phosphorus), its ability to make Hittorf-chains and novel structures when mixed with lithium.[1]

P21/c Li1Si1

The creation of new materials is both difficult and expensive. It is very difficult to “see” the structure of these materials over the length scales that they work and very expensive to create prototype materials to test. Whilst experimental physics can use x-rays, high-energy electrons or neutrons to infer the structure of these materials, this inference is made LiP double helix infinite chainmuch more robust when combined with theoretical predictions of the kinds of structures that can be formed. In a computer we use quantum mechanical calculations to simulate the results of these kinds of experiments, helping to understand materials and suggest new materials with the kind of properties desired. The theoretical prediction of even very simple structures has, until recently, been out of bounds due to the large number of possible atomic arrangements. The ab initio random structure searching method (AIRSS) [2] uses a stochastic approach to suggest different structural configurations of atoms within a material. By searching over a range of stoichiometries it is possible to model how a battery is charged as we can predict the structural changes of the electrodes as a charging potential is applied. Using this knowledge it is possible to perform further theoretical analysis such as NMR and EELS (electron energy loss) spectroscopies and predict charge and discharge voltages.

[1] “Inorganic Double Helix Structures of Unusually Simple Li-P Species”, Alexander S. Ivanov, Andrew J. Morris, Konstantin V. Bozhenko, Chris J. Pickard and Alexander I. Boldyrev, [Angewandte Chemie 51 33, 8330-8333 (2012)]

[2] C.J. Pickard, and R.J. Needs , Phys. Rev. Lett. 97, 045504 (2006)

Quantum Monte Carlo simulations of solid hydrogen


Jonathan Lloyd-Williams, Bartomeu Monserrat, Pablo López Ríos, Neil Drummond, and Richard Needs

Hydrogen has the simplest electronic structure of any atom but its bulk properties are surprisingly complex. As hydrogen can only form a single covalent bond, it is expected to remain molecular to very high pressures. Several solid phases of hydrogen have been observed; the low-pressure phase I which consists of a quantum crystal of rotating molecules on a hexagonal close-packed lattice transforms at pressures of about 110 GPa to the broken symmetry phase II, in which the mean molecular orientations are ordered, and then to phase III at about 150 GPa. A new phase IV was recently observed in room temperature experiments and has generated great excitement within the field. Despite many years of intensive study, the arrangements of the molecules in phases II and III remain undetermined. The structure of phase IV is also unknown.

Experimental measurements have only been able to provide limited information about the molecular orientations in phases II and III of hydrogen. The information obtained from X-ray diffraction is limited because of the weak scattering of hydrogen atoms. So far the most important experimental data have come from infra-red and Raman vibrational spectroscopy. Theoretical identification of the structures is also difficult because of the need to consider the quantum statistics of the nuclei, the large zero-point vibrational energies, and the small enthalpy differences between competing phases.

To make progress in understanding the high pressure phase diagram of hydrogen, we will perform diffusion Monte Carlo (DMC) calculations for the candidate phases (P21/c-24, C2/c-24, Cmca-12, Cmca-4, and Pc-48) at several different densities which will allow us to evaluate the enthalpy. DMC is a zero-temperature method but estimating the small room temperature corrections should be straightforward. We will first perform a set of static lattice DMC calculations and add zero-point motion effects from density functional theory to allow some preliminary comparisons with experiment. We will then perform a second set of calculations including the full zero-point motion effects within DMC.

Project begun: January 2013

Expected completion: December 2014