本文整理汇总了Python中vasp.Vasp.vasprc方法的典型用法代码示例。如果您正苦于以下问题:Python Vasp.vasprc方法的具体用法?Python Vasp.vasprc怎么用?Python Vasp.vasprc使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类vasp.Vasp
的用法示例。
在下文中一共展示了Vasp.vasprc方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: Atoms
# 需要导入模块: from vasp import Vasp [as 别名]
# 或者: from vasp.Vasp import vasprc [as 别名]
from ase import Atoms, Atom
from vasp import Vasp
Vasp.vasprc(mode=None)
#Vasp.log.setLevel(10)
import matplotlib.pyplot as plt
import numpy as np
from ase.dft import DOS
import pylab as plt
a = 3.9 # approximate lattice constant
b = a / 2.
bulk = Atoms([Atom('Pd', (0.0, 0.0, 0.0))],
cell=[(0, b, b),
(b, 0, b),
(b, b, 0)])
kpts = [8, 10, 12, 14, 16, 18, 20]
calcs = [Vasp('bulk/pd-dos-k{0}-ismear-5'.format(k),
encut=300,
xc='PBE',
kpts=[k, k, k],
atoms=bulk) for k in kpts]
Vasp.wait(abort=True)
for calc in calcs:
# this runs the calculation
if calc.potential_energy is not None:
dos = DOS(calc, width=0.2)
d = dos.get_dos() + k / 4.0
e = dos.get_energies()
plt.plot(e, d, label='k={0}'.format(k))
else:
pass
plt.xlabel('energy (eV)')