Towards the high-throughput prediction of finite-temperature properties using the quasi-harmonic approximation.

Lattice dynamics calculations within the quasi-harmonic approximation (QHA) provide an infrastructure for modelling the finite-temperature properties of periodic solids at a modest computational cost. With the recent widespread interest in materials discovery by data mining, a database of computed finite-temperature properties would be highly desirable. In this work we provide a first step toward this goal with a comparative study of the accuracy of five exchange-correlation functionals, spanning the local density approximation (LDA), generalised-gradient approximation (GGA) and meta-GGA levels of theory, for predicting the properties of ten Group 1, 2 and 12 binary metal oxides. We find that the predictions are bounded by the LDA, which tends to underestimate lattice parameters and cell volumes relative to experiments, but yields the most accurate results for bulk moduli, expansion coefficients and Grüneisen parameters, and the PBE GGA, which shows the opposite behaviour. The PBEsol GGA gives the best overall predictions of the lattice parameters and volumes whilst also giving relatively reliable results for other properties. Our results demonstrate that, given a suitable choice of functional, a variety of finite-temperature properties can be predicted with useful accuracy, and hence that high-throughout QHA calculations are technically feasible.

Breathing Behaviour Modification of Gallium MIL-53 Metal-Organic Frameworks Induced by the Bridging Framework Inorganic Anion.

Controlling aspects of the µ2 -X- bridging anion in the metal-organic framework Ga-MIL-53 [GaX(bdc)] (X- =(OH)- or F- , bdc=1, 4-benzenedicarboxylate) is shown to direct the temperature at which thermally induced breathing transitions of this framework occur. In situ single crystal X-ray diffraction studies reveal that substituting 20 % of (OH)- in [Ga(OH)(bdc)] (1) for F- to produce [Ga(OH)0.8 F0.2 (bdc)] (2) stabilises the large pore (lp) form relative to the narrow pore (np) form, causing a well-defined decrease in the onset of the lp to np transition at higher temperatures, and the adsorption/desorption of nitrogen at lower temperatures through np to lp to intermediate (int) pore transitions. These in situ diffraction studies have also yielded a more plausible crystal structure of the int-[GaX(bdc)] ⋅ H2 O phases and shown that increasing the heating rate to a flash heating regime can enable the int-[GaX(bdc)] ⋅ H2 O to lp-[GaX(bdc)] transition to occur at a lower temperature than np-[GaX(bdc)] via an unreported pathway.

Phase stability of the tin monochalcogenides SnS and SnSe: a quasi-harmonic lattice-dynamics study.

The tin monochalcogenides SnS and SnSe adopt four different crystal structures, viz. orthorhombic Pnma and Cmcm and cubic rocksalt and π-cubic (P213) phases, each of which has optimal properties for a range of potential applications. This rich phase space makes it challenging to identify the conditions under which the different phases are obtained. We have performed first-principles quasi-harmonic lattice-dynamics calculations to assess the relative stabilities of the four phases of SnS and SnSe. We investigate dynamical stability through the presence or absence of imaginary modes in the phonon dispersion curves, and we compute Helmholtz and Gibbs free energies to evaluate the thermodynamic stability. We also consider applied pressures up to 15 GPa to obtain simulated temperature-pressure phase diagrams. Finally, the relationships between the orthorhombic crystal phases are investigated by explicitly mapping the potential-energy surfaces along the imaginary harmonic phonon modes in the Cmcm phase, and the relationships between the cubic phases are found by transition-state modelling using the climbing-image nudged elastic-band method.