| 5 | 1/1 | 返回列表 |
| 查看: 1641 | 回復(fù): 6 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
zhangguangping木蟲 (著名寫手)
|
[交流]
【求助】有哪位高手能夠在這兒給大家講解一下siesta中如何通過simplex去優(yōu)化基組! 已有3人參與
|
|||
|
如題,在siesta的/Util/Optimizer目錄下有一個(gè)代碼叫simplex,它能夠?qū)w系優(yōu)化基組,哪位高手在此給大家詳細(xì)講解一下他的原理和使用方法啊。里面自帶的文檔比較少,自己看起來比較費(fèi)勁。 在此謝謝高手了! |
第一性原理和電化學(xué) |

木蟲 (著名寫手)
|
正如下面的帖子: http://www.gaoyang168.com/bbs/viewthread.php?tid=2485045&fpage=1 會(huì)發(fā)現(xiàn),使用的基組會(huì)使結(jié)果影響較大。 使用程序給出的DZP參數(shù)產(chǎn)生的基組計(jì)算的結(jié)果不令人滿意的話,就需要自己手動(dòng)調(diào)整基組。調(diào)整基組的zeta數(shù)目和rc的大小。所以應(yīng)該得到一個(gè)最優(yōu)化的數(shù)值。但是自己憑空設(shè)想是不可能的。而程序自己產(chǎn)生的那個(gè)參數(shù)肯定不是一個(gè)最優(yōu)化的參數(shù)。所以出來這個(gè)問題。 怎么樣使用simplex來優(yōu)化這個(gè)基組。 %block PAO.Basis H 1 n=1 0 2 P 1 5.471 3.953 1.000 1.000 %endblock PAO.Basis 上面的這個(gè)rc是使用DZP參數(shù)程序給出的結(jié)果。 我想在這個(gè)附近對(duì)這個(gè)四個(gè)參數(shù)進(jìn)行優(yōu)化一下,得到最好的結(jié)果。 想知道simplex使用的什么方法做這件事情的,還有就是具體怎么做去找到最優(yōu)化的rc. 謝謝你的討論。 [ Last edited by zhangguangping on 2010-10-14 at 11:58 ] |

銀蟲 (小有名氣)
銀蟲 (小有名氣)
|
我知道你的問題了,因?yàn)槲以?jīng)也想針對(duì)不同的體系,優(yōu)化具體原子軌道的rc 就我所知道的,以及看到simplex文件夾下給出的例子,參考其中的TEMPLATE輸入文件,將其中要優(yōu)化的rc,改成參數(shù)的形式$s_z1_rc_O之類的,然后再看另一個(gè)文件VARS,設(shè)置對(duì)應(yīng)參數(shù)的取值范圍,第三個(gè)輸入文件就是PARAMS,這里面給出的應(yīng)該是優(yōu)化方法的參量,類似步長,判斷收斂的標(biāo)準(zhǔn)之類,設(shè)置好輸入文件后,修改run_script.sh腳本,SIESTA=$HOME/bin/siesta,類似這種執(zhí)行文件的位置,name=$1應(yīng)該是你執(zhí)行該腳本時(shí)輸入的任務(wù)名字吧,設(shè)置好后,直接執(zhí)行simplex應(yīng)該就可以運(yùn)行了,當(dāng)然如果是queuing-system,你需要寫個(gè)提交腳本 至此,就是我的理解 |
木蟲 (著名寫手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0854AI CV方向招收調(diào)劑 +3 | 章小魚567 2026-03-23 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
[考研]
|
pk3725069 2026-03-19 | 17/850 |
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|
|
[考博] 申博26年 +4 | 八6八68 2026-03-19 | 4/200 |
|
|
[考研] 307求調(diào)劑 +5 | 超級(jí)伊昂大王 2026-03-24 | 5/250 |
|
|
[考研] 305分求調(diào)劑(食品工程) +5 | Sxy112 2026-03-21 | 7/350 |
|
|
[考研] 材料與化工328分調(diào)劑 +4 | 。,。,。,。i 2026-03-23 | 4/200 |
|
|
[考研] 材料專碩英一數(shù)二306 +8 | z1z2z3879 2026-03-18 | 8/400 |
|
|
[考研] 生物學(xué)一志愿985,分?jǐn)?shù)349求調(diào)劑 +6 | zxts12 2026-03-21 | 9/450 |
|
|
[考研]
|
2117205181 2026-03-21 | 8/400 |
|
|
[考研] 324求調(diào)劑 +6 | lucky呀呀呀鴨 2026-03-20 | 6/300 |
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
[考研]
|
Grand777 2026-03-21 | 3/150 |
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
[考研] 22 350 本科985求調(diào)劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研]
|
然11 2026-03-19 | 4/200 |
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|