| 查看: 1622 | 回复: 6 | ||||
zhangguangping木虫 (著名写手)
|
[交流]
【求助】有哪位高手能够在这儿给大家讲解一下siesta中如何通过simplex去优化基组！ 已有3人参与
|
|
如题，在siesta的/Util/Optimizer目录下有一个代码叫simplex,它能够对体系优化基组，哪位高手在此给大家详细讲解一下他的原理和使用方法啊。里面自带的文档比较少，自己看起来比较费劲。 在此谢谢高手了! |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理和电化学 |
» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Siesta 安装的问题。。。求高手指点
已经有4人回复
- 求助siesta收敛问题
已经有3人回复
- 【求助】安装siesta3.0出错：segmentation violation signal raised。请高手指教
已经有8人回复
- 【求助】谁能具体说一下各个基组的优缺点
已经有5人回复
corbin
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 457.1
- 散金: 4
- 帖子: 113
- 在线: 93.5小时
- 虫号: 1045382
- 注册: 2010-06-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2010-10-14 08:40:36
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
|
正如下面的帖子： 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
corbin
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 457.1
- 散金: 4
- 帖子: 113
- 在线: 93.5小时
- 虫号: 1045382
- 注册: 2010-06-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2):谢谢讨论 2010-10-14 10:07:09
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2):谢谢讨论 2010-10-14 10:07:09
|
我知道你的问题了，因为我曾经也想针对不同的体系，优化具体原子轨道的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
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
5楼2010-10-14 19:23:05
yangjun102
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 486.5
- 散金: 24
- 帖子: 248
- 在线: 351.2小时
- 虫号: 763027
- 注册: 2009-05-04
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
6楼2010-12-09 18:40:20
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
7楼2010-12-12 21:33:47
