24小时热门版块排行榜    

查看: 3807  |  回复: 18
【奖励】 本帖被评价4次,作者zhang668增加金币 3.2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zhang668

木虫 (正式写手)


[资源] 【原创】写了个集群下应用测试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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

微6凉

铁虫 (小有名气)


引用回帖:
17楼: Originally posted by zhang668 at 2011-01-05 10:00:02
这个是pwscf的脚本。。。。
至于VASP,测试K点的话,就不需要KPOINTS, 同样测试截断能,不需要INCAR。因为这些文件都会在脚本里生成。直接提交脚本就行了。

楼主您好,我想请问一下提交脚本用什么命令呢
19楼2015-09-02 16:37:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhang668 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见