应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】关于UGBS和ANO-RCC用在DKH计算的精度
71/1返回列表
查看: 1422  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

exabyss916

新虫 (小有名气)

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

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

两种方法能量相差几百个Hatree（计算的能量为10000个Hatree左右）。
请问到底哪种方法是合理的呢？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-05-17 12:08:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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一下，我的建议是你不要只比较能量，看看两者对其它性质的描述有区别没？
赞一下(2人) 回复此楼
好好学习，天天向上！嘿嘿
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轨道占据是否一致。还可能发生了基组依赖
赞一下(2人) 回复此楼
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
4楼2010-05-18 21:33:09
已阅   回复此楼   关注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的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+5):谢谢提示 2010-05-19 15:10:37
UHF/UGBS能量
Rb: -2978.86122825
Cs: -7780.59443740
总和：-10759.45566565，和你算的RHF能量差不多。

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

如果忽略相关能，应该和UGBS更接近。所以你的ANO-RCC结果是错的。

Gaussian使用一般收缩基组时效率非常低，例如Rb作ANO收缩和不收缩，UHF计算时间之比大约20:1，而molpro、molcas等程序为1:3。ANO-RCC改用非收缩形式，也许能得到和UGBS接近的结果。

还有个可能存在的问题。Cs ANO-RCC的s高斯函数指数远远大于UGBS，Gaussian在做积分时可能有问题。
赞一下(2人) 回复此楼
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
6楼2010-05-19 13:42:44
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 exabyss916 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然函评费到账 +16 羊腰板 2026-08-21 17/850 2026-08-23 23:02 by tianxiaochun
[基金申请] 让我中一个面上吧！ +12 大萍1987 2026-08-20 15/750 2026-08-23 20:44 by 大萍1987
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[基金申请] 什么时候开奖？ +9 CrisMessi 2026-08-18 10/500 2026-08-23 12:03 by 丶昵称占用
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 范进中举一文的中心思想 +3 炎黄贵胄 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 时间戳变了，能看出什么问题？ +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭