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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 【求助】H2 Binding Energy@SIESTA
141/2返回列表12下一页
查看: 2026  |  回复: 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):给个红包，谢谢回帖交流
引用回帖:
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的回帖

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的回帖

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的回帖

zhangguangping

木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by minmin_0082003 at 2011-03-02 02:23:37:
考虑自旋的情况下，认为你要计算两个H原子，这个我也试过的。
还有一种可能就是因为你的Basis，不知道你优化了没有？

请问，怎么优化basis？请指教一二！
赞一下(1人) 回复此楼
弘德明志博学笃行
6楼2011-03-02 10:28:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by zhangguangping at 2011-03-02 10:28:52:
请问，怎么优化basis？请指教一二！

不好意思，生病了，没来实验室。那个basis的优化，我优化过一次。
具体细节记不太清楚了，但是是在siesta/Util/Optimizer 可以参考里面的Example。
你可以根据你需要的不同来选择优化的参数。比如你希望你的basis 可以得到你想要的O2键长，它会生成不同的basis后，去优化O2结构分子，一直到你想要的键长为止
赞一下(1人) 回复此楼
7楼2011-03-06 16:53:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

houfeifan

金虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by minmin_0082003 at 2011-03-06 16:53:22:
不好意思，生病了，没来实验室。那个basis的优化，我优化过一次。
具体细节记不太清楚了，但是是在siesta/Util/Optimizer 可以参考里面的Example。
你可以根据你需要的不同来选择优化的参数。比如你希望你的b ...

请问你对MS里面的CASTEP熟悉么，我遇到类似的结合能计算的问题，想请教一下你！
赞一下(1人) 回复此楼
8楼2011-03-06 18:10:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1): 谢谢交流 2011-03-06 22:36:17
引用回帖:
Originally posted by minmin_0082003 at 2011-03-06 09:53:22:
不好意思，生病了，没来实验室。那个basis的优化，我优化过一次。
具体细节记不太清楚了，但是是在siesta/Util/Optimizer 可以参考里面的Example。
你可以根据你需要的不同来选择优化的参数。比如你希望你的b ...

你的意思是说根据你要的物理量来找basis了。这样话，对于键长是好的basis对别的物理量有可能就不是好的。
赞一下(2人) 回复此楼
弘德明志博学笃行
9楼2011-03-06 18:26:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by zhangguangping at 2011-03-06 18:26:39:
你的意思是说根据你要的物理量来找basis了。这样话，对于键长是好的basis对别的物理量有可能就不是好的。

对的，我没有同时优化两个参数。我是优化完一个参数后， 就去试试别的参数，不会太离谱就用了
赞一下(1人) 回复此楼
10楼2011-03-07 11:00:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 leila3409 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 311求调剂 +5 冬十三 2026-03-15 5/250 2026-03-15 18:38 by 无际的草原
[考研] 0703化学调剂 290分有科研经历，论文在投 +3 腻腻gk 2026-03-14 3/150 2026-03-15 17:28 by 小物理化学
[考研] 274求调剂 +4 时间点 2026-03-13 4/200 2026-03-15 15:29 by Rambo13
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 331求调剂（0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 求调剂，药学 +3 归零lbm 2026-03-09 5/250 2026-03-14 02:21 by JourneyLucky
[考研] 一志愿郑大070303，338分，求调剂 +4 dadawaf 2026-03-10 5/250 2026-03-14 01:20 by lsw010101
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 26考研调剂 +3 ying123. 2026-03-10 3/150 2026-03-14 00:18 by JourneyLucky
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 工科调剂 +4 Jiang191123！ 2026-03-11 4/200 2026-03-13 15:15 by Miko19
[考研] 285求调剂 +4 ytter 2026-03-12 4/200 2026-03-13 14:48 by jxchenghu
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[考研] 279求调剂 +3 莫xiao 2026-03-10 4/200 2026-03-11 08:06 by 斩魂滴兔子！
[考研] 调剂 +5 呵唔哦豁 2026-03-10 5/250 2026-03-10 22:00 by 28375m
[考研] 085602化工求调剂 +7 董boxing 2026-03-10 7/350 2026-03-10 17:07 by BruceLiu320
信息提示
关闭
请填处理意见
关闭