本文介绍ovito分析rdf的方法以及参数设置。

专栏已经介绍过多种rdf的计算方法。

如：在in文件中计算rdf：

lammps案例：液体平衡态rdf计算模拟练习

编程后处理计算rdf：

lammps后处理：python ovito编程计算不同原子对的rdf

补充计算rdf：

lammps数据后处理：模拟完成后补充计算rdf方法2

还有一种方法：在ovito中计算rdf。以水分子的模拟为例，把dump.xyz文件导入到ovito中，添加“Coordination analysis”

ovito会自动计算所有原子对之间的rdf，同时在下方给出曲线。

cutoff radius：决定了曲线图横坐标的范围

Number of histogram bins：计算rdf时分的块数为100上图中的参数和下面的lammps命令是等效的

compute gofr all rdf 100  1 1  2 2  1 2  cutoff 10.0# O-O, H-H, O-H

如果需要导出rdf的数值，可以点击“show in data inspector”，然后点击保存图标，即可把数据保存到文本文件，从而使用其他绘图软件绘制rdf曲线。

默认情况下，ovito会给出所有原子对的rdf，如果仅需要个别原子的rdf，可以先使用“select type”命令选择特定原子，然后再选择“coordination analysis”如下图仅计算了1-1的rdf

ovito计算rdf比较方便，但也有个缺点，只能输出某一帧的rdf值，无法对一定步数内rdf值进行平均（pro版本可以），如需平均值，可选择之前介绍的计算方法。

---

公众号：lammps加油站