24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wh_xy

金虫 (正式写手)

应助: 1 (幼儿园)
金币: 1022.9
散金: 15
红花: 2
帖子: 880
在线: 178.3小时
虫号: 499004
注册: 2008-02-05
专业: 凝聚态物性 II ：电子结构

[求助] 关于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
~~~~~~~~~~

回复此楼

» 猜你喜欢

请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有6人回复
孩子确诊有中度注意力缺陷已经有14人回复
三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

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

Wien2k运算过程lapw1总是出现错误已经有7人回复
wien2k的w2web的一点问题已经有29人回复
[有志者分享]毕业前收集的学习资料——Wien2k篇已经有55人回复
wien2k初始化出现的问题已经有3人回复
第一次使用wien2k的入门级问题已经有8人回复
wien2k编译出错已经有6人回复
wien2k在终端下计算能带，改怎么做？已经有6人回复
WIEN2k的能带数据文件case.spaghetti_ene。你们都是怎么分析提取数据？和我的看法。已经有8人回复
并行WIEN2K安装失败，如何解决，谢谢！已经有9人回复
在vasp计算里，有没有好的只优化原子位置和晶胞体积的方法？已经有8人回复
VASP优化结构遇到的问题。已经有9人回复
利用CASTEP对建立的晶胞进行几何优化的问题（急！！！）已经有13人回复
【求助】Charge Density of WIEN2K 已经有8人回复
【求助】安装wien2k提示“Variable name must begin with a letter” 已经有13人回复
【求助】wien2k中K-point个数已经有9人回复
【求助】简单问题----关于vasp晶格优化已经有7人回复

1楼2011-07-31 14:14:07

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

慕容无敌11

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 170.3
帖子: 7
在线: 4.6小时
虫号: 2925515
注册: 2014-01-10
专业: 凝聚态物性 II ：电子结构

请问楼主，计算物质存在反演时，进行结构优化，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只这样修改这两个地方对不对？谢谢