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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【原创】ecutwfc and ecutrho
201/1返回列表
查看: 3640  |  回复: 19
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【奖励】 本帖被评价10次，作者xirainbow增加金币 6.25 个

xirainbow

木虫 (正式写手)

[资源] 【原创】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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2009-07-16 11:30:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

xirainbow

木虫 (正式写手)
★ ★ ★
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。而这个十倍是需要手动设置的
所以把自己的教训写在这里：）
赞一下(17人) 回复此楼
3楼2009-07-17 10:08:30
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

wuchenwf

荣誉版主 (职业作家)
  貌似楼主发的是自己解决的一个经验
赞一下 回复此楼
2楼2009-07-16 22:32:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjr

铁杆木虫 (正式写手)

★★★ 三星级,支持鼓励

你是牛人呀！！
回复此楼
6楼2009-09-17 10:20:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)

★★★ 三星级,支持鼓励

嗯,嗯, 需要注意一下使用的是什么PP,如果是USPP,那么ecutrho=(8~10)*ecutwfc
赞一下 回复此楼
8楼2009-09-22 23:21:42
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaoguoying1981

新虫 (小有名气)
★ ★
zzy870720z(金币+2):感谢交流 2010-05-13 08:39:50
对于USPS，ecutrho必需大于4*ecutwfc。具体应该是多少是需要测的吧
就像测ecut和kmesh一样
赞一下(1人) 回复此楼
9楼2009-12-03 19:05:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enola

捐助贵宾 (正式写手)

★★★ 三星级,支持鼓励

楼主辛苦呢，我们基本4
赞一下 回复此楼
10楼2010-05-12 15:18:20
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★
zzy870720z(金币+1):感谢交流 2010-05-13 08:40:41
USPP一般8被足够了！
你不设置一般来说
也不会错，但是结果的精度可能会差一点。
赞一下(1人) 回复此楼
11楼2010-05-12 17:50:59
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

★★★ 三星级,支持鼓励

我有一个笔记，记了USPP情况下的ecutwfc 和 ecutrho的设置测试方法，但没有记录出处：
首先按照默认的 ecutrho=4* ecutwfc ，并不断增大ecutwfc，直至收敛，然后把由此定下来的ecutrho保持住，再不断减小ecutwfc，找到可以让结果收敛的最小ecutwfc值。

另外一个常用的设置方法是让ecutrho=（8~10）* ecutwfc，一般8倍就差不多了。
赞一下(1人) 回复此楼
12楼2010-05-12 18:11:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)
引用回帖:
Originally posted by 老虎大王 at 2010-05-12 18:11:38:
我有一个笔记，记了USPP情况下的ecutwfc 和 ecutrho的设置测试方法，但没有记录出处：
首先按照默认的 ecutrho=4* ecutwfc ，并不断增大ecutwfc，直至收敛，然后把由此定下来的ecutrho保持住，再不断减小ecutwfc ...

“首先按照默认的 ecutrho=4* ecutwfc ，并不断增大ecutwfc，直至收敛，然后把由此定下来的ecutrho保持住，再不断减小ecutwfc，找到可以让结果收敛的最小ecutwfc值。”
这句话有意思，头一次听说：）
赞一下 回复此楼
13楼2010-05-12 20:54:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★
zzy870720z(金币+1):感谢iu 2010-05-13 08:41:17
呵呵。我也觉得有点麻烦。也许我表达的不够准确。总之看起来的意思就是，如果使用USPP，当ecutrho足够大的时候，其四分之一大小的ecutwfc就显得浪费了，这个时候降低ecutwfc，也可以满足收敛。这个意思其实与“ecutrho必须大于ecutwfc的四倍”是倒过来理解的。不过我还没有真正试过。供参考。这是以前记在本子上的，我尽量再找找出处。

[ Last edited by 老虎大王 on 2010-5-12 at 22:18 ]
赞一下(1人) 回复此楼
14楼2010-05-12 22:12:42
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)
引用回帖:
Originally posted by 老虎大王 at 2010-05-12 22:12:42:
呵呵。我也觉得有点麻烦。也许我表达的不够准确。总之看起来的意思就是，如果使用USPP，当ecutrho足够大的时候，其四分之一大小的ecutwfc就显得浪费了，这个时候降低ecutwfc，也可以满足收敛。这个意思其实与“ec ...

解释的非常好：P
谢谢啦：）
赞一下 回复此楼
15楼2010-05-12 23:09:12
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★
zzy870720z(金币+1):感谢交流 2010-05-13 12:14:22
呵呵，上述说法的出处我找到了，原来是在Pw_form里。
http://www.democritos.it/piperma ... ecember/003396.html
我原来看到，就随手记在本子上了。不过我自己没有试过，仅供参考啊。
赞一下(1人) 回复此楼
16楼2010-05-13 09:08:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wutianhuan

铁杆木虫 (著名写手)

★★★★★ 五星级,优秀推荐

要注意的地方真多啊，谢谢了
赞一下 回复此楼
17楼2010-05-13 15:26:28
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

157806400

金虫 (初入文坛)
★
cenwanglai(金币+1):谢谢参与讨论 2010-05-30 21:51:25
超软赝势必须要4倍以上吗？我发过一篇文章都用的6倍左右，我也没作测试，不会是错误的吧？不过我增大ecutrho测试收敛了。用的是6倍的比例。
赞一下(1人) 回复此楼
18楼2010-05-30 21:24:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
lvjian85964楼
2009-07-17 11:26   回复  
 
freshgirl5楼
2009-07-17 14:13   回复  
 
sandongyuwang7楼
2009-09-18 20:22   回复  
 
茶味果冻19楼
2014-10-22 19:05   回复  
一般  顶一下，感谢分享！
吕小强强20楼
2018-01-09 11:29   回复  
三星好评  顶一下！～
相关版块跳转 我要订阅楼主 xirainbow 的主题更新
201/1返回列表
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭