| 查看: 1513 | 回复: 6 | ||||
zhangguangping木虫 (著名写手)
|
[交流]
【求助】有哪位高手能够在这儿给大家讲解一下siesta中如何通过simplex去优化基组！ 已有3人参与
|
|
如题，在siesta的/Util/Optimizer目录下有一个代码叫simplex,它能够对体系优化基组，哪位高手在此给大家详细讲解一下他的原理和使用方法啊。里面自带的文档比较少，自己看起来比较费劲。 在此谢谢高手了! |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理和电化学 |
» 猜你喜欢
职称评审没过，求安慰
已经有41人回复
-
回收溶剂求助
已经有7人回复
-
硝基苯如何除去
已经有3人回复
-
A期刊撤稿
已经有4人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
求助文献
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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