本文介绍非晶sio2的建模方法。

sio2的建模方法有很多，可以下载cif文件使用atomsk转为data文件，也可以直接在ms中导入sio2结构，这些方法建的sio2都是标准的sio2晶体。

如果要建立非晶状态的sio2，就是把排列规则的sio2变成无序排列的sio2。

非晶建模方法有很多，在lammps中，可以通过高温高压的方法打乱sio2的排列实现非晶建模。

初始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

模拟得到的非晶结构如下图所示：