应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【原创】写了个集群下应用测试K点或截断能的脚本
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 4036  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【奖励】 本帖被评价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
回复此楼

» 收录本帖的淘帖专辑推荐

材料计算模拟实用技巧 VASP

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-10-24 18:28:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fan159

金虫 (正式写手)
因为我也不懂VASP
赞一下 回复此楼
6楼2010-10-24 18:54:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答

fan159

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
看不懂
是不是多加点注释比较好啊？
赞一下(1人) 回复此楼
2楼2010-10-24 18:32:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhang668

木虫 (正式写手)
得到的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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhang668

木虫 (正式写手)
引用回帖:
Originally posted by fan159 at 2010-10-24 18:32:43:
看不懂
是不是多加点注释比较好啊？

这个我觉得还好吧，VASP中有类似的脚本，哪里不好懂？
赞一下 回复此楼
4楼2010-10-24 18:43:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[材料工程] 一志愿C9材料与化工专业总分300求调剂 +8 曼111 2026-03-24 9/450 2026-03-28 07:58 by YYYYX1234
[考研] 265求调剂 +8 小木虫085600 2026-03-27 8/400 2026-03-27 22:16 by 无际的草原
[考研] 265求调剂11408 +3 刘小鹿lu 2026-03-27 3/150 2026-03-27 20:53 by nihaoar
[考研] 287求调剂 +10 land xuxu 2026-03-26 10/500 2026-03-27 15:33 by 帕尔马拉特
[考研] 22408 359分调剂 +3 Qshers 2026-03-27 3/150 2026-03-27 12:22 by wxiongid
[考研] 333求调剂 +3 question挽风 2026-03-23 3/150 2026-03-27 11:29 by 不吃魚的貓
[考研] 求调剂 +3 刘柯@ 2026-03-24 4/200 2026-03-27 11:28 by shangxh
[考研] 求调剂，一志愿 南京航空航天大学大学 ，080500材料科学与工程学硕 +4 @taotao 2026-03-26 5/250 2026-03-27 08:10 by hypershenger
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-26 3/150 2026-03-27 07:58 by chemisry
[考研] 343求调剂 +4 赠我一本书 2026-03-23 4/200 2026-03-27 00:40 by wxiongid
[考研] 325求调剂 +3 Aoyijiang 2026-03-23 3/150 2026-03-26 20:46 by 不吃魚的貓
[考研] 机械学硕310分，数一英一，一志愿211本科双非找调剂信息 +3 @357 2026-03-25 3/150 2026-03-26 16:34 by by.MENG
[考研] 291 求调剂 +7 化工2026届毕业� 2026-03-21 8/400 2026-03-26 11:25 by AlenQIN.
[考研] 化学调剂一志愿上海交通大学336分-本科上海211 +4 小鱼爱有机 2026-03-25 4/200 2026-03-26 10:19 by aa331100
[考研] 一志愿天津大学339材料与化工求调剂 +3 江往卖鱼 2026-03-26 3/150 2026-03-26 09:42 by 王小欠i
[考研] 334分 一志愿武理-080500 材料求调剂 +4 李李不服输 2026-03-25 4/200 2026-03-25 21:26 by 星空星月
[考研] 考研一志愿苏州大学初始315（英一）求调剂 +3 sbdksD 2026-03-24 4/200 2026-03-25 18:16 by xcjcqu
[考研] 0854人工智能方向招收调剂 +4 章小鱼567 2026-03-24 4/200 2026-03-25 13:29 by 2177681040
[有机交流] 有机合成求助 20+3 FENGSHUJEI 2026-03-23 5/250 2026-03-24 19:31 by 88817753
[考研] 305分求调剂（食品工程） +5 Sxy112 2026-03-21 7/350 2026-03-24 12:27 by 544594351
信息提示
关闭
请填处理意见
关闭