Dr Joseph Prentice

Abstract:

We demonstrate the capability of embedded mean-field theory (EMFT) within the linear-scaling density-functional-theory code ONETEP, which enables DFT-in-DFT quantum embedding calculations on systems containing thousands of atoms at a fraction of the cost of a full calculation. We perform simulations on a wide range of systems from molecules to complex nanostructures to demonstrate the performance of our implementation with respect to accuracy and efficiency. This work paves the way for the application of this class of quantum embedding method to large-scale systems that are beyond the reach of existing implementations.

Abstract:

Calcium silicate perovskite (CaSiO3) is one of the major mineral components of the lower mantle, but has been the subject of relatively little work compared to the more abundant Mg-based materials. One of the major problems related to CaSiO3 that is still the subject of research is its crystal structure under lower mantle conditions – a cubic Pm¯3m structure is accepted in general, but some have suggested that lower-symmetry structures may be relevant. In this work, we use a fully first principles vibrational self-consistent field (VSCF) method to perform high accuracy anharmonic vibrational calculations on several candidate structures at a variety of points along the geotherm near the base of the lower mantle, in order to investigate the stability of the cubic structure and related distorted structures. Our results show that the cubic structure is the most stable throughout the lower mantle, and that this result is robust against the effects of thermal expansion.

Abstract:

High-level vibrational calculations have been used to investigate anharmonicity in a wide variety of materials using density-functional-theory (DFT) methods. We have developed a new and efficient approach for describing strongly-anharmonic systems using a vibrational self-consistent-field (VSCF) method. By far the most computationally expensive part of the calculations is the mapping of an accurate Born-Oppenheimer (BO) energy surface within the region of interest. Here we present an improved method which reduces the computational cost of the mapping. In this approach we use data from a set of energy calculations for different vibrational distortions of the materials and the corresponding forces on the atoms. Results using both energies and forces are presented for the test cases of the hydrogen molecule, solid hydrogen under high pressure including mapping of two-dimensional subspaces of the BO surface, and the bcc phases of the metals Li and Zr. The use of forces data speeds up the anharmonic calculations by up to 40%.