24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1801  |  回复: 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的回帖

学至老

铜虫 (初入文坛)

【答案】应助回帖

★ ★
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的回帖

慕容无敌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的回帖

ljw4010

荣誉版主 (职业作家)

小木虫从头派教主

你的文件没问题,run_lapw -ec 0.0001需要加-p!你这五个点没出来,是不是因为你的设置值计算了一个点就结束了,还是报错了,产看一下show dayfile
不要说话
6楼2014-12-14 13:30:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

l1009568589

新虫 (初入文坛)

请问楼主这个问题解决了吗? 我也碰到同样的问题,请楼主帮忙,多谢!
7楼2015-11-24 23:47:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wh_xy 的主题更新
信息提示
请填处理意见