| 查看: 3561 | 回复: 5 | ||
dwd0826铁杆木虫 (著名写手)
|
[求助]
VASP计算H原子的能量的一点困惑
|
|
刚在网上看到一篇帖子，用VASP计算H原子的能量，氢原子的能量参考值为-13.6eV INCAR文件如下： PREC=ACCURATE ISPIN=2 NELMDL = 5 ISMEAR = 0; SIGMA=0.05 采用如下的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，并采用GGA的赝势，得到的结果： Free energy of the ion-electron system (eV) --------------------------------------------------- alpha Z PSCENC = 0.00621465 Ewald energy TEWEN = -1.38027565 -1/2 Hartree DENC = -6.91107031 -V(xc)+E(xc) XCENC = 2.32601856 PAW double counting = 0.70286470 -0.71256934 entropy T*S EENTRO = 0.00434843 eigenvalues EBANDS = -7.71257941 atomic energy EATOM = 12.52153358 --------------------------------------------------- free energy TOTEN = -1.15551480 eV energy without entropy = -1.15986323 energy(sigma->0) = -1.15696428 通过TOTEN-EATOM计算单个原子的能量，氢原子的能量约等于-13.67eV =================================================== 现有点困惑：结合能=晶体中每个原子的能量-孤立的单个原子的能量 这里的“孤立原子的能量”指得是上面TOTEN-EATOM所计算的单个原子的能量？还是OUTCAR中"energy without entropy=”的取值？ 我对上面TOTEN-EATOM概念比较模糊，有谁能帮忙具体解释下其中的意义吗？ |
回复此楼» 本帖已获得的红花（最新10朵）
0Reliance» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有59人回复
-
filecode=后面第一个是大写字母
已经有5人回复
-
时间戳又变了8-15
已经有21人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有10人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- VASP算同一物质超晶胞Band与原胞的差别很大
已经有12人回复
- vasp小白一枚，求大神指教，单原子能量计算问题，内附计算文件
已经有5人回复
- 单层原子态密度VASP计算设置问题
已经有9人回复
- vasp做表面吸附如何确定要固定的原子坐标
已经有3人回复
- vasp 计算单分子/原子 时晶胞大小的选取
已经有6人回复
- vasp计算含有磷原子结构的态密度时磷原子为什么还有d态
已经有14人回复
- 如何用VASP计算H原子的振动频率
已经有16人回复
- vasp计算界面结构时原子弛豫的问题
已经有12人回复
- vasp计算单个氧原子和分子的能量
已经有20人回复
- 在vasp计算里，有没有好的只优化原子位置和晶胞体积的方法？
已经有8人回复
- vasp进行固定某些原子的某些坐标的结构优化
已经有5人回复
- vasp可以固定单个原子的磁矩吗？
已经有5人回复
- 如何用VASP计算体系中每一类原子的DOS？
已经有10人回复
- vasp如何优化指定的几个原子
已经有9人回复
- vasp计算磁矩的困惑
已经有10人回复
- 【求助】关于锂原子vasp计算后的1s2s2p分波态密度处理
已经有12人回复
- 【求助】计算单个Fe原子，无法收敛，请教
已经有19人回复
uuv2010
荣誉版主 (职业作家)
- 1ST强帖: 6
- 应助: 160 (高中生)
- 贵宾: 3.14
- 金币: 3677.7
- 散金: 5807
- 红花: 79
- 沙发: 35
- 帖子: 4564
- 在线: 916.6小时
- 虫号: 966216
- 注册: 2010-03-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 第一性原理
2楼2011-12-07 17:29:15
stractor
金虫 (著名写手)
- 1ST强帖: 6
- 应助: 91 (初中生)
- 金币: 1487.6
- 散金: 666
- 红花: 30
- 沙发: 1
- 帖子: 2224
- 在线: 426.8小时
- 虫号: 394087
- 注册: 2007-06-06
- 专业: 半导体材料
3楼2011-12-07 17:56:13
4楼2013-10-29 17:38:56
5楼2013-10-29 17:41:49
0Reliance
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2474.8
- 散金: 100
- 红花: 2
- 帖子: 378
- 在线: 84.5小时
- 虫号: 4113441
- 注册: 2015-09-30
- 性别: MM
- 专业: 原子和分子物理
6楼2016-03-24 09:18:48
