| 查看: 3238 | 回复: 9 | ||||||
| 本帖产生 1 个 QC强帖 ，点击这里进行查看 | ||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||||
ggdh专家顾问 (正式写手)
|
[交流]
【ggdh个人文集】【benchmark】 ORCA2.8 vs gaussian09 ——opt 已有8人参与
|
|||||
|
cpu: Intel(R) Core(TM)2 CPU 6700 @ 2.66GHz 内存: 2G 分子:  任务: 优化结构 初始结构: gauss viewer5.0 画好后clean orca2.8输入文件： ! PAL2 opt B3LYP def2-SVP def2-SVP/J TightSCF RIJCOSX * xyzfile 0 1 cdgtest.xyz gaussian 09 输入文件： %chk=test.chk %nprocshared=2 %mem=100mw # b3lyp/6-31g(d) opt nosymm Title Card Required 0 1 。。。。。。 结果： orca： FINAL SINGLE POINT ENERGY -952.120481689 FINAL SINGLE POINT ENERGY -952.169767216 FINAL SINGLE POINT ENERGY -952.190790616 FINAL SINGLE POINT ENERGY -952.195998095 FINAL SINGLE POINT ENERGY -952.198594914 FINAL SINGLE POINT ENERGY -952.199839212 FINAL SINGLE POINT ENERGY -952.200335262 FINAL SINGLE POINT ENERGY -952.200665774 FINAL SINGLE POINT ENERGY -952.200855001 FINAL SINGLE POINT ENERGY -952.200928599 FINAL SINGLE POINT ENERGY -952.200933194 FINAL SINGLE POINT ENERGY -952.200937961 FINAL SINGLE POINT ENERGY -952.200942761 *** FINAL ENERGY EVALUATION AT THE STATIONARY POINT *** FINAL SINGLE POINT ENERGY -952.200942440 104.384 min （程序内显示时间和用time命令显示的时间基本一致） 两台机器并行：（4cpu，TCP/IP) 66min 去掉rijcosx （2cpu） 15步收敛 FINAL SINGLE POINT ENERGY -952.200296973 与rijcosx的能量差 0.405（kcal/mol) 162min 这样看来用orca 不用rij 基本没有什么意义。 改成b3lyp/G 6-31g(d) (客场作战）21步收敛 FINAL SINGLE POINT ENERGY -953.388962626 208min gaussian 09 SCF Done: E(RB3LYP) = -953.329057275 A.U. after 15 cycles SCF Done: E(RB3LYP) = -953.374871272 A.U. after 13 cycles SCF Done: E(RB3LYP) = -953.393322126 A.U. after 13 cycles SCF Done: E(RB3LYP) = -953.403290144 A.U. after 13 cycles SCF Done: E(RB3LYP) = -953.404332759 A.U. after 12 cycles SCF Done: E(RB3LYP) = -953.404813043 A.U. after 12 cycles SCF Done: E(RB3LYP) = -953.405060848 A.U. after 12 cycles SCF Done: E(RB3LYP) = -953.405100591 A.U. after 11 cycles SCF Done: E(RB3LYP) = -953.405147977 A.U. after 11 cycles SCF Done: E(RB3LYP) = -953.405174150 A.U. after 10 cycles SCF Done: E(RB3LYP) = -953.405199212 A.U. after 11 cycles SCF Done: E(RB3LYP) = -953.405211913 A.U. after 10 cycles SCF Done: E(RB3LYP) = -953.405233140 A.U. after 11 cycles SCF Done: E(RB3LYP) = -953.405250086 A.U. after 10 cycles SCF Done: E(RB3LYP) = -953.405263444 A.U. after 10 cycles SCF Done: E(RB3LYP) = -953.405267855 A.U. after 10 cycles SCF Done: E(RB3LYP) = -953.405269166 A.U. after 9 cycles SCF Done: E(RB3LYP) = -953.405269230 A.U. after 8 cycles 90 min gaussian03 SCF Done: E(RB+HF-LYP) = -953.329057270 A.U. after 15 cycles SCF Done: E(RB+HF-LYP) = -953.374871273 A.U. after 13 cycles SCF Done: E(RB+HF-LYP) = -953.394104107 A.U. after 13 cycles SCF Done: E(RB+HF-LYP) = -953.402606978 A.U. after 12 cycles SCF Done: E(RB+HF-LYP) = -953.404254818 A.U. after 12 cycles SCF Done: E(RB+HF-LYP) = -953.404749169 A.U. after 13 cycles SCF Done: E(RB+HF-LYP) = -953.405048977 A.U. after 12 cycles SCF Done: E(RB+HF-LYP) = -953.405110503 A.U. after 11 cycles SCF Done: E(RB+HF-LYP) = -953.405146670 A.U. after 11 cycles SCF Done: E(RB+HF-LYP) = -953.405170930 A.U. after 10 cycles SCF Done: E(RB+HF-LYP) = -953.405191140 A.U. after 10 cycles SCF Done: E(RB+HF-LYP) = -953.405209730 A.U. after 11 cycles SCF Done: E(RB+HF-LYP) = -953.405225304 A.U. after 11 cycles SCF Done: E(RB+HF-LYP) = -953.405237410 A.U. after 10 cycles SCF Done: E(RB+HF-LYP) = -953.405252850 A.U. after 10 cycles SCF Done: E(RB+HF-LYP) = -953.405261436 A.U. after 10 cycles SCF Done: E(RB+HF-LYP) = -953.405267561 A.U. after 9 cycles SCF Done: E(RB+HF-LYP) = -953.405268741 A.U. after 8 cycles SCF Done: E(RB+HF-LYP) = -953.405269102 A.U. after 8 cycles SCF Done: E(RB+HF-LYP) = -953.405269273 A.U. after 8 cycles SCF Done: E(RB+HF-LYP) = -953.405269280 A.U. after 7 cycles SCF Done: E(RB+HF-LYP) = -953.405269270 A.U. after 7 cycles SCF Done: E(RB+HF-LYP) = -953.405269276 A.U. after 7 cycles SCF Done: E(RB+HF-LYP) = -953.405269273 A.U. after 7 cycles 130min 基函数不一样，因为orca 推荐aldrich基组 而gaussian推荐pople机组，所以一般用orca优化就会选择def2-svp（特别是你要用rijcosx的话），而用gaussian优化就会选择6-31g（d）所以暂时没有统一基组。 gaussian： 322 basis functions, 604 primitive gaussians orca：# of primitive gaussian shells ... 295 # of primitive gaussian functions ... 557 谢谢beefly提醒，确实忘了对称性这码事了。重新用gaussian09计算。计算结果和没加nosymm一样，看来在优化过程中gaussian没有找到对称性 谢谢solarman的意见，soscf默认是打开的， SOSCF CNVSOSCF .... on Start iteration SOSCFMaxIt .... 150 Startup grad/error SOSCFStart .... 0.003300 我把soscfstart 调到0.01测试中。。。时间一样长。可能是基组小的原因，看不出差别。 scf的收敛指标如下： orca2.8 Convergence Tolerance: Convergence Check Mode ConvCheckMode .... Total+1el-Energy Energy Change TolE .... 1.000e-08 Eh 1-El. energy change .... 1.000e-05 Eh Orbital Gradient TolG .... 1.000e-05 Orbital Rotation angle TolX .... 1.000e-05 DIIS Error TolErr .... 5.000e-07 gaussian09 Requested convergence on RMS density matrix=1.00D-08 within 128 cycles. Requested convergence on MAX density matrix=1.00D-06. Requested convergence on energy=1.00D-06. gaussian09中虽然有densityfit，但是目前只能用于纯泛函。 似乎gaussian完胜啊。orca即使用了rijcosx也不及gaussian快，且orca的收敛精度比gaussian低，几何优化步数比gaussian少。不过这个分子比较小就是了。后续benchmark进行中。。。 [ Last edited by yjcmwgk on 2010-9-26 at 21:12 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 Gaussian |
» 猜你喜欢
Comsol 6.0 桶内肥料产热产湿过程模拟求助
已经有1人回复
-
分子动力学 Plumed教程求助
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有122人回复
-
研究生学术新锐奖｜宁波微纳材料研讨会参会通知
已经有2人回复
-
申请2027年材料/化学类博士
已经有6人回复
-
求相关合成条件
已经有14人回复
-
2027春季申博自荐
已经有6人回复
-
青基的材料青椒们，可以把硕士联培生送我这里，要研二男生
已经有68人回复
-
大湾区大学刘天辉课题组招收2027级博士生（与中山大学物理学院联合培养）
已经有74人回复
-
申博
已经有2人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 2012湿地农业可持续发展——洪山论坛——“湿地农业与涝渍灾害防治”国际研讨会
已经有19人回复
- J Opt. A:Pure Appl. Opt 已经停刊了吗？
已经有13人回复
- Gaussian计算得到的前线分子轨道图能说明什么？
已经有7人回复
- 从刀笔小吏到第一重臣——李鸿章全传【转载】
已经有215人回复
- 问个基础问题——树脂是热固性还是热塑性与是否结晶有关吗？
已经有8人回复
- Gaussian09中文用户参考手册（PDF，DOC）
已经有724人回复
- 【求助】同一个结构，G09和G03 opt计算的能量差别大吗？
已经有7人回复
lihb734
铁杆木虫 (职业作家)
站在计算化学入门的门槛上
- QC强帖: 2
- 应助: 132 (高中生)
- 金币: 7495.4
- 散金: 9378
- 红花: 74
- 帖子: 3107
- 在线: 1177.3小时
- 虫号: 448403
- 注册: 2007-11-01
- 性别: GG
- 专业: 理论和计算化学
3楼2010-09-26 08:52:05
ikea1984
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1074.8
- 散金: 132
- 红花: 2
- 帖子: 130
- 在线: 152.6小时
- 虫号: 1083533
- 注册: 2010-08-27
- 专业: 凝聚态物性 II ：电子结构
2楼2010-09-26 06:29:17
sculhf
荣誉版主 (知名作家)
断翼天使
-
专家经验: +1 - 应助: 243 (大学生)
- 贵宾: 4.018
- 金币: 43355.9
- 散金: 1638
- 红花: 64
- 沙发: 16
- 帖子: 7688
- 在线: 2008.8小时
- 虫号: 497659
- 注册: 2008-02-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 基金申请
4楼2010-09-26 09:23:17
shalene
至尊木虫 (著名写手)
- 应助: 4 (幼儿园)
- 金币: 10199.1
- 散金: 2
- 帖子: 2953
- 在线: 459.9小时
- 虫号: 3657
- 注册: 2003-03-18
- 性别: GG
- 专业: 化学动力学
5楼2010-09-26 09:32:29