No Arabic abstract
Peierls-Boltzmann transport equation, coupled with third-order anharmonic lattice dynamics calculations, has been widely used to model lattice thermal conductivity ($kappa_{l}$) in bulk crystals. However, its application to materials with structural phase transition at relatively high temperature is fundamentally challenged by the presence of lattice instabilities (imaginary phonon modes). Additionally, its accuracy suffers from the absence of higher-than-third-order phonon scattering processes, which are important near/above the Debye temperature. In this letter, we present an effective scheme that combines temperature-induced anharmonic phonon renormalization and four-phonon scattering to resolve these two theoretical challenges. We apply this scheme to investigate the lattice dynamics and thermal transport properties of GeTe, which undergoes a second-order ferroelectric phase transition from rhombohedral $alpha$-GeTe to rocksalt $beta$-GeTe at about 700~K. Our results on the high-temperature phase $beta$-GeTe at 800~K confirm the stabilization of $beta$-GeTe by temperature effects. We find that considering only three-phonon scattering leads to significantly overestimated $kappa_{l}$ of 3.8~W/mK at 800~K, whereas including four-phonon scattering reduces $kappa_{l}$ to 1.7~W/mK, a value comparable with experiments. To explore the possibility to further suppress $kappa_{l}$, we show that alloying $beta$-GeTe with heavy cations such as Pb and Bi can effectively reduce $kappa_{l}$ to about 1.0~W/mK, whereas sample size needs to be around 10nm through nanostructuring to achieve a comparable reduction of $kappa_{l}$.
We extend recent textit{ab initio} calculations of the electronic band structure and the phonon dispersion relations of rhombohedral GeTe to calculations of the density of phonon states and the temperature dependent specific heat. The results are compared with measurements of the specific heat. It is discovered that the specific heat depends on hole concentration, not only in the very low temperature region (Sommerfeld term) but also at the maximum of $C_p/T^3$ (around 16 K). To explain this phenomenon, we have performed textit{ab initio} lattice dynamical calculations for GeTe rendered metallic through the presence of a heavy hole concentration ($p$ $sim$ 2$times$ 10$^{21}$ cm$^{-3}$). They account for the increase observed in the maximum of $C_p/T^3$.
Lattice dynamical methods used to predict phase-transformations in crystals typically evaluate the harmonic phonon spectra and therefore do not work in frequent and important situations where the crystal structure is unstable in the harmonic approximation, such as the $beta$ structure when it appears as a high-temperature phase of the shape memory alloy (SMA) NiTi. Here it is shown by self consistent {it ab initio} lattice dynamical calculations (SCAILD) that the critical temperature for the pre-martensitic $R$ to $beta$ phase-transformation in NiTi can be effectively calculated with good accuracy, and that the $beta$-phase is a result primarily of the stabilizing interaction between different lattice vibrations.
Ferroelectric domain walls are boundaries between regions with different polarization orientations in a ferroelectric material. Using first principles calculations, we characterize all different types of domain walls forming on ($11bar{1}$), ($111$) and ($1bar{1}0$) crystallographic planes in thermoelectric GeTe. We find large structural distortions in the vicinity of most of these domain walls, which are driven by polarization variations. We show that such strong strain-order parameter coupling will considerably reduce the lattice thermal conductivity of GeTe samples containing domain walls with respect to single crystal. Our results thus suggest that domain engineering is a promising path for enhancing the thermoelectric figure of merit of GeTe.
An accurate and easily extendable method to deal with lattice dynamics of solids is offered. It is based on first-principles molecular dynamics simulations and provides a consistent way to extract the best possible harmonic - or higher order - potential energy surface at finite temperatures. It is designed to work even for strongly anharmonic systems where the traditional quasiharmonic approximation fails. The accuracy and convergence of the method are controlled in a straightforward way. Excellent agreement of the calculated phonon dispersion relations at finite temperature with experimental results for bcc Li and bcc Zr is demonstrated.
Lattice vibrations of point defects are essential for understanding non-radiative electron and hole capture in semiconductors as they govern properties including persistent photoconductivity and Shockley-Read-Hall recombination rate. Although the harmonic approximation is sufficient to describe a defect with small lattice relaxation, for cases of large lattice relaxation it is likely to break down. We describe a first-principles procedure to account for anharmonic carrier capture and apply it to the important case of the textit{DX} center in GaAs. This is a system where the harmonic approximation grossly fails. Our treatment of the anharmonic Morse-like potentials accurately describes the observed electron capture barrier, predicting the absence of quantum tunnelling at low temperature, and a high hole capture rate that is independent of temperature. The model also explains the origin of the composition-invariant electron emission barrier. These results highlight an important shortcoming of the standard approach for describing point defect ionization that is accompanied by large lattice relaxation, where charge transfer occurs far from the equilibrium configuration.