vasp计算脚本放在服务器的位置,让搬砖变得轻松——vasp提交系列任务用的脚本...

在上adv. metal课的时候就曾经遇到过这个问题,当然那个时候还是全靠手动去做nuclei move的锅,不过现在也会遇到这样的情况,就是我要提交一系列的任务,而这一系列的任务是有一定的规律的,目前从我遇到过的case来看,主要有两种

Case 1.

这种情况每个任务里面不一样的地方只有其中一个参数,参数是单一的或者多个参数之间有某种数学关系

e.g. 要做一个convergence test,提交相同cut off energy不同k-point的一系列任务,或者POSCAR里面成比例地修改键长

用改POSCAR这个case来说,需要完成两个script:

生成POSCAR的script

首先理清思路,在我们的一系列的任务中,POSCAR是一直在变的,所以在最开始的输入文件中,我们不事先提供POSCAR这个文件,而是用每个文件夹中的script去生成跟任务文件夹对应的POSCAR文件,然后运行vasp程序->得到结果->收集结果

写一个生成POSCAR的script pos.py,这里用python写了一个,用c也可以,完全看习惯

import sys

import math

a = float(sys.argv[1]) # lattice constant

fout = open('POSCAR', 'w') # make a file named as POSCAR

fout.write(

'''G_monolayer

1.0

%.10f 0.0000000000 0.0000000000

%.10f %.10f 0.0000000000

0.0000000000 0.0000000000 3.0000000000

C

2

Cartesian

0.000000000 0.000000000 2.171259403

1.233863017 0.712369702 2.171259403

'''%(a, -a/2, math.sqrt(3)*a/2))

# construct a primitive cell with constant c direction but variable a b directions

fout.close()

这里我们是为了改变一个graphene的primitive cell的lattice

  • 0
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值