24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1066  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

exabyss916

新虫 (小有名气)

[交流] 【求助】关于UGBS和ANO-RCC用在DKH计算的精度已有3人参与

我用两种方法分别计算了RbCs系统,由于计算量太大只计算了单点能量,但既然方法都一样只是Basis set不一样即使单点能量应该也可以直接比较。

两种方法能量相差几百个Hatree(计算的能量为10000个Hatree左右)。
请问到底哪种方法是合理的呢?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

exabyss916

新虫 (小有名气)

引用回帖:
Originally posted by beefly at 2010-05-19 13:42:44:
UHF/UGBS能量
Rb: -2978.86122825
Cs: -7780.59443740
总和:-10759.45566565,和你算的RHF能量差不多。

ANO-RCC主页上的CASPT2能量:
Rb: -2979.02
Cs: -7780.95

如果忽略相关能,应该和UGBS更接近。 ...

非常感谢详细的回复和相关数据,我也问了gaussian技术支持,他们也说用应该用非收缩基组计算,那我以后就都用UGBS系列计算。
7楼2010-05-19 14:30:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

exabyss916

新虫 (小有名气)

更多计算细节:

基组1:
# MP2/gen Integral=DKH
这里gen产生的是ano-rcc基组,结果如下:
1\1\GINC-NODE04\SP\RMP2-FC\Gen\Cs1Rb1\name\08-May-2010\0\\# MP2/gen
Integral=DKH\\\0,1\Rb\Cs,1,4.3\\Version=EM64L-G09RevA.01\State=1-SG\HF
=-10488.3784426\MP2=-10488.3888615\RMSD=7.597e-09\PG=C*V [C*(Rb1Cs1)]\

基组2:
# MP2/UGBS guess(indo) Integral=DKH
结果如下:
1\1\GINC-NODE05\SP\RMP2-FC\UGBS\Cs1Rb1\name\10-May-2010\0\\# MP2/UGB
S guess(indo) Integral=DKH\\\0,1\Rb\Cs,1,4.3\\Version=EM64L-G09RevA.01
\State=1-SG\HF=-10759.4475783\MP2=-10759.4590548\RMSD=2.655e-09\PG=C*V  [C*(Rb1Cs1)]\\@
2楼2010-05-18 15:04:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天空空

木虫 (正式写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
aylayl08(金币+3):谢谢指点 2010-05-18 16:29:46
一般而言,Primitive (uncontracted) basis set 是最精确的。 UGBS就是Primitive (uncontracted) basis set。但是你最好check一下,我的建议是你不要只比较能量,看看两者对其它性质的描述有区别没?
好好学习,天天向上!嘿嘿
3楼2010-05-18 16:01:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
aylayl08(金币+2):谢谢提示 2010-05-19 15:10:27
不可能差这么大,即使和最小基组比较,最多只能差(1 hartree*重原子数)。检查一下HF轨道占据是否一致。还可能发生了基组依赖
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
4楼2010-05-18 21:33:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见