24小时热门版块排行榜    

查看: 4339  |  回复: 19
【奖励】 本帖被评价10次,作者xirainbow增加金币 6.25 个
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

[资源] 【原创】ecutwfc and ecutrho

以下内容来自我的博客: blog.sina.com.cn/nkasir

     When I use PWgui to edit input file of pw.x, a warning always jumps out:"some variable are not set: List of unset variables: ecutrho".

     In my mind, in general plane wave code, only the energy cutoff for wavefunctions is adequate. Therefore, I never set ecutrho manually.

     This is the instruction of ecutrho:"

   Type:           REAL
   Default:        4 * ecutwfc
   Description:    kinetic energy cutoff (Ry) for charge density and potential
                   May be larger ( for ultrasoft PP ) or somewhat smaller
                   ( but not much smaller ) than the default value. Note that
                   if you have norm-conserving PP only, setting it to a larger
                   value than the default is a waste of time."



     But now, I know that it is only true for norm-conserving pseudopotentials, not for utral-soft ones.

     I find the following statement at: http://encina.northwestern.edu/i ... nd_PWSCF_in_general

"We are going to use ultrasoft pseudopotentials here. In the case of norm-conserving pseudopotentials from Lab 2, ecutrho (the charge density cutoff) is automatically determined by 4*ecutwfc. However, in the case of ultrasoft pseudopotentials, we need an augmented charge around the ion core, so ecutrho should be higher than 4*ecutwfc. The usual value is 25-35 Ry for ecutwfc and 200-300 Ry for ecutrho. You might want to do a convergence check. Keep in mind that the value you should look at is the energy difference or force, not the absolute value of the energy (the energy will not converge unless you use very, very high ecutwfc and ecutrho). "

[ Last edited by xirainbow on 2009-7-16 at 11:31 ]
回复此楼

» 猜你喜欢

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

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

157806400

金虫 (初入文坛)


★
cenwanglai(金币+1):谢谢参与讨论 2010-05-30 21:51:25
超软赝势必须要4倍以上吗?我发过一篇文章都用的6倍左右,我也没作测试,不会是错误的吧?不过我增大ecutrho测试收敛了。用的是6倍的比例。
18楼2010-05-30 21:24:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答

wuchenwf

荣誉版主 (职业作家)


  貌似楼主发的是自己解决的一个经验
2楼2009-07-16 22:32:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★
freshgirl(金币+1,VIP+0):谢谢参与~ 7-17 14:13
wuchenwf(金币+2,VIP+0):谢谢分享 12-3 21:57
wuchenwf ,您好。
上面是在我学习ESPRESSO时被忽略的一个问题。
每次使用USP,我都采用了默认的ecutrho(既四倍的ecutwfc)
虽然GUI会提醒ecutrho没有设置,但是我没有做任何改动,直接计算。

后来才发现,norm-conserving时,ecutrho=4*ecutwfc是正确的设置
而USPS下,通常ecutrho=10*ecutwfc。而这个十倍是需要手动设置的
所以把自己的教训写在这里:)
3楼2009-07-17 10:08:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)


★★★ 三星级,支持鼓励

嗯,嗯, 需要注意一下使用的是什么PP,如果是USPP,那么ecutrho=(8~10)*ecutwfc
8楼2009-09-22 23:21:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
2009-09-18 20:22   回复  
 
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:34 by ZvyPK8n6Nfki
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:18 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-28 23:19 by b8fCuTHqckEt
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:45 by ez6fGg9abYaj
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +3 tianxiaoxian 2026-09-28 6/300 2026-09-28 22:44 by tianxiaoxian
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:28 by ez6fGg9abYaj
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-28 22:21 by ez6fGg9abYaj
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 17:30 by ZYXdhzDAy9ZX
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +7 zju2000 2026-09-26 8/400 2026-09-28 16:42 by beefly
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 15:32 by ZYXdhzDAy9ZX
信息提示
请填处理意见