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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 计算能带时的三个问题，求大神指点
111/2返回列表12下一页
查看: 2704  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

MilchR

捐助贵宾 (小有名气)

[求助] 计算能带时的三个问题，求大神指点

1.计算能带时，步骤relax-scf-nscf-band。其中，在scf和nscf计算的&system中要不要加入nband项？
2.在scf和nscf中，K_POINTS可不可以用automatic？还是必须用tpiba？这个“tpiba”是什么意思？
3. 我在木虫的一个算能带的贴子里看到&system中有一项 exxdiv_treatment = 'gygi-baldereschi' ，请问这是什么项？
希望能解释的明白一些～
回复此楼

» 猜你喜欢
1楼 2015-12-14 22:25:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)
1. 算能带只需要scf -> bands两步，再用bands.x处理下数据；
2. tpiba是说明下面的K点的坐标是笛卡尔坐标，一般用不到；
3. 这个参数是杂化泛函用的。
赞一下 回复此楼

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

MilchR
2楼2015-12-15 11:34:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

MilchR

捐助贵宾 (小有名气)
送红花一朵
引用回帖:
2楼: Originally posted by liqizuiyang at 2015-12-15 11:34:39
1. 算能带只需要scf -> bands两步，再用bands.x处理下数据；
2. tpiba是说明下面的K点的坐标是笛卡尔坐标，一般用不到；
3. 这个参数是杂化泛函用的。

我看了QE里对于tpiba的描述，确实是笛卡尔坐标。不过我有疑问:对于不同的格子，tpiba下的k点取法是不一样的吧?那么应该怎么取？
我进行nscf时使用了automatic的k，但是前期不能收敛，换用tpiba很快就算完了。不知道这是为何？您能解答吗？

发自小木虫Android客户端
赞一下 回复此楼
3楼2015-12-15 12:26:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)
引用回帖:
3楼: Originally posted by MilchR at 2015-12-15 12:26:11
我看了QE里对于tpiba的描述，确实是笛卡尔坐标。不过我有疑问:对于不同的格子，tpiba下的k点取法是不一样的吧?那么应该怎么取？
我进行nscf时使用了automatic的k，但是前期不能收敛，换用tpiba很快就算完了。不知道 ...

用tpiba要换算坐标，很麻烦。

nscf最好用automatic，如果不能收敛，建议加展宽。用tpiba除非能确定自己能确保高对称点坐标和对应的权重无误。
赞一下 回复此楼

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

MilchR
4楼2015-12-15 14:29:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

MilchR

捐助贵宾 (小有名气)
送红花一朵
引用回帖:
4楼: Originally posted by liqizuiyang at 2015-12-15 14:29:51
用tpiba要换算坐标，很麻烦。

nscf最好用automatic，如果不能收敛，建议加展宽。用tpiba除非能确定自己能确保高对称点坐标和对应的权重无误。...

我使用了展宽0.02，band structure calculation 开始之后连续十多个不能收敛，因为是本子算的怕过热所以关掉了...

发自小木虫Android客户端
赞一下 回复此楼
5楼2015-12-15 16:15:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

MilchR

捐助贵宾 (小有名气)
引用回帖:
4楼: Originally posted by liqizuiyang at 2015-12-15 14:29:51
用tpiba要换算坐标，很麻烦。

nscf最好用automatic，如果不能收敛，建议加展宽。用tpiba除非能确定自己能确保高对称点坐标和对应的权重无误。...

换用tpiba后很快收敛了，这可以说明这套tpiba可以用吗？

发自小木虫Android客户端
赞一下 回复此楼
6楼2015-12-15 16:16:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)
引用回帖:
6楼: Originally posted by MilchR at 2015-12-15 16:16:37
换用tpiba后很快收敛了，这可以说明这套tpiba可以用吗？
...

发下K_POINTS设置。
赞一下 回复此楼
7楼2015-12-15 18:59:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

MilchR

捐助贵宾 (小有名气)
引用回帖:
7楼: Originally posted by liqizuiyang at 2015-12-15 18:59:54
发下K_POINTS设置。...

K_POINTS (tpiba)
36
0.500000000    0.500000000    0.500000000      1.000000000
0.400000000    0.400000000    0.400000000      2.000000000
0.300000000    0.300000000    0.300000000      3.000000000
0.200000000    0.200000000    0.200000000      4.000000000
0.100000000    0.100000000    0.100000000      5.000000000
0.000000000    0.000000000    0.000000000      6.000000000
0.000000000    0.000000000    0.100000000      7.000000000
0.000000000    0.000000000    0.200000000      8.000000000
0.000000000    0.000000000    0.300000000      9.000000000
0.000000000    0.000000000    0.400000000     10.000000000
0.000000000    0.000000000    0.500000000     11.000000000
0.000000000    0.000000000    0.600000000     12.000000000
0.000000000    0.000000000    0.700000000     13.000000000
0.000000000    0.000000000    0.800000000     14.000000000
0.000000000    0.000000000    0.900000000     15.000000000
0.000000000    0.000000000    1.000000000     16.000000000
0.000000000    0.100000000    1.000000000     17.000000000
0.000000000    0.200000000    1.000000000     18.000000000
0.000000000    0.300000000    1.000000000     19.000000000
0.000000000    0.400000000    1.000000000     20.000000000
0.000000000    0.500000000    1.000000000     21.000000000
0.000000000    0.600000000    1.000000000     22.000000000
0.000000000    0.700000000    1.000000000     23.000000000
0.000000000    0.800000000    1.000000000     24.000000000
0.000000000    0.900000000    1.000000000     25.000000000
0.000000000    1.000000000    1.000000000     26.000000000
0.000000000    0.900000000    0.900000000     27.000000000
0.000000000    0.800000000    0.800000000     28.000000000
0.000000000    0.700000000    0.700000000     29.000000000
0.000000000    0.600000000    0.600000000     30.000000000
0.000000000    0.500000000    0.500000000     31.000000000
0.000000000    0.400000000    0.400000000     32.000000000
0.000000000    0.300000000    0.300000000     33.000000000
0.000000000    0.200000000    0.200000000     34.000000000
0.000000000    0.100000000    0.100000000     35.000000000
0.000000000    0.000000000    0.000000000     36.000000000
赞一下 回复此楼
8楼2015-12-16 07:58:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)
引用回帖:
8楼: Originally posted by MilchR at 2015-12-16 07:58:01
K_POINTS (tpiba)
36
0.500000000    0.500000000    0.500000000      1.000000000
0.400000000    0.400000000    0.400000000      2.000000000
0.300000000    0.300000000    0.300000000      3.0000000 ...

别这样设置。

直接用automatic算一遍scf，再用crystal_b算bands就可以了。
赞一下 回复此楼
9楼2015-12-16 09:23:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

MilchR

捐助贵宾 (小有名气)
引用回帖:
9楼: Originally posted by liqizuiyang at 2015-12-16 09:23:31
别这样设置。

直接用automatic算一遍scf，再用crystal_b算bands就可以了。...

我发现一开始不能自洽是因为celldm错了...改了之后就正常了。
太抱歉了，犯这种低级错误→_→
请问怎样给金币，我把悬赏都给您吧

发自小木虫Android客户端
赞一下 回复此楼
10楼2015-12-16 10:15:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 MilchR 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 好奇怪的filecode +3 布布和一二 2026-08-08 4/200 2026-08-08 17:46 by zhanghaozhu
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:22 by oEVWOejN9taj
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:19 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[基金申请] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 16:27 by oEVWOejN9taj
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:22 by oEVWOejN9taj
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 15:59 by oEVWOejN9taj
[论文投稿] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 15:07 by oEVWOejN9taj
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 KXLV3nuBVBY7 2026-08-07 6/300 2026-08-08 14:01 by oEVWOejN9taj
[基金申请] 关于代码变化问题，想知道的进来 +11 且听虎啸 2026-08-07 14/700 2026-08-08 12:50 by 布布和一二
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 8月时间戳变的，举个手。玩一下，释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 有没有H口的？有收到消息的吗？ +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
信息提示
关闭
请填处理意见
关闭