| 20 | 1/1 | 返回列表 |
| 查看: 3809 | 回復(fù): 19 | ||
| 【獎勵】 本帖被評價10次,作者xirainbow增加金幣 6.25 個 | ||
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時被忽略的一個問題。 每次使用USP,我都采用了默認(rèn)的ecutrho(既四倍的ecutwfc) 雖然GUI會提醒ecutrho沒有設(shè)置,但是我沒有做任何改動,直接計(jì)算。 后來才發(fā)現(xiàn),norm-conserving時,ecutrho=4*ecutwfc是正確的設(shè)置 而USPS下,通常ecutrho=10*ecutwfc。而這個十倍是需要手動設(shè)置的 所以把自己的教訓(xùn)寫在這里:) |
金蟲 (著名寫手)
木蟲 (正式寫手)
|
呵呵。我也覺得有點(diǎn)麻煩。也許我表達(dá)的不夠準(zhǔn)確?傊雌饋淼囊馑季褪,如果使用USPP,當(dāng)ecutrho足夠大的時候,其四分之一大小的ecutwfc就顯得浪費(fèi)了,這個時候降低ecutwfc,也可以滿足收斂。這個意思其實(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ā)表 | |
|---|---|---|---|---|
|
[考研] 319求調(diào)劑 +4 | 小力氣珂珂 2026-03-20 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +3 | 苑豆豆 2026-03-20 | 3/150 |
|
|
[考研] 考研調(diào)劑 +4 | 來好運(yùn)來來來 2026-03-21 | 4/200 |
|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +5 | 生物工程調(diào)劑 2026-03-17 | 9/450 |
|
|
[考研] 333求調(diào)劑 +5 | 87639 2026-03-21 | 7/350 |
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[考研] 工科0856求調(diào)劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
[考研] 一志愿中國石油大學(xué)(華東) 本科齊魯工業(yè)大學(xué) +3 | 石能偉 2026-03-17 | 3/150 |
|
|
[考研] 307求調(diào)劑 +10 | 冷笙123 2026-03-17 | 10/500 |
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
[考研] 261求B區(qū)調(diào)劑,科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
[考研] 08工學(xué)調(diào)劑 +5 | 用戶573181 2026-03-20 | 5/250 |
|