Source author record

Jonathan M. Skelton

Jonathan M. Skelton appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
1topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2022arXiv

The Physical Significance of Imaginary Phonon Modes in Crystals

The lattice vibrations (phonon modes) of crystals underpin a large number of material properties. The harmonic phonon spectrum of a solid is the simplest description of its structural dynamics and can be straightforwardly derived from the Hellman-Feynman forces obtained in a ground-state electronic structure calculation. The presence of imaginary harmonic modes in the spectrum indicates that a structure is a local maximum on the athermal structural potential-energy surface and can yield important insight into the fundamental nature and physical properties of a material. In this review article, we discuss the physical significance of imaginary harmonic modes and distinguish between cases where imaginary modes are indicative of such phenomena, and those where they reflect technical problems in the calculations. We outline basic approaches for exploring and renormalising imaginary modes, and demonstrate their utility through a set of three case studies in the materials sciences.

preprint2020arXiv

Accuracy of hybrid functionals with non-self-consistent Kohn-Sham orbitals for predicting the properties of semiconductors

Accurately modeling the electronic structure of materials is a persistent challenge to high-throughput screening. A promising means of balancing accuracy against computational cost are non-self-consistent calculations with hybrid density-functional theory, where the electronic band energies are evaluated using a hybrid functional from orbitals obtained with a less demanding (semi-)local functional. We have quantified the performance of this technique for predicting the physical properties of sixteen tetrahedral semiconductors with bandgaps from 0.2-5.5 eV. Provided the base functional predicts a non-metallic electronic structure, bandgaps within 5 % of the PBE0 and HSE06 gaps can be obtained with an order of magnitude reduction in computing time. The positions of the valence and conduction band extrema and the Fermi level are well reproduced, further enabling calculation of the band dispersion, density of states, and dielectric properties using Fermi's Golden Rule. While the error in the non-self-consistent total energies is ~50 meV atom-1, the energy-volume curves are reproduced accurately enough to obtain the equilibrium volume and bulk modulus with minimal error. We also test the dielectric-dependent scPBE0 functional and obtain bandgaps and dielectric constants to within 2.5 % of the self-consistent results, which amount to a significant improvement over self-consistent PBE0 for a similar computational cost. We identify cases where the non-self-consistent approach is expected to perform poorly, and demonstrate that partial self-consistency provides a practical and efficient workaround. Finally, we perform proof-of-concept calculations on CoO and NiO to demonstrate the applicability of the technique to strongly-correlated open-shell transition-metal oxides.

preprint2020arXiv

Lattice dynamics of $\textit{Pnma}$ Sn(S$_{1-x}$Se$_{x}$) solid solutions: energetics, phonon spectra and thermal transport

Alloying is widely used as a means to fine-tune the properties of thermoelectric materials by reducing the lattice thermal conductivity. However, the effects of compositional variation on the lattice dynamics of alloy systems are not well understood, due in part to the difficulty of building realistic first-principles models of structurally-complex solid solutions. This work builds on our previous study of Sn$_{n}$(S$_{1-x}$Se$_{x}$)$_{m}$ solid solutions [Gunn $\textit{et al.}$, $\textit{Chem. Mater.}$ $\textbf{31}$, $\textit{10}$, 3672, $\textbf{2019}$] to explore the lattice dynamics of the $\textit{Pnma}$ Sn(S$_{1-x}$Se$_{x}$) system, which has been widely studied for potential thermoelectric applications. We find that the vibrational internal energy and entropy have a large quantitative impact on the mixing free energy and are likely to be particularly important in alloy systems with competing phases. The thermodynamically-averaged phonon dispersions and density of states curves show that alloying preserves the structure of the low-frequency bands of modes associated with the Sn sublattice but broadens the high-frequency chalcogen bands into a near-continuous spectrum at the 50/50 mixed composition. This results in a general reduction in the phonon mode group velocities and an increase in the number of energy-conserving scattering channels for heat-carrying low-frequency modes, which is consistent with the decrease in thermal conductivity observed in experimental measurements. Finally, we discuss some of the limitations of our first-principles modelling approach and propose methods to address these in future studies.

preprint2016arXiv

Anharmonicity in the high-temperature Cmcm phase of SnSe: soft modes and three-phonon interactions

The layered semiconductor SnSe is one of the highest-performing thermoelectric materials known. We demonstrate, through a first-principles lattice-dynamics study, that the high-temperature Cmcm phase is a dynamic average over lower-symmetry minima separated by very small energetic barriers. Compared to the low-temperature Pnma phase, the Cmcm phase displays a phonon softening and enhanced three-phonon scattering, leading to an anharmonic damping of the low-frequency modes and hence the thermal transport. We develop a renormalisation scheme to quantify the effect of the soft modes on the calculated properties, and confirm that the anharmonicity is an inherent feature of the Cmcm phase. These results suggest a design concept for thermal insulators and thermoelectric materials, based on displacive instabilities, and highlight the power of lattice-dynamics calculations for materials characterization.

preprint2016arXiv

Direct Observation of Dynamic Symmetry Breaking above Room Temperature in Methylammonium Lead Iodide Perovskite

Lead halide perovskites such as methylammonium lead triiodide (MAPI) have outstanding optical and electronic properties for photovoltaic applications, yet a full understanding of how this solution processable material works so well is currently missing. Previous research has revealed that MAPI possesses multiple forms of static disorder regardless of preparation method, which is surprising in light of its excellent performance. Using high energy resolution inelastic X-ray (HERIX) scattering, we measure phonon dispersions in MAPI and find direct evidence for another form of disorder in single crystals: large amplitude anharmonic zone-edge rotational instabilities of the PbI_6 octahedra that persist to room temperature and above, left over from structural phase transitions that take place tens to hundreds of degrees below. Phonon calculations show that the orientations of the methylammonium couple strongly and cooperatively to these modes. The result is a non-centrosymmetric, instantaneous local structure, which we observe in atomic pair distribution function (PDF) measurements. This local symmetry breaking is unobservable by Bragg diffraction, but can explain key material properties such as the structural phase sequence, ultra low thermal transport, and large minority charge carrier lifetimes despite moderate carrier mobility.

preprint2016arXiv

Free energy of ligand removal in the metal-organic framework UiO-66

We report an investigation of the "missing-linker phenomenon" in the Zr-based metal-organic framework UiO-66 using atomistic forcefield and quantum chemical methods. For a vacant benzene dicarboxylate ligand, the lowest energy charge capping mechanism involves acetic acid or Cl-/H2O. The calculated defect free energy of formation is remarkably low, consistent with the high defect concentrations reported experimentally. A dynamic structural instability is identified for certain higher defect concentrations. In addition to the changes in material properties upon defect formation, we assess the formation of molecular aggregates, which provide an additional driving force for ligand loss. These results are expected to be of relevance to a wide range of metal-organic frameworks.

preprint2016arXiv

Suppression of lattice thermal conductivity by mass-conserving cation mutation in multi-component semiconductors

In semiconductors almost all heat is conducted by phonons (lattice vibrations), which is limited by their quasi-particle lifetimes. Phonon-phonon interactions represent scattering mechanisms that produce thermal resistance. In thermoelectric materials, this resistance due to anharmonicity should be maximised for optimal performance. We use a first-principles lattice-dynamics approach to explore the changes in lattice dynamics across an isostructural series where the average atomic mass is conserved: ZnS to CuGaS$_2$ to Cu$_2$ZnGeS$_4$. Our results demonstrate an enhancement of phonon interactions in the multernary materials, and confirm that lattice thermal conductivity can be controlled independently of the average mass and local coordination environments.

preprint2015arXiv

Crystal structure optimisation using an auxiliary equation of state

Standard procedures for local crystal-structure optimisation involve numerous energy and force calculations. It is common to calculate an energy-volume curve, fitting an equation of state around the equilibrium cell volume. This is a computationally intensive process, in particular for low-symmetry crystal structures where each isochoric optimisation involves energy minimisation over many degrees of freedom. Such procedures can be prohibitive for non-local exchange-correlation functionals or other 'beyond' density functional theory electronic structure techniques, particularly where analytical gradients are not available. We present a simple approach for efficient optimisation of crystal structures based on a known equation of state. The equilibrium volume can be predicted from one single-point calculation, and refined with successive calculations if required. The approach is validated for PbS, PbTe, ZnS and ZnTe using nine density functionals, and applied to the quaternary semiconductor Cu$_{2}$ZnSnS$_{4}$ and the magnetic metal-organic framework HKUST-1.

preprint2015arXiv

Lattice dynamics and vibrational spectra of the orthorhombic, tetragonal and cubic phases of methylammonium lead iodide

The hybrid halide perovskite CH3NH3PbI3 exhibits a complex structural behaviour, with successive transitions between orthorhombic, tetragonal and cubic polymorphs at ca. 165 K and 327 K. Herein we report first-principles lattice dynamics (phonon spectrum) for each phase of CH3NH3PbI3. The equilibrium structures compare well to solutions of temperature-dependent powder neutron diffraction. By following the normal modes we calculate infrared and Raman intensities of the vibrations, and compare them to the measurement of a single crystal where the Raman laser is controlled to avoid degradation of the sample. Despite a clear separation in energy between low frequency modes associated with the inorganic PbI3 network and high-frequency modes of the organic CH3NH3+ cation, significant coupling between them is found, which emphasises the interplay between molecular orientation and the corner-sharing octahedral networks in the structural transformations. Soft modes are found at the boundary of the Brillouin zone of the cubic phase, consistent with displacive instabilities and anharmonicity involving tilting of the PbI6 octahedra around room temperature.

preprint2015arXiv

Solid-state chemistry of glassy antimony oxides

The amorphous phases of antimony oxide at three different oxidation levels have been investigated by {\it ab initio} molecular-dynamics calculations using a simulated `melt-quench' approach. The atomic and electronic structure of the amorphous phases are analyzed and compared to their crystalline counterparts. The amorphous structure of the trioxide ({\it a\,}--Sb$_2$O$_3$) resembles the $β$-phase of the crystalline form (valentinite), in agreement with previous observations. In the relaxed athermal structure, however, more senarmontite-like structural features are apparent. The phase diagram of the amorphous phases with respect to oxygen chemical potential has been calculated within the framework of {\it ab initio} thermodynamics. At elevated oxidation levels, the resulting tetraoxide and pentoxide materials show distinct variation in electronic structure compared to the trioxide, with the electronic density of states indicating a narrowing of the electronic gap with increasing oxidation level.