| 查看: 2185 | 回复: 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 ~~~~~~~~~~ |
回复此楼» 猜你喜欢
朋友圈看到的
已经有9人回复
-
能否退出参与的面上项目解除限项
已经有13人回复
-
让我中一个面上吧！
已经有16人回复
-
放榜前的不淡定
已经有17人回复
-
科研孤儿太难了
已经有19人回复
-
明天应该可查了！？
已经有4人回复
-
什么时候开奖？
已经有11人回复
-
范进中举一文的中心思想
已经有6人回复
-
建议基金发布提前给出明确的时间点
已经有14人回复
-
2026国自然函评费到账
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
★
fzx2008(金币+1): 谢谢指教 2011-11-29 19:20:17
fzx2008(金币+1): 谢谢指教 2011-11-29 19:20:17
|
本帖内容被屏蔽 |
3楼2011-11-09 20:52:00
2楼2011-08-11 19:47:46
4楼2011-11-29 17:55:36
5楼2014-06-24 10:32:04
20