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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 【求助】H2 Binding Energy@SIESTA
51/1返回列表
查看: 2032  |  回复: 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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 311求调剂 +13 冬十三 2026-03-15 14/700 2026-03-21 22:10 by peike
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +6 嘉年新程 2026-03-15 6/300 2026-03-21 17:07 by Dream007008
[考研] 268求调剂 +9 简单点0 2026-03-17 9/450 2026-03-21 15:37 by lature00
[考研] 一志愿华中科技大学071000，求调剂 +3 沿岸有贝壳6 2026-03-21 3/150 2026-03-21 10:35 by 暮云清寒
[考研] 306求调剂 +4 chuanzhu川烛 2026-03-18 4/200 2026-03-21 08:25 by laoshidan
[考研] 316求调剂 +6 梁茜雯 2026-03-19 6/300 2026-03-21 06:32 by Ecowxq666！
[考研] 307求调剂 +3 wyyyqx 2026-03-17 3/150 2026-03-21 03:20 by JourneyLucky
[考研] 一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。 +3 晨昏线与星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 22408 344分 求调剂 一志愿 华电计算机技术 +4 solanXXX 2026-03-20 4/200 2026-03-20 23:49 by alg094825
[考研] 287求调剂 +7 晨昏线与星海 2026-03-19 8/400 2026-03-20 22:19 by JourneyLucky
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 319求调剂 +3 小力气珂珂 2026-03-20 3/150 2026-03-20 19:47 by JourneyLucky
[考研] 085410人工智能专硕317求调剂（0854都可以） +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 【同济软件】软件（085405）考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
信息提示
关闭
请填处理意见
关闭