lammps教程：非晶sio2建模方法

大家好，我是小马老师。
本文介绍非晶sio2的建模方法。
sio2的建模方法有很多，可以下载cif文件使用atomsk转为data文件，也可以直接在ms中导入sio2结构，这些方法建的sio2都是标准的sio2晶体。
如果要建立非晶状态的sio2，就是把排列规则的sio2变成无序排列的sio2。
非晶建模方法有很多，在lammps中，可以通过高温高压的方法打乱sio2的排列实现非晶建模。
代码中的温度和步数可根据需要自己调整。
下面给出in文件案例代码。

units           metal
atom_style      atomic
boundary        p p p
timestep        0.001
neighbor        2 bin 
neigh_modify    every 1 delay 0 check yes
#读取data文件
read_data       sio2.data
#势函数设置
pair_style      tersoff
pair_coeff      * * SiO.tersoff O Si
#温度初始化
velocity all create 7000 98989
#高温弛豫
dump 1 all atom 1000 dump.xyz
fix 1 all npt temp 7000 7000 0.1 iso 100 100 1
run 50000
unfix 1
#降温
fix 1 all npt temp 7000 300 0.1 iso 1 1 1
run 50000
unfix 1
#300K弛豫
fix 1 all nvt temp 300 300 0.1 
run 50000
#保存非晶文件
write_data sio2_a.data

结果：

本公众号号正式推出一对一咨询辅导服务，根据课题方向不同详细讲解对应的in文件编写方法、模拟关键技术、数据后处理，经一对一咨询辅导后能够独立编写出in文件。

具体内容可查看公众号一对一培训说明：lammps一对一培训

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

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