【LAMMPS翻译系列】pair_write命令
pair_write命令将原子对间所定义的势函数,以距离作为自变量,将对应的能量和受力写入到文件中。
使用语法
pair_write itype jtype N style inner outer file keyword Qi Qj
- itype,jtype = 2 atom types
- N = # of values
- style = r or rsq or bitmap
- inner,outer = inner and outer cutoff (distance units)
- file = name of file to write values to
- keyword = section name in file for this set of tabulated values
- Qi,Qj = 2 atom charges (charge units) (optional)
使用举例
pair_write 1 3 500 r 1.0 10.0 table.txt LJ pair_write 1 1 1000 rsq 2.0 8.0 table.txt Yukawa_1_1 -0.5 0.5
使用介绍
该命令将原子对间所定义的势函数,以距离作为自变量,将对应的能量和受力写入到文件中。如果你想绘制势函数曲线或调试势函数,这个命令会比较有用。在写文件的时候,如果文件已经存在了,那么数据会添加到文件的末尾,这样就可以将多个能量-受力数据块写入到同一文件中。
在命令中,itype和jtype指定要计算能量和力的原子类型,inner和outer为要计算的距离区间,计算中所用到的其他参数由相应的pair_coeff命令指定。如果这里设置的style是r,那么会输出N个在r上均匀分布的能量-受力;如果style是rsq,那么会输出N个距离在r^2上均匀分布的能量-受力。
举个例子,如果N=7,style=r,inner=1.0,outer=4.0,那么会计算的值有:r = 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0。
如果style是bitmap[位图],那么会有2^N个数据按着特定的格式写入到文件中,其顺序与pair_coeff命令读入pair_style table所指定文件的顺序是一致的。为了能够保证位图数据表的精度,需要设置:N>=12,inner小于最近的两个原子间的距离,outer小于等于势函数的截断。
如果要计算势函数的
是有电荷的,那么可以为这对原子指定电荷量。如果没有指定,默认会使用Qi=Qj=1.0。
所写出文件的格式是与pair_style命令中table类型所需要的输入文件是一样的,其实keyword指定的是数据段的名称。写入到文件中的每一行都是由序号(1-N)、原子距离(距离单位)、能量(能量单位)、受力(力单位)构成的。
使用限制
所有原子对间的力场参数和其他类型的相互作用都必须在使用该命令之前已经被设置好了。
受限于程序中原子对间力的计算方式,inner必须设置为大于0的数,即便在r=0所对应的相互作用是有限值。
对于EAM势函数而言,该命令只能输出其对势项,而不能输出嵌入项。