Predicting defect stability and annealing kinetics in two-dimensional PtSe2using steepest entropy ascent quantum thermodynamics

Aimen Younis; Fazel Baniasadi; Michael R von Spakovsky; William T Reynolds Jr

doi:10.1088/1361-648X/aca3f1

Predicting defect stability and annealing kinetics in two-dimensional PtSe₂using steepest entropy ascent quantum thermodynamics

J Phys Condens Matter. 2022 Dec 15;35(7). doi: 10.1088/1361-648X/aca3f1.

Authors

Aimen Younis¹, Fazel Baniasadi², Michael R von Spakovsky³, William T Reynolds Jr²

Affiliations

¹ Mechanical Engineering Department, Virginia Tech, Blacksburg, VA 24061, United States of America.
² Materials Science Engineering Department, Virginia Tech, Blacksburg, VA 24061, United States of America.
³ Center for Energy Systems Research, Mechanical Engineering Department, Virginia Tech, Blacksburg, VA 24061, United States of America.

PMID: 36395516
DOI: 10.1088/1361-648X/aca3f1

Abstract

The steepest-entropy-ascent quantum thermodynamic (SEAQT) framework was used to calculate the stability of a collection of point defects in 2D PtSe₂and predict the kinetics with which defects rearrange during thermal annealing. The framework provides a non-equilibrium, ensemble-based framework with a self-consistent link between mechanics (both quantum and classical) and thermodynamics. It employs an equation of motion derived from the principle of steepest entropy ascent (maximum entropy production) to predict the time evolution of a set of occupation probabilities that define the states of a system undergoing a non-equilibrium process. The system is described by a degenerate energy landscape of eigenvalues, and the entropy is found from the occupation probabilities and the eigenlevel degeneracies. Scanning tunneling microscopy was used to identify the structure and distribution of point defects observed experimentally in a 2D PtSe₂film. A catalog of observed defects includes six unique point defects (vacancies and anti-site defects on Pt and Se sublattices) and twenty combinations of multiple point defects in close proximity. The defect energies were estimated with density functional theory, while the degeneracies, or density of states, for the 2D film with all possible combinations or arrangements of cataloged defects was constructed using a non-Markovian Monte-Carlo approach (i.e. the Replica-Exchange-Wang-Landau algorithm (Vogelet al2013Phys. Rev. Lett.110210603)) with a q-state Potts model. The energy landscape and associated degeneracies were determined for a 2D PtSe₂film two molecules thick and30×30unit cells in area (total of 5400 atoms). The SEAQT equation of motion was applied to the energy landscape to determine how an arbitrary density and arrangement of the six defect types evolve during annealing. Two annealing processes were modeled: heating from 77 K (-196 ∘C) to 523 K (250 ^∘C) and isothermal annealing at 523 K. The SEAQT framework predicted defect configurations, which were consistent with experimental STM images.

Keywords: 2D materials; PtSe2; annealing kinetics; defect configuration; defect density; energy landscape; steepest entropy ascent.