24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

dreamboat

木虫 (著名写手)

应助: 8 (幼儿园)
金币: 3090.3
帖子: 2746
在线: 593.4小时
虫号: 442693

[交流] pwscf 收敛性测试问题

在测试Ecut和k点时，只看总能吗？
最后这个总的应力要不要看啊
有时候算出来很大的应力怎么解决，谢谢

total stress  (Ry/bohr**3)                (GPa)    P= 28.53
0.00193923 0.00000000 0.00000000       28.527    0.000    0.000
0.00000000 0.00193923 0.00000000       0.000 28.527    0.000
0.00000000 0.00000000 0.00193923       0.000    0.000 28.527

[ Last edited by dreamboat on 2011-4-7 at 19:54 ]

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

难道每次计算前都咬进行收敛性测试吗？已经有14人回复
PWscf计算homo lumo或者valence conduction band，问题已解决已经有14人回复
【求助成功】castep 收敛性测试已经有22人回复
【求助成功】vasp测试K点收敛性的脚本问题已经有20人回复
【求助】pwscf计算中的warning要紧吗？如何消除呢？已经有9人回复
【已解决】【求助】为何PWscf不能得到类似于CASTEP那样的收敛测试曲线呢？已经有14人回复
【求助】pwscf并行已经有7人回复
【求助】pwscf 能带计算已经有16人回复
【求助】收敛性测试的问题已经有40人回复
【求助】castep收敛性测试问题已经有20人回复
【求助】PWSCF中Ti的赝势是不是有问题呢？已经有17人回复
【求助】PWSCF计算介电极化时为何出现如下提示？已经有19人回复
【求助】请教PWSCF晶体建模和自洽计算，不收敛已经有21人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-04-02 16:49:12

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

★ ★ ★
dreamboat(金币+1): xiexie 2011-04-09 00:28:17
sunyang1988(金币+3): 谢谢交流 2011-04-09 12:44:41

引用回帖:

Originally posted by souledge at 2011-04-08 15:30:14:
一般精度的计算，仅仅能量收敛就足够了，考虑能量和力的关系为
F = - d(E) / d(r)
所以当力为0时，能量一定收敛，反过来就说不准了～因为不能完全连续取点，最多只能是拟合二项式求最低，理论上可能还是有偏差的 ...

也对也不对。
力的精度通常只是能量精度的开方，因此通常在算力的时候至少应该保持力的收敛性在10-3之内。关于这个问题，已经有了比较不错的进展，可以参考文献：
http://iopscience.iop.org/0295-5075/88/2/20001
http://prl.aps.org/abstract/PRL/v104/i24/e245001
其基本思想就是在保证能量足够收敛的情况下，比如10-4，以大的力的收敛误差来计算分子动力学，这样可以大大提高计算效率。

但是压强或者说应力的收敛与力的收敛不是一回事。一般只要在一定精度下能量收敛了，力就可以说收敛了。但是压强则不同，压强的计算需要更多的平面波基展开，以保证其收敛性。这是有计算方法决定的。

赞一下(1人)

回复此楼

9楼2011-04-08 19:19:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 17 个回答

souledge

专家顾问 (著名写手)

★ ★
sunyang1988(金币+2): 谢谢指点 2011-04-02 19:04:24
dreamboat(金币+2): 谢谢，用哪个fortran好，每次都装成功了，还是找不到 2011-04-02 19:34:20

请楼主先说明编译不成功的报错信息。
打个不恰当的比方，你能知道下一秒地球上会去世的人是会怎么去世么？
PWscf的编译比VASP简单的多，就是一个configure和make all的过程，至于中件在库和编译器上了什么差错显然依赖于每个人系统不同的软件安装和设置。就算是SISSA那边，也没法直接就跟你说一步一步怎么做啊，尤其在没有给任何信息的情况下

赞一下(1人)

回复此楼

2楼2011-04-02 19:03:26

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖