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

木蟲 (著名寫手)

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

銀蟲 (小有名氣)
|
我知道你的問題了,因為我曾經(jīng)也想針對不同的體系,優(yōu)化具體原子軌道的rc 就我所知道的,以及看到simplex文件夾下給出的例子,參考其中的TEMPLATE輸入文件,將其中要優(yōu)化的rc,改成參數(shù)的形式$s_z1_rc_O之類的,然后再看另一個文件VARS,設(shè)置對應(yīng)參數(shù)的取值范圍,第三個輸入文件就是PARAMS,這里面給出的應(yīng)該是優(yōu)化方法的參量,類似步長,判斷收斂的標準之類,設(shè)置好輸入文件后,修改run_script.sh腳本,SIESTA=$HOME/bin/siesta,類似這種執(zhí)行文件的位置,name=$1應(yīng)該是你執(zhí)行該腳本時輸入的任務(wù)名字吧,設(shè)置好后,直接執(zhí)行simplex應(yīng)該就可以運行了,當然如果是queuing-system,你需要寫個提交腳本 至此,就是我的理解 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 南京大學化學調(diào)劑 +9 | 景隨風 2026-03-29 | 14/700 |
|
|---|---|---|---|---|
|
[考博] 26申博自薦 +6 | whh869393 2026-03-24 | 6/300 |
|
|
[考研] 材料科學與工程求調(diào)劑 +6 | 深V宿舍吧 2026-03-29 | 6/300 |
|
|
[考研] 296求調(diào)劑 +6 | 彼岸t 2026-03-29 | 6/300 |
|
|
[考研] 調(diào)劑310 +12 | 溫柔的晚安 2026-03-25 | 13/650 |
|
|
[考研] 一志愿北京化工大學材料與化工(085600)296求調(diào)劑 +17 | 稻妻小編 2026-03-26 | 17/850 |
|
|
[考研] 322求調(diào)劑 +7 | 宋明欣 2026-03-27 | 7/350 |
|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +12 | 路癡小琪 2026-03-27 | 12/600 |
|
|
[考研] 312,生物學求調(diào)劑 +3 | 小譯同學abc 2026-03-28 | 3/150 |
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
[考研] 286求調(diào)劑 +4 | 丟掉懶惰 2026-03-27 | 7/350 |
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學-環(huán)境工程專碩-298 +4 | 手機用戶 2026-03-26 | 4/200 |
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 07化學280分求調(diào)劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
[考研] 一志愿陜師大生物學071000,298分,求調(diào)劑 +5 | SYA! 2026-03-23 | 5/250 |
|
|
[考研] 325求調(diào)劑 +5 | 李嘉圖·S·路 2026-03-23 | 5/250 |
|
|
[考研] 0703化學求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|