24小时热门版块排行榜    

查看: 1637  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zhangguangping

木虫 (著名写手)

[交流] 【求助】有哪位高手能够在这儿给大家讲解一下siesta中如何通过simplex去优化基组! 已有3人参与

如题,在siesta的/Util/Optimizer目录下有一个代码叫simplex,它能够对体系优化基组,哪位高手在此给大家详细讲解一下他的原理和使用方法啊。里面自带的文档比较少,自己看起来比较费劲。
在此谢谢高手了!
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理和电化学

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

弘德明志博学笃行
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

引用回帖:
Originally posted by corbin at 2010-10-14 03:03:25:


我知道你的问题了,因为我曾经也想针对不同的体系,优化具体原子轨道的rc

就我所知道的,以及看到simplex文件夹下给出的例子,参考其中的TEMPLATE输入文件,将其中要优化的rc,改成参数的形式$s_z1_rc_O之 ...

谢谢你的指点!
弘德明志博学笃行
5楼2010-10-14 19:23:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

corbin

银虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
曾经试图用过,觉得就是在那几个输入文件中设置想要优化的参数,但是对于如何提交到cluster上,没具体用过,你可以说说具体问题,咱探讨探讨
2楼2010-10-14 08:40:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

引用回帖:
Originally posted by corbin at 2010-10-14 01:40:36:
曾经试图用过,觉得就是在那几个输入文件中设置想要优化的参数,但是对于如何提交到cluster上,没具体用过,你可以说说具体问题,咱探讨探讨

正如下面的帖子:
http://muchong.com/bbs/viewthread.php?tid=2485045&fpage=1
会发现,使用的基组会使结果影响较大。
使用程序给出的DZP参数产生的基组计算的结果不令人满意的话,就需要自己手动调整基组。调整基组的zeta数目和rc的大小。所以应该得到一个最优化的数值。但是自己凭空设想是不可能的。而程序自己产生的那个参数肯定不是一个最优化的参数。所以出来这个问题。
怎么样使用simplex来优化这个基组。
%block PAO.Basis               
H                1                  
n=1   0   2 P   1                  
   5.471      3.953   
   1.000      1.000   
%endblock PAO.Basis
上面的这个rc是使用DZP参数程序给出的结果。
我想在这个附近对这个四个参数进行优化一下,得到最好的结果。
想知道simplex使用的什么方法做这件事情的,还有就是具体怎么做去找到最优化的rc.
谢谢你的讨论。

[ Last edited by zhangguangping on 2010-10-14 at 11:58 ]
弘德明志博学笃行
3楼2010-10-14 09:39:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

corbin

银虫 (小有名气)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+2):谢谢讨论 2010-10-14 10:07:09
引用回帖:
Originally posted by zhangguangping at 2010-10-14 09:39:25:

正如下面的帖子:
http://muchong.com/bbs/viewthread.php?tid=2485045&fpage=1
会发现,使用的基组会使结果影响较大。
使用程序给出的DZP参数产生的基组计算的结果不令人满意的话,就需要自己 ...

我知道你的问题了,因为我曾经也想针对不同的体系,优化具体原子轨道的rc

就我所知道的,以及看到simplex文件夹下给出的例子,参考其中的TEMPLATE输入文件,将其中要优化的rc,改成参数的形式$s_z1_rc_O之类的,然后再看另一个文件VARS,设置对应参数的取值范围,第三个输入文件就是PARAMS,这里面给出的应该是优化方法的参量,类似步长,判断收敛的标准之类,设置好输入文件后,修改run_script.sh脚本,SIESTA=$HOME/bin/siesta,类似这种执行文件的位置,name=$1应该是你执行该脚本时输入的任务名字吧,设置好后,直接执行simplex应该就可以运行了,当然如果是queuing-system,你需要写个提交脚本

至此,就是我的理解
4楼2010-10-14 10:03:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0805 316求调剂 +3 大雪深藏 2026-03-18 3/150 2026-03-21 18:55 by 学员8dgXkO
[考研] 332求调剂 +3 凤凰院丁真 2026-03-20 3/150 2026-03-21 10:27 by luoyongfeng
[考研] 0856材料专硕353求调剂 +3 NIFFFfff 2026-03-20 3/150 2026-03-21 10:23 by luoyongfeng
[考研] 材料工程(专)一志愿985 初试335求调剂 +3 hiloiy 2026-03-17 4/200 2026-03-21 03:04 by JourneyLucky
[考研] 二本跨考郑大材料306英一数二 +3 z1z2z3879 2026-03-17 3/150 2026-03-21 02:29 by JourneyLucky
[考研] 324分 085600材料化工求调剂 +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[考研] 一志愿重庆大学085700资源与环境专硕,总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 材料专业求调剂 +6 hanamiko 2026-03-18 6/300 2026-03-21 00:24 by JourneyLucky
[考研] 296求调剂 +6 www_q 2026-03-18 10/500 2026-03-20 23:56 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 5/250 2026-03-20 22:11 by 云游重阳
[考研] 260求调剂 +3 朱芷琳 2026-03-20 3/150 2026-03-20 20:35 by 学员8dgXkO
[考研] 环境工程调剂 +9 大可digkids 2026-03-16 9/450 2026-03-20 17:38 by 醉在风里
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
信息提示
请填处理意见