| 查看: 2150 | 回复: 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 |
回复此楼» 猜你喜欢
药学硕士，第一、第二作者已发表6 篇 SCI，药理方向及相关方向2026年/2027年博士申请
已经有4人回复
-
26年博士申请自荐-电催化
已经有3人回复
-
中国地质大学（北京）博士招生补录，数理学院材料科学与工程专业和材料与化工专业
已经有6人回复
-
收到国自然专家邀请后几年才会有本子送过来评
已经有4人回复
-
考博
已经有5人回复
-
26年申博自荐-计算机视觉
已经有4人回复
-
药化及相关博士的申请
已经有3人回复
-
一篇MDPI论文改变了学习工作和生活
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
fengshiquan
铁杆木虫 (正式写手)
- 应助: 10 (幼儿园)
- 金币: 5830.3
- 散金: 80
- 红花: 5
- 帖子: 599
- 在线: 1569.1小时
- 虫号: 884725
- 注册: 2009-10-26
- 专业: 原子和分子物理
12楼2011-03-07 11:40:30
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