24小时热门版块排行榜    

查看: 2689  |  回复: 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. 最后对于赝势一般是用软件提供的赝势,这里也能够优化来提高准确度?需要怎么做?
先谢谢了!!!
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lqh91

木虫 (正式写手)

引用回帖:
4楼: Originally posted by c_chi at 2014-09-11 14:45:11
问题还是没解决,我主要是想优化MnO2的基组,我看了下例子不是很明白,优化参数和原子轨道有什么关系吗?能解释下下面的内容:
1.有个如下的VASRS文件,应该是优化的阈值,但是怎么取?有什么规则还不是很明白?为 ...

请问楼主弄清楚这些值代表的意思了吗?最近我也在用siesta算氧化铁,需要优化基组,但是完全不懂这些值代表的意思,请问楼主可否赐教?
6楼2019-02-28 22:55:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

liqizuiyang

木虫 (著名写手)

★ ★ ★
franch: 金币+3, 谢谢回帖交流。。。 2014-09-22 22:29:11
推荐楼主去看下openmx的帮助文档和项目组发的文章,我记得有篇PRB文章专门讲基组优化的。虽然方法不同,但都是LCAO,应该有相通的部分。

LCAO使用原子轨道作为基组,原子轨道可分为径向函数和球谐函数两部分。球谐函数有确定的表达式,径向函数有slater型,gauss型和数值型。数值型就是把径向函数分成一系列格点,存储格点上的函数值,计算时再用插值法还原。基组优化实际上就是调整格点上的值。
2楼2014-09-09 16:01:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c_chi

木虫 (正式写手)

引用回帖:
2楼: Originally posted by liqizuiyang at 2014-09-09 16:01:22
推荐楼主去看下openmx的帮助文档和项目组发的文章,我记得有篇PRB文章专门讲基组优化的。虽然方法不同,但都是LCAO,应该有相通的部分。

LCAO使用原子轨道作为基组,原子轨道可分为径向函数和球谐函数两部分。球 ...

谢谢,我去研究下opemx的文档和文章。
3楼2014-09-10 08:30:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c_chi

木虫 (正式写手)

问题还是没解决,我主要是想优化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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-28 05:13 by winsaint
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +11 晴天加油 2026-08-26 12/600 2026-08-28 01:22 by li09465290@
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +9 Kittylucky 2026-08-27 9/450 2026-08-28 01:08 by jurkat.1640
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 XLGfUIdOmEAO 2026-08-27 4/200 2026-08-27 22:26 by 3RAFyYXYSwGS
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见