24小时热门版块排行榜    

CyRhmU.jpeg
查看: 516  |  回复: 0

yongning88

铜虫 (小有名气)

[交流] 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
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 yongning88 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见