| 3 | 1/1 | 返回列表 |
| 查看: 1162 | 回復(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)問這個(gè)是怎么回事?為什么很難收斂? 難道是贗勢(shì)里Ba為非線性nsp贗勢(shì)原因還是什么?按理說USPP贗勢(shì)截?cái)嗄軕?yīng)該很小啊 [ Last edited by zhang668 on 2010-10-27 at 23:10 ] |
木蟲 (正式寫手)
木蟲 (小有名氣)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
|
然11 2026-03-19 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +3 | 77chaselx 2026-03-20 | 3/150 |
|
|
[考研]
|
簡(jiǎn)木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 3/150 |
|
|
[考研] 299求調(diào)劑 +3 | shxchem 2026-03-20 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +3 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 330求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
[考研] 【同濟(jì)軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
[考研] 303求調(diào)劑 +4 | 睿08 2026-03-17 | 6/300 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|