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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 计算能带时的三个问题，求大神指点
111/2返回列表12下一页
查看: 2581  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 289 分105500药学专硕求调剂(找B区学校) +4 白云123456789 2026-04-13 4/200 2026-04-16 00:18 by 粉沁若尘
[考研] 291分调剂 +11 上岸小莹加油 2026-04-09 12/600 2026-04-15 21:55 by noqvsozv
[考研] 22专硕求调剂 +9 haoyun上岸 2026-04-11 11/550 2026-04-15 14:47 by elainzaizai
[考研] 211本科材料化工求调剂 +19 YHLAH 2026-04-11 23/1150 2026-04-14 22:25 by fenglj492
[考研] 272分材料子求调剂 +41 Loy0361 2026-04-10 54/2700 2026-04-14 18:00 by lhj2009
[考研] 26药学专硕105500求调剂 +6 喽哈加油 2026-04-13 6/300 2026-04-14 16:40 by zhouxiaoyu
[考研] 279学硕食品专业求调剂院校 20+6 孤独的狼爱吃羊 2026-04-12 28/1400 2026-04-14 15:44 by zs92450
[考研] 279求调剂 +12 张番茄不炒蛋 2026-04-11 12/600 2026-04-14 15:38 by zs92450
[基金申请] RY：中国产出的科学垃圾论文，绝对数量和比例都世界第一 +6 zju2000 2026-04-14 17/850 2026-04-14 14:34 by jurkat.1640
[考研] 085600材料与化工329分求调剂 +24 叶zilin 2026-04-13 25/1250 2026-04-14 09:20 by 试管破裂
[考研] 305求调剂 +8 玛卡巴卡boom 2026-04-11 8/400 2026-04-14 09:04 by pengliang8036
[考研] 085404 298分求调剂 +11 呼啦呼啦呼呼呼 2026-04-10 12/600 2026-04-14 08:38 by wfj257
[考研] 本科西工大 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-13 23:08 by pies112
[考研] 一志愿中南大学 0855 机械 286 求调剂 +11 不会吃肉 2026-04-12 11/550 2026-04-13 21:59 by bljnqdcc
[考研] 293求调剂 +16 我爱高数高数爱� 2026-04-12 18/900 2026-04-13 21:47 by 学员JpLReM
[考研] 一志愿085802 323分求调剂 +13 drizzle_9 2026-04-12 14/700 2026-04-13 10:26 by Faiz5552
[考研] 331求调剂 +5 王国帅 2026-04-11 5/250 2026-04-11 22:56 by 溪涧流水
[考研] 一志愿厦大0856，306求调剂 +15 Bblinging 2026-04-11 15/750 2026-04-11 22:53 by 314126402
[考研] 机械专硕270求调剂，接受跨专业 +12 老师看看我吧aba 2026-04-09 14/700 2026-04-11 10:21 by laoshidan
[考研] 22408 327分求调剂 +4 韵风kon 2026-04-10 4/200 2026-04-11 09:51 by 猪会飞
信息提示
关闭
请填处理意见
关闭