-
Theoretical study of the stability of the tetradymite-like phases of Sb$_2$S$_3$, Bi$_2$S$_3$, and Sb$_2$Se$_3$
Authors:
E. Lora da Silva,
J. M. Skelton,
P. Rodríguez-Hernández,
A. Muñoz,
M. C. Santos,
D. Martínez-García,
R. Vilaplana,
F. J. Manjón
Abstract:
We report a comparative theoretical study of the \textit{Pnma} and \textit{R-3m} phases of Sb$_2$S$_3$, Bi$_2$S$_3$, and Sb$_2$Se$_3$ close to ambient pressure.
Our enthalpy calculations at 0 K show that at ambient pressure the \textit{R-3m} (tetradymite-like) phase of Sb$_2$Se$_3$ is energetically more stable than the \textit{Pnma} phase, contrary to what is observed for Sb$_2$S$_3$ and Bi$_2$S…
▽ More
We report a comparative theoretical study of the \textit{Pnma} and \textit{R-3m} phases of Sb$_2$S$_3$, Bi$_2$S$_3$, and Sb$_2$Se$_3$ close to ambient pressure.
Our enthalpy calculations at 0 K show that at ambient pressure the \textit{R-3m} (tetradymite-like) phase of Sb$_2$Se$_3$ is energetically more stable than the \textit{Pnma} phase, contrary to what is observed for Sb$_2$S$_3$ and Bi$_2$S$_3$, and irrespective of the exchange-correlation functional employed in the calculations.
The result for Sb$_2$Se$_3$ is in contradiction to experiments where all three compounds are usually grown in the \textit{Pnma} phase.
This result is further confirmed by free-energy calculations taking into account the temperature dependence of the unit-cell volumes and phonon frequencies.
Lattice dynamics and elastic tensor calculations further show that both \textit{Pnma} and \textit{R-3m} phases of Sb$_2$Se$_3$ are dynamically and mechanically stable at zero applied pressure.
Since these results suggest that the formation of the \textit{R-3m} phase for Sb$_2$Se$_3$ should be feasible at close to ambient conditions, we provide a theoretical crystal structure and simulated Raman and infrared spectra to help in its identification. We also discuss the results of the two published works that have claimed to have synthesized tetradymite-like Sb$_2$Se$_3$.
Finally, the stability of the \textit{R-3m} phase across the three group-15 \textit{A$_2$X$_3$} sesquichalcogenides is analysed based on their van der Waals gap and X-X in-plane geometry.
△ Less
Submitted 2 June, 2022;
originally announced June 2022.
-
The Physical Significance of Imaginary Phonon Modes in Crystals
Authors:
Ioanna Pallikara,
Prakriti Kayastha,
Jonathan M. Skelton,
Lucy D. Whalley
Abstract:
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 str…
▽ More
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.
△ Less
Submitted 2 March, 2022;
originally announced March 2022.
-
Accuracy of hybrid functionals with non-self-consistent Kohn-Sham orbitals for predicting the properties of semiconductors
Authors:
Jonathan M. Skelton,
David S. D. Gunn,
Sebastian Metz,
Stephen C. Parker
Abstract:
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…
▽ More
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.
△ Less
Submitted 1 May, 2020;
originally announced May 2020.
-
Lattice dynamics of $\textit{Pnma}$ Sn(S$_{1-x}$Se$_{x}$) solid solutions: energetics, phonon spectra and thermal transport
Authors:
Jonathan M. Skelton
Abstract:
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…
▽ More
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.
△ Less
Submitted 30 January, 2020;
originally announced January 2020.
-
Theoretical study of phase transitions in Sb2S3, Bi2S3 and Sb2Se3 under compression
Authors:
E. Lora da Silva,
J. M. Skelton,
P. Rodríguez-Hernández,
A. Muñoz,
D. Martínez-García,
F. J. Manjón
Abstract:
We report a theoretical study of Sb2S3, Sb2Se3 and Bi2S3 sesquichalcogenides at hydrostatic pressures up to 60 GPa. We explore the possibility that the R-3m, C2/m, C2/c and Im-3m phases observed in sesquichalcogenides with heavier cations, viz. Bi2Se3, Bi2Te3 and Sb2Te3, could also be formed in the former compounds as suggested by recent experiments. Our calculations show that the C2/m and C2/c ph…
▽ More
We report a theoretical study of Sb2S3, Sb2Se3 and Bi2S3 sesquichalcogenides at hydrostatic pressures up to 60 GPa. We explore the possibility that the R-3m, C2/m, C2/c and Im-3m phases observed in sesquichalcogenides with heavier cations, viz. Bi2Se3, Bi2Te3 and Sb2Te3, could also be formed in the former compounds as suggested by recent experiments. Our calculations show that the C2/m and C2/c phases are energetically unstable in any of the three compounds and over the entire range of pressures examined. In contrast, the disordered bcc-like Im-3m phase is energetically stable at high pressures, but only the Sb2Se3 phase does not show dynamical instabilities below 60 GPa. Our calculations show that the Pnma structure is the most energetically stable Sb2S3 and Bi2S3 phase at ambient pressure whereas, surprisingly, the R-3m phase of Sb2Se3 is predicted to have the lowest enthalpy energy at 0 GPa, in contradiction to experimental evidence. Our calculations show that both Pnma and R-3m phases are dynamical and mechanically stable at room pressure. These results suggest that the formation of the R-3m phase in Sb2Se3 could be feasible at close to ambient conditions. To aid its identification, we provide a theoretical crystal structure (lattice and atomic parameters) and complete infrared and Raman spectra for this phase.
△ Less
Submitted 16 November, 2019; v1 submitted 19 April, 2019;
originally announced April 2019.
-
Electronic and Phonon Instabilities in Bilayer Graphene under Applied External Bias
Authors:
E. Lora da Silva,
M. C. Santos,
J. M. Skelton,
Tao Yang,
T. Santos,
S. C. Parker,
A. Walsh
Abstract:
We have performed electronic-structure and lattice-dynamics calculations on the AB and AA structures of bilayer graphene. We study the effect of external electric fields and compare results obtained with different levels of theory to existing theoretical and experimental results. Application of an external field to the AB bilayer alters the electronic spectrum, with the bands changing under bias f…
▽ More
We have performed electronic-structure and lattice-dynamics calculations on the AB and AA structures of bilayer graphene. We study the effect of external electric fields and compare results obtained with different levels of theory to existing theoretical and experimental results. Application of an external field to the AB bilayer alters the electronic spectrum, with the bands changing under bias from a parabolic to a "Mexican hat" double-well structure. This results in a semi-metal-to-semiconductor phase transition, with the size of the induced electronic band-gap being tuneable through the field strength. A reduction of continuous symmetry from a hexagonal to a triangular lattice is also evidenced through in-plane electronic charge inhomogeneities between the sublattices. When spin-orbit coupling is turned on for the AB system, we find that the bulk gap decreases, gradually increasing for larger intensities of the bias. Under large bias the energy dispersion recovers the Mexican hat structure, since the energy interaction between the layers balances the coupling interaction. We find that external bias perturbs the harmonic phonon spectra and leads to anomalous behaviour of the out-of-plane flexural ZA and layer-breathing ZO' modes. For the AA system, the electronic and phonon dispersions both remain stable under bias, but the phonon spectrum exhibits zone-center imaginary modes due to layer-sliding dynamical instabilities.
△ Less
Submitted 18 September, 2018;
originally announced September 2018.
-
Acoustic Phonon Lifetimes Limit Thermal Transport in Methylammonium Lead Iodide
Authors:
Aryeh Gold-Parker,
Peter M. Gehring,
Jonathan M. Skelton,
Ian C. Smith,
Dan Parshall,
Jarvist M. Frost,
Hemamala I. Karunadasa,
Aron Walsh,
Michael F. Toney
Abstract:
Hybrid organic-inorganic perovskites (HOIPs) have become an important class of semiconductors for solar cells and other optoelectronic applications. Electron-phonon coupling plays a critical role in all optoelectronic devices, and although the lattice dynamics and phonon frequencies of HOIPs have been well studied, little attention has been given to phonon lifetimes. We report the first high-preci…
▽ More
Hybrid organic-inorganic perovskites (HOIPs) have become an important class of semiconductors for solar cells and other optoelectronic applications. Electron-phonon coupling plays a critical role in all optoelectronic devices, and although the lattice dynamics and phonon frequencies of HOIPs have been well studied, little attention has been given to phonon lifetimes. We report the first high-precision measurements of acoustic phonon lifetimes in the hybrid perovskite methylammonium lead iodide (MAPI), using inelastic neutron spectroscopy to provide high energy resolution and fully deuterated single crystals to reduce incoherent scattering from hydrogen. Our measurements reveal extremely short lifetimes on the order of picoseconds, corresponding to nanometer mean free paths and demonstrating that acoustic phonons are unable to dissipate heat efficiently. Lattice-dynamics calculations using ab-initio third-order perturbation theory indicate that the short lifetimes stem from strong three-phonon interactions and a high density of low-energy optical phonon modes related to the degrees of freedom of the organic cation. Such short lifetimes have significant implications for electron-phonon coupling in MAPI and other HOIPs, with direct impacts on optoelectronic devices both in the cooling of hot carriers and in the transport and recombination of band edge carriers. These findings illustrate a fundamental difference between HOIPs and conventional photovoltaic semiconductors and demonstrate the importance of understanding lattice dynamics in the effort to develop metal halide perovskite optoelectronic devices.
△ Less
Submitted 17 July, 2018;
originally announced July 2018.
-
Spontaneous Octahedral Tilting in the Cubic Inorganic Caesium Halide Perovskites CsSnX$_3$ and CsPbX$_3$ (X = F, Cl, Br, I)
Authors:
Ruo Xi Yang,
Jonathan M. Skelton,
E. Lora da Silva,
Jarvist M. Frost,
Aron Walsh
Abstract:
The local crystal structures of many perovskite-structured materials deviate from the average space group symmetry. We demonstrate, from lattice-dynamics calculations based on quantum chemical force constants, that all the caesium-lead and caesium-tin halide perovskites exhibit vibrational instabilities associated with octahedral titling in their high-temperature cubic phase. Anharmonic double-wel…
▽ More
The local crystal structures of many perovskite-structured materials deviate from the average space group symmetry. We demonstrate, from lattice-dynamics calculations based on quantum chemical force constants, that all the caesium-lead and caesium-tin halide perovskites exhibit vibrational instabilities associated with octahedral titling in their high-temperature cubic phase. Anharmonic double-well potentials are found for zone-boundary phonon modes in all compounds with barriers ranging from 108 to 512 meV. The well depth is correlated with the tolerance factor and the chemistry of the composition, but is not proportional to the imaginary harmonic phonon frequency. We provide quantitative insights into the thermodynamic driving forces and distinguish between dynamic and static disorder based on the potential-energy landscape. A positive band gap deformation (spectral blueshift) accompanies the structural distortion, with implications for understanding the performance of these materials in applications areas including solar cells and light-emitting diodes.
△ Less
Submitted 14 September, 2017; v1 submitted 1 August, 2017;
originally announced August 2017.
-
Chemical and Lattice Stability of the Tin Sulfides
Authors:
Jonathan M. Skelton,
Lee A. Burton,
Fumiyasu Oba,
Aron Walsh
Abstract:
The tin sulfides represent a materials platform for earth-abundant semiconductor technologies. We present a first-principles study of the five known and proposed phases of SnS together with SnS2 and Sn2S3. Lattice-dynamics techniques are used to evaluate the dynamical stability and temperature-dependent thermodynamic free energy, and we also consider the effect of dispersion forces on the energeti…
▽ More
The tin sulfides represent a materials platform for earth-abundant semiconductor technologies. We present a first-principles study of the five known and proposed phases of SnS together with SnS2 and Sn2S3. Lattice-dynamics techniques are used to evaluate the dynamical stability and temperature-dependent thermodynamic free energy, and we also consider the effect of dispersion forces on the energetics. The recently identified π-cubic phase of SnS is found to be metastable with respect to the well-known orthorhombic Pnma/Cmcm equilibrium. The Cmcm phase is a low-lying saddle point between Pnma minima on the potential-energy surface, and is observed as an average structure at high temperatures. Bulk rocksalt and zincblende phases are found to be dynamically unstable, and we show that whereas rocksalt SnS can potentially be stabilised under a reduction of the lattice constant, the hypothetical zincblende phase proposed in several earlier studies is extremely unlikely to form. We also investigate the stability of Sn2S3 with respect to SnS and SnS2, and find that both dispersion forces and vibrational contributions to the free energy are required to explain its experimentally-observed resistance to decomposition.
△ Less
Submitted 16 March, 2017; v1 submitted 1 March, 2017;
originally announced March 2017.
-
Phonon anharmonicity, lifetimes and thermal transport in CH$_3$NH$_3$PbI$_3$ from many-body perturbation theory
Authors:
Lucy D. Whalley,
Jonathan M. Skelton,
Jarvist M. Frost,
Aron Walsh
Abstract:
Lattice vibrations in CH$_3$NH$_3$PbI$_3$ are strongly interacting, with double well instabilities present at the Brillouin zone boundary. Analysis within a first-principles lattice dynamics framework reveals anharmonic potentials with short phonon quasi-particle lifetimes and mean-free paths. The phonon behaviour is distinct from the inorganic semiconductors GaAs and CdTe where three-phonon inter…
▽ More
Lattice vibrations in CH$_3$NH$_3$PbI$_3$ are strongly interacting, with double well instabilities present at the Brillouin zone boundary. Analysis within a first-principles lattice dynamics framework reveals anharmonic potentials with short phonon quasi-particle lifetimes and mean-free paths. The phonon behaviour is distinct from the inorganic semiconductors GaAs and CdTe where three-phonon interaction strengths are three orders of magnitude smaller. The implications for the applications of hybrid halide perovskites arising from thermal conductivity, band-gap deformation, and charge-carrier scattering through electron-phonon coupling, are presented.
△ Less
Submitted 12 October, 2016; v1 submitted 3 September, 2016;
originally announced September 2016.
-
A general forcefield for accurate phonon properties of metal-organic frameworks
Authors:
J. K. Bristow,
J. M. Skelton,
K. L. Svane,
A. Walsh,
J. D. Gale
Abstract:
We report the development of a forcefield capable of reproducing accurate lattice dynamics of metal-organic frameworks. Phonon spectra, thermodynamic and mechanical properties, such as free energies, heat capacities and bulk moduli, are calculated using the quasi-harmonic approximation to account for anharmonic behaviour due to thermal expansion. Comparison to density functional theory calculation…
▽ More
We report the development of a forcefield capable of reproducing accurate lattice dynamics of metal-organic frameworks. Phonon spectra, thermodynamic and mechanical properties, such as free energies, heat capacities and bulk moduli, are calculated using the quasi-harmonic approximation to account for anharmonic behaviour due to thermal expansion. Comparison to density functional theory calculations of properties such as Gr$\mathrm{\ddot{u}}$neisen parameters, bulk moduli and thermal expansion supports the accuracy of the derived forcefield model. Material properties are also reported in a full analysis of the lattice dynamics of an initial subset of structures including: MOF-5, IRMOF-10, UiO-66, UiO-67, NOTT-300, MIL-125, MOF-74 and MOF-650.
△ Less
Submitted 25 July, 2016;
originally announced July 2016.
-
Direct Observation of Dynamic Symmetry Breaking above Room Temperature in Methylammonium Lead Iodide Perovskite
Authors:
Alexander N. Beecher,
Octavi E. Semonin,
Jonathan M. Skelton,
Jarvist M. Frost,
Maxwell W. Terban,
Haowei Zhai,
Ahmet Alatas,
Jonathan S. Owen,
Aron Walsh,
Simon J. L. Billinge
Abstract:
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 ligh…
▽ More
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.
△ Less
Submitted 5 October, 2016; v1 submitted 29 June, 2016;
originally announced June 2016.
-
Suppression of lattice thermal conductivity by mass-conserving cation mutation in multi-component semiconductors
Authors:
Taizo Shibuya,
Jonathan M. Skelton,
Adam J. Jackson,
Kenju Yasuoka,
Atsushi Togo,
Isao Tanaka,
Aron Walsh
Abstract:
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 th…
▽ More
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.
△ Less
Submitted 15 June, 2016;
originally announced June 2016.
-
Free energy of ligand removal in the metal-organic framework UiO-66
Authors:
Jessica K. Bristow,
Katrine L. Svane,
Davide Tiana,
Jonathan M. Skelton,
Julian D. Gale,
Aron Walsh
Abstract:
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 re…
▽ More
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.
△ Less
Submitted 19 February, 2016;
originally announced February 2016.
-
Anharmonicity in the high-temperature Cmcm phase of SnSe: soft modes and three-phonon interactions
Authors:
Jonathan M. Skelton,
Lee A. Burton,
Stephen C. Parker,
Aron Walsh,
Chang-Eun Kim,
Aloysius Soon,
John Buckeridge,
Alexey A. Sokol,
C. Richard A. Catlow,
Atsushi Togo,
Isao Tanaka
Abstract:
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-ph…
▽ More
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.
△ Less
Submitted 30 July, 2016; v1 submitted 11 February, 2016;
originally announced February 2016.
-
Solid-state chemistry of glassy antimony oxides
Authors:
Chang-Eun Kim,
Jonathan M. Skelton,
Aron Walsh,
Aloysius Soon
Abstract:
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…
▽ More
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.
△ Less
Submitted 28 September, 2015;
originally announced September 2015.
-
Crystal structure optimisation using an auxiliary equation of state
Authors:
Adam J. Jackson,
Jonathan M. Skelton,
Christopher H. Hendon,
Keith T. Butler,
Aron Walsh
Abstract:
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…
▽ More
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.
△ Less
Submitted 12 October, 2015; v1 submitted 31 July, 2015;
originally announced July 2015.
-
Lattice dynamics and vibrational spectra of the orthorhombic, tetragonal and cubic phases of methylammonium lead iodide
Authors:
Federico Brivio,
Jarvist M. Frost,
Jonathan M. Skelton,
Adam J. Jackson,
Oliver J. Weber,
Mark T. Weller,
Alejandro R. Goni,
Aurelien M. A. Leguy,
Piers R. F. Barnes,
Aron Walsh
Abstract:
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 follo…
▽ More
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.
△ Less
Submitted 23 September, 2015; v1 submitted 28 April, 2015;
originally announced April 2015.
-
Thermal physics of the lead chalcogenides PbS, PbSe, and PbTe from first principles
Authors:
Jonathan Michael Skelton,
Stephen C. Parker,
Atsushi Togo,
Isao Tanaka,
Aron Walsh
Abstract:
The lead chalcogenides represent an important family of functional materials, in particular due to the benchmark high-temperature thermoelectric performance of PbTe. A number of recent investigations, experimental and theoretical, have aimed to gather insight into their unique lattice dynamics and electronic structure. However, the majority of first-principles modelling has been performed at fixed…
▽ More
The lead chalcogenides represent an important family of functional materials, in particular due to the benchmark high-temperature thermoelectric performance of PbTe. A number of recent investigations, experimental and theoretical, have aimed to gather insight into their unique lattice dynamics and electronic structure. However, the majority of first-principles modelling has been performed at fixed temperatures, and there has been no comprehensive and systematic computational study of the effect of temperature on the material properties. We report a comparative lattice-dynamics study of the temperature dependence of the properties of PbS, PbSe and PbTe, focussing particularly on those relevant to thermoelectric performance, viz. phonon frequencies, lattice thermal conductivity, and electronic band structure. Calculations are performed within the quasi-harmonic approximation, with the inclusion of phonon-phonon interactions from many-body perturbation theory, which are used to compute phonon lifetimes and predict the lattice thermal conductivity. The results are critically compared against experimental data and other calculations, and add new insight to on-going research on the PbX compounds in relation to the off-centring of Pb at high temperatures, which is shown to be related to phonon softening. The agreement with experiment suggests that this method could serve as a straightforward, powerful and generally-applicable means of investigating the temperature dependence of material properties from first principles.
△ Less
Submitted 24 May, 2014;
originally announced May 2014.