| 查看: 1927 | 回复: 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问题 |
» 猜你喜欢
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有7人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有7人回复
-
磺酰氟产物，毕不了业了！
已经有5人回复
-
论文终于录用啦！满足毕业条件了
已经有25人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有3人回复
-
自荐读博
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有5人回复
-
不自信的我
已经有5人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】FORTRAN运行错误，请大家帮助
已经有1人回复
- 【求助】请教高人：过热蒸气在1atm下不发生分解的最高温度。￥送50金￥
已经有10人回复
3楼2011-03-01 09:25:54
13楼2011-03-07 12:37:43
14楼2011-03-07 12:40:12