本文整理汇总了Python中vasp.Vasp.get_potential_energy方法的典型用法代码示例。如果您正苦于以下问题:Python Vasp.get_potential_energy方法的具体用法?Python Vasp.get_potential_energy怎么用?Python Vasp.get_potential_energy使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类vasp.Vasp
的用法示例。
在下文中一共展示了Vasp.get_potential_energy方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: Vasp
# 需要导入模块: from vasp import Vasp [as 别名]
# 或者: from vasp.Vasp import get_potential_energy [as 别名]
from vasp import Vasp
from ase import Atom, Atoms
# parent metals
atoms = Vasp('bulk/alloy/cu').get_atoms()
cu = atoms.get_potential_energy() / len(atoms)
atoms = Vasp('bulk/alloy/pd').get_atoms()
pd = atoms.get_potential_energy() / len(atoms)
atoms = Atoms([Atom('Cu', [-3.672, 3.672, 3.672]),
Atom('Cu', [0.000, 0.000, 0.000]),
Atom('Cu', [-10.821, 10.821, 10.821]),
Atom('Pd', [-7.246, 7.246, 7.246])],
cell=[[-5.464, 3.565, 5.464],
[-3.565, 5.464, 5.464],
[-5.464, 5.464, 3.565]])
calc = Vasp('bulk/alloy/cu3pd-1',
xc='PBE',
encut=350,
kpts=[8, 8, 8],
nbands=34,
ibrion=2,
isif=3,
nsw=10,
atoms=atoms)
e3 = atoms.get_potential_energy()
Vasp.wait(abort=True)
for atom in atoms:
if atom.symbol == 'Cu':
e3 -= cu
else:
e3 -= pd
e3 /= len(atoms)