24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1068  |  回复: 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-18 21:33:09:
不可能差这么大,即使和最小基组比较,最多只能差(1 hartree*重原子数)。检查一下HF轨道占据是否一致。还可能发生了基组依赖

确实不太一样,轨道能级也不一样:
基组1:
       Occupied  (SG) (SG) (SG) (SG) (PI) (PI) (SG) (SG) (SG) (PI)
                 (PI) (SG) (PI) (PI) (SG) (PI) (PI) (DLTA) (DLTA)
                 (SG) (SG) (SG) (PI) (PI) (DLTA) (DLTA) (PI) (PI)
                 (SG) (SG) (PI) (PI) (DLTA) (DLTA) (PI) (PI) (SG)
                 (SG) (SG) (PI) (PI) (SG) (SG) (PI) (PI) (SG)
Alpha  occ. eigenvalues -- **********-589.04971-212.16515-190.33495-190.33491
Alpha  occ. eigenvalues -- -190.33491 -89.74429 -45.77127 -38.89039 -38.87668
Alpha  occ. eigenvalues --  -38.87668 -38.54419 -38.54412 -38.54412 -27.83521
Alpha  occ. eigenvalues --  -27.83519 -27.83519 -27.83516 -27.83516 -15.37358
Alpha  occ. eigenvalues --   -9.34989  -7.61111  -7.61066  -7.61066  -7.60985
Alpha  occ. eigenvalues --   -7.60985  -6.93297  -6.93297  -6.93286  -5.70744
Alpha  occ. eigenvalues --   -5.70351  -5.70351  -3.27885  -3.27885  -3.27865
Alpha  occ. eigenvalues --   -3.27865  -3.27861  -2.11817  -1.33481  -0.70938
Alpha  occ. eigenvalues --   -0.70938  -0.70614  -0.49512  -0.49075  -0.49075
Alpha  occ. eigenvalues --   -0.15553


基组2:
       Occupied  (SG) (SG) (SG) (SG) (PI) (PI) (SG) (SG) (PI) (PI)
                 (SG) (SG) (PI) (PI) (SG) (PI) (PI) (DLTA) (DLTA)
                 (SG) (SG) (PI) (PI) (SG) (SG) (PI) (PI) (SG) (PI)
                 (PI) (DLTA) (DLTA) (SG) (PI) (PI) (DLTA) (DLTA)
                 (SG) (SG) (SG) (PI) (PI) (SG) (PI) (PI) (SG)
Alpha  occ. eigenvalues -- **********-561.64098-212.13438-190.30259-190.30255
Alpha  occ. eigenvalues -- -190.30255 -77.44271 -68.46967 -68.46962 -68.46962
Alpha  occ. eigenvalues --  -45.73965 -38.51214 -38.51201 -38.51201 -27.80330
Alpha  occ. eigenvalues --  -27.80327 -27.80327 -27.80318 -27.80318 -12.54071
Alpha  occ. eigenvalues --   -9.56621  -9.56605  -9.56605  -9.31841  -6.90121
Alpha  occ. eigenvalues --   -6.90103  -6.90103  -4.62803  -4.62798  -4.62798
Alpha  occ. eigenvalues --   -4.62785  -4.62785  -3.24712  -3.24705  -3.24705
Alpha  occ. eigenvalues --   -3.24690  -3.24690  -1.55366  -1.30469  -0.79716
Alpha  occ. eigenvalues --   -0.79562  -0.79562  -0.67816  -0.67732  -0.67732
Alpha  occ. eigenvalues --   -0.11975
5楼2010-05-19 08:44:17
已阅   回复此楼   关注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的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见