24小时热门版块排行榜    

查看: 1442  |  回复: 13

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
回复此楼
已阅   回复此楼   关注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),这个值与实验值差距很 ...

能不能告诉我,结合能是如何计算呢?从哪里看出来结合能的大小?
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 并不是你要计算体系的总能量?
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,不知道你优化了没有?
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?请指教一二!
弘德明志博学笃行
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结构分子,一直到你想要的键长为止
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熟悉么,我遇到类似的结合能计算的问题,想请教一下你!
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对别的物理量有可能就不是好的。
弘德明志博学笃行
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对别的物理量有可能就不是好的。

对的,我没有同时优化两个参数。我是优化完一个参数后, 就去试试别的参数,不会太离谱就用了
10楼2011-03-07 11:00:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 leila3409 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 关于2023的收入 +33 小龙虾2008 2024-06-10 34/1700 2024-06-15 23:01 by zeolitess
[教师之家] 每次骚扰女学生的都是院系领导,而不是普通教师,小编们要注意措辞正确 +7 zju2000 2024-06-15 9/450 2024-06-15 22:54 by zeolitess
[考博] 上海交大招收材料化学方向科研助理/“申请考核”博士(请勿回复帖子或站内投条) +3 灵梦and紫 2024-06-12 4/200 2024-06-15 20:58 by 1822836277
[教师之家] 请问事业编制和年薪制冲突吗? +6 ZHONGWU_U 2024-06-14 6/300 2024-06-15 20:16 by Ermito
[基金申请] 博后基金,博管会会提前知道消息吗? +3 yuyiang 2024-06-13 3/150 2024-06-15 19:52 by Lynn212
[基金申请] 为什么我的博后基金还在流动站审核中?不会是学院给我卡了吧? +14 王凯12 2024-06-13 26/1300 2024-06-15 15:22 by 好人与坏人
[考博] 2025考博 +4 自强不息a?a 2024-06-15 6/300 2024-06-15 12:36 by 自强不息a?a
[基金申请] 有没有机械的前辈分享一下评上海优都是什么成果啊 +7 wulala800 2024-06-10 7/350 2024-06-15 09:33 by 晓目崇
[基金申请] 博后基金,以往的结果点不开,怎么回事呢?最后一次机会了,两次都没中前面。 +7 kyukitu 2024-06-14 13/650 2024-06-15 06:46 by 我是王小帅
[基金申请] 博士后基金需要结题吗? +8 zhouchuck 2024-06-13 8/400 2024-06-14 17:27 by liuyupu132
[论文投稿] ACS Nano投稿后分配到副编辑手里12天了,能确定送审了吗? +5 潇洒怡惜 2024-06-12 10/500 2024-06-14 09:56 by 潇洒怡惜
[论文投稿] 最近写了一篇控制优化领域的文章,可以投哪里啊?有没有水一些的期刊推荐 +7 香瓜木香 2024-06-12 13/650 2024-06-14 07:05 by 香瓜木香
[论文投稿] 文章proof要求使用机构的邮箱 5+3 不可不信缘 2024-06-11 11/550 2024-06-14 07:00 by 3001160025
[有机交流] 旋蒸请教 +6 办公室摸鱼 2024-06-09 9/450 2024-06-13 21:31 by wuxiaominggo
[考博] 博导选择 +3 bing85977 2024-06-12 3/150 2024-06-13 15:34 by 我是邱尧
[基金申请] C口食品 +6 好好不过分 2024-06-09 19/950 2024-06-13 12:26 by zydemon525
[论文投稿] with editor日期变更 +3 慎独的小花卷 2024-06-12 8/400 2024-06-13 11:00 by 慎独的小花卷
[找工作] 成都产品质量检测研究院 200+3 鲸鱼663 2024-06-11 8/400 2024-06-13 08:10 by 加纳居士
[硕博家园] 申博 +3 悦悦小小鱼 2024-06-12 3/150 2024-06-12 15:11 by chen5805
[教师之家] 公办双非,学生论文升学就业都不行。一本大学的论文升学就业没问题吧? +5 河西夜郎 2024-06-09 5/250 2024-06-10 17:59 by yyallen2003
信息提示
请填处理意见