| 查看: 2323 | 回复: 7 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
c_chi木虫 (正式写手)
|
[求助]
siesta基组优化问题 已有1人参与
|
|
|
和PW方法相比，LCAO方法有时候的计算准确度确实有所不够，特别是对于过渡金属和贵金属感觉不太容易计算，之前用siesta优化一个Beta相的MnO2，发现结果基本是错的。问了些人，也在论坛里翻了其他人的帖子，说是可以通过优化基组来达到PW的准确度（工具：/Util/Optimizer目录下的simplex），但是还是不怎么明白。 http://muchong.com/bbs/viewthread.php?tid=2486484&fpage=1&target=blank http://muchong.com/bbs/viewthread.php?tid=2485045&fpage=1&target=self&page=1 http://muchong.com/bbs/viewthread.php?tid=3146982 也看了高手提供的一篇这方面的文章，看了下，没能全看懂。 Optimal strictly localized basis sets for noble metal surfaces.（PHYSICAL REVIEW B79, 075441 2009） 这里希望有人能够解一下疑惑： 1详细.解释下基组，大家经验来看什么情况下使用什么基组可以达到可行的精度？ 2.解释下siesta基组优化的具体方法和技巧，或者提供这方面的资料？或者可不可以参考高斯的基组来设置siesta的基组？ 3.我在siesta计算过渡金属氧化物时，感觉scf比较难以收敛，如果优化基组后可能更难，计算量更大，我一般的做法是加大DM.NumberPulay并减小DM.MixingWeight，如下。请问我的做法合理吗？还有什么其他方法？ DM.NumberPulay 8 DM.MixingWeight 0.1 4. 最后对于赝势一般是用软件提供的赝势，这里也能够优化来提高准确度？需要怎么做？ 先谢谢了！！！ |
回复此楼» 猜你喜欢
材料调剂
已经有4人回复
-
266分，求材料相关专业调剂
已经有13人回复
-
315求调剂
已经有6人回复
-
材料专硕 调剂
已经有10人回复
-
一志愿华北电力大学（北京），材料科学与工程学硕265，求调剂
已经有10人回复
-
一志愿0817化学工程与技术，求调剂
已经有13人回复
-
265求调剂
已经有17人回复
-
08工科，295，接受跨专业调剂
已经有8人回复
-
274求调剂
已经有8人回复
-
320分，材料与化工专业，求调剂
已经有19人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 基组or结构问题，B3lyp优化成功，B3PW91/WB97XD优化全出错
已经有16人回复
- 请教如何加快siesta的收敛速度
已经有9人回复
- Siesta从官网下载赝势计算出错
已经有7人回复
- 求助siesta收敛问题
已经有3人回复
- 赝试基组输入出错，求教
已经有5人回复
- siesta结构优化收敛标准
已经有5人回复
- SIESTA计算能级不准怎么解决
已经有10人回复
- 求助siesta优化问题
已经有3人回复
- siesta 计算电子占据数的问题
已经有3人回复
- 求助赝势基组的问题
已经有9人回复
- 关于SIESTA中如何计算功函数的一种做法探讨
已经有8人回复
- 手动写入原子轨道基组的原则和技巧
已经有5人回复
- SIESTA几何优化不收敛是不是就说明了这个结构是不稳定的?
已经有4人回复
- siesta计算弱相互作用VDW
已经有17人回复
- 【求助】siesta 优化晶胞
已经有3人回复
- 【分享】SIESTA中Bader analysis
已经有43人回复
- 【求助】siesta 计算DOS应该注意什么问题
已经有3人回复
- 【讨论】关于如何优化大体系的问题
已经有9人回复
- 【求助】siesta多个/种分子优化的问题
已经有7人回复
- 【求助】siesta输入文件中赝势文件的作用？
已经有8人回复
- 【求助】谁能具体说一下各个基组的优缺点
已经有5人回复
- 【求助】siesta中关于几何结构优化的问题【完结】
已经有17人回复
- 【求助】求gamess 输入使用赝势基组与分裂基组的例子
已经有6人回复
c_chi
木虫 (正式写手)
- 应助: 17 (小学生)
- 金币: 2449.7
- 散金: 228
- 红花: 2
- 帖子: 551
- 在线: 493.6小时
- 虫号: 1253512
- 注册: 2011-04-02
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
5楼2014-09-17 14:04:48
liqizuiyang
木虫 (著名写手)
- 应助: 59 (初中生)
- 金币: 6859.5
- 散金: 261
- 红花: 94
- 帖子: 2093
- 在线: 708.9小时
- 虫号: 2366686
- 注册: 2013-03-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2014-09-09 16:01:22
c_chi
木虫 (正式写手)
- 应助: 17 (小学生)
- 金币: 2449.7
- 散金: 228
- 红花: 2
- 帖子: 551
- 在线: 493.6小时
- 虫号: 1253512
- 注册: 2011-04-02
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
3楼2014-09-10 08:30:24
c_chi
木虫 (正式写手)
- 应助: 17 (小学生)
- 金币: 2449.7
- 散金: 228
- 红花: 2
- 帖子: 551
- 在线: 493.6小时
- 虫号: 1253512
- 注册: 2011-04-02
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
|
问题还是没解决，我主要是想优化MnO2的基组，我看了下例子不是很明白，优化参数和原子轨道有什么关系吗？能解释下下面的内容： 1.有个如下的VASRS文件，应该是优化的阈值，但是怎么取？有什么规则还不是很明白？为什么要这样取这些轨道？ s_z1_rc_O 2.5 8.0 spl_s_O 0.05 0.5 0.15 p_z1_rc_O 2.5 8.0 spl_p_O 0.05 0.5 d_z1_rc_O 2.5 8.0 s_z1_rc_H 2.0 8.0 spl_s_H 0.05 0.8 p_z1_rc_H 2.0 8.0 2.同时运行的输入文件有个基组的block，这个和上面的优化阈值有什么关系吗？谢谢了！！！ %block PAO.Basis # Define Basis set O 3 n=2 0 2 S $spl_s_O $s_z1_rc_O 0.0 1.000 1.000 n=2 1 2 S $spl_p_O $p_z1_rc_O 0.0 1.000 1.000 n=3 2 1 $d_z1_rc_O 1.000 H 2 n=1 0 2 S $spl_s_H $s_z1_rc_H 0.0 1.000 1.000 n=2 1 1 $p_z1_rc_H 1.000 %endblock PAO.Basis |
4楼2014-09-11 14:45:11
10