Molecular Dynamics Simulation on Orientation-Dependent Mechanical Behaviors of ZnO Single Crystals Under Nanoindentation

作者：Xiaolin Zhu, Jijun Li, Shiting Yang, Weiguang Zhang, Xiuxia Li, Hui Tang, Fengchao Lang, Lin Lin, Xiaohu Hou, Xueping Zhao, Jiayi Chen · 发表于：Materials · 年份：2025 · DOI：10.3390/ma18163905 · 被引用次数：46 · 研究领域：Metal and Thin Film Mechanics、ZnO doping and properties、Advanced ceramic materials synthesis

The present study aims to investigate the orientation-dependent mechanical behaviors of ZnO single crystals under nanoindentation by molecular dynamics simulation. The load-indentation depth curves, atomic displacement, shear strain and dislocations for the c-plane, m-plane and a-plane ZnO single crystals were analyzed in detail. The simulation results showed that the elastic deformation stage of the loading curves for the three oriented ZnO single crystals can be described well by the Herz elastic contact model. The Young modulus values for the c-plane, m-plane and a-plane ZnO were calculated to be 122.5 GPa, 158.3 GPa and 170.5 GPa, respectively. The onset of plastic deformation occurred first in a-plane ZnO, then in m-plane ZnO, and lastly in c-planeZnO. The atomic displacement vectors in the three oriented ZnO single crystals were in good agreement with the primary activated slip systems predicted by the maximum Schmid factor. For the c-plane ZnO, the activated pyramidal {112¯2}<112¯3> slip system led to a complex dislocation pattern surrounding the indenter. A U-shaped prismatic half-loop was formed in the [211¯0] direction, confirming the activation of the prismatic {101¯0}<112¯0> slip system. For the m-plane ZnO, the activated prismatic {101¯0}<112¯0> slip system led to the preferential nucleation of dislocations along the 11¯20 and [2¯110] directions. A prismatic loop was formed and emitted along the [2¯110] direction, governed by a confined glide on {101¯0} planes. F...