In this theoretical study, we investigate the origins of the very low thermal conductivity of tin selenide (SnSe) using ab-initio calculations. We obtained high-temperature lattice thermal conductivity values that are close to those of amorphous compounds. We also found a strong anisotropy between the three crystallographic axes: one of the in-plane directions conducts heat much more easily than the other. Our results are compatible with most of the experimental literature on SnSe, and differ markedly from the more isotropic values reported by a recent study.The ongoing quest for improved thermoelectric materials has spanned the past six decades. 1,2 Good thermoelectrics are characterized by their high figures of merit ZT = P T /κ. This requires a low thermal conductivity κ, and a high power factor P = σS 2 (where σ is the electrical conductivity, and S the thermopower) across their intended range of working temperatures. Nanotechnology has been able to greatly improve on earlier single-crystal thermoelectric materials 3 via hindering phonon flow while keeping favorable electric properties, even when starting from a poor performer such as bulk silicon. 4 Thus, it was somewhat surprising when single-crystal SnSe recently emerged as the best thermoelectric material ever measured. Experimental measurements on monocrystalline samples 5 point to a record figure of merit of ZT = 2.6 at T = 923 K. This is due mainly to its very low lattice thermal conductivity κ . Intense interest in SnSe was spurred by these results, and to date two independent sets of measurements on polycrystals have been published. 6,7 The two series are compatible with each other, yet they both display higher values of κ than those reported in Ref. 5 for the single crystal. This contrasts with the expectation that grain boundaries should decrease the thermal conductivity below that of the single crystal. Remarkably, preexisting studies of single crystals 8 reported even higher values of κ, with a directional maximum of around 1.8 W m −1 K −1 at room temperature.To address this controversy, we studied phonon transport in SnSe from first principles based on the Boltzmann transport equation (BTE) formalism. Such modeling studies offer substantially more detailed information than is readily available from the experiments, and have demonstrated both their wide range of applicability and their predictive power. 9-13 For the present work, we focused on the Pnma-symmetric phase of SnSe. This is the stable phase from room temperature up to ∼ 750 − 800 K. 5 Hence, it is the pertinent phase for most of the aforementioned measurements. We started with an atomistic description of Pnma SnSe extracted from the AFLOWLIB.org consortium repository [auid=aflow:d188eb9b8df60e58]. 14,15 We relaxed the structure without any geometrical constraint using the density functional theory (DFT) software package VASP. 16 All technical details are presented in the supplementary material. We obtained unit cell lengths of a = 11.72Å, b = 4.20Å and c = 4.55Å, which are in reasonab...