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

返回列表

下一頁(yè)

leila3409

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

應(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)題

» 猜你喜歡

求調(diào)劑院校信息已經(jīng)有4人回復(fù)
085600材料與化工306 已經(jīng)有4人回復(fù)
286求調(diào)劑已經(jīng)有10人回復(fù)
328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷已經(jīng)有9人回復(fù)
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑已經(jīng)有4人回復(fù)
286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！已經(jīng)有8人回復(fù)
資源與環(huán)境調(diào)劑申請(qǐng)(333分) 已經(jīng)有5人回復(fù)
280求調(diào)劑已經(jīng)有12人回復(fù)
269專碩求調(diào)劑已經(jīng)有5人回復(fù)
材料學(xué)碩301分求調(diào)劑已經(jīng)有7人回復(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的回帖

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

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

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

minmin_0082003

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

莫愁莫愁

1ST強(qiáng)帖: 1
應(yīng)助: 2 (幼兒園)
貴賓: 0.08
金幣: 921.3
散金: 149
帖子: 910
在線: 218.4小時(shí)
蟲(chóng)號(hào): 366989
注冊(cè): 2007-05-10
性別: MM
專業(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的回帖

zhangguangping

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

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è)紅包，謝謝回帖交流

引用回帖:

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

請(qǐng)問(wèn)，怎么優(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

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

莫愁莫愁

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

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

引用回帖:

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

不好意思，生病了，沒(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ù)。比如你希望你的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

金蟲(chóng) (小有名氣)

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

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

引用回帖:

Originally posted by minmin_0082003 at 2011-03-06 16: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 ...

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

贊一下(1人)

回復(fù)此樓

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

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

zhangguangping

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

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

minmin_0082003

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

莫愁莫愁

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

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

引用回帖:

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

對(duì)的，我沒(méi)有同時(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 的主題更新

返回列表

下一頁(yè)

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

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

[考研] 299求調(diào)劑 +4	某某某某位 2026-03-21	4/200	2026-03-21 16:30 by barlinike
[考研] 085601調(diào)劑 358分 +3	zzzzggh 2026-03-20	4/200	2026-03-21 10:21 by luoyongfeng
[考研] 306求調(diào)劑 +4	chuanzhu川燭 2026-03-18	4/200	2026-03-21 08:25 by laoshidan
[考研] 265求調(diào)劑 +3	Jack?k?y 2026-03-17	3/150	2026-03-21 03:17 by JourneyLucky
[考研] 材料專業(yè)求調(diào)劑 +6	hanamiko 2026-03-18	6/300	2026-03-21 00:24 by JourneyLucky
[考研] 324求調(diào)劑 +5	lucky呀呀呀鴨 2026-03-20	5/250	2026-03-20 22:30 by 促天成
[考研] 317求調(diào)劑 +5	申子申申 2026-03-19	9/450	2026-03-20 22:26 by JourneyLucky
[考研] 085600材料與化工 +8	安全上岸！ 2026-03-16	8/400	2026-03-20 22:13 by luoyongfeng
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +3	困困困困坤坤 2026-03-20	3/150	2026-03-20 20:38 by 學(xué)員8dgXkO
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +5	Charlieyq 2026-03-19	5/250	2026-03-20 20:35 by JourneyLucky
[考研] 材料學(xué)求調(diào)劑 +4	Stella_Yao 2026-03-20	4/200	2026-03-20 20:28 by ms629
[考研] 材料與化工專碩調(diào)劑 +7	heming3743 2026-03-16	7/350	2026-03-20 19:31 by zhukairuo
[考研] 0856調(diào)劑，是學(xué)校就去 +8	sllhht 2026-03-19	9/450	2026-03-20 14:25 by 無(wú)懈可擊111
[考研] 求調(diào)劑 +3	暗涌afhb 2026-03-16	3/150	2026-03-20 00:28 by 河南大學(xué)校友
[論文投稿] 申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變，給編輯又寫(xiě)郵件了，沒(méi)回復(fù) 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 0856求調(diào)劑 +3	劉夢(mèng)微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015

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

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

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

» 猜你喜歡

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

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

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