[Lammps] Lammps模拟方英石降温过程中NTP系综晶格剧烈膨胀

[复制链接 Copy URL]

Epikur

20 帖子	0 威望	147 eV

积分: 167

Level 3 能力者

发消息 Send PM

电梯直达 Go to

楼主

发表于 Post on 2021-8-9 10:19:08 | 只看该作者 Only view this author |只看大图 View big picture 回帖奖励

|倒序浏览 Reverse view |阅读模式 Reading model

各位老师好，我是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做，发现完成模拟的结构中出现孤立的氧原子，一些硅氧四面体结构被破坏掉了
还请各位老师不吝赐教，稍加指点

result.png (220.44 KB, 下载次数 Times of downloads: 22)