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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 计算能带时的三个问题，求大神指点
111/2返回列表12下一页
查看: 2545  |  回复: 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下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭