大家好，我是小马老师。
本文介绍atomsk孪晶多晶建模方法。
atomsk多晶建模的原理是先建立一个晶胞，然后编写一个polycrystal.txt文件，设定最终模型的尺寸和晶粒的个数。
atomsk按照晶粒个数随机生成对应的节点node坐标，根据Voronoi算法生成每个晶粒空间，并把最初设定的晶胞填充到晶粒空间内，从而得到多晶文件。
如果把最初的晶胞换成孪晶，最终得到的多晶就是孪晶多晶结构。
如Al的孪晶多晶建模代码为：
atomsk --create fcc 4.02 Al orient [11-2] [111] [-110] -duplicate 1 8 1 Al_cell.xsf
atomsk Al_cell.xsf -mirror 0 Y -wrap Al_mirror.xsf
atomsk --merge Y 2 Al_cell.xsf Al_mirror.xsf Al_final.cfg
atomsk --polycrystal Al_final.cfg polycrystal.txt Al_twin_polyX.lmp
对应的polycrystal.txt为：
box 400 200 100
random 30
最终的结构：
根据以上原理，可以建成任意孪晶多晶结构。