计算化学公社

标题: 求助：径向分布函数问题 [打印本页]

作者
Author: fdy哈哈哈 时间: 2022-5-27 17:41
标题: 求助：径向分布函数问题
本帖最后由 fdy哈哈哈于 2022-5-27 18:05 编辑

请教各位老师，首先在听取了sob老师的建议增加采样帧数后，做出的曲线明显要平滑很多了。其次，对于我计算出来的结果，我有地方不明白，想向各位老师请教。图1是我算出来的水体系的O-H径向分布函数，其中大约2、3.5埃的峰对应形成的氢键。随着r的增大，图1的g（r）趋向于1，但是图2的取趋向于0，这是为什么呢？
补充说明：
图1：SPC水模型晶胞参数24.4埃 cubic rdf截断15埃
图2：将水分子放入MOF单晶胞内进行的MD模拟。MOF晶胞参数16.991埃 cubic，手动放进6个水分子进行结构弛豫，然后进行动力学计算，系综NPT 温控NHL 压控Berendse 步长1fs，时长500ps 步数500000 每200步输出一次结果 rdf截断15埃

作者
Author: sobereva 时间: 2022-5-28 01:48
搞清楚rdf的定义、计算公式，自然就知道怎么回事
图2的情况都不是像一般水溶液体系那样水作为溶剂那样占体系大部分区域，本来就没理由收敛到1左右

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