| 查看: 3565 | 回复: 6 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】VASP 计算分子解离能问题，为什么总是出错 已有4人参与
|
|||
|
最近算了几个分子的解离能，包括F2 ，Cl2, H2分子，只有H的结果与实验值吻合较好。我的计算公式是Edis= E(molecule) -2 E(atom)。 其中，H2分子以及F，Cl，H计算时加了自旋，Cl2 和F2分子没有加自旋。结果Cl2和F2的结果比实验值大了1点几个eV, H2解离能算出来4.6左右，与实验值4.4左右还算对应不错。 请问有人遇到过类似问题吗？ 以下是从别人的帖子转过来的，我也用的类似的方法： ------------------------------------------------------------------------------------ 氢气分子的解离能，也就是结合能,根据资料中给出的是约4.48eV。(G. Kresse & J. Hafner, Surface Sci. 459 (2000) 287) 为此，首先要计算一个氢气分子的孤立能量，再减去两个孤立氢原子能量，将得到氢气分子的结合能。 先计算单个原子能量，选取PAW_PBE文件夹下的H下面的POTCAR 用到的其他输入文件如下： INCAR: SYSTEM = H atom in a box ISMEAR = 0 ! Gaussian smearing SIGMA = 0.01 ENCUT = 350.0 KPOINTS: Automatic mesh 0 Monkhorst Pack 1 1 1 0. 0. 0. POSCAR: H atom in a box 1.0 ! universal scaling parameters 7.0 0.0 0.0 ! lattice vector a(1) 0.0 8.0 0.0 ! lattice vector a(2) 0.0 0.0 9.0 ! lattice vector a(3) 1 ! number of atoms cart ! positions in cartesian coordinates 0 0 0 用上述文件计算得到 TOTEN(H) = +0.000854eV (H atom)，用PAW_GGA得到类似结果。 然后计算氢分子能量，用类似的输入文件： INCAR: SYSTEM = H2 dimer in a box ISMEAR = 0 ! Gaussian smearing NSW = 5 ! 5 ionic steps IBRION = 2 ! use the conjugate gradient algorithm ENCUT = 350.0 POTIM = 0.1 KPOINTS: Automatic mesh 0 Monkhorst Pack 1 1 1 0. 0. 0. POSCAR: H2 molecule in a box 1.0 ! universal scaling parameters 8.0 0.0 0.0 ! lattice vector a(1) 0.0 8.0 0.0 ! lattice vector a(2) 0.0 0.0 8.0 ! lattice vector a(3) 2 ! number of atoms cart ! positions in cartesian coordinates 0 0 0 ! first atom 0 0 0.5 ! second atom 根据定义，E=-(TOTEN(H2)-2*TOTEN(H)) = 6.68eV，键长 = 0.75102 A。比文献中的结合能大了不少。 出现上述结果的原因是在计算单个H原子能量的时候没有指定基态为spin polarized state。为了得到正确的解离能，计算单个H原子能量的时候需要指定ISPIN=2。加上以上的tag后，TOTEN(H atom [spin-polarized])=-1.10351 eV，用公式重新计算解离能： E_binding=-(TOTEN(H2)-2*TOTEN(H [spin-polarized ])) = 4.48eV |
回复此楼» 猜你喜欢
极限拉扯，如何抉择
已经有15人回复
-
申请科研项目纯靠实力能成功吗？
已经有8人回复
-
中科院海洋所招收海洋材料相关领域博士研究生及研究助理
已经有4人回复
-
40+中年人失业后，还能找到好工作吗？
已经有3人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp计算过渡态问题
已经有13人回复
- vasp自洽计算K点不匹配问题，急！！
已经有15人回复
- VASP计算：LORBIT设置出现问题
已经有5人回复
- vasp计算在进行静态计算时，出错。
已经有3人回复
- VASP计算态密度出现问题
已经有12人回复
- 【求助】vasp计算出错，不知道问题出在哪里
已经有5人回复
- 【求助】vasp计算过渡态结果出错
已经有17人回复
- 【求助】vasp5 光学计算结果问题
已经有4人回复
- 【求助】新手求助vasp计算氧分子结合能的问题
已经有12人回复
- 【求助】VASP计算cluster出错！
已经有5人回复
- 【转帖】用VASP计算氢分子的能量
已经有8人回复
stractor
金虫 (著名写手)
- 1ST强帖: 6
- 应助: 91 (初中生)
- 金币: 1487.6
- 散金: 666
- 红花: 30
- 沙发: 1
- 帖子: 2224
- 在线: 426.8小时
- 虫号: 394087
- 注册: 2007-06-06
- 专业: 半导体材料
7楼2010-11-02 23:27:58
stractor
金虫 (著名写手)
- 1ST强帖: 6
- 应助: 91 (初中生)
- 金币: 1487.6
- 散金: 666
- 红花: 30
- 沙发: 1
- 帖子: 2224
- 在线: 426.8小时
- 虫号: 394087
- 注册: 2007-06-06
- 专业: 半导体材料
2楼2010-11-02 16:16:53
3楼2010-11-02 20:07:58
yidiandian2703
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 77
- 散金: 15
- 帖子: 51
- 在线: 31.4小时
- 虫号: 621019
- 注册: 2008-10-09
- 专业: 凝聚态物性I:结构、力学和
4楼2010-11-02 20:29:28
200