如果不涉及相变,可以用晶体密度,忽略固液体积差别。 也可以先用NPT充分弛豫,平衡之后得到一个平均体积(密度),将之用于NVT模拟。 也可以使用实验的液体密度。 三种方法中以第二种较为常用,因为它提供的是与你的势函数有关的密度参数。虽然这样得到的数据通常比实验值为小(体积偏大),但它保证了计算的自洽性,并使得计算中的原子应力(或体系的应力张量)降低。不过,除非你专门研究密度或体积问题,一般不会有人认真关心模拟中的密度参数,只要它与实验值的偏离不是很离谱。 |
欢迎光临 研发埠 (http://bbs.yanfabu.com/) | Powered by Discuz! X3.2 |