| 查看: 1215 | 回复: 6 | ||||
[交流]
【求助】关于UGBS和ANO-RCC用在DKH计算的精度 已有3人参与
|
|
我用两种方法分别计算了RbCs系统，由于计算量太大只计算了单点能量，但既然方法都一样只是Basis set不一样即使单点能量应该也可以直接比较。 两种方法能量相差几百个Hatree（计算的能量为10000个Hatree左右）。 请问到底哪种方法是合理的呢？ |
回复此楼» 猜你喜欢
阴离子交换膜电解二氧化碳还原
已经有0人回复
-
分子动力学模拟合作需求
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有124人回复
-
九江学院2026年最新高层次人才招聘公告
已经有0人回复
-
深圳信息职业技术学院2026年教师招聘公告（最新）
已经有0人回复
-
求ENDNOTE破解软件安装包
已经有4人回复
-
研究生招生
已经有0人回复
-
河海大学 （211.双一流高校）钙钛矿光电实验室丁勇团队欢迎有志青年！
已经有3人回复
-
黑龙江省原子能研究院化学课题组招收2026级化学专业或化学相关专业硕士生
已经有0人回复
-
山西大学资源与环境工程研究所化工分离/环境功能材料小组招收2026年博士生
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 使用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
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13001.7
- 散金: 27307
- 红花: 300
- 帖子: 4672
- 在线: 926.9小时
- 虫号: 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
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13001.7
- 散金: 27307
- 红花: 300
- 帖子: 4672
- 在线: 926.9小时
- 虫号: 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
