| 5 | 1/1 | 返回列表 |
| 查看: 2059 | 回復(fù): 13 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
[交流]
【求助】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 |
siesta問(wèn)題 |
銅蟲 (初入文壇)
銅蟲 (初入文壇)
金蟲 (正式寫手)
莫愁莫愁
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 359求調(diào)劑 +6 | 鄧鄧鄧書書 2026-03-25 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 一志愿西安交大材料學(xué)碩(英一數(shù)二)347,求調(diào)劑到高分子/材料相關(guān)專業(yè) +3 | zju51 2026-03-31 | 5/250 |
|
|
[考研] 299求調(diào)劑 +8 | 嗯嗯嗯嗯2 2026-03-27 | 8/400 |
|
|
[考研] 308求調(diào)劑 +11 | 墨墨漠 2026-03-25 | 11/550 |
|
|
[考研] 考研生物與醫(yī)藥調(diào)劑 +6 | 鐵憨憨123425 2026-03-31 | 6/300 |
|
|
[考研] 英一數(shù)一總分334求調(diào)劑 +4 | 陳陽(yáng)坤 2026-03-31 | 4/200 |
|
|
[考研] 272求調(diào)劑,接受跨專業(yè)調(diào)劑! +3 | 閑魚盧 2026-03-31 | 3/150 |
|
|
[考研] 0856調(diào)劑 +6 | 曲聽筠 2026-03-30 | 6/300 |
|
|
[考研] 370求調(diào)劑 +3 | 080700調(diào)劑 2026-03-30 | 3/150 |
|
|
[考研] 本科211總分289,08工學(xué)真心求調(diào)劑 +3 | utopiaE 2026-03-30 | 3/150 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +9 | hanamiko 2026-03-27 | 9/450 |
|
|
[考研] 抱歉 +3 | 田洪有 2026-03-30 | 3/150 |
|
|
[有機(jī)交流] 考研調(diào)劑 +8 | watb 2026-03-26 | 8/400 |
|
|
[考研] 275求調(diào)劑 +15 | Micky11223 2026-03-25 | 20/1000 |
|
|
[考研] 340求調(diào)劑 +6 | Amber00 2026-03-26 | 6/300 |
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|