| 查看: 1862 | 回复: 13 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[交流]
【求助】H2 Binding Energy@SIESTA已有6人参与
|
||||
|
因为我需要计算H在Fe中的性质，所以首先想确保H势函数的精确性，我算了一下孤立空间中的H2分子的性质，键长为0.745 Ang，还比较正常。（实验值0.741）可是结合能是7.17 eV（per molecular），这个值与实验值差距很大（实验值4.75eV）。 我换了势函数（有LDA），计算结果与GGA的还是比较类似，我又换了基函数，结果还是不理想，所以想请问有没有高手能指点一下，可能是哪里的错误？我的输入文件已经帖在了下面。 另外，我还用了VASP做测试，不考虑自旋时，距离为0.75 Ang，结合能4.55 eV。与实验值非常稳合，可是当考虑了自旋后，H－H距离变为3.48 Ang，结合能小到0.019 eV。不知道有没有高手遇到过这种情况，谢谢！ SIESTA input ----------------------------------------------------- SystemName H2 SystemLabel H NumberOfSpecies 1 %block ChemicalSpeciesLabel 1 1 H %endblock ChemicalSpeciesLabel # Basis #PAO.basis-size DZP PAO.energy-shift 0.015 Ry #300 meV %block PAO.Basis H 2 0.22 n=1 0 2 E 2.07 0.00 4.971 1.771 1.000 1.000 n=2 1 1 E 0.89 0.01 4.988 1.000 %endblock PAO.Basis Diag.DivideAndConquer F UseSaveData T MD.TypeOfRun cg MD.NumCGsteps 100 MaxSCFIterations 400 SolutionMethod diagon xc.functional GGA # Exchange-correlation functional xc.authors PBE # Exchange-correlation version WriteCoorXmol .true. WriteForces .true. DM.MixingWeight 0.04 # New DM amount for next SCF cycle DM.NumberPulay 7 %block kgrid_Monkhorst_Pack 5 0 0 0.5 0 5 0 0.5 0 0 5 0.5 %endblock kgrid_Monkhorst_Pack NumberOfAtoms 2 LatticeConstant 10. Ang %block LatticeVectors 1.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 1.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 1.0000000000000000 %endblock LatticeVectors # Atomic coordinates AtomicCoordinatesFormat Fractional %block AtomicCoordinatesAndAtomicSpecies 0. 0. 0. 1 0. 0. 0.075 1 %endblock AtomicCoordinatesAndAtomicSpecies |
回复此楼» 收录本帖的淘帖专辑推荐
 siesta问题 |
» 猜你喜欢
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有9人回复
-
心脉受损
已经有6人回复
-
博士申请都是内定的吗？
已经有8人回复
-
博士读完未来一定会好吗
已经有32人回复
-
读博
已经有5人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
Springer期刊投稿求助
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
【求助】FORTRAN运行错误，请大家帮助
已经有1人回复- 【求助】请教高人：过热蒸气在1atm下不发生分解的最高温度。￥送50金￥
已经有10人回复
minmin_0082003
金虫 (正式写手)
莫愁莫愁
- 1ST强帖: 1
- 应助: 2 (幼儿园)
- 贵宾: 0.08
- 金币: 921.3
- 散金: 149
- 帖子: 910
- 在线: 218.4小时
- 虫号: 366989
- 注册: 2007-05-10
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
5楼2011-03-02 09:23:37
minmin_0082003
金虫 (正式写手)
莫愁莫愁
- 1ST强帖: 1
- 应助: 2 (幼儿园)
- 贵宾: 0.08
- 金币: 921.3
- 散金: 149
- 帖子: 910
- 在线: 218.4小时
- 虫号: 366989
- 注册: 2007-05-10
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
7楼2011-03-06 16:53:22
minmin_0082003
金虫 (正式写手)
莫愁莫愁
- 1ST强帖: 1
- 应助: 2 (幼儿园)
- 贵宾: 0.08
- 金币: 921.3
- 散金: 149
- 帖子: 910
- 在线: 218.4小时
- 虫号: 366989
- 注册: 2007-05-10
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
10楼2011-03-07 11:00:14
minmin_0082003
金虫 (正式写手)
莫愁莫愁
- 1ST强帖: 1
- 应助: 2 (幼儿园)
- 贵宾: 0.08
- 金币: 921.3
- 散金: 149
- 帖子: 910
- 在线: 218.4小时
- 虫号: 366989
- 注册: 2007-05-10
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
11楼2011-03-07 11:01:06