应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 对费米能级下所有态的积分parchg与CHGCAR差别很大
111/2返回列表12下一页
查看: 2173  |  回复: 10
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

kgdu

金虫 (正式写手)

[求助] 对费米能级下所有态的积分parchg与CHGCAR差别很大

对费米能级下所有态的积分parchg与CHGCAR（静态计算得到的）差别很大，计算时加了SOC的，而且parchg的电子数目比CHGCAR少很多，不知道大家有没有遇到过这样的问题！计算的体系是绝缘体。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-01-06 20:05:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

09njtpw

铁杆木虫 (职业作家)
请问楼主对费米能级下所有态的积分parchg是怎么操作的？我可以试试我的数据看看有什么问题
赞一下 回复此楼
2楼2012-01-06 22:24:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kgdu

金虫 (正式写手)
★ ★
fzx2008(金币+2): 鼓励交流 2012-01-06 23:11:51
引用回帖:
2楼: Originally posted by 09njtpw at 2012-01-06 22:24:15:
请问楼主对费米能级下所有态的积分parchg是怎么操作的？我可以试试我的数据看看有什么问题

我是计算的 Band decomposed chargedensity，见 http://hi.baidu.com/datavault/bl ... 678cc537d3ca50.html
对费米能级下所有的态积分，
LPARD=.TRUE.
NBMOD=-3
EINT=-60.0 0.0
LSEPB=.FALSE.
LSEPK=.FALSE.
计算的体系最低能级在-55.0eV（相对费米能级），这样应该就包括了所有的态。对于不加SOC的情况，这样就算的结果和CHGCAR相差很小，但是加SOC后的PARCHG，与CHGCAR相差较大。我现在就不知道这个PARCHG到底准不准了。欢迎讨论！
赞一下(1人) 回复此楼
3楼2012-01-06 22:47:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

09njtpw

铁杆木虫 (职业作家)
引用回帖:
3楼: Originally posted by kgdu at 2012-01-06 22:47:51:
我是计算的 Band decomposed chargedensity，见 http://hi.baidu.com/datavault/bl ... 678cc537d3ca50.html
对费米能级下所有的态积分，
LPARD=.TRUE.
NBMOD=-3
EINT=-60.0 0.0
L ...

计算的时候一定要保存WAVECAR文件么？我之前算的东西都没有输出这个文件，因为太大了
本来想试试，估计帮不上忙了，有其他方法么？
赞一下 回复此楼

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

kgdu
4楼2012-01-07 10:21:07
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kgdu

金虫 (正式写手)
送鲜花一朵
引用回帖:
4楼: Originally posted by 09njtpw at 2012-01-07 10:21:07:
计算的时候一定要保存WAVECAR文件么？我之前算的东西都没有输出这个文件，因为太大了
本来想试试，估计帮不上忙了，有其他方法么？

谢谢你了，我自己再看看吧！
赞一下 回复此楼
5楼2012-01-07 15:19:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kjsun

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
WDD880227(金币+2): 感谢交流 2012-01-08 12:28:32
kgdu(金币+3): ★有帮助 2012-01-08 19:14:03
两种可能：
1.-60eV以下有能级，你没有包括进去，测试的方法是算一个 -500 到 -60 区间，看看parchg有没有值。
2.就是EINT=-60.0 0.0 这里的0.0错了， 一般情况下，0.0不费米能级的位置。解决方法一.你search一下outcar里E-fermi在哪里，比如在-2.2然后用-60.0 -2.2 重新算一下。
解决方法二： EINT=-60.0 重新算，系统自动补全费米能级位置。
赞一下(1人) 回复此楼
6楼2012-01-08 10:52:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kgdu

金虫 (正式写手)
引用回帖:
6楼: Originally posted by kjsun at 2012-01-08 10:52:54:
两种可能：
1.-60eV以下有能级，你没有包括进去，测试的方法是算一个 -500 到 -60 区间，看看parchg有没有值。
2.就是EINT=-60.0 0.0 这里的0.0错了， 一般情况下，0.0不费米能级的位置。解决方法一.你search一 ...

我先看了态密度的，最低能级在-55eV，我也算过更低的能量态的parchg 没有能带在里面，全是0. 我设置的NBMOD=-3就是相对于fermi能级的，我看了输出选择的能量也是fermi能量下的所有的。我现在怀疑是vasp的bug，对于noncollinear感觉PARCHG不对！
赞一下 回复此楼
7楼2012-01-08 19:19:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kjsun

木虫 (正式写手)
★ ★ ★
franch(金币+3): 谢谢回帖交流 2012-01-09 12:16:16
引用回帖:
7楼: Originally posted by kgdu at 2012-01-08 19:19:38:
我先看了态密度的，最低能级在-55eV，我也算过更低的能量态的parchg 没有能带在里面，全是0. 我设置的NBMOD=-3就是相对于fermi能级的，我看了输出选择的能量也是fermi能量下的所有的。我现在怀疑是vasp的bug，对 ...

你计算的态密度是也是要设置范围的，你不设置系统就是默认，默认值未必会计算所有能级。
一般上，DFT的深层能级都作为原子实（也就是定值）来处理的，但是我个人觉得vasp里面肯定是包括了深层能级的贡献的，（证据是当你计算芯位移的时候，你就会发现这些深层能级是随着环境发生了变化的，说明vasp对这些深层能级进行了处理，而不是真当作不变的原子实),但是一般上，我们不会注意这些深层能级，vasp也一般放在默认范围之外。
我猜可能是原因。这些深层能级一般能量很低，大约在负几百eV的范围内，你计算更低的能量态最低是多少啊？
赞一下(1人) 回复此楼

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

kgdu
8楼2012-01-09 12:04:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kgdu

金虫 (正式写手)
送鲜花一朵
引用回帖:
楼: Originally posted by kjsun at 2012-01-09 12:04:37:
你计算的态密度是也是要设置范围的，你不设置系统就是默认，默认值未必会计算所有能级。
一般上，DFT的深层能级都作为原子实（也就是定值）来处理的，但是我个人觉得vasp里面肯定是包括了深层能级的贡献的，（证 ...

谢谢指导
回复此楼
9楼2012-01-09 20:56:07
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kgdu

金虫 (正式写手)
引用回帖:
8楼: Originally posted by kjsun at 2012-01-09 12:04:37:
你计算的态密度是也是要设置范围的，你不设置系统就是默认，默认值未必会计算所有能级。
一般上，DFT的深层能级都作为原子实（也就是定值）来处理的，但是我个人觉得vasp里面肯定是包括了深层能级的贡献的，（证 ...

谢谢你的回复。奇怪的是我对同样的体系，做collinear的计算，再计算fermi能级下的parchg，和chgcar几乎相同，但是加上SOC后，静态计算，再次计算fermi能级下的parchg，这是就和我静态计算的chgcar不一样了，差距很大。感觉vasp有问题。
下面是DOS
赞一下 回复此楼
10楼2012-01-09 21:07:07
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kgdu 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭