电脑桌面
添加小米粒文库到电脑桌面
安装后可以在桌面快捷访问

VASP几个计算实例VIP免费

VASP几个计算实例_第1页
1/18
VASP几个计算实例_第2页
2/18
VASP几个计算实例_第3页
3/18
用VASP 计算H原子的能量 氢原子的能量为 。在这一节中,我们用VASP 计算H 原子的能量。对于原子计算,我们可以采用如下的INCAR 文件 PREC=ACCURATE NELMDL = 5 make five delays till charge mixing ISMEAR = 0; SIGMA=0.05 use smearing method 采用如下的KPOINTS文件。由于增加K 点的数目只能改进描述原子间的相互作用,而在单原子计算中并不需要。所以我们只需要一个K 点。 Monkhorst Pack 0 Monkhorst Pack 1 1 1 0 0 0 采用如下的POSCAR 文件 atom 1 15.00000 .00000 .00000 .00000 15.00000 .00000 .00000 .00000 15.00000 1 cart 0 0 0 采用标准的H 的POTCAR 得到结果如下: k-point 1 : 0.0000 0.0000 0.0000 band No. band energies occupation 1 -6.3145 1.00000 2 -0.0527 0.00000 3 0.4829 0.00000 4 0.4829 0.00000 我们可以看到,电子的能级不为 。 Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 0.00060791 Ewald energy TEWEN = -1.36188267 -1/2 Hartree DENC = -6.27429270 -V(xc)+E(xc) XCENC = 1.90099128 PAW double counting = 0.00000000 0.00000000 entropy T*S EENTRO = -0.02820948 eigenvalues EBANDS = -6.31447362 atomic energy EATOM = 12.04670449 --------------------------------------------------- free energy TOTEN = -0.03055478 eV energy without entropy = -0.00234530 energy(sigma->0) = -0.01645004 我们可以看到 也不等于 。 在上面的计算中有个问题,就是 H 原子有 spin,而在上面的计算中我们并没有考虑到spin。所以如果我们改用 LSDA 近似,在 INCAR 中用 ISPIN=2 的 tag,则得到如下结果: k-point 1 : 0.0000 0.0000 0.0000 band No. band energies occupation 1 -7.2736 1.00000 2 -0.1229 0.00000 3 0.4562 0.00000 4 0.4562 0.00000 5 0.4562 0.00000 spin component 2 k-point 1 : 0.0000 0.0000 0.0000 band No. band energies occupation 1 -2.4140 0.00000 2 -0.0701 0.00000 3 0.5179 0.00000 4 0.5179 0.00000 5 0.5179 0.00000 Free energy of the ion-electron system ...

1、当您付费下载文档后,您只拥有了使用权限,并不意味着购买了版权,文档只能用于自身使用,不得用于其他商业用途(如 [转卖]进行直接盈利或[编辑后售卖]进行间接盈利)。
2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。
3、如文档内容存在违规,或者侵犯商业秘密、侵犯著作权等,请点击“违规举报”。

碎片内容

VASP几个计算实例

确认删除?
VIP
微信客服
  • 扫码咨询
会员Q群
  • 会员专属群点击这里加入QQ群
客服邮箱
回到顶部