Temperature-dependent phonon spectrum of transition metal dichalcogenides calculated from the spectral energy density: Lattice thermal conductivity as an application

Arash Mobaraki, Cem Sevik, Haluk Yapicioglu, Deniz Çakır, and Oğuz Gülseren
Phys. Rev. B 100, 035402 – Published 2 July 2019
PDFHTMLExport Citation

Abstract

Predicting the mechanical and thermal properties of quasi-two-dimensional (2D) transition metal dichalcogenides (TMDs) is an essential task necessary for their implementation in device applications. Although rigorous density-functional-theory–based calculations are able to predict mechanical and electronic properties, mostly they are limited to zero temperature. Classical molecular dynamics facilitates the investigation of temperature-dependent properties, but its performance highly depends on the potential used for defining interactions between the atoms. In this study, we calculated temperature-dependent phonon properties of single-layer TMDs, namely, MoS2, MoSe2, WS2, and WSe2, by utilizing Stillinger-Weber–type potentials with optimized sets of parameters with respect to the first-principles results. The phonon lifetimes and contribution of each phonon mode in thermal conductivities in these monolayer crystals are systematically investigated by means of the spectral-energy-density method based on molecular dynamics simulations. The obtained results from this approach are in good agreement with previously available results from the Green-Kubo method. Moreover, detailed analysis of lattice thermal conductivity, including temperature-dependent mode decomposition through the entire Brillouin zone, shed more light on the thermal properties of these 2D crystals. The LA and TA acoustic branches contribute most to the lattice thermal conductivity, while ZA mode contribution is less because of the quadratic dispersion around the Brillouin zone center, particularly in MoSe2 due to the phonon anharmonicity, evident from the redshift, especially in optical modes, by increasing temperature. For all the considered 2D crystals, the phonon lifetime values are compelled by transition metal atoms, whereas the group velocity spectrum is dictated by chalcogen atoms. Overall, the lattice thermal conductivity is linearly proportional with inverse temperature.

  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Received 16 December 2018

DOI:https://doi.org/10.1103/PhysRevB.100.035402

©2019 American Physical Society

Physics Subject Headings (PhySH)

Condensed Matter, Materials & Applied Physics

Authors & Affiliations

Arash Mobaraki1, Cem Sevik2, Haluk Yapicioglu3, Deniz Çakır4, and Oğuz Gülseren1,*

  • 1Department of Physics, Bilkent University, Ankara 06800, Turkey
  • 2Department of Mechanical Engineering, Faculty of Engineering, Eskisehir Technical University, Eskisehir, TR 26555, Turkey
  • 3Department of Industrial Engineering, Faculty of Engineering, Eskisehir Technical University, Eskisehir, TR 26555, Turkey
  • 4Department of Physics and Astrophysics, University of North Dakota, Grand Forks, North Dakota 58202, USA

  • *gulseren@fen.bilkent.edu.tr

Article Text (Subscription Required)

Click to Expand

Supplemental Material (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand
Issue

Vol. 100, Iss. 3 — 15 July 2019

Reuse & Permissions
Access Options
CHORUS

Article Available via CHORUS

Download Accepted Manuscript
Author publication services for translation and copyediting assistance advertisement

Authorization Required


×
×

Images

×

Sign up to receive regular email alerts from Physical Review B

Log In

Cancel
×

Search


Article Lookup

Paste a citation or DOI

Enter a citation
×