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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】关于UGBS和ANO-RCC用在DKH计算的精度
71/1返回列表
查看: 1104  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭