版塊導(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í)熱門版塊排行榜

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

返回列表

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

leila3409

捐助貴賓 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 616.2
散金: 406
紅花: 7
帖子: 310
在線: 71.5小時(shí)
蟲(chóng)號(hào): 985885
注冊(cè): 2010-03-30
性別: MM
專業(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的還是比較類似，我又換了基函數(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ù)此樓

» 收錄本帖的淘帖專輯推薦

siesta問(wèn)題

» 猜你喜歡

299求調(diào)劑已經(jīng)有8人回復(fù)
一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑已經(jīng)有6人回復(fù)
300求調(diào)劑，材料科學(xué)英一數(shù)二已經(jīng)有8人回復(fù)
招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑已經(jīng)有5人回復(fù)
070305高分子化學(xué)與物理 304分求調(diào)劑已經(jīng)有7人回復(fù)
289求調(diào)劑已經(jīng)有13人回復(fù)
一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分已經(jīng)有9人回復(fù)
081200-11408-276學(xué)碩求調(diào)劑已經(jīng)有6人回復(fù)
調(diào)劑求院校招收已經(jīng)有5人回復(fù)
調(diào)劑310 已經(jīng)有8人回復(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的回帖

zhangguangping

木蟲(chóng) (著名寫手)

1ST強(qiáng)帖: 23
應(yīng)助: 71 (初中生)
貴賓: 0.031
金幣: 1973.1
散金: 9059
紅花: 79
帖子: 2681
在線: 2128.8小時(shí)
蟲(chóng)號(hào): 529624
注冊(cè): 2008-03-20
性別: GG
專業(yè): 原子和分子物理

★ ★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流
zzy870720z(金幣+1): 謝謝交流 2011-03-06 22:36:17

引用回帖:

Originally posted by minmin_0082003 at 2011-03-06 09:53:22:
不好意思，生病了，沒(méi)來(lái)實(shí)驗(yàn)室。那個(gè)basis的優(yōu)化，我優(yōu)化過(guò)一次。
具體細(xì)節(jié)記不太清楚了，但是是在siesta/Util/Optimizer 可以參考里面的Example。
你可以根據(jù)你需要的不同來(lái)選擇優(yōu)化的參數(shù)。比如你希望你的b ...

你的意思是說(shuō)根據(jù)你要的物理量來(lái)找basis了。這樣話，對(duì)于鍵長(zhǎng)是好的basis對(duì)別的物理量有可能就不是好的。

贊一下(2人)

回復(fù)此樓

弘德明志博學(xué)篤行

9樓2011-03-06 18:26:39

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

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

wangsong1016

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

1ST強(qiáng)帖: 2
應(yīng)助: 0 (幼兒園)
金幣: 169
散金: 4
帖子: 45
在線: 26.3小時(shí)
蟲(chóng)號(hào): 408613
注冊(cè): 2007-06-21
性別: GG
專業(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的回帖

leila3409

捐助貴賓 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 616.2
散金: 406
紅花: 7
帖子: 310
在線: 71.5小時(shí)
蟲(chóng)號(hào): 985885
注冊(cè): 2010-03-30
性別: MM
專業(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的回帖

wangsong1016

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

1ST強(qiáng)帖: 2
應(yīng)助: 0 (幼兒園)
金幣: 169
散金: 4
帖子: 45
在線: 26.3小時(shí)
蟲(chóng)號(hào): 408613
注冊(cè): 2007-06-21
性別: GG
專業(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的回帖

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

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

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

[考研] 11408軟件工程求調(diào)劑 +3	Qiu學(xué)ing 2026-03-28	3/150	2026-03-28 21:50 by zhq0425
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過(guò) +7	幸運(yùn)哩哩 2026-03-22	11/550	2026-03-28 20:27 by 唐沐兒
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +4	小張zxy張 2026-03-26	8/400	2026-03-28 19:25 by lbsjt
[考研] 一志愿廈門大學(xué)化學(xué)學(xué)碩307求調(diào)劑 +10	y7czhao 2026-03-26	10/500	2026-03-28 14:23 by 唐沐兒
[考研] 0703化學(xué)求調(diào)劑 +9	奶油草莓. 2026-03-22	10/500	2026-03-28 13:30 by 唐沐兒
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8	曼111 2026-03-24	9/450	2026-03-28 07:58 by YYYYX1234
[考研] 0703化學(xué)求調(diào)劑，各位老師看看我！�。� +5	祁祺祺 2026-03-25	5/250	2026-03-27 21:44 by 東方豬豬
[考研] 272求調(diào)劑 +7	腳滑的守法公民 2026-03-27	7/350	2026-03-27 17:23 by laoshidan
[考研] 0856調(diào)劑 +5	求求讓我有書(shū)讀?/a> 2026-03-26	6/300	2026-03-27 15:12 by caszguilin
[考研] 化學(xué)308分求調(diào)劑 +8	你好明天你好 2026-03-23	9/450	2026-03-27 14:01 by 楊光于青云
[考研] 314求調(diào)劑 +3	溪云珂 2026-03-26	3/150	2026-03-27 11:20 by sanrepian
[考研] 考研調(diào)劑 +9	小蠟新筆 2026-03-26	9/450	2026-03-27 11:10 by 不吃魚(yú)的貓
[考研] 281求調(diào)劑 +3	亞克西good 2026-03-26	5/250	2026-03-26 19:48 by 不吃魚(yú)的貓
[考研] 352求調(diào)劑 +4	大米飯！ 2026-03-22	4/200	2026-03-26 16:40 by 不吃魚(yú)的貓
[考研] 289求調(diào)劑 +17	碩星赴 2026-03-23	17/850	2026-03-26 16:18 by 不吃魚(yú)的貓
[考研] 材料與化工328分調(diào)劑 +6	。，。，。，。i 2026-03-23	6/300	2026-03-25 22:30 by 418490947
[考研] 296求調(diào)劑 +4	汪�。�！ 2026-03-25	7/350	2026-03-25 16:41 by 汪�。�！
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7	小羊睡著了? 2026-03-23	10/500	2026-03-25 02:24 by 清風(fēng)拂揚(yáng)。 m
[考研] 一志愿南航材料專317分求調(diào)劑 +5	炸呀炸呀炸薯?xiàng)l 2026-03-23	5/250	2026-03-24 16:52 by 星空星月
[考研] 求老師收我 +3	zzh16938784 2026-03-23	3/150	2026-03-23 12:56 by ztnimte