|
|
本帖最后由 chieko 于 2023-3-10 13:35 编辑
在做一个氧化铝固液界面的模拟，现在先只做羟基化氧化铝的模型。力场用CLAYFF，使用Sobtop生成拓扑文件，再跑em，nvt和npt平衡。最后产生相也是npt。em和nvt目前看都正常，但是npt总是报错，提示表面羟基的O-H键bonds that rotated more than 30 degrees。具体相关信息如下：
1、模型：在MS中切割氧化铝0001表面，暴露O，加H，上下加真空层，用castep优化下结构。
2、用Sobtop基于CLAYFF力场建立拓扑。原子类型、各种力场参数我都反复检查了，没有问题。
3、em：不约束hbonds，emtol = 50.0，完成后发现固体扁了约2 A；（不知道正常与否）
4、NVT：约束hbonds，步长1 fs，v-rescale热浴，经过100 ps从0 K退火到298K，保持100 ps，全程200 ps。轨迹正常，固体形态维持，羟基H规律振动。
5、NPxyLzT：xy方向berendsen控压1 bar，z方向维持盒子尺寸（这是文献中的做法，我也试了各向同性控压，结果一样），步长1 fs，约束hbonds（这里和nvt我也试了不约束，因为步长不大，但是也一样报错），ref_p=5.0。接着nvt的继续算，报错。提示表面羟基的O-H键bonds that rotated more than 30 degrees。查看轨迹，模拟几乎刚开始（0.1 ps左右），部分表面的O-H键就开始快速旋转，并且这些羟基上O原子的位置也不能很好保持。
我反复确认过力场参数，应该没有问题（我也算过不加羟基的氧化铝表面，npt没有任何问题）。nvt也没有问题，npt下各向异性控压的mdp设置我觉得也没有问题 也给盒子填满水计算，还是不行。
这么排查了一周多，始终无法解决，希望Sob老师和大家给点意见，谢谢！
所有结构、拓扑、mdp和运行产生的文件都打包在压缩包中了，terminal的报错在Typescript中。再次感谢各位！
|
|
发表于 Post on 2023-3-9 21:08:17
|
只看该作者 Only view this author
|只看大图 View big picture
|倒序浏览 Reverse view
|阅读模式 Reading model
收藏 Add to favorites
楼主 Author