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

返回列表

leila3409

捐助貴賓 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 616.2
散金: 406
紅花: 7
帖子: 310
在線: 71.5小時(shí)
蟲號(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)問有沒有高手能指點(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。不知道有沒有高手遇到過這種情況，謝謝！

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問題

» 猜你喜歡

求調(diào)劑已經(jīng)有7人回復(fù)
275求調(diào)劑已經(jīng)有6人回復(fù)
招08考數(shù)學(xué) 已經(jīng)有7人回復(fù)
303求調(diào)劑已經(jīng)有6人回復(fù)
306求調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑已經(jīng)有7人回復(fù)
考研調(diào)劑已經(jīng)有4人回復(fù)
285求調(diào)劑已經(jīng)有6人回復(fù)
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑已經(jīng)有7人回復(fù)
318求調(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)教高人：過熱蒸氣在1atm下不發(fā)生分解的最高溫度。￥送50金￥已經(jīng)有10人回復(fù)

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

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

wangsong1016

銅蟲 (初入文壇)

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

★
小木蟲(金幣+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ì)算呢？從哪里看出來結(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í)
蟲號(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ì)算呢？從哪里看出來結(jié)合能的大小？

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

贊一下

回復(fù)此樓

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

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

wangsong1016

銅蟲 (初入文壇)

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

★ ★
小木蟲(金幣+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左右。和你說的7.7eV差的比較多。是不是在輸出中選擇了不同的能量，得到了你的結(jié)果，而我選Total 并不是你要計(jì)算體系的總能量？

贊一下(2人)

回復(fù)此樓

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

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

minmin_0082003

金蟲 (正式寫手)

莫愁莫愁

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

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

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

贊一下(2人)

回復(fù)此樓

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

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

zhangguangping

木蟲 (著名寫手)

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

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

引用回帖:

Originally posted by minmin_0082003 at 2011-03-02 02:23:37:
考慮自旋的情況下，認(rèn)為你要計(jì)算兩個(gè)H原子，這個(gè)我也試過的。
還有一種可能就是因?yàn)槟愕腂asis，不知道你優(yōu)化了沒有？

請(qǐng)問，怎么優(yōu)化basis？請(qǐng)指教一二！

贊一下(1人)

回復(fù)此樓

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

6樓2011-03-02 10:28:52

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

minmin_0082003

金蟲 (正式寫手)

莫愁莫愁

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

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

引用回帖:

Originally posted by zhangguangping at 2011-03-02 10:28:52:
請(qǐng)問，怎么優(yōu)化basis？請(qǐng)指教一二！

不好意思，生病了，沒來實(shí)驗(yàn)室。那個(gè)basis的優(yōu)化，我優(yōu)化過一次。
具體細(xì)節(jié)記不太清楚了，但是是在siesta/Util/Optimizer 可以參考里面的Example。
你可以根據(jù)你需要的不同來選擇優(yōu)化的參數(shù)。比如你希望你的basis 可以得到你想要的O2鍵長(zhǎng)，它會(huì)生成不同的basis后，去優(yōu)化O2結(jié)構(gòu)分子，一直到你想要的鍵長(zhǎng)為止

贊一下(1人)

回復(fù)此樓

7樓2011-03-06 16:53:22

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

houfeifan

金蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 959.4
散金: 23
帖子: 144
在線: 198.8小時(shí)
蟲號(hào): 916103
注冊(cè): 2009-11-29
性別: GG
專業(yè): 集成電路設(shè)計(jì)與測(cè)試

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

引用回帖:

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

請(qǐng)問你對(duì)MS里面的CASTEP熟悉么，我遇到類似的結(jié)合能計(jì)算的問題，想請(qǐng)教一下你！

贊一下(1人)

回復(fù)此樓

8樓2011-03-06 18:10:06

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

zhangguangping

木蟲 (著名寫手)

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

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

引用回帖:

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

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

贊一下(2人)

回復(fù)此樓

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

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

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

minmin_0082003

金蟲 (正式寫手)

莫愁莫愁

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

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

引用回帖:

Originally posted by zhangguangping at 2011-03-06 18:26:39:
你的意思是說根據(jù)你要的物理量來找basis了。這樣話，對(duì)于鍵長(zhǎng)是好的basis對(duì)別的物理量有可能就不是好的。

對(duì)的，我沒有同時(shí)優(yōu)化兩個(gè)參數(shù)。我是優(yōu)化完一個(gè)參數(shù)后，就去試試別的參數(shù)，不會(huì)太離譜就用了

贊一下(1人)

回復(fù)此樓

10樓2011-03-07 11:00:14

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 leila3409 的主題更新

返回列表

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

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

[考研] 一志愿070300浙大化學(xué)358分，求調(diào)劑！ +3	酥酥魚.. 2026-03-21	3/150	2026-03-22 11:31 by 楊楊楊紫
[考研] 0856材料專碩353求調(diào)劑 +4	NIFFFfff 2026-03-20	4/200	2026-03-22 09:49 by 2026paper
[考研] 317求調(diào)劑 +9	申子申申 2026-03-19	15/750	2026-03-21 17:31 by 學(xué)員8dgXkO
[考研] 265求調(diào)劑 +12	梁梁校校 2026-03-19	14/700	2026-03-21 13:38 by lature00
[考研] 306求調(diào)劑 +4	chuanzhu川燭 2026-03-18	4/200	2026-03-21 08:25 by laoshidan
[考研] 材料工程（專）一志愿985 初試335求調(diào)劑 +3	hiloiy 2026-03-17	4/200	2026-03-21 03:04 by JourneyLucky
[考研] 085700資源與環(huán)境308求調(diào)劑 +12	墨墨漠 2026-03-18	13/650	2026-03-21 01:42 by JourneyLucky
[考研] 274求調(diào)劑 +10	S.H1 2026-03-18	10/500	2026-03-20 23:51 by JourneyLucky
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3	hisfailed 2026-03-19	6/300	2026-03-20 23:43 by hisfailed
[考研] 考研調(diào)劑求學(xué)校推薦 +3	伯樂29 2026-03-18	5/250	2026-03-20 22:59 by JourneyLucky
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9	Doleres 2026-03-19	9/450	2026-03-20 22:36 by JourneyLucky
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4	@taotao 2026-03-19	4/200	2026-03-20 22:14 by JourneyLucky
[考研] 085600材料與化工 +8	安全上岸！ 2026-03-16	8/400	2026-03-20 22:13 by luoyongfeng
[考研] 一志愿武理材料工程348求調(diào)劑 +3	￣^￣゜汗 2026-03-19	4/200	2026-03-20 21:01 by zhukairuo
[考研] 工科材料085601 279求調(diào)劑 +7	困于星晨 2026-03-17	9/450	2026-03-20 17:38 by 無懈可擊111
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3	小木蟲200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm

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

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

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

» 猜你喜歡

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

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

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