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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】关于UGBS和ANO-RCC用在DKH计算的精度
51/1返回列表
查看: 1557  |  回复: 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

新虫 (小有名气)
引用回帖:
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一下，我的建议是你不要只比较能量，看看两者对其它性质的描述有区别没？
赞一下(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的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 X5bHxOULxFqb 2026-10-30 26/1300 2026-11-01 10:22 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 07:01 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 07:01 by YXM1FDCSE8Bn
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 06:53 by YXM1FDCSE8Bn
[基金申请] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[基金申请] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:43 by YXM1FDCSE8Bn
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 05:41 by YXM1FDCSE8Bn
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 04:04 by YXM1FDCSE8Bn
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
[找工作] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 03:56 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 10/500 2026-11-01 03:51 by YXM1FDCSE8Bn
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:25 by YXM1FDCSE8Bn
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-11-01 00:41 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 X5bHxOULxFqb 2026-10-30 12/600 2026-11-01 00:15 by YXM1FDCSE8Bn
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 X5bHxOULxFqb 2026-10-30 13/650 2026-10-31 17:12 by hKe7U1kad5Sg
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +8 zuesli 2026-10-27 14/700 2026-10-30 22:44 by 好运连连公主
[基金申请] 论文的挂名怎么评判贡献？基金申请中，排名第一的挂名文章 +3 biily 2026-10-25 3/150 2026-10-30 18:33 by 旁观2020
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
信息提示
关闭
请填处理意见
关闭