| 20 | 1/1 | 返回列表 |
| 查看: 3813 | 回復(fù): 19 | ||
| 【獎(jiǎng)勵(lì)】 本帖被評(píng)價(jià)10次,作者xirainbow增加金幣 6.25 個(gè) | ||
xirainbow木蟲 (正式寫手)
|
[資源]
【原創(chuàng)】ecutwfc and ecutrho
|
|
|
以下內(nèi)容來自我的博客: 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 ] |
木蟲 (正式寫手)
|
wuchenwf ,您好。 上面是在我學(xué)習(xí)ESPRESSO時(shí)被忽略的一個(gè)問題。 每次使用USP,我都采用了默認(rèn)的ecutrho(既四倍的ecutwfc) 雖然GUI會(huì)提醒ecutrho沒有設(shè)置,但是我沒有做任何改動(dòng),直接計(jì)算。 后來才發(fā)現(xiàn),norm-conserving時(shí),ecutrho=4*ecutwfc是正確的設(shè)置 而USPS下,通常ecutrho=10*ecutwfc。而這個(gè)十倍是需要手動(dòng)設(shè)置的 所以把自己的教訓(xùn)寫在這里:) |
榮譽(yù)版主 (職業(yè)作家)
捐助貴賓 (正式寫手)
金蟲 (著名寫手)
木蟲 (正式寫手)
|
呵呵。我也覺得有點(diǎn)麻煩。也許我表達(dá)的不夠準(zhǔn)確。總之看起來的意思就是,如果使用USPP,當(dāng)ecutrho足夠大的時(shí)候,其四分之一大小的ecutwfc就顯得浪費(fèi)了,這個(gè)時(shí)候降低ecutwfc,也可以滿足收斂。這個(gè)意思其實(shí)與“ecutrho必須大于ecutwfc的四倍”是倒過來理解的。不過我還沒有真正試過。供參考。這是以前記在本子上的,我盡量再找找出處。 [ Last edited by 老虎大王 on 2010-5-12 at 22:18 ] |
木蟲 (正式寫手)
|
呵呵,上述說法的出處我找到了,原來是在Pw_form里。 http://www.democritos.it/piperma ... ecember/003396.html 我原來看到,就隨手記在本子上了。不過我自己沒有試過,僅供參考啊。 |










| 20 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +6 | tcxiaoxx 2026-03-20 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研]
|
酥酥魚.. 2026-03-21 | 3/150 |
|
|
[考研] 280求調(diào)劑 +11 | 咕嚕曉曉 2026-03-18 | 12/600 |
|
|
[考研] 求助 +5 | 夢(mèng)里的無言 2026-03-21 | 6/300 |
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
[考研] 321求調(diào)劑 +9 | 何潤采123 2026-03-18 | 11/550 |
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 材料與化工專碩調(diào)劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
[考研] 281求調(diào)劑(0805) +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|