| 7 | 1/1 | 返回列表 |
| 查看: 1940 | 回復(fù): 6 | ||
wh_xy金蟲 (正式寫手)
|
[求助]
關(guān)于Wien2k的體積優(yōu)化疑問求助。
|
|
剛學(xué)習(xí)使用wien2k,按照UG的指導(dǎo)進(jìn)行體積優(yōu)化時出現(xiàn)問題: 我的optimize.job文件如下: 運(yùn)行后僅僅得到一個數(shù)據(jù)點(diǎn),為什么? 正常的是應(yīng)該有5個點(diǎn),并連成曲線。 請幫助!謝謝! ~~~~~~~~~~ #!/bin/csh -f # Modify this script according to your needs: # Uncomment one of the lines ... to adjust # starting electron density: either use # clmextrapol # or a clmsum file from a previous run (with smaller k-mesh,...) # convergence criteria, # spin-polarization (change run_lapw to runsp_lapw) # activate min_lapw # modify the save_lapw command if (-e TiC.clmsum && ! -z TiC.clmsum) then x dstart -super endif if (-e TiC.clmup && ! -z TiC.clmup) then x dstart -super -up x dstart -super -dn endif foreach i ( \ TiC_vol_-10.0 \ TiC_vol__-5.0 \ TiC_vol___0.0 \ TiC_vol___5.0 \ TiC_vol__10.0 \ ) rm TiC.struct # NFS-bug cp $i.struct TiC.struct # Please uncomment and adapt any of the lines below according to your needs # if you have a previous optimize-run: # cp $i.clmsum TiC.clmsum # cp $i.clmup TiC.clmup # cp $i.clmdn TiC.clmdn # if you want to start with dstart: # x dstart # -c # x dstart -up # -c # x dstart -dn # -c # recommended option: use charge extrapolation clmextrapol_lapw if (-e TiC.clmup && ! -z TiC.clmup) then clmextrapol_lapw -up clmextrapol_lapw -dn endif run_lapw -ec 0.0001 # -p -it # runsp_lapw -ec 0.0001 # min -I -j "run_lapw -I -fc 1.0 -i 40 " set stat = $status if ($stat) then echo "ERROR status in" $i exit 1 endif save_lapw ${i} # save_lapw -f -d XXX $i end ~~~~~~~~~~ 網(wǎng)站上給出的optimize.job文件如下: ~~~~~~~~~~ #!/bin/csh -f # Modify this script according to your needs: # Uncomment one of the lines ... # Change run_lapw to runsp_lapw or use different convergence criterium # Change save_lapw -d XXX foreach i ( \ TiC_vol_-10.0 \ TiC_vol__-5.0 \ TiC_vol___0.0 \ TiC_vol___5.0 \ TiC_vol__10.0 \ ) cp $i.struct TiC.struct # cp $i.clmsum TiC.clmsum x dstart # run_lapw -ec 0.0001 -in1new 3 -in1orig -renorm run_lapw -ec 0.0001 set stat = $status if ($stat) then echo "ERROR status in" $i exit 1 endif save_lapw $i # save_lapw -f -d XXX $i end ~~~~~~~~~~ |
銀蟲 (初入文壇)
|
本帖內(nèi)容被屏蔽 |
銀蟲 (初入文壇)
新蟲 (初入文壇)
|
請問樓主,計算物質(zhì)存在反演時,進(jìn)行結(jié)構(gòu)優(yōu)化,optimize.job文件修改只是修改這兩個地方,#if you want to start with dstart : x dstart -c ,還有一個地方就是run_lapw -cc 0.0001 -p 。但是,為什么回出現(xiàn)下面這個錯誤呢?error command: ~/lapw1cpara -c lapw1.def failed。而且是有的結(jié)構(gòu)會出現(xiàn)這個錯誤有的不會。樓主有沒有遇到這個問題,或者請告知,體積優(yōu)化時,optimize.job只這樣修改這兩個地方對不對?謝謝 |
榮譽(yù)版主 (職業(yè)作家)
小木蟲從頭派教主

新蟲 (初入文壇)
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 295求調(diào)劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿西南交大,求調(diào)劑 +4 | 材化逐夢人 2026-03-18 | 4/200 |
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +16 | 不知道叫什么! 2026-03-15 | 18/900 |
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 0703一志愿211 285分求調(diào)劑 +5 | ly3471z 2026-03-13 | 5/250 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[考研] 265求調(diào)劑 +4 | 威化餅07 2026-03-12 | 4/200 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +7 | 綠豆芹菜湯 2026-03-12 | 7/350 |
|
|
[考研] 307求調(diào)劑 +5 | 超級伊昂大王 2026-03-12 | 5/250 |
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生,2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|