| 查看: 2688 | 回复: 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 |
» 猜你喜欢
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有10人回复
-
夏威夷大学生物胶体与界面实验室诚招全奖博士生和博士后
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有261人回复
-
求助立方晶体In(OH)3的CIF卡片
已经有2人回复
-
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340
已经有2人回复
-
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料
已经有0人回复
-
七嗪类物质合成求助
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
-
[电子教材]天津大学物理化学（第七版）上、下册
已经有1人回复
-
阴离子交换膜电解二氧化碳还原
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 2012湿地农业可持续发展——洪山论坛——“湿地农业与涝渍灾害防治”国际研讨会
已经有19人回复
- J Opt. A:Pure Appl. Opt 已经停刊了吗？
已经有13人回复
- Gaussian计算得到的前线分子轨道图能说明什么？
已经有7人回复
- 从刀笔小吏到第一重臣——李鸿章全传【转载】
已经有215人回复
- 问个基础问题——树脂是热固性还是热塑性与是否结晶有关吗？
已经有8人回复
- Gaussian09中文用户参考手册（PDF，DOC）
已经有724人回复
- 【求助】同一个结构，G09和G03 opt计算的能量差别大吗？
已经有7人回复
shalene
至尊木虫 (著名写手)
- 应助: 4 (幼儿园)
- 金币: 10199.1
- 散金: 2
- 帖子: 2953
- 在线: 459.9小时
- 虫号: 3657
- 注册: 2003-03-18
- 性别: GG
- 专业: 化学动力学
5楼2010-09-26 09:32:29
ikea1984
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1074.8
- 散金: 132
- 红花: 2
- 帖子: 130
- 在线: 152.6小时
- 虫号: 1083533
- 注册: 2010-08-27
- 专业: 凝聚态物性 II ：电子结构
2楼2010-09-26 06:29:17
lihb734
铁杆木虫 (职业作家)
站在计算化学入门的门槛上
- QC强帖: 2
- 应助: 132 (高中生)
- 金币: 7495.4
- 散金: 9378
- 红花: 74
- 帖子: 3107
- 在线: 1177.3小时
- 虫号: 448403
- 注册: 2007-11-01
- 性别: GG
- 专业: 理论和计算化学
3楼2010-09-26 08:52:05
sculhf
荣誉版主 (知名作家)
断翼天使
-
专家经验: +1 - 应助: 243 (大学生)
- 贵宾: 4.018
- 金币: 43334.1
- 散金: 1574
- 红花: 64
- 沙发: 16
- 帖子: 7659
- 在线: 2006.6小时
- 虫号: 497659
- 注册: 2008-02-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 基金申请
4楼2010-09-26 09:23:17