def load_default_vasp(vasp: Vasp, structure=None):
vasp.has_nlep = False
vasp.maxiter = 100
if not vasp.kpoints:
if structure == None:
vasp.kpoints = "\n0\nAuto\n14"
else:
kpoints_density = 4000
vasp.kpoints = pylada.gen_kpts(structure, kpoints_density)
vasp.program = '$VASP_PYLADA' # default vasp program
try:
vasp.npar = int(os.environ['PBS_NUM_NODES'])
except:
vasp.npar = int(os.environ['SLURM_JOB_NUM_NODES'])
vasp.prec = "accurate"
vasp.ediff = 1.0e-6 # total, not per atom
vasp.ediffg = -0.02
vasp.addgrid = True
vasp.ismear = 0
vasp.sigma = 0.01
vasp.lwave = True
vasp.lcharg = True
vasp.lmaxmix = 4
vasp.convergence = 1.0e-5
vasp.algo = "Normal"
vasp.lorbit = 11
vasp.nsw = 5000
return (vasp)
def set_441(vasp: Vasp, structure=None):
x = 4
y = 4
z = 1
packing = 'Gamma'
vasp.kpoints = "Gamma_Mesh\n0\n{}\n{} {} {}".format(packing, x, y, z)
return vasp
def set_333(vasp: Vasp, structure=None):
x = 3
y = 3
z = 3
packing = 'Gamma'
vasp.kpoints = "Gamma_Mesh\n0\n{}\n{} {} {}".format(packing, x, y, z)
vasp.add_keyword('auto_gamma', False)
return vasp
def set_gamma(vasp: Vasp, structure=None):
x = 1
y = 1
z = 1
packing = 'Gamma'
vasp.kpoints = "Gamma_Mesh\n0\n{}\n{} {} {}".format(packing, x, y, z)
vasp.ismear = 0
vasp.add_keyword('kpar', 1)
return vasp
def set_kpoint(vasp: Vasp, structure):
lengths = [
sum([x**2 for x in structure.cell.transpose()[i]])**(1 / 2)
for i in range(3)
] # using distance formula to get vector lengths
kpoints = [
math.ceil(min(lengths) / x * kpoint) for x in lengths
] # scaling number of kpoints for shorter vectors
packing = 'Gamma'
vasp.kpoints = "Gamma_Mesh\n0\n{}\n{} {} {}".format(
packing, kpoints[0], kpoints[1], kpoints[2])
return vasp
def set_kpoints_auto_20(vasp: Vasp, structure=None):
x = 20
packing = 'Auto'
vasp.kpoints = "Automatic\n0\n{}\n{}".format(packing, x)
kpoint_str = vasp.kpoints.split('\n')[3]
kpoints = [int(x) for x in kpoint_str.split()]
if len(kpoints) > 1:
num_kpoints = kpoints[0] * kpoints[1] * kpoints[2]
else:
density = kpoints[0]
num_kpoints = math.ceil(
density / np.linalg.norm(structure.cell[:, 0])) * math.ceil(
density / np.linalg.norm(structure.cell[:, 1])) * math.ceil(
density / np.linalg.norm(structure.cell[:, 2]))
if num_kpoints < 4:
vasp.ismear = 0
return vasp
def bulk_standard(vasp: Vasp, structure):
# Electronic
vasp.add_keyword('GGA', 'PS')
vasp.ismear = -5
vasp.prec = "Accurate"
vasp.nelm = 60
vasp.ediff = 1e-5
vasp.nelmdl = 0
# Ionic
vasp.nsw = 5000
vasp.ediffg = -0.02
# Output
vasp.lwave = True
vasp.lcharg = True
#TODO: Get these automatically
x=2; y=2 ; z=2
packing = 'Gamma'
vasp.kpoints = "Gamma_Mesh\n0\n{}\n{} {} {}".format(packing, x, y, z)
return vasp
def mnte_standard(vasp: Vasp, structure):
# Electronic
vasp.add_keyword('METAGGA', 'SCAN')
vasp.ismear = -5
vasp.prec = "Accurate"
vasp.nelm = 200
vasp.ediff = 1e-5
vasp.nelmdl = 0
vasp.add_keyword('KPAR', 2)
vasp.ispin = 2
vasp.encut=800
# Ionic
vasp.nsw = 5000
vasp.ediffg = -0.02
# Output
vasp.lwave = True
vasp.lcharg = True
#TODO: Get these automatically
x=22 ; y='' ; z=''
packing = 'Auto'
vasp.kpoints = "Automatic\n0\n{}\n{} {} {}".format(packing, x, y, z)
return vasp
def set_kpoint(vasp: Vasp, structure):
vasp.kpoints = kpoints
return vasp
def set_kpoints_1704(vasp: Vasp, structure=None):
kpts = Kpoints.automatic_density(pyl_to_pmg(structure), 1704)
vasp.kpoints = str(kpts)
return vasp