24小时热门版块排行榜    

查看: 1943  |  回复: 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的回帖

学至老

银虫 (初入文坛)


fzx2008(金币+1): 谢谢交流 2011-11-29 19:20:32
重新初始化,编辑结构文件
4楼2011-11-29 17:55:36
已阅   回复此楼   关注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的回帖

慕容无敌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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工调剂 324分 +8 llllkkkhh 2026-03-18 8/400 2026-03-18 21:01 by Catalysis25
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-16 10/500 2026-03-18 20:41 by Wangjingyue
[考研] 材料专业求调剂 +5 hanamiko 2026-03-18 5/250 2026-03-18 20:19 by 星空星月
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 295求调剂 +3 一志愿京区211 2026-03-18 5/250 2026-03-18 17:03 by zhaoqian0518
[考研] 0703化学调剂 +4 pupcoco 2026-03-17 7/350 2026-03-18 12:14 by djl2006
[考研] 280求调剂 +6 咕噜晓晓 2026-03-18 7/350 2026-03-18 11:25 by 无际的草原
[考研] 301求调剂 +9 yy要上岸呀 2026-03-17 9/450 2026-03-18 08:58 by 无际的草原
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 268求调剂 +7 好运连绵不绝 2026-03-12 8/400 2026-03-17 20:28 by xilongliang
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[考研] 278求调剂 +3 Yy7400 2026-03-13 3/150 2026-03-17 08:24 by laoshidan
[考研] 一志愿,福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 328求调剂 +3 5201314Lsy! 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
信息提示
请填处理意见