24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3059  |  回复: 16
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

dreamboat

木虫 (著名写手)


[交流] 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 ]
回复此楼

» 猜你喜欢

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

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

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★
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,以大的力的收敛误差来计算分子动力学,这样可以大大提高计算效率。

但是压强或者说应力的收敛与力的收敛不是一回事。一般只要在一定精度下能量收敛了,力就可以说收敛了。但是压强则不同,压强的计算需要更多的平面波基展开,以保证其收敛性。这是有计算方法决定的。
9楼2011-04-08 19:19:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答
★ ★
sunyang1988(金币+2): 谢谢指点 2011-04-02 19:04:24
dreamboat(金币+2): 谢谢,用哪个fortran好,每次都装成功了,还是找不到 2011-04-02 19:34:20
请楼主先说明编译不成功的报错信息。
打个不恰当的比方,你能知道下一秒地球上会去世的人是会怎么去世么?
PWscf的编译比VASP简单的多,就是一个configure和make all的过程,至于中件在库和编译器上了什么差错显然依赖于每个人系统不同的软件安装和设置。就算是SISSA那边,也没法直接就跟你说一步一步怎么做啊,尤其在没有给任何信息的情况下
2楼2011-04-02 19:03:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hn19870519

木虫 (正式写手)


dreamboat(金币+1): 2011-04-07 19:55:21
楼上正解。。。pwscf编译确实比vasp容易的很多。。。
3楼2011-04-02 20:22:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988(金币+1): 谢谢交流 2011-04-02 22:55:53
dreamboat(金币+2): xiexie 2011-04-07 19:56:01
fortran的话就用gfortran就足够了啊~在终端中直接sudo apt-get install gfortran就可以了~
4楼2011-04-02 22:49:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见