24小時(shí)熱門版塊排行榜

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

zhang668

木蟲(chóng) (正式寫(xiě)手)

1ST強(qiáng)帖: 3
應(yīng)助: 11 (小學(xué)生)
貴賓: 0.061
金幣: 1735
散金: 1294
紅花: 11
帖子: 843
在線: 324.9小時(shí)
蟲(chóng)號(hào): 733242
注冊(cè): 2009-03-28
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[交流] 【求助】截?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 ]

回復(fù)此樓

» 猜你喜歡

303求調(diào)劑已經(jīng)有5人回復(fù)
一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑已經(jīng)有7人回復(fù)
材料學(xué)求調(diào)劑已經(jīng)有5人回復(fù)
考研調(diào)劑已經(jīng)有4人回復(fù)
281求調(diào)劑已經(jīng)有4人回復(fù)
0805 316求調(diào)劑已經(jīng)有6人回復(fù)
085601求調(diào)劑總分293英一數(shù)二已經(jīng)有3人回復(fù)
08工學(xué)調(diào)劑已經(jīng)有17人回復(fù)
340求調(diào)劑已經(jīng)有4人回復(fù)
311求調(diào)劑已經(jīng)有3人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

關(guān)于收斂性測(cè)試問(wèn)題已經(jīng)有6人回復(fù)
收斂測(cè)試，這種情況我該繼續(xù)測(cè)試么,簡(jiǎn)直就是跳躍式的嘛已經(jīng)有24人回復(fù)
關(guān)于dmol3的收斂性測(cè)試已經(jīng)有12人回復(fù)
難道每次計(jì)算前都咬進(jìn)行收斂性測(cè)試嗎？已經(jīng)有14人回復(fù)
【求助】K點(diǎn)取得越大越不收斂，什么原因！已經(jīng)有20人回復(fù)
【求助成功】castep 收斂性測(cè)試已經(jīng)有22人回復(fù)
【已解決】【求助】為何PWscf不能得到類似于CASTEP那樣的收斂測(cè)試曲線呢？已經(jīng)有14人回復(fù)
【其他】關(guān)于收斂測(cè)試的問(wèn)題已經(jīng)有22人回復(fù)
【求助】收斂測(cè)試怎么做？？？已經(jīng)有17人回復(fù)
【求助】請(qǐng)問(wèn)收斂測(cè)試的標(biāo)準(zhǔn) 已經(jīng)有15人回復(fù)
【求助】收斂性測(cè)試的問(wèn)題已經(jīng)有40人回復(fù)
【求助】castep收斂性測(cè)試問(wèn)題已經(jīng)有20人回復(fù)
【求助】請(qǐng)教PWSCF晶體建模和自洽計(jì)算，不收斂已經(jīng)有21人回復(fù)
【求助】是不是什么計(jì)算都得進(jìn)行收斂測(cè)試？已經(jīng)有8人回復(fù)

1樓 2010-10-27 11:49:52

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ylli

木蟲(chóng) (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1558.8
紅花: 3
帖子: 267
在線: 41.6小時(shí)
蟲(chóng)號(hào): 277021
注冊(cè): 2006-09-03
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

BaTiO3是金屬???

★ ★
zhang668(金幣+1):謝謝參與
aylayl08(金幣+1):謝謝提點(diǎn) 2010-10-27 20:08:45
zhang668(金幣+1): 2010-10-27 23:08:09

注意BaTiO3是金屬嗎? occupations = 'smearing', smearing = 'gaussian', degauss = 0.01
這幾個(gè)參數(shù)只對(duì)金屬使用!!

贊一下(2人)

回復(fù)此樓

3樓2010-10-27 12:33:36

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 3 個(gè)回答

zhang668

木蟲(chóng) (正式寫(xiě)手)

1ST強(qiáng)帖: 3
應(yīng)助: 11 (小學(xué)生)
貴賓: 0.061
金幣: 1735
散金: 1294
紅花: 11
帖子: 843
在線: 324.9小時(shí)
蟲(chóng)號(hào): 733242
注冊(cè): 2009-03-28
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

贊一下

回復(fù)此樓

2樓2010-10-27 12:02:14

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 3 個(gè)回答

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 資源與環(huán)境調(diào)劑申請(qǐng)(333分) +7	holy J 2026-03-21	7/350	2026-03-24 17:24 by xiaohai104
[考研] 材料學(xué)碩，求調(diào)劑 6+3	糖葫蘆888ll 2026-03-22	7/350	2026-03-24 17:11 by hello七七
[考研] 材料考研調(diào)劑生 +3	黃粱一夢(mèng)千年 2026-03-24	3/150	2026-03-24 17:00 by barlinike
[考研] 292求調(diào)劑 +4	鵝鵝鵝額額額額?/a> 2026-03-24	4/200	2026-03-24 16:41 by peike
[考博] 26申博自薦 +3	whh869393 2026-03-24	3/150	2026-03-24 09:55 by 21018060
[考研] 335求調(diào)劑 +4	yuyu宇 2026-03-23	5/250	2026-03-23 23:49 by Txy@872106
[考研] 材料專業(yè)求調(diào)劑 +11	hanamiko 2026-03-18	11/550	2026-03-23 23:12 by peike
[考研] 工科0856求調(diào)劑 +5	沐析汀汀 2026-03-21	5/250	2026-03-23 17:56 by 海瑟薇-
[考研] 291 求調(diào)劑 +4	化工2026屆畢業(yè)?/a> 2026-03-21	5/250	2026-03-23 16:46 by 化工2026屆畢業(yè)?/a>
[考研] 352求調(diào)劑 +3	大米飯！ 2026-03-22	3/150	2026-03-22 23:28 by king123！
[考研] 275求調(diào)劑 +6	shansx 2026-03-22	8/400	2026-03-22 15:27 by barlinike
[考研] 0703化學(xué)調(diào)劑 +4	妮妮ninicgb 2026-03-21	4/200	2026-03-21 18:39 by 學(xué)員8dgXkO
[考研] 材料學(xué)碩333求調(diào)劑 +3	北道巷 2026-03-18	3/150	2026-03-21 18:17 by 學(xué)員8dgXkO
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7	墨墨漠 2026-03-20	7/350	2026-03-21 16:36 by barlinike
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9	Doleres 2026-03-19	9/450	2026-03-20 22:36 by JourneyLucky
[考研] 一志愿西南交通專碩材料355 本科雙非求調(diào)劑 +5	西南交通專材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 353求調(diào)劑 +3	拉鉤不許變 2026-03-20	3/150	2026-03-20 19:56 by JourneyLucky
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3	葵梓衛(wèi)隊(duì) 2026-03-20	3/150	2026-03-20 19:28 by zhukairuo
[考研] 招收調(diào)劑碩士 +4	lidianxing 2026-03-19	12/600	2026-03-20 12:25 by lidianxing
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3	小木蟲(chóng)200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu