|
|
[求助]
LAMMPS中势函数电子密度具体是什么?
最近在研究LAMMPS中的势函数文件,有几个问题想请教下;
假设Nrho=500, drho=0.001, Nr=500, dr=0.002, cutoff=1
对于对势,我知道可以从0到1之间以0.002为间隔取500个值,计算出相应的对势,再乘以相应的r值就是势函数文件中的对势值。
对于电子密度就不是很明白了,势函数文件中的“electron density”是指公式A2的值还是公式A4的值?还有对于f(r),可以和对势一样的方法求得,但对于rho就不知道怎么求了,假设有5个原子,按照公式A2,
ρ_1=f_2 (r_12 )+f_3 (r_13 )+f_4 (r_14 )+f_5 (r_15 )
ρ_2=f_1 (r_21 )+f_3 (r_23 )+f_4 (r_24 )+f_5 (r_25 )
.......
这里面r_12, r_13......是怎么取值的?
![LAMMPS中势函数电子密度具体是什么?]()
1.PNG
![LAMMPS中势函数电子密度具体是什么?-1]()
2.PNG
![LAMMPS中势函数电子密度具体是什么?-2]()
3.PNG |
» 猜你喜欢
关于filecode,很负责任的告诉大家
已经有5人回复
关于豆爷回答的JTJC与%2F数量
已经有7人回复
2026国自然放榜时间
已经有3人回复
国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率
已经有7人回复
好奇怪的filecode
已经有4人回复
售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
售SCI文章,我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
|