24小时热门版块排行榜    

查看: 2125  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 jRl3mE6ddZGq 2026-08-03 7/350 2026-08-04 03:35 by k7OM8YghWbkC
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 03:34 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 5/250 2026-08-04 03:30 by k7OM8YghWbkC
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 03:27 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 4/200 2026-08-04 03:25 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 12/600 2026-08-04 03:17 by k7OM8YghWbkC
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 02:10 by Ocx3OKOzGU4q
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 8/400 2026-08-03 16:02 by yanghuping
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 16/800 2026-08-03 12:23 by fuweiguochen
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 12/600 2026-08-03 08:22 by syl200707
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
请填处理意见