计算化学公社

标题: 求助,为什么Ih冰在II相的条件下,并没有产生相变 [打印本页]

作者
Author:
gkg    时间: 2022-12-30 17:32
标题: 求助,为什么Ih冰在II相的条件下,并没有产生相变
本帖最后由 gkg 于 2022-12-30 17:49 编辑

模拟方法:全原子力场,tip4p/ice水模型,Ih冰全周期,npt系综温度200k,z方向压强1bar,xy压强为0,采用deform压缩模拟,压缩速度与这篇文章相同 (, 下载次数 Times of downloads: 9) 。模型如右图 (, 下载次数 Times of downloads: 13)


从tip4p/ice模型预测的相图来看,

(, 下载次数 Times of downloads: 12)

在400mpa时,应变为II相,但模拟中达到400mpa时为什么只是晶格有了压缩,但并没有产生相变呢?

(, 下载次数 Times of downloads: 17) (, 下载次数 Times of downloads: 11)

即使在400mpa下npt平衡5ns也几乎没有变化。在什么条件下才能产生Ih到二相的相变呢?希望得到各位老师指点。





作者
Author:
sobereva    时间: 2022-12-31 00:41
从TIP4P/ice模拟的相图来看,400mpa 200K恰好在边界位置。模拟设置别取得这么临界,尝试更高的温度和更大的压强




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3