| 查看: 2047 | 回复: 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问题 |
» 猜你喜欢
306求调剂
已经有5人回复
-
284求调剂
已经有4人回复
-
寻找调剂
已经有4人回复
-
一志愿中南化学（0703）总分337求调剂
已经有10人回复
-
324求调剂
已经有6人回复
-
一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。
已经有3人回复
-
289求调剂
已经有7人回复
-
319求调剂
已经有4人回复
-
298求调剂一志愿211
已经有3人回复
-
求调剂一志愿海大，0703化学学硕304分，有大创项目，四级已过
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】FORTRAN运行错误，请大家帮助
已经有1人回复
- 【求助】请教高人：过热蒸气在1atm下不发生分解的最高温度。￥送50金￥
已经有10人回复
houfeifan
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 959.4
- 散金: 23
- 帖子: 144
- 在线: 198.8小时
- 虫号: 916103
- 注册: 2009-11-29
- 性别: GG
- 专业: 集成电路设计与测试
8楼2011-03-06 18:10:06
wangsong1016
铜虫 (初入文坛)
- 1ST强帖: 2
- 应助: 0 (幼儿园)
- 金币: 169
- 散金: 4
- 帖子: 45
- 在线: 26.3小时
- 虫号: 408613
- 注册: 2007-06-21
- 性别: GG
- 专业: 物理
2楼2011-02-28 21:48:53
3楼2011-03-01 09:25:54
wangsong1016
铜虫 (初入文坛)
- 1ST强帖: 2
- 应助: 0 (幼儿园)
- 金币: 169
- 散金: 4
- 帖子: 45
- 在线: 26.3小时
- 虫号: 408613
- 注册: 2007-06-21
- 性别: GG
- 专业: 物理
4楼2011-03-01 16:44:13
