| 5 | 1/1 | 返回列表 |
| 查看: 2311 | 回復(fù): 7 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
c_chi木蟲 (正式寫手)
|
[求助]
siesta基組優(yōu)化問題 已有1人參與
|
||
|
和PW方法相比,LCAO方法有時(shí)候的計(jì)算準(zhǔn)確度確實(shí)有所不夠,特別是對(duì)于過渡金屬和貴金屬感覺不太容易計(jì)算,之前用siesta優(yōu)化一個(gè)Beta相的MnO2,發(fā)現(xiàn)結(jié)果基本是錯(cuò)的。問了些人,也在論壇里翻了其他人的帖子,說(shuō)是可以通過優(yōu)化基組來(lái)達(dá)到PW的準(zhǔn)確度(工具:/Util/Optimizer目錄下的simplex),但是還是不怎么明白。 http://www.gaoyang168.com/bbs/viewthread.php?tid=2486484&fpage=1&target=blank http://www.gaoyang168.com/bbs/viewthread.php?tid=2485045&fpage=1&target=self&page=1 http://www.gaoyang168.com/bbs/viewthread.php?tid=3146982 也看了高手提供的一篇這方面的文章,看了下,沒能全看懂。 Optimal strictly localized basis sets for noble metal surfaces.(PHYSICAL REVIEW B79, 075441 2009) 這里希望有人能夠解一下疑惑: 1詳細(xì).解釋下基組,大家經(jīng)驗(yàn)來(lái)看什么情況下使用什么基組可以達(dá)到可行的精度? 2.解釋下siesta基組優(yōu)化的具體方法和技巧,或者提供這方面的資料?或者可不可以參考高斯的基組來(lái)設(shè)置siesta的基組? 3.我在siesta計(jì)算過渡金屬氧化物時(shí),感覺scf比較難以收斂,如果優(yōu)化基組后可能更難,計(jì)算量更大,我一般的做法是加大DM.NumberPulay并減小DM.MixingWeight,如下。請(qǐng)問我的做法合理嗎?還有什么其他方法? DM.NumberPulay 8 DM.MixingWeight 0.1 4. 最后對(duì)于贗勢(shì)一般是用軟件提供的贗勢(shì),這里也能夠優(yōu)化來(lái)提高準(zhǔn)確度?需要怎么做? 先謝謝了。! |
木蟲 (正式寫手)
|
問題還是沒解決,我主要是想優(yōu)化MnO2的基組,我看了下例子不是很明白,優(yōu)化參數(shù)和原子軌道有什么關(guān)系嗎?能解釋下下面的內(nèi)容: 1.有個(gè)如下的VASRS文件,應(yīng)該是優(yōu)化的閾值,但是怎么?有什么規(guī)則還不是很明白?為什么要這樣取這些軌道? s_z1_rc_O 2.5 8.0 spl_s_O 0.05 0.5 0.15 p_z1_rc_O 2.5 8.0 spl_p_O 0.05 0.5 d_z1_rc_O 2.5 8.0 s_z1_rc_H 2.0 8.0 spl_s_H 0.05 0.8 p_z1_rc_H 2.0 8.0 2.同時(shí)運(yùn)行的輸入文件有個(gè)基組的block,這個(gè)和上面的優(yōu)化閾值有什么關(guān)系嗎?謝謝了。! %block PAO.Basis # Define Basis set O 3 n=2 0 2 S $spl_s_O $s_z1_rc_O 0.0 1.000 1.000 n=2 1 2 S $spl_p_O $p_z1_rc_O 0.0 1.000 1.000 n=3 2 1 $d_z1_rc_O 1.000 H 2 n=1 0 2 S $spl_s_H $s_z1_rc_H 0.0 1.000 1.000 n=2 1 1 $p_z1_rc_H 1.000 %endblock PAO.Basis |
木蟲 (著名寫手)
|
推薦樓主去看下openmx的幫助文檔和項(xiàng)目組發(fā)的文章,我記得有篇PRB文章專門講基組優(yōu)化的。雖然方法不同,但都是LCAO,應(yīng)該有相通的部分。 LCAO使用原子軌道作為基組,原子軌道可分為徑向函數(shù)和球諧函數(shù)兩部分。球諧函數(shù)有確定的表達(dá)式,徑向函數(shù)有slater型,gauss型和數(shù)值型。數(shù)值型就是把徑向函數(shù)分成一系列格點(diǎn),存儲(chǔ)格點(diǎn)上的函數(shù)值,計(jì)算時(shí)再用插值法還原。基組優(yōu)化實(shí)際上就是調(diào)整格點(diǎn)上的值。 |
木蟲 (正式寫手)
木蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 359求調(diào)劑 +5 | 王了個(gè)楠 2026-03-25 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 085600 286分 材料求調(diào)劑 +11 | 麻辣魷魚 2026-03-27 | 12/600 |
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +18 | Mr. Z 2026-03-25 | 19/950 |
|
|
[考研] 070300化學(xué)354求調(diào)劑 +7 | 101次希望 2026-03-28 | 7/350 |
|
|
[考研] 求調(diào)劑 +7 | 青春裁為三截 2026-03-29 | 7/350 |
|
|
[考研] 求調(diào)劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷 +8 | cccchenso 2026-03-29 | 8/400 |
|
|
[考研]
|
nnnnnnn5 2026-03-25 | 11/550 |
|
|
[考研] 299求調(diào)劑 +10 | 15188958825 2026-03-25 | 10/500 |
|
|
[考研] 329求調(diào)劑 +10 | 鈕恩雪 2026-03-25 | 10/500 |
|
|
[考研] 332求92調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
[考研] 一志愿太原理工安全工程300分,求調(diào)劑 +5 | 0857求調(diào)劑. 2026-03-24 | 6/300 |
|
|
[考研] 求調(diào)劑 +7 | 爭(zhēng)取九點(diǎn)睡 2026-03-28 | 8/400 |
|
|
[考研] 322求調(diào)劑 +5 | 舊吢 2026-03-24 | 5/250 |
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲085600 2026-03-27 | 8/400 |
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|