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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 对费米能级下所有态的积分parchg与CHGCAR差别很大
111/2返回列表12下一页
查看: 2410  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 有没有学校材料专业收跨调(一志愿085410) +6 momo(上岸版) 2026-04-06 9/450 2026-04-11 22:38 by wj165256
[考研] 295分求调剂 +9 ?要上岸? 2026-04-10 9/450 2026-04-11 10:48 by qingpingzhu
[考研] 268分085602化学工程调剂 +27 月照花林。 2026-04-09 27/1350 2026-04-11 10:40 by maddjdld
[考研] 生物学调剂 可调剂到生物与医药 +8 李政莹 2026-04-06 9/450 2026-04-11 10:36 by wwj2530616
[考研] 0831生医工第一轮调剂失败求助 +10 小熊睿睿_s 2026-04-11 13/650 2026-04-11 10:04 by maddjdld
[考研] 263能源动力专硕求调剂 +4 加大号饭盒袋 2026-04-10 4/200 2026-04-10 20:52 by gong120082
[考研] 22408 366分，本科211，一志愿西工大 +4 Rubt 2026-04-09 4/200 2026-04-10 19:51 by chemisry
[考研] 一志愿京区985，085401，与本科专业一致，电子信息工程， +4 阳光开朗的男孩 2026-04-10 4/200 2026-04-10 18:27 by shenrf
[考研] 08工学 309分求调剂 +6 Yin DY 2026-04-08 6/300 2026-04-10 09:18 by Delta2012
[考研] 284求调剂 +7 让我上岸吧阿西 2026-04-09 7/350 2026-04-09 18:59 by haironglove
[考研] 一志愿中国科学院上海有机所，有机化学356分找调剂 +11 Nadiums 2026-04-09 11/550 2026-04-09 18:04 by lijunpoly
[考研] 085400电子信息类（川大控制工程）求调剂可跨专业 求老师联系 +3 626776879 2026-04-08 3/150 2026-04-09 16:05 by 猪会飞
[考研] 328求调剂 +17 lftmya 2026-04-07 18/900 2026-04-09 08:05 by 5268321
[考研] 327求调剂 +10 Xxjc1107. 2026-04-06 11/550 2026-04-09 01:21 by lature00
[考研] 材料与化工专硕306分找合适调剂 +27 沧海轻舟e 2026-04-06 28/1400 2026-04-08 22:06 by wdyheheeh
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 22:00 by hemengdong
[考研] 一志愿南科大生物学297分，求调剂推荐 +8 Y-yyusx 2026-04-06 9/450 2026-04-07 19:38 by biomichael
[考研] 328求调剂 +4 ghhh88888 2026-04-06 5/250 2026-04-07 14:45 by ghhh88888
[考研] 22408 331分求调剂 +4 y__1 2026-04-06 4/200 2026-04-06 17:26 by 土木硕士招生
[考研] 377求调剂 +6 by.ovo 2026-04-05 6/300 2026-04-05 22:18 by dongzh2009
信息提示
关闭
请填处理意见
关闭