| 查看: 4506 | 回复: 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 |
» 猜你喜欢
建议基金发布提前给出明确的时间点
已经有13人回复
-
2026国自然函评费到账
已经有7人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有20人回复
-
科研孤儿太难了
已经有18人回复
-
看来今天不会放榜了？
已经有11人回复
-
时间戳又变了
已经有19人回复
-
filecode，4个jtjc了
已经有14人回复
-
让我中一个面上吧！
已经有13人回复
-
人气不行了
已经有7人回复
-
感觉是下周放榜了
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- BET测试中回滞环变宽,什么原因？
已经有10人回复
- 中科院上海应用物理研究所核能项目招聘熔盐质量控制人员（熔盐理化组）
已经有7人回复
- 请教Windows系统下perl脚本的ARGV应用问题，谢谢
已经有5人回复
- 氨基树脂应用方面的求助
已经有23人回复
- 求助怎样写论文综述
已经有7人回复
- 写国基过程中，这样的问题你们有遇到吗？我前几年没遇到。
已经有3人回复
- 原子荧光测试遇到问题，怎么办？
已经有5人回复
- 【活动】Show 脚本，赢奖励
已经有23人回复
- UL 黄卡申请 阻燃、RTI、紫外老化的测试方法及标准
已经有131人回复
- 请教ORCA的并行及集群上的脚本文件
已经有14人回复
- 【求助成功】vasp测试K点收敛性的脚本问题
已经有20人回复
8楼2010-10-24 22:42:00
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
