24小时热门版块排行榜    

查看: 3778  |  回复: 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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳又变了 +11 wuchongjun 2026-08-20 14/700 2026-08-20 16:26 by 虫友是什么虫
[基金申请] 科研孤儿太难了 +4 我4大白菜 2026-08-20 4/200 2026-08-20 15:58 by 陈晨晨陈啊
[基金申请] 今天系统多次维护,明天很可能放榜! +9 zju2000 2026-08-16 10/500 2026-08-20 15:56 by Leogzhya
[基金申请] 放榜前的不淡定 40+3 snowwithsea 2026-08-19 12/600 2026-08-20 15:56 by lfy8008
[基金申请] 言之凿凿,明天官网会通知放榜? +14 医学老男孩 2026-08-20 16/800 2026-08-20 15:47 by LNP@mRNA
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +16 Ldrop2023 2026-08-13 20/1000 2026-08-20 11:07 by xskun
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 21/1050 2026-08-20 10:37 by Vivilian
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
信息提示
请填处理意见