lammps不规则界面建模技巧

大家好，我是小马老师。

本文介绍使用Atomsk软件建立不规则界面的方法。

在一些lammps模拟中会涉及到不规则界面对模拟结果的影响，如凹凸界面对摩擦力的影响、凹凸界面对纳米流体的影响等等。

Atomsk可根据设定条件重新调整原子的位置，从而生成不同的结构。

下面以Al板为例介绍不规则界面的建模代码。

（1）生成2层原子的Al板

atomsk --create fcc 4.046 Al -duplicate 20 20 1 Al_supercell.cfg

结果如图所示：

（2）新建disp.txt文件，并输入以下内容：

# Atom displacements
displacement function
uz = 10*cos(4*pi*x/Hx)

第一句是注释，运行时不起作用，仅说明文件功能。

第二句说明对原子施加displacement（移位）操作，原子位移受函数function控制

第三句设定原子z坐标调整为“10cos(4pi*x/Hx)”

（3）生成不规则界面

atomsk Al_supercell.cfg -properties disp.txt final.lmp

运行结果如下：

可以看出，原子z坐标与x坐标有关，截面为一条余弦线，公式中的10和4分别控制幅度和周期。

disp.txt更改为：

# Atom displacements
displacement function
uz = 2*cos(4*pi*x/Hx) + 2*cos(4*pi*y/Hy)

运行结果如下图所示：

更换不同的公式，可得到不同的界面结构，大家有兴趣的话可以尝试一下。

公众号案例代码以及学习交流，请加QQ群：754749935

扫描关注微信公众号：lammps加油站。