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

24小时热门版块排行榜    

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

天空空

木虫 (正式写手)
★ ★ ★ ★
小木虫(金币+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的回帖
查看全部 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的回帖

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的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 初试 317 +4 半拉月丙 2026-03-20 4/200 2026-03-21 16:41 by macy2011
[考研] 297求调剂 +11 戏精丹丹丹 2026-03-17 11/550 2026-03-21 15:43 by ColorlessPI
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +5 嘉年新程 2026-03-15 5/250 2026-03-21 15:20 by 巴黎材料
[考研] 306求调剂 +4 chuanzhu川烛 2026-03-18 4/200 2026-03-21 08:25 by laoshidan
[考研] 316求调剂 +6 梁茜雯 2026-03-19 6/300 2026-03-21 06:32 by Ecowxq666！
[考研] 070300化学319求调剂 +7 锦鲤0909 2026-03-17 7/350 2026-03-21 03:46 by JourneyLucky
[考研] 一志愿天津大学化学工艺专业（081702）315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 一志愿重庆大学085700资源与环境专硕，总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 一志愿西南交大，求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 323求调剂 +3 洼小桶 2026-03-18 3/150 2026-03-20 22:54 by JourneyLucky
[考研] 一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 350求调剂 +5 weudhdk 2026-03-19 5/250 2026-03-20 22:04 by luoyongfeng
[考研] 本人考085602 化学工程 专硕 +19 不知道叫什么！ 2026-03-15 21/1050 2026-03-20 20:48 by zhukairuo
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 一志愿985，本科211，0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 085601专硕，总分342求调剂，地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 0854可跨调剂，一作一项核心论文五项专利，省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
信息提示
关闭
请填处理意见
关闭