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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 对费米能级下所有态的积分parchg与CHGCAR差别很大
111/2返回列表12下一页
查看: 2413  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 335求调剂 +5 想上岸呀！！ 2026-04-12 5/250 2026-04-12 07:41 by laoshidan
[考研] 305求调剂 +7 玛卡巴卡boom 2026-04-11 7/350 2026-04-12 07:35 by zhouxiaoyu
[考研] 化学工程调剂289 +44 yang婷 2026-04-07 50/2500 2026-04-12 02:36 by 秋豆菜芽
[考研] 303求调剂 +14 SereinQ 2026-04-10 15/750 2026-04-11 20:43 by 蓝云思雨
[考研] 290调剂生物0860 +17 哇哈哈，。 2026-04-11 19/950 2026-04-11 20:20 by dongdian1
[考研] 271求调剂 +21 2261744733 2026-04-11 21/1050 2026-04-11 20:19 by Angchemist
[考研] 求调剂，一志愿大连理工大学354分 +5 雨声余生 2026-04-11 6/300 2026-04-11 16:12 by 雨声余生
[考研] 296求调剂 +13 汪！？！ 2026-04-10 15/750 2026-04-11 10:31 by 逆水乘风
[考研] 材料工程085601，270求调剂 +30 @ASDF1234 2026-04-08 32/1600 2026-04-11 10:30 by Delta2012
[考研] 材料与化工调剂 +12 否极泰来2026 2026-04-10 13/650 2026-04-11 00:28 by wangjihu
[考研] 344求调剂 +7 丶风雪夜归人丶 2026-04-09 7/350 2026-04-10 12:05 by pengliang8036
[考研] 274求调剂 +5 山阿蔓 2026-04-07 5/250 2026-04-09 15:28 by 18828373951
[考研] 308求调剂 +17 墨墨漠 2026-04-06 17/850 2026-04-09 09:25 by 壹往無前
[考研] 286求调剂 +19 Faune 2026-04-08 20/1000 2026-04-09 08:36 by 哦哦123
[考研] 353求调剂 +8 晴空万里air 2026-04-07 8/400 2026-04-09 00:18 by GouQ
[考研] 307求调剂 +14 超级伊昂大王 2026-04-06 14/700 2026-04-08 07:03 by 无际的草原
[考研] 338求调剂 +8 wxygxsaaaaa 2026-04-06 8/400 2026-04-08 06:58 by 无际的草原
[考研] 277、学硕，求调剂 数一104， +11 瓶子PZ 2026-04-07 12/600 2026-04-07 23:30 by 一只好果子?
[考研] 工科 22408 267求推荐 +4 wanwan00 2026-04-05 5/250 2026-04-06 22:47 by chenzhimin
[考研] 324求调剂 +3 k可乐 2026-04-05 4/200 2026-04-06 09:54 by 蓝云思雨
信息提示
关闭
请填处理意见
关闭