24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2324  |  回复: 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的回帖

c_chi

木虫 (正式写手)

没有虫子对这个有研究啊?
5楼2014-09-17 14:04:48
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 英一数一408,总分284,二战真诚求调剂 +13 12.27 2026-03-30 15/750 2026-04-03 14:41 by 氮气气气
[考研] 0854求调剂 +3 assdll 2026-04-03 3/150 2026-04-03 14:32 by fxue1114
[考研] 一志愿华北电力大学(北京),材料科学与工程学硕265,求调剂 +9 yelck 2026-04-03 10/500 2026-04-03 13:20 by yelck
[考研] 材料与化工306分找调剂 +21 沧海轻舟e 2026-04-02 22/1100 2026-04-03 12:47 by xingguangj
[考研] 309求调剂 +4 刘刘刘1231 2026-04-02 5/250 2026-04-03 12:04 by 1753564080
[考研] 初试成绩337找调剂 +3 ??? ?. ? 2026-04-03 3/150 2026-04-03 11:43 by 土木硕士招生
[考研] 279求调剂 +6 qazplm0852 2026-04-02 6/300 2026-04-03 10:03 by 蓝云思雨
[考研] 一志愿安徽大学0817化学工程与技术,求调剂 +14 我不是只因 2026-04-02 15/750 2026-04-03 09:49 by 蓝云思雨
[考研] 312求调剂 +4 赊月色 2026-04-02 5/250 2026-04-03 08:21 by fangshan711
[基金申请] esi高被引论文是不是能对中标有所加分和帮助呢 +4 redcom 2026-04-01 5/250 2026-04-02 20:48 by redcom
[考研] 279求调剂 +5 傅文秋 2026-04-02 5/250 2026-04-02 18:10 by 笔落锦州
[考研] 求调剂 302分初试 0854 +5 伶可乐 2026-04-02 5/250 2026-04-02 17:53 by 笔落锦州
[考研] 318求调剂 +3 笃行致远. 2026-03-31 4/200 2026-04-02 15:56 by Jaylen.
[考研] 282求调剂 +13 呼吸都是减肥 2026-04-01 13/650 2026-04-02 14:10 by baoball
[考研] 266求调剂 +4 学员97LZgn 2026-04-02 4/200 2026-04-02 13:03 by yulian1987
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +10 哇呼哼呼哼 2026-04-01 11/550 2026-04-02 11:31 by lnilvy
[考研] 311求调剂 +16 zchqwer 2026-04-01 18/900 2026-04-02 10:15 by 不吃魚的貓
[考研] 301求调剂 +13 A_JiXing 2026-04-01 13/650 2026-04-02 09:01 by sanrepian
[考研] 求调剂 生物学 377分 +6 zzll03 2026-03-31 6/300 2026-03-31 17:33 by 唐沐儿
[考研] 085601 329分调剂 +6 yzsa12 2026-03-31 6/300 2026-03-31 15:23 by yanflower7133
信息提示
请填处理意见