| 查看: 4521 | 回复: 18 | ||||
| 【奖励】 本帖被评价4次，作者zhang668增加金币 3.2 个 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[资源]
【原创】写了个集群下应用测试K点或截断能的脚本 已有4人参与
|
||||
|
由于我用的是学校服务器，scf生成的总能等信息都写在了任务提交系统指定的结果文件output..+（变化ID） 文件中，提交任务时的指定输出没用，因此写了个在指定任务提交系统下的脚本测试截断能和K点。 #!/bin/sh for i in 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200 do cat > scf.in < &control calculation= 'scf' restart_mode='from_scratch' pseudo_dir = './' outdir ='./' prefix = 'BiAlO3' tprnfor = .true. / &system ibrav = 1, celldm(1)=7.09, nat=5, ntyp = 3, ecutwfc = $i, ecutrho = 300, / &electrons diagonalization = 'david' conv_thr = 1.0e-8 mixing_beta = 0.7 / ATOMIC_SPECIES Bi 208.98 Bi.pbe-d-mt.UPF Al 26.98 Al.pbe-sp-van.UPF O 15.9994 O.pbe-van_ak.UPF ATOMIC_POSITIONS Bi 0 0 0 Al 0.5 0.5 0.5 O 0.5 0.5 0 O 0.5 0 0.5 O 0 0.5 0.5 K_POINTS AUTOMATIC 5 5 5 0 0 0 ! pw.x -inp scf.in E=`grep 'Harris-Foulkes estimate ' output.* |tail -1 |awk '{printf "%7.6f \n", $4}'` echo $i $E >>Evi.dat done |
回复此楼» 收录本帖的淘帖专辑推荐
 材料计算模拟实用技巧 | VASP |
» 猜你喜欢
2026国自然函评费到账
已经有17人回复
-
让我中一个面上吧！
已经有15人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
什么时候开奖？
已经有10人回复
-
今天放榜吗？
已经有16人回复
-
93BebMhtakh前后11位开头都是大写
已经有8人回复
-
范进中举一文的中心思想
已经有4人回复
-
只有每年这种时候来逛逛小木虫
已经有26人回复
-
filecode，4个jtjc了
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- BET测试中回滞环变宽,什么原因？
已经有10人回复
- 中科院上海应用物理研究所核能项目招聘熔盐质量控制人员（熔盐理化组）
已经有7人回复
- 请教Windows系统下perl脚本的ARGV应用问题，谢谢
已经有5人回复
- 氨基树脂应用方面的求助
已经有23人回复
- 求助怎样写论文综述
已经有7人回复
- 写国基过程中，这样的问题你们有遇到吗？我前几年没遇到。
已经有3人回复
- 原子荧光测试遇到问题，怎么办？
已经有5人回复
- 【活动】Show 脚本，赢奖励
已经有23人回复
- UL 黄卡申请 阻燃、RTI、紫外老化的测试方法及标准
已经有131人回复
- 请教ORCA的并行及集群上的脚本文件
已经有14人回复
- 【求助成功】vasp测试K点收敛性的脚本问题
已经有20人回复
★ ★
zzy870720z(金币+2):谢谢交流 2011-01-05 14:18:12
zzy870720z(金币+2):谢谢交流 2011-01-05 14:18:12
|
个人测试截断能和晶格常数两个参数，采用与CASTEP下相同的方法，结果居然随着晶格常数增加，总能一直在降低……相当无语……提高截断能和K点都不成，但是理论上总能应该有个最低点，以表示收敛啊…… 不过值得一提的是，同样的晶体结构，同样的截断能（把Ry换算为eV的话）、K点设置，居然PWscf的速度比CASTEP快得有点难以接受…… 顺便，建议在每个单独的计算后，rm -rf $DIR_TEMP/*，以防止临时文件对后续计算有影响，至少CASTEP下偶尔就有这个毛病，相当烦人…… [ Last edited by souledge on 2011-1-1 at 20:35 ] |
13楼2011-01-01 20:33:28
2楼2010-10-24 18:32:43
|
得到的Evi.dat文件中，同一组的i，E数据重复出现次数和并行节点数一样多，这个貌似没办法解决。下面是我测试的结果：（第一列是截断能，单位Ry） 60 -397.677201 70 -397.680237 80 -397.685099 90 -397.689479 100 -397.691504 110 -397.692147 120 -397.692890 130 -397.694475 140 -397.696578 150 -397.698359 160 -397.699449 170 -397.700197 180 -397.701290 190 -397.703283 200 -397.706246 发现一个奇怪的现象，能量变化幅度不是随之截断能的增加而逐渐减小的，而是在120Ry达到最小后，就开始增加了。不知道大家有没有遇到过这种情况？是什么原因呢？ [ Last edited by zhang668 on 2010-10-24 at 18:38 ] |
3楼2010-10-24 18:36:54
4楼2010-10-24 18:43:31
40