24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1651  |  回复: 2
【奖励】 本帖被评价2次,作者maomao7910增加金币 1.6
本帖产生 1 个 1ST强帖 ,点击这里进行查看

maomao7910

金虫 (小有名气)


[资源] 自发极化设置LCALCPOL时存在的问题

如题,用LBERRY = .TRUE. 的方法候博已经给出详细例子,
http://valenhou.blog.edu.cn/2008/133226.html ,我不再赘述。vasp5.2后加入了LCALCPOL=.TRUE.开关,可以一次性算出电子的极化贡献。不过目前这样设置后在OUTCAR中并没有出来如手册上所讲的
Ionic dipole moment: p[ion]=(     2.25510     2.25510     2.25510 ) electrons Angst

Total electronic dipole moment: p[elc]=(     0.00000     0.00000     0.00000 ) electrons Angst
两项,而只是给出
CALCP:  cpu time  334.26: real time  341.92

Total macroscopic dipole moment: p[mac]=(     0.00000     0.00000     0.19288 ) electrons Angst
这一项,是电子极化的贡献,那么离子极化的贡献如何解决?给两个办法,如上面候博讲的求原子移动和未移动时的berry phase,但这里只求一次,即 IGPAR = 1 or2 or3就可以了。另一种方法。如果你看关于berry phase 的手册 http://cms.mpi.univie.ac.at/vasp/Berry_phase/node4.html
你会发现离子极化项其实就是 (价电荷*原子位移)所以,自己手动一算就出来了。先在势函数里面找到你要移动的原子的价电子数  ZVAL*(原子移动前的坐标-原子移动后的坐标)就等于离子极化的贡献。
欢迎指正!

[ Last edited by youzhizhe on 2011-7-6 at 12:59 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

shagua8352

银虫 (小有名气)


★★★★★ 五星级,优秀推荐

高见!
2楼2012-01-07 19:20:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

李天帮

木虫 (正式写手)


★★★★★ 五星级,优秀推荐

楼主给力啊,很巧我也遇到了这个问题,受教了感谢楼主分享经验。
3楼2014-05-29 09:55:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 maomao7910 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复(可上传附件)
信息提示
请填处理意见