24小时热门版块排行榜    

查看: 1141  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 04:47 by LH5NkK5Ud4bh
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-16 04:11 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 03:59 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-16 03:02 by DgNGHc3h5tPl
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-16 02:50 by DgNGHc3h5tPl
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-16 02:14 by DgNGHc3h5tPl
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 19:59 by CgyNCDVNhGVg
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-15 19:23 by CgyNCDVNhGVg
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-15 19:11 by CgyNCDVNhGVg
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:59 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:13 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 09:38 by VMDqgnlKimtZ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:26 by C79jjtjAKjEn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:49 by 5BDX0d0WFp7t
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
信息提示
请填处理意见