24小时热门版块排行榜    

查看: 2193  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +14 爱看书的可乐 2026-08-26 16/800 2026-08-26 16:47 by skyjudy
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 基金不中,共勉 +3 eulota 2026-08-26 3/150 2026-08-26 16:17 by 路yyyy
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +6 process2012 2026-08-23 9/450 2026-08-26 15:33 by process2012
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 只有每年这种时候来逛逛小木虫 +26 yaoyewhu2008 2026-08-20 28/1400 2026-08-26 13:49 by zzuzxg
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见