| 7 | 1/1 | 返回列表 |
| 查看: 1634 | 回復(fù): 6 | ||||
zhangguangping木蟲(chóng) (著名寫(xiě)手)
|
[交流]
【求助】有哪位高手能夠在這兒給大家講解一下siesta中如何通過(guò)simplex去優(yōu)化基組! 已有3人參與
|
|
如題,在siesta的/Util/Optimizer目錄下有一個(gè)代碼叫simplex,它能夠?qū)w系優(yōu)化基組,哪位高手在此給大家詳細(xì)講解一下他的原理和使用方法啊。里面自帶的文檔比較少,自己看起來(lái)比較費(fèi)勁。 在此謝謝高手了! |
第一性原理和電化學(xué) |

銀蟲(chóng) (小有名氣)
木蟲(chóng) (著名寫(xiě)手)
|
正如下面的帖子: 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ù)。所以出來(lái)這個(gè)問(wèn)題。 怎么樣使用simplex來(lái)優(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 ] |

銀蟲(chóng) (小有名氣)
|
我知道你的問(wèn)題了,因?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)化方法的參量,類似步長(zhǎng),判斷收斂的標(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,你需要寫(xiě)個(gè)提交腳本 至此,就是我的理解 |
木蟲(chóng) (著名寫(xiě)手)

銀蟲(chóng) (小有名氣)
木蟲(chóng) (著名寫(xiě)手)

| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +6 | 墨墨漠 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑 +10 | 妮妮ninicgb 2026-03-15 | 14/700 |
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +8 | zlingli 2026-03-13 | 8/400 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 一志愿天大材料與化工(085600)總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚(yú) 2026-03-18 | 6/300 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風(fēng) 2026-03-17 | 5/250 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛(ài)生物生物愛(ài)?/a> 2026-03-17 | 3/150 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[基金申請(qǐng)]
今年的國(guó)基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化,一志愿:哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|