| 查看: 1067 | 回复: 6 | ||||
[交流]
【求助】关于UGBS和ANO-RCC用在DKH计算的精度已有3人参与
|
|
我用两种方法分别计算了RbCs系统，由于计算量太大只计算了单点能量，但既然方法都一样只是Basis set不一样即使单点能量应该也可以直接比较。 两种方法能量相差几百个Hatree（计算的能量为10000个Hatree左右）。 请问到底哪种方法是合理的呢？ |
回复此楼» 猜你喜欢
超导转变温度NaN
已经有4人回复
-
华南师范大学先进光电子研究院电子纸团队诚聘2026级博士研究生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有55人回复
-
QE利用声子计算的一系列lambda和Tc如何选
已经有8人回复
-
理论计算合作
已经有2人回复
-
小木虫的论文辅导靠谱吗？有没有用过的同学？
已经有1人回复
-
QE计算声子ph.out的Γ点出现虚频（-30cm-1）vasp计算没有
已经有0人回复
-
40-50万/年，中国散裂中子源诚聘计算模拟方向博士后
已经有81人回复
-
2026博士申请—化学能源方向，本人熟练分子动力学模拟的使用
已经有4人回复
-
书籍求助：汽车市场营销理论与实务（电子版）——章小平
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
使用Multiwfn绘制原子轨道图形、研究原子壳层结构及相对论效应的影响
已经有17人回复
|
更多计算细节： 基组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
3楼2010-05-18 16:01:41
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +457 - QC强帖: 12
- 应助: 86 (初中生)
- 金币: 12717.9
- 散金: 27307
- 红花: 300
- 帖子: 4663
- 在线: 925.3小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
4楼2010-05-18 21:33:09
|
确实不太一样，轨道能级也不一样： 基组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
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +457 - QC强帖: 12
- 应助: 86 (初中生)
- 金币: 12717.9
- 散金: 27307
- 红花: 300
- 帖子: 4663
- 在线: 925.3小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+5):谢谢提示 2010-05-19 15:10:37
小木虫(金币+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在做积分时可能有问题。 |
6楼2010-05-19 13:42:44
7楼2010-05-19 14:30:51