New ab initio based pair potential for accurate simulation of phase transitions in ZnO

Journal article (2014)

Authors

S. Wang External organisation
Z. Fan Engineering Thermodynamics - Mechanical, Maritime and Materials Engineering
RS Koster External organisation
T.J.H. Vlugt Engineering Thermodynamics - Mechanical, Maritime and Materials Engineering
Research Group
Engineering Thermodynamics (Mechanical, Maritime and Materials Engineering) (TU Delft)
More Info
expand_more

Published Date

2014

Language

English

Reuse Rights

Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.

Faculty

Mechanical, Maritime and Materials Engineering

Department

Process and Energy

Research Group

Engineering Thermodynamics

Abstract

A set of interatomic pair potentials is developed for ZnO based on the partially charged rigid ion model (PCRIM). The derivation of the potentials combines lattice inversion, empirical fitting, and ab initio energy surface fitting. We show that, 0 despite the low number of parameters in this model (8), a wide range of physical properties is accurately reproduced using the new potential model. The calculated lattice parameters and elastic constants of ZnO in the wurtzite (WZ) phase, as well as the lattice parameters and stabilities of ZnO in other high-pressure and metastable phases, agree well with experiments and with density functional theory (DFT) calculations. The calculated transition pressure of the wurtzite-to-rocksalt (WZ-to-RS) transition is 12.3 GPa. A wurtzite-to-honeycomb (WZ-to-HC) phase transition induced by uniaxial pressure along the c-axis is simulated by means of molecular dynamics (MD) simulations. The WZ-to-HC transition takes place at an uniaxial pressure of 8.8 GPa while the reverse transition takes place at 2.9 GPa, which is consistent with DFT calculations. Other physical properties, including phonon dispersion, phonon density of states, and melting point calculated using our ZnO potential model are in good agreement with experimental and DFT data. Limitations of the novel ZnO potential model are also discussed.