版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話(huà)題的所有回帖

leila3409

捐助貴賓 (正式寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 616.2
散金: 406
紅花: 7
帖子: 310
在線(xiàn): 71.5小時(shí)
蟲(chóng)號(hào): 985885
注冊(cè): 2010-03-30
性別: MM
專(zhuān)業(yè): 爆炸與沖擊動(dòng)力學(xué)

[交流] 【求助】H2 Binding Energy@SIESTA 已有6人參與

因?yàn)槲倚枰?jì)算H在Fe中的性質(zhì)，所以首先想確保H勢(shì)函數(shù)的精確性，我算了一下孤立空間中的H2分子的性質(zhì)，鍵長(zhǎng)為0.745 Ang，還比較正常。（實(shí)驗(yàn)值0.741）可是結(jié)合能是7.17 eV（per molecular），這個(gè)值與實(shí)驗(yàn)值差距很大（實(shí)驗(yàn)值4.75eV）。
我換了勢(shì)函數(shù)（有LDA），計(jì)算結(jié)果與GGA的還是比較類(lèi)似，我又換了基函數(shù)，結(jié)果還是不理想，所以想請(qǐng)問(wèn)有沒(méi)有高手能指點(diǎn)一下，可能是哪里的錯(cuò)誤？我的輸入文件已經(jīng)帖在了下面。

另外，我還用了VASP做測(cè)試，不考慮自旋時(shí)，距離為0.75 Ang，結(jié)合能4.55 eV。與實(shí)驗(yàn)值非常穩(wěn)合，可是當(dāng)考慮了自旋后，H－H距離變?yōu)?.48 Ang，結(jié)合能小到0.019　eV。不知道有沒(méi)有高手遇到過(guò)這種情況，謝謝！

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

回復(fù)此樓

» 收錄本帖的淘帖專(zhuān)輯推薦

siesta問(wèn)題

» 猜你喜歡

環(huán)境工程 085701，267求調(diào)劑已經(jīng)有9人回復(fù)
085601一志愿西北工業(yè)大學(xué)初試346 已經(jīng)有4人回復(fù)
材料與化工（0856）304求B區(qū)調(diào)劑已經(jīng)有10人回復(fù)
本2一志愿C9-333分，材料科學(xué)與工程，求調(diào)劑已經(jīng)有4人回復(fù)
一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑已經(jīng)有6人回復(fù)
070300化學(xué)279求調(diào)劑已經(jīng)有7人回復(fù)
085601一志愿中山大學(xué)深圳材料工程330求調(diào)劑已經(jīng)有5人回復(fù)
293分求調(diào)劑，外語(yǔ)為俄語(yǔ) 已經(jīng)有4人回復(fù)
0856調(diào)劑已經(jīng)有3人回復(fù)
一志愿大連理工大學(xué)材料求調(diào)劑已經(jīng)有6人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【求助】FORTRAN運(yùn)行錯(cuò)誤，請(qǐng)大家?guī)椭?/a> 已經(jīng)有1人回復(fù)

【求助】請(qǐng)教高人：過(guò)熱蒸氣在1atm下不發(fā)生分解的最高溫度。￥送50金￥已經(jīng)有10人回復(fù)

1樓 2011-02-28 11:11:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

leila3409

捐助貴賓 (正式寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 616.2
散金: 406
紅花: 7
帖子: 310
在線(xiàn): 71.5小時(shí)
蟲(chóng)號(hào): 985885
注冊(cè): 2010-03-30
性別: MM
專(zhuān)業(yè): 爆炸與沖擊動(dòng)力學(xué)

引用回帖:

Originally posted by wangsong1016 at 2011-02-28 21:48:53:
能不能告訴我，結(jié)合能是如何計(jì)算呢？從哪里看出來(lái)結(jié)合能的大�。�

E(H2)-2E(isolate H atom)

贊一下

回復(fù)此樓

3樓2011-03-01 09:25:54

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 14 個(gè)回答

wangsong1016

銅蟲(chóng) (初入文壇)

1ST強(qiáng)帖: 2
應(yīng)助: 0 (幼兒園)
金幣: 169
散金: 4
帖子: 45
在線(xiàn): 26.3小時(shí)
蟲(chóng)號(hào): 408613
注冊(cè): 2007-06-21
性別: GG
專(zhuān)業(yè): 物理

★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流

引用回帖:

Originally posted by leila3409 at 2011-02-28 11:11:01:
因?yàn)槲倚枰?jì)算H在Fe中的性質(zhì)，所以首先想確保H勢(shì)函數(shù)的精確性，我算了一下孤立空間中的H2分子的性質(zhì)，鍵長(zhǎng)為0.745 Ang，還比較正常。（實(shí)驗(yàn)值0.741）可是結(jié)合能是7.17 eV（per molecular），這個(gè)值與實(shí)驗(yàn)值差距很 ...

能不能告訴我，結(jié)合能是如何計(jì)算呢？從哪里看出來(lái)結(jié)合能的大��？

贊一下(1人)

回復(fù)此樓

2樓2011-02-28 21:48:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wangsong1016

銅蟲(chóng) (初入文壇)

1ST強(qiáng)帖: 2
應(yīng)助: 0 (幼兒園)
金幣: 169
散金: 4
帖子: 45
在線(xiàn): 26.3小時(shí)
蟲(chóng)號(hào): 408613
注冊(cè): 2007-06-21
性別: GG
專(zhuān)業(yè): 物理

★ ★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流
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)

我用你的文件計(jì)算了一下，得到的輸出文件中Total =-31.341350eV，氫原子的基態(tài)能級(jí)是-13.6eV，這樣相減，得到的能量差是4.14eV左右。和你說(shuō)的7.7eV差的比較多。是不是在輸出中選擇了不同的能量，得到了你的結(jié)果，而我選Total 并不是你要計(jì)算體系的總能量？

贊一下(2人)

回復(fù)此樓

4樓2011-03-01 16:44:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

minmin_0082003

金蟲(chóng) (正式寫(xiě)手)

莫愁莫愁

1ST強(qiáng)帖: 1
應(yīng)助: 2 (幼兒園)
貴賓: 0.08
金幣: 921.3
散金: 149
帖子: 910
在線(xiàn): 218.4小時(shí)
蟲(chóng)號(hào): 366989
注冊(cè): 2007-05-10
性別: MM
專(zhuān)業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

★ ★ ★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流
zzy870720z(金幣+2): 謝謝指點(diǎn) 2011-03-02 16:17:07

考慮自旋的情況下，認(rèn)為你要計(jì)算兩個(gè)H原子，這個(gè)我也試過(guò)的。
還有一種可能就是因?yàn)槟愕腂asis，不知道你優(yōu)化了沒(méi)有？

贊一下(2人)

回復(fù)此樓

5樓2011-03-02 09:23:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 14 個(gè)回答

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 材料與化工調(diào)劑一志愿大連海事085600，349 +6	吃的不少 2026-03-30	6/300	2026-03-31 04:17 by fmesaito
[考研] 材料科學(xué)與工程調(diào)劑 +10	深V宿舍吧 2026-03-30	11/550	2026-03-31 04:12 by fmesaito
[考研] 一志愿085600中科院寧波所276分求調(diào)劑 +4	材料學(xué)257求調(diào)劑 2026-03-28	5/250	2026-03-31 03:57 by fmesaito
[考研] 哈爾濱工業(yè)大學(xué)材料與化工專(zhuān)碩378求調(diào)劑 +3	塔比烏斯 2026-03-30	3/150	2026-03-30 22:55 by 無(wú)際的草原
[考研] 08工科求調(diào)劑286 +4	tgs_001 2026-03-28	4/200	2026-03-30 21:03 by dophin1985
[考研] 285求調(diào)劑 +6	AZMK 2026-03-29	9/450	2026-03-30 21:02 by dophin1985
[考研] 346求調(diào)劑一志愿070303有機(jī)化學(xué) +7	蘿卜燉青菜 2026-03-28	8/400	2026-03-30 19:59 by 無(wú)際的草原
[考研] 287求調(diào)劑 +14	land xuxu 2026-03-26	14/700	2026-03-30 18:38 by 544594351
[考研] 0856求調(diào)劑 +8	楒桉 2026-03-28	8/400	2026-03-30 10:00 by wzy-lxz
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +8	柳溪邊 2026-03-26	8/400	2026-03-29 20:42 by 唐沐兒
[考研] 321求調(diào)劑 +7	璞玉~~ 2026-03-25	8/400	2026-03-29 06:41 by 544594351
[考研] 343求調(diào)劑 +5	愛(ài)羈絆 2026-03-28	5/250	2026-03-28 20:53 by 唐沐兒
[考研] 283求調(diào)劑 +3	A child 2026-03-28	3/150	2026-03-28 15:41 by ms629
[考研] 340求調(diào)劑 +5	jhx777 2026-03-27	5/250	2026-03-28 04:18 by fmesaito
[考研] 352分化工與材料 +5	海納百川Ly 2026-03-27	5/250	2026-03-28 03:39 by fmesaito
[考研] 一志愿鄭大085600，310分求調(diào)劑 +5	李瀟可 2026-03-26	5/250	2026-03-27 11:14 by 不吃魚(yú)的貓
[論文投稿] Journal of Mechanical Science and Technology +3	Russ_ss 2026-03-25	5/250	2026-03-27 10:49 by 陸小果畫(huà)大餅
[考研] 求調(diào)劑一志愿本科北科大化學(xué) 343 +6	13831862839 2026-03-24	7/350	2026-03-26 22:57 by 不吃魚(yú)的貓
[考研] 334分一志愿武理-080500 材料求調(diào)劑 +4	李李不服輸 2026-03-25	4/200	2026-03-25 21:26 by 星空星月
[考研] 考研一志愿蘇州大學(xué)初始315（英一）求調(diào)劑 +3	sbdksD 2026-03-24	4/200	2026-03-25 18:16 by xcjcqu