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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » pwscf 收敛性测试问题
171/1返回列表
查看: 3062  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭