Lammps模拟——铜纳米流体

1初始模型

      该in文件模拟水体系中铜纳米颗粒。水采用SPC/E，铜采用EAM势能。铜颗粒个数为25个。模型的初始构型：

2弛豫过程

       在正式模拟该纳米流体的性质之前需要对模型采用npt系统进行弛豫。   

      首先采用fix move命令直接将纳米颗粒随机移动分散至水体系中。根据molecule关键词对25个铜颗粒分组，然后以随机的速度对25个铜颗粒在z方向上进行移动，并跑12000步。然后unfix掉所有的fix move命令，继续弛豫500000步，得到弛豫后的模型，如下图。图中为方便显示铜纳米颗粒的位置水分子中氧原子和氢原子设置为了小点。在in文件末尾写出弛豫后的data文件，然后基于新的data文件继续进行纳米流体粘度和热导率等性质的模拟。

3弛豫过程in文件

     初始data文件（before_relax.data）和弛豫完成后的data文件࿰