应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 【求助】H2 Binding Energy@SIESTA
51/1返回列表
查看: 2406  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

leila3409

捐助贵宾 (正式写手)

[交流] 【求助】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问题

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-02-28 11:11:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangsong1016

铜虫 (初入文坛)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
mazuju028(金币+1): 谢谢交流 2011-03-07 08:07:09
引用回帖:
Originally posted by leila3409 at 2011-03-01 09:25:54:
E(H2)-2E(isolate H atom)

我用你的文件计算了一下，得到的输出文件中Total =-31.341350eV，氢原子的基态能级是-13.6eV，这样相减，得到的能量差是4.14eV左右。和你说的7.7eV差的比较多。是不是在输出中选择了不同的能量，得到了你的结果，而我选Total 并不是你要计算体系的总能量？
赞一下(2人) 回复此楼
4楼2011-03-01 16:44:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

wangsong1016

铜虫 (初入文坛)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by leila3409 at 2011-02-28 11:11:01:
因为我需要计算H在Fe中的性质，所以首先想确保H势函数的精确性，我算了一下孤立空间中的H2分子的性质，键长为0.745 Ang，还比较正常。（实验值0.741）可是结合能是7.17 eV（per molecular），这个值与实验值差距很 ...

能不能告诉我，结合能是如何计算呢？从哪里看出来结合能的大小？
赞一下(1人) 回复此楼
2楼2011-02-28 21:48:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leila3409

捐助贵宾 (正式写手)
引用回帖:
Originally posted by wangsong1016 at 2011-02-28 21:48:53:
能不能告诉我，结合能是如何计算呢？从哪里看出来结合能的大小？

E(H2)-2E(isolate H atom)
赞一下 回复此楼
3楼2011-03-01 09:25:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2): 谢谢指点 2011-03-02 16:17:07
考虑自旋的情况下，认为你要计算两个H原子，这个我也试过的。
还有一种可能就是因为你的Basis，不知道你优化了没有？
赞一下(2人) 回复此楼
5楼2011-03-02 09:23:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +31 医学老男孩 2026-08-13 69/3450 2026-08-16 00:16 by mo_mo_mo_mo
[基金申请] 时间戳又变了8-15 +12 archvillain 2026-08-15 23/1150 2026-08-15 20:47 by foyo404
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
信息提示
关闭
请填处理意见
关闭