24小时热门版块排行榜    

查看: 1952  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wh_xy

金虫 (正式写手)

[求助] 关于Wien2k的体积优化疑问求助!!

刚学习使用wien2k,按照UG的指导进行体积优化时出现问题:
我的optimize.job文件如下:
运行后仅仅得到一个数据点,为什么?
正常的是应该有5个点,并连成曲线。
请帮助!谢谢!
~~~~~~~~~~
#!/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
~~~~~~~~~~

网站上给出的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
~~~~~~~~~~
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

慕容无敌11

新虫 (初入文坛)

请问楼主,计算物质存在反演时,进行结构优化,optimize.job文件修改只是修改这两个地方,#if you want to start with dstart :
x dstart -c ,还有一个地方就是run_lapw -cc 0.0001 -p 。但是,为什么回出现下面这个错误呢?error command:  ~/lapw1cpara  -c  lapw1.def  failed。而且是有的结构会出现这个错误有的不会。楼主有没有遇到这个问题,或者请告知,体积优化时,optimize.job只这样修改这两个地方对不对?谢谢
5楼2014-06-24 10:32:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

学至老

银虫 (初入文坛)

【答案】应助回帖

★ ★
wh_xy(金币+10): 不错!学习了!谢谢! 2011-08-11 20:16:00
gzqdyouxia(金币+2): 谢谢指导 2011-08-13 11:02:12
命令行前面 # 标志为注释标志,表示其号的命令行不起作用,需要根据所需将符号去掉。另外结构文件中的原子球半径也需要降低,否则无法进行优化。
2楼2011-08-11 19:47:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leithron

禁虫 (正式写手)


fzx2008(金币+1): 谢谢指教 2011-11-29 19:20:17
本帖内容被屏蔽

3楼2011-11-09 20:52:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

学至老

银虫 (初入文坛)


fzx2008(金币+1): 谢谢交流 2011-11-29 19:20:32
重新初始化,编辑结构文件
4楼2011-11-29 17:55:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 354求调剂 +6 Tyoumou 2026-03-18 9/450 2026-03-21 20:47 by lbsjt
[考研] 材料与化工(0856)304求B区调剂 +3 邱gl 2026-03-20 7/350 2026-03-21 19:05 by 15709483992
[考研] 0805 316求调剂 +3 大雪深藏 2026-03-18 3/150 2026-03-21 18:55 by 学员8dgXkO
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-18 3/150 2026-03-21 18:17 by 学员8dgXkO
[考研] 317求调剂 +9 申子申申 2026-03-19 15/750 2026-03-21 17:31 by 学员8dgXkO
[考研] 材料与化工(0856)304求 B区 调剂 +3 邱gl 2026-03-21 3/150 2026-03-21 13:47 by lature00
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 南昌大学材料专硕311分求调剂 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[考研] 324分 085600材料化工求调剂 +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 085410人工智能专硕317求调剂(0854都可以) +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考研] 生物学调剂招人!!! +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 288求调剂,一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
信息提示
请填处理意见