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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » pwscf 收敛性测试问题
171/1返回列表
查看: 3770  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

» 猜你喜欢

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

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

查看全部散金贴
1楼 2011-04-02 16:49:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

hn19870519

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

souledge

专家顾问 (著名写手)
★
sunyang1988(金币+1): 谢谢交流 2011-04-02 22:55:53
dreamboat(金币+2): xiexie 2011-04-07 19:56:01
fortran的话就用gfortran就足够了啊~在终端中直接sudo apt-get install gfortran就可以了~
赞一下(1人) 回复此楼
4楼2011-04-02 22:49:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feierxide

铜虫 (小有名气)
★
dreamboat(金币+1): xiexie 2011-04-08 10:24:51
zhang668(金币+1): 多谢交流 2011-04-08 10:36:19
我是一般看总能做收敛測试,但正确的做法如何也不太清楚
赞一下(1人) 回复此楼
5楼2011-04-08 10:12:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhang668

木虫 (正式写手)
★
sunyang1988(金币+1): 谢谢指点 2011-04-08 11:53:58
dreamboat(金币+1): xiexie 我测试K点时 总能在小K点时低 然后一下子跳跃很高 然后再慢慢降低 这个是不是要取后面趋于不变的K点 2011-04-08 14:21:01
一般只用看总能收敛即可，我测试过pwscf Encut和K点，有时候会出现总能收敛不稳定情况
赞一下(1人) 回复此楼
6楼2011-04-08 10:37:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★
sunyang1988(金币+1): 谢谢指点 2011-04-08 11:54:11
dreamboat(金币+1): xiexie 2011-04-08 14:19:09
这要看你需要计算什么物理性质了。
一般来说，电子压强的收敛比能量的收敛更慢更难，所以如果你要算压强的话，必须以压强收敛为准！
赞一下(1人) 回复此楼
7楼2011-04-08 11:46:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
★
franch(金币+1): 谢谢回帖交流 2011-04-08 17:36:40
dreamboat(金币+2): xiexie 2011-04-09 00:24:52
一般精度的计算，仅仅能量收敛就足够了，考虑能量和力的关系为
F = - d(E) / d(r)
所以当力为0时，能量一定收敛，反过来就说不准了～因为不能完全连续取点，最多只能是拟合二项式求最低，理论上可能还是有偏差的……
如果进行精度要求高，或者对晶格常数、截断能比较敏感的性能进行计算，那就还是用力，甚至stress标准来判断收敛吧……

[ Last edited by souledge on 2011-4-8 at 15:31 ]
赞一下(1人) 回复此楼
8楼2011-04-08 15:30:14
已阅   回复此楼   关注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的回帖

souledge

专家顾问 (著名写手)
★
sunyang1988(金币+1): 谢谢交流 2011-04-09 12:44:51
dreamboat(金币+1): xiexie 2011-04-26 14:54:44
引用回帖:
Originally posted by goldenfisher at 2011-04-08 19:19:14:
也对也不对。
力的精度通常只是能量精度的开方，因此通常在算力的时候至少应该保持力的收敛性在10-3之内。关于这个问题，已经有了比较不错的进展，可以参考文献：
[url]http://iopscience.iop.org/0295-5075/ ...

压强……突然想起来，在PWscf的scf或者velax计算中，能加入压强么？似乎只看到添加电场的玩法……
顺便，这两个文献都是……嘿嘿～

[ Last edited by souledge on 2011-4-9 at 11:43 ]
赞一下(1人) 回复此楼
10楼2011-04-09 11:40:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feierxide

铜虫 (小有名气)
dreamboat(金币+1): xiexie 2011-04-26 14:54:53
不知道这个对楼主有帮助么 http://muchong.com/bbs/viewthread.php?tid=2512497&fpage=1
赞一下 回复此楼
11楼2011-04-12 11:00:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
dreamboat(金币+1): xiexie 2011-04-26 14:55:01
引用回帖:
Originally posted by souledge at 2011-04-09 11:40:33:
压强……突然想起来，在PWscf的scf或者velax计算中，能加入压强么？似乎只看到添加电场的玩法……
顺便，这两个文献都是……嘿嘿～

[ Last edited by souledge on 2011-4-9 at 11:43 ]

加压强那不是最平常的事情吗
赞一下 回复此楼
12楼2011-04-12 11:07:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
dreamboat(金币+1): xiexie 2011-04-26 14:55:17
引用回帖:
Originally posted by goldenfisher at 2011-04-12 11:07:38:
加压强那不是最平常的事情吗

俺愣是没找到在PW中从哪儿加……
赞一下 回复此楼
13楼2011-04-12 11:30:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
dreamboat(金币+1): xiexie 2011-04-26 14:55:28
引用回帖:
Originally posted by souledge at 2011-04-12 11:30:34:
俺愣是没找到在PW中从哪儿加……

我晕。

calculation='vc-relax'
赞一下 回复此楼
14楼2011-04-12 22:12:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
dreamboat(金币+1): xiexie 2011-04-26 14:55:43
引用回帖:
Originally posted by goldenfisher at 2011-04-12 22:12:59:
我晕。

calculation='vc-relax'

&control
calculation = 'vc-relax' ,
/
&cell
press = 0.D0 ,
/
是这个么？PW只能给静水压么……
赞一下 回复此楼
15楼2011-04-12 23:07:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★
zzy870720z(金币+1): 谢谢建议 2011-04-13 22:32:22
dreamboat(金币+1): xiexie 2011-04-26 14:55:54
引用回帖:
Originally posted by souledge at 2011-04-12 23:07:56:
&control
calculation = 'vc-relax' ,
/
&cell
press = 0.D0 ,
/
是这个么？PW只能给静水压么……

动力学压强也可以啊，NPT系综，calculation='vc-md'
不过就比较trick了
赞一下(1人) 回复此楼
16楼2011-04-13 17:08:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
Originally posted by goldenfisher at 2011-04-13 17:08:34:
动力学压强也可以啊，NPT系综，calculation='vc-md'
不过就比较trick了

那听起来，就很想知道CASTEP中加压强是不是比较trick了~哈哈~
赞一下(1人) 回复此楼
17楼2011-04-26 17:37:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dreamboat 的主题更新
171/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +55 医学老男孩 2026-08-13 123/6150 2026-08-19 23:44 by stan_et
[基金申请] 今天基金会出结果吗？20260819 +13 kkkl_v 2026-08-19 13/650 2026-08-19 22:37 by 宝贝虫子
[基金申请] 时间戳变了，能看出什么问题？ +12 基诺咪客 2026-08-17 16/800 2026-08-19 22:04 by kk 的小木虫
[基金申请] 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
[基金申请] 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
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖？ +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 重要消息，中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 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
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
关闭
请填处理意见
关闭