24小时热门版块排行榜    

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

zhangguangping

木虫 (著名写手)

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

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

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

第一性原理和电化学

» 猜你喜欢

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

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

zhangguangping

木虫 (著名写手)

引用回帖:
Originally posted by yangjun102 at 2010-12-09 11:40:20:



楼主是否使用了Optimizer优化了基组,我看了一下比较复杂,不知如何下手,楼主能否分享一下使用的方法?

不前还不会。呵呵呵
弘德明志博学笃行
7楼2010-12-12 21:33:47
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +4 1孙悟空 2026-03-17 4/200 2026-03-18 17:59 by fivewind
[考研] 一志愿武理材料305分求调剂 +5 想上岸的鲤鱼 2026-03-18 6/300 2026-03-18 17:53 by 无际的草原
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +8 Ncdx123456 2026-03-13 9/450 2026-03-18 14:40 by haxia
[考研] 085601材料工程专硕求调剂 +6 慕寒mio 2026-03-16 6/300 2026-03-18 14:26 by 007_lilei
[考研] 280求调剂 +6 咕噜晓晓 2026-03-18 7/350 2026-03-18 11:25 by 无际的草原
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 332求调剂 +6 Zz版 2026-03-13 6/300 2026-03-17 17:03 by ruiyingmiao
[考研] 290求调剂 +6 孔志浩 2026-03-12 11/550 2026-03-17 14:41 by 周舟舟77
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[基金申请] 国自科面上基金字体 +6 iwuli 2026-03-12 7/350 2026-03-16 21:18 by sculhf
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 一志愿华中师范071000,325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 328化工专硕求调剂 +4 。,。,。,。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考研] 321求调剂(食品/专硕) +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
信息提示
请填处理意见