# Quantitative predictions of the thermal conductivity in transition metal dichalcogenides: The impact of point defects in MoS_{2} and WS_{2} monolayers

S. Mahendran,
J. Carrete,
A. Isacsson,
G. K. H. Madsen,
and
P. Erhart

arXiv:2310.09405

doi: 10.48550/arXiv.2310.09405

Download PDF

Transition metal dichalcogenides are investigated for various applications at the nanoscale thanks to their unique combination of properties and dimensionality. For many of the anticipated applications, heat conduction plays an important role. At the same time, these materials often contain relatively large amounts of point defects. Here, we provide a systematic analysis of the impact of intrinsic and selected extrinsic defects on the lattice thermal conductivity of MoS_{2} and WS_{2} monolayers. We combine Boltzmann transport theory and the Green’ s function-based ** T**-matrix approach for the calculation of scattering rates. The force constants for the defect configurations are obtained from density functional theory calculations via a regression approach, which allows us to sample a rather large number of defects at a moderate computational cost and to systematically enforce both the translational and rotational acoustic sum rules. The calculated lattice thermal conductivity is in quantitative agreement with experimental data for heat transport and defect concentrations for both MoS

_{2}and WS

_{2}. Crucially, this demonstrates that the strong deviation from a 1/

*T*-temperature dependence of the lattice thermal conductivity observed experimentally, can be fully explained by the presence of point defects. We furthermore predict the scattering strengths of the intrinsic defects to decrease in the sequence V

_{Mo}≈ V

_{2S}

^{=}> V

_{2S}

^{⊥}> V

_{S}> S

_{ad}in both materials, while the scattering rates for the extrinsic (adatom) defects decrease with increasing mass such that Li

_{ad}> Na

_{ad}> K

_{ad}. Compared to earlier work, we find that both intrinsic and extrinsic adatoms are relatively weak scatterers. We attribute this difference to the treatment of the translational and rotational acoustic sum rules, which if not enforced can lead to spurious contributions in the zero-frequency limit.