| 7 | 1/1 | 返回列表 |
| 查看: 1635 | 回復(fù): 6 | ||||
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)化方法的參量,類似步長(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,你需要寫個(gè)提交腳本 至此,就是我的理解 |
木蟲 (著名寫手)

銀蟲 (小有名氣)
木蟲 (著名寫手)

| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 274求調(diào)劑 +7 | S.H1 2026-03-18 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考! +3 | lemonzzn 2026-03-17 | 4/200 |
|
|
[考研]
|
.6lL 2026-03-18 | 7/350 |
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +3 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 5/250 |
|
|
[考研] 321求調(diào)劑 +8 | 何潤(rùn)采123 2026-03-18 | 10/500 |
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
[考研] 一志愿985,本科211,0817化學(xué)工程與技術(shù)319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院,一志愿華科070200學(xué)碩,346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽(yáng)花天天開心 2026-03-16 | 4/200 |
|
|
[考研] 一志愿,福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭(zhēng) 2026-03-15 | 3/150 |
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|