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

S Wang, Z Fan, RS Koster, C Fang, MA van Huis, AO Yalcin, FD Tichelaar, HW Zandbergen, TJH Vlugt

Research output: Contribution to journalArticleScientificpeer-review

30 Citations (Scopus)


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.
Original languageEnglish
Pages (from-to)11050-11061
Number of pages12
JournalThe Journal of Physical Chemistry C
Issue number20
Publication statusPublished - 2014


Dive into the research topics of 'New ab initio based pair potential for accurate simulation of phase transitions in ZnO'. Together they form a unique fingerprint.

Cite this