24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

cenwanglai

荣誉版主 (知名作家)

[交流] 【求助完毕】bader计算和ELF计算的流程和精度

希望大家先帮看看我的计算流程原则上是不是正确的，然后说说计算精度的问题。

我的计算流程： geo-opt---->scf---->elf and bader

scf 后, mkdir elf bader

cp scf/* elf/.

INCAR做如下修改：
ICHARG =11
#MAGMOM = * * *
LELF = T
LWAVE = F

其它不变。k-mesh不变

cp scf/* bader/.

INCAR做如下修改：
ICHARG = 1
#MAGMOM=* * *
LAECHG = T
LWAVE = F

其他不变。k-mesh不变。

以上两个性质计算中，elf是做非自洽计算的，我的理解是，elf仅仅作为初始电荷密度的数值处理，不需要做自洽计算。而bader需要原子价层电子和内层电子，一般vasp的PAW方法不提供内层电子，所以需要以上一步scf电荷密度做输入，来补全电荷密度。原则上这两个任务设置有问题没？

精度方面，elf计算为了提高精度（图的饱满平滑），需要调高NGX(Y,Z)，然后重新做静态自洽计算。这个不必提高k采样密度，保持不变就好。

bader计算说要适当调高NGX(Y,Z)，使得总电子数等于实际电子数。这个总电子数是不是为bader分析结果文件ACF.dat的最后一行给出的电荷数？

[ Last edited by cenwanglai on 2010-12-16 at 16:14 ]

回复此楼

» 收录本帖的淘帖专辑推荐

VASP	vasp计算资料	计算化学-辅助工具	第一性
第一性

» 本帖已获得的红花（最新10朵）

wujiandongde

yzhang2010

» 猜你喜欢

球磨粉体时遇到了大的问题，请指教！已经有8人回复
江汉大学解明教授课题组招博士研究生/博士后已经有3人回复

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

matlab 积分计算算不出结果，求帮助！！！已经有10人回复
分享一篇ATK计算金纳米管（分子修饰）的输运计算已经有95人回复
发一个shake91 一维土层地震反应的计算程序（伯克利出品带说明书）已经有59人回复
vasp计算DOS,以及如何判断参数的设置是合理的，达到要求的精度已经有7人回复
【求助完毕】集群上用vasp计算ELF有ELFCAR而没有DOSCAR是啥问题？已经有15人回复
【求助】请问下FORTRAN的计算精度能不能达到MATLAB的计算精度高呢？已经有3人回复
请教药典关于精确度及计算的问题请大家赐教多谢了！！已经有5人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2010-12-14 13:39:14

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

first_yg

银虫 (小有名气)

应助: 1 (幼儿园)
金币: 463.5
帖子: 132
在线: 81.8小时
虫号: 374955

★
zsl_321123(金币+1): 鼓励交流 2011-05-01 10:46:15

引用回帖:

Originally posted by cenwanglai at 2010-12-14 13:39:14:
希望大家先帮看看我的计算流程原则上是不是正确的，然后说说计算精度的问题。

我的计算流程： geo-opt---->scf---->elf and bader

scf 后, mkdir elf bader
[box=white]cp scf/* elf/.

INCAR做 ...

bader计算说要适当调高NGX(Y,Z)，使得总电子数等于实际电子数。
本文来自: 小木虫论坛 http://muchong.com/bbs/viewthread.php?tid=2699399&fpage=2

怎么能把NGX(Y,Z)设置的正好使总电子数等于实际电子数呢？