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

wh_xy

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

應(yīng)助: 1 (幼兒園)
金幣: 1022.9
散金: 15
紅花: 2
帖子: 880
在線(xiàn): 178.3小時(shí)
蟲(chóng)號(hào): 499004
注冊(cè): 2008-02-05
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] 關(guān)于Wien2k的體積優(yōu)化疑問(wèn)求助�。�

剛學(xué)習(xí)使用wien2k，按照UG的指導(dǎo)進(jìn)行體積優(yōu)化時(shí)出現(xiàn)問(wèn)題：
我的optimize.job文件如下：
運(yùn)行后僅僅得到一個(gè)數(shù)據(jù)點(diǎn)，為什么？
正常的是應(yīng)該有5個(gè)點(diǎn)，并連成曲線(xiàn)。
請(qǐng)幫助！謝謝！
～～～～～～～～～～
#!/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
~~~~~~~~~~

回復(fù)此樓

» 猜你喜歡

求調(diào)劑已經(jīng)有3人回復(fù)
265求調(diào)劑已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有6人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有12人回復(fù)
286分人工智能專(zhuān)業(yè)請(qǐng)求調(diào)劑愿意跨考！已經(jīng)有3人回復(fù)
329求調(diào)劑已經(jīng)有5人回復(fù)
申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變，給編輯又寫(xiě)郵件了，沒(méi)回復(fù) 已經(jīng)有4人回復(fù)
材料學(xué)碩318求調(diào)劑已經(jīng)有5人回復(fù)
一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑已經(jīng)有6人回復(fù)
081700化工學(xué)碩調(diào)劑已經(jīng)有3人回復(fù)

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

Wien2k運(yùn)算過(guò)程lapw1總是出現(xiàn)錯(cuò)誤已經(jīng)有7人回復(fù)
wien2k的w2web的一點(diǎn)問(wèn)題已經(jīng)有29人回復(fù)
[有志者分享]畢業(yè)前收集的學(xué)習(xí)資料——Wien2k篇已經(jīng)有55人回復(fù)
wien2k初始化出現(xiàn)的問(wèn)題已經(jīng)有3人回復(fù)
第一次使用wien2k的入門(mén)級(jí)問(wèn)題已經(jīng)有8人回復(fù)
wien2k編譯出錯(cuò) 已經(jīng)有6人回復(fù)
wien2k在終端下計(jì)算能帶，改怎么做？已經(jīng)有6人回復(fù)
WIEN2k的能帶數(shù)據(jù)文件case.spaghetti_ene。你們都是怎么分析提取數(shù)據(jù)？和我的看法。已經(jīng)有8人回復(fù)
并行WIEN2K安裝失敗，如何解決，謝謝！已經(jīng)有9人回復(fù)
在vasp計(jì)算里，有沒(méi)有好的只優(yōu)化原子位置和晶胞體積的方法？已經(jīng)有8人回復(fù)
VASP優(yōu)化結(jié)構(gòu)遇到的問(wèn)題。已經(jīng)有9人回復(fù)
利用CASTEP對(duì)建立的晶胞進(jìn)行幾何優(yōu)化的問(wèn)題（急�。。。� 已經(jīng)有13人回復(fù)
【求助】Charge Density of WIEN2K 已經(jīng)有8人回復(fù)
【求助】安裝wien2k提示“Variable name must begin with a letter” 已經(jīng)有13人回復(fù)
【求助】wien2k中K-point個(gè)數(shù) 已經(jīng)有9人回復(fù)
【求助】簡(jiǎn)單問(wèn)題----關(guān)于vasp晶格優(yōu)化已經(jīng)有7人回復(fù)

1樓 2011-07-31 14:14:07

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

leithron

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

★
fzx2008(金幣+1): 謝謝指教 2011-11-29 19:20:17

本帖內(nèi)容被屏蔽

3樓2011-11-09 20:52:00

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

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

學(xué)至老

銀蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 233.8
帖子: 33
在線(xiàn): 36.5小時(shí)
蟲(chóng)號(hào): 752477
注冊(cè): 2009-04-19
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

【答案】應(yīng)助回帖

★ ★
wh_xy(金幣+10): 不錯(cuò)！學(xué)習(xí)了！謝謝！ 2011-08-11 20:16:00
gzqdyouxia(金幣+2): 謝謝指導(dǎo) 2011-08-13 11:02:12

命令行前面 # 標(biāo)志為注釋標(biāo)志，表示其號(hào)的命令行不起作用，需要根據(jù)所需將符號(hào)去掉。另外結(jié)構(gòu)文件中的原子球半徑也需要降低，否則無(wú)法進(jìn)行優(yōu)化。

贊一下(1人)

回復(fù)此樓

2樓2011-08-11 19:47:46

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

學(xué)至老

銀蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 233.8
帖子: 33
在線(xiàn): 36.5小時(shí)
蟲(chóng)號(hào): 752477
注冊(cè): 2009-04-19
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
fzx2008(金幣+1): 謝謝交流 2011-11-29 19:20:32

重新初始化，編輯結(jié)構(gòu)文件

贊一下(1人)

回復(fù)此樓

4樓2011-11-29 17:55:36

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

慕容無(wú)敵11

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 170.3
帖子: 7
在線(xiàn): 4.6小時(shí)
蟲(chóng)號(hào): 2925515
注冊(cè): 2014-01-10
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

請(qǐng)問(wèn)樓主，計(jì)算物質(zhì)存在反演時(shí)，進(jìn)行結(jié)構(gòu)優(yōu)化，optimize.job文件修改只是修改這兩個(gè)地方，#if you want to start with dstart :
x dstart -c ,還有一個(gè)地方就是run_lapw -cc 0.0001 -p 。但是，為什么回出現(xiàn)下面這個(gè)錯(cuò)誤呢？error command： ~/lapw1cpara -c lapw1.def failed。而且是有的結(jié)構(gòu)會(huì)出現(xiàn)這個(gè)錯(cuò)誤有的不會(huì)。樓主有沒(méi)有遇到這個(gè)問(wèn)題，或者請(qǐng)告知，體積優(yōu)化時(shí)，optimize.job只這樣修改這兩個(gè)地方對(duì)不對(duì)？謝謝

贊一下

回復(fù)此樓

5樓2014-06-24 10:32:04

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

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

[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +9	Ymlll 2026-03-18	12/600	2026-03-20 00:01 by 23Postgrad
[考研] 307求調(diào)劑 +9	冷笙123 2026-03-17	9/450	2026-03-19 22:44 by 學(xué)員8dgXkO
[考研] 求調(diào)劑 +3	Mqqqqqq 2026-03-19	3/150	2026-03-19 14:11 by peike
[考研] 材料考研調(diào)劑 +3	xwt。 2026-03-19	3/150	2026-03-19 11:22 by w沐陽(yáng)w
[考研] 材料專(zhuān)碩英一數(shù)二306 +5	z1z2z3879 2026-03-18	5/250	2026-03-19 07:43 by BruceLiu320
[考研] 化學(xué)工程321分求調(diào)劑 +15	大米飯！ 2026-03-15	18/900	2026-03-18 14:52 by haxia
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專(zhuān)利，省、國(guó)級(jí)證書(shū)40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +7	手機(jī)用戶(hù) 2026-03-17	7/350	2026-03-18 14:34 by vgtyfty
[考研] 297求調(diào)劑 +8	戲精丹丹丹 2026-03-17	8/400	2026-03-18 14:30 by laoshidan
[考研] 299求調(diào)劑 +5	△小透明* 2026-03-17	5/250	2026-03-18 11:49 by 盡舜堯1
[考研] 303求調(diào)劑 +4	睿08 2026-03-17	6/300	2026-03-18 11:01 by Iveryant
[考研] 334求調(diào)劑 +3	志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[基金申請(qǐng)] 今年的國(guó)基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤(rùn)色
[考研] 0703一志愿211 285分求調(diào)劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6	太想進(jìn)步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 中科大材料專(zhuān)碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu
[碩博家園] 085600 260分求調(diào)劑 +3	天空還下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空還下雨么
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu