大家好，我是俭朴的烤鸡老师。
本文介绍如何使用lammps计算径向分布函数g®。
径向分布函数（Radial distribution function）是指给定某个粒子的坐标，其他粒子在空间的分布几率。
在lammps中对应的命令为：compute rdf。
语法格式为：
compute ID group-ID rdf Nbin itype1 jtype1 itype2 jtype2 ... keyword/value ...
其中，Nbin为分片数量，一般选择数值为100-500之间。
itype1表示中心原子，jtype1为分布原子。
itype1和jtype1应该成对设置，表示计算第i种原子周围出现第j种原子的概率。
如果不设置itype1、jtype1表示全部原子周围出现其它原子的概率。
计算结果可配合fix ave/time命令输出。
例如：
compute myRDF all rdf 50fix 1 all ave/time 100 1 100 c_myRDF[*] file tmp.rdf mode vector
以c_myRDF[1]为横坐标，c_myRDF[2]为纵坐标，即可绘制g®曲线。
下面给出一个具体案例：
unitsrealatom_styleatomic# 读取data文件read_dataKr.data# 力场设置pair_style lj/cut 10.0pair_coeff 1 1 0.34319013128 3.591neighbor 1.0 binneigh_modify every 1 check yes#输出设置thermo 1000thermo_style custom step temp pe ke etotaldump 2 all custom 5000 md.lammpstrj id type x y z#计算rdftimestep 1.0compute myRDF all rdf 1000 1 1fix 1 all nvefix 2 all ave/time 1 100000 100000 c_myPDF[*] file rdf.dat mode vectorrun 100000unfix 2unfix 1
更多lammps教程案例请关注微信公号：
lammps加油站