24小时热门版块排行榜    

查看: 839  |  回复: 3

nigelallan

银虫 (初入文坛)

[求助] 对于MgNi团簇使用什么高精度的方法和基组比较合适

我是用b3lyp/lanl2dz来算的,审稿人说选用高精度的计算方法和基组的计算MgNi, Ni2,和Mg2;将计算结果与之比较。我应该选用什么方法和基组比较合适呢?
回复此楼

» 猜你喜欢

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

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

冬天里的骄阳

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★
gmy1990(金币+2): 2011-07-07 23:15:05
跟实验值键长,电离能,亲和能等等对比,符合实验值的就是合适的。真要提高精度的话,考虑下修正后的lanl2dz,cc-pVTZ-pp这样的
2楼2011-07-07 23:06:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

【答案】应助回帖

nigelallan(金币+5): 感谢大牛指点! 2011-07-08 09:19:01
这些分子都属于强关联体系,不能用常规方法去算。刚去世不久的大牛Roos为了算好含Ni的小体系,花费了20多年和别人写了CASSCF和CASPT2程序。需要用非常大的活性空间,结合至少4zeta级别的CC-PVnZ或ANO基组,才能得到正确结果:J. Chem. Phys. 101, 4893, 1994.

而你用B3LYP/DZ基组去算这类小分子,这种计算水平放到2000年以前还说得过去,今天肯定是不行的。

一点建议:

1。最好的结果,用CASPT2+相对论效应+4zeta基组。
好处:不仅得到基态和一大堆激发态;还可以考虑旋轨耦合,每个分子的结果都够发一篇文章。
缺点:CASPT2对活性空间的设置极其敏感,弄不好会得到毫无意义的结果。另外,计算量很大。即使熟悉这种计算,没有一个月的机时也算不完一个分子。

2。退而求其次,用CCSD(T)+相对论效应+4zeta基组。
最好换程序。Gaussian、GAMESS慢10至100倍,甚至根本算不动。

3。最容易实现的办法:用杂化泛函+相对论效应。因为DFT对基组质量的依赖性不如Post-HF方法大,基组质量可以降到triple-zeta。
只要别遇到对DFT很挑剔的审稿人,就能蒙混过关。否则还要对不同杂化泛函进行测试,因为强关联体系对杂化泛函的精确交换因子的依赖性很大,B3LYP未必是最佳选择。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
3楼2011-07-08 00:15:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mchen10

金虫 (正式写手)

补充一下楼上大牛的:

双原子体系MRCI / CASPT2计算量还算适中,含Mg原子的计算量很小,Ni2稍稍有点大(20 e in 12 orbital?),但绝对可以算
MR-CCSD(T) (我没用过)应该也适用

CCSD(T)方法算的话,小体系可以做CBS的外推 (NZ, N=2,3,4,5或N=2,3,4)
更痛苦的是梦没醒路已经不能走了
4楼2011-07-08 06:04:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 nigelallan 的主题更新
信息提示
请填处理意见