计算化学公社

标题: Lammps模拟方英石降温过程中NTP系综晶格剧烈膨胀 [打印本页]

作者
Author: Epikur 时间: 2021-8-9 10:19
标题: Lammps模拟方英石降温过程中NTP系综晶格剧烈膨胀
各位老师好，我是lammps新手，想用Lammps模拟beta-critobalite在5000K降温到300K，得到无定形氧化硅结构，用来进行后续VASP的计算结构是3x3x2的超晶胞，力场是Morse，NPT系综先在5000 K模拟16 ps，NVT用16 ps降温至4000 K，最后NPT用59.2 ps降温到300 K
结果发现在npt降温过程中，晶格膨胀的非常厉害（5000 K模拟的时候从3 nm增加到了3.4 nm，降温到室温膨胀到4.x nm了，结构上SiO2“缩成了一团”，完全不对了
如果用NVT做，发现完成模拟的结构中出现孤立的氧原子，一些硅氧四面体结构被破坏掉了
还请各位老师不吝赐教，稍加指点

欢迎光临计算化学公社 (http://bbs.keinsci.com/)