| 3 | 1/1 | 返回列表 |
| 查看: 1163 | 回復(fù): 2 | |||
zhang668木蟲 (正式寫手)
|
[交流]
【求助】截?cái)嗄軠y(cè)試不收斂是怎么回事 已有1人參與
|
|
如題,我用pwscf測(cè)試了BaTiO3的截?cái)嗄,結(jié)果發(fā)現(xiàn)截?cái)嗄芎茈y收斂,基本要到450Ry才能達(dá)到1.0e-4Ry精度,我同樣測(cè)試了其它參數(shù)的影響,發(fā)現(xiàn)不論怎么換其他參數(shù)都是這個(gè)結(jié)果。 測(cè)試腳本: #!/bin/sh for i in 60 70 80 90 100 110 120 130 140 ...... 570 do cat > scf.in < &control calculation= 'scf' restart_mode='from_scratch' pseudo_dir = './' outdir ='./' prefix = 'BaTiO3' tprnfor = .true. / &system ibrav = 1, celldm(1)=7.57, nat=5, ntyp = 3, ecutwfc = $i, ecutrho = 600,occupations = 'smearing', smearing = 'gaussian', degauss = 0.01 / &electrons diagonalization = 'david' conv_thr = 1.0e-8 mixing_beta = 0.3 / ATOMIC_SPECIES Ba 137.327 Ba.pbe-nsp-van.UPF Ti 47.867 Ti.pbe-sp-van_ak.UPF O 15.9994 O.pbe-van_ak.UPF ATOMIC_POSITIONS Ba 0 0 0 Ti 0.5 0.5 0.5 O 0.5 0.5 0 O 0.5 0 0.5 O 0 0.5 0.5 K_POINTS AUTOMATIC 7 7 7 0 0 0 ! pw.x -inp scf.in E=`grep 'Harris-Foulkes estimate ' output.* |tail -1 |awk '{printf "%7.6f \n", $4}'` echo $i $E >>Evk.dat done 截?cái)嗄茉俅,?80已算不下去了,請(qǐng)問(wèn)這個(gè)是怎么回事?為什么很難收斂? 難道是贗勢(shì)里Ba為非線性nsp贗勢(shì)原因還是什么?按理說(shuō)USPP贗勢(shì)截?cái)嗄軕?yīng)該很小啊 [ Last edited by zhang668 on 2010-10-27 at 23:10 ] |
木蟲 (正式寫手)
木蟲 (小有名氣)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 6/300 |
|
|
[考研] 324求調(diào)劑 +5 | lucky呀呀呀鴨 2026-03-20 | 5/250 |
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 261求B區(qū)調(diào)劑,科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
[考研] 環(huán)境工程調(diào)劑 +9 | 大可digkids 2026-03-16 | 9/450 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 266求調(diào)劑 +5 | 陽(yáng)陽(yáng)哇塞 2026-03-14 | 10/500 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 303求調(diào)劑 +4 | 睿08 2026-03-17 | 6/300 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|