| 查看: 3236 | 回复: 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 |
» 猜你喜欢
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
Postdoc position, materials simulations/machine learning, U. of South Carolina
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有169人回复
-
祈福青基
已经有34人回复
-
牛来！米来！面来！
已经有8人回复
-
2027申博咨询
已经有11人回复
-
寻求CrC2Al的cif文件，感谢
已经有0人回复
-
ORCA 纯新手向入门指南0.7版
已经有0人回复
-
哈尔滨工业大学（深圳）绿色合成化学团队招收27年春季或秋季入学的博士生
已经有5人回复
-
【接测试】TPD-TPR、TPO、脉冲吸附、原位红外，电化学原位红外，站内联系
已经有0人回复
-
吉林大学化学学院理论化学研究所刘靖尧教授课题组招收2027级申请考核制博士研究生
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 2012湿地农业可持续发展——洪山论坛——“湿地农业与涝渍灾害防治”国际研讨会
已经有19人回复
- J Opt. A:Pure Appl. Opt 已经停刊了吗？
已经有13人回复
- Gaussian计算得到的前线分子轨道图能说明什么？
已经有7人回复
- 从刀笔小吏到第一重臣——李鸿章全传【转载】
已经有215人回复
- 问个基础问题——树脂是热固性还是热塑性与是否结晶有关吗？
已经有8人回复
- Gaussian09中文用户参考手册（PDF，DOC）
已经有724人回复
- 【求助】同一个结构，G09和G03 opt计算的能量差别大吗？
已经有7人回复
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
- 金币: 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
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +459 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13403
- 散金: 27307
- 红花: 302
- 帖子: 4686
- 在线: 929小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
6楼2010-09-26 11:18:11
solarman
新虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 1210.9
- 散金: 492
- 红花: 2
- 帖子: 458
- 在线: 210.5小时
- 虫号: 599917
- 注册: 2008-09-11
- 专业: 凝聚态物性 II ：电子结构
7楼2010-09-26 11:50:35
chrinide
木虫 (正式写手)
- QC强帖: 1
- 应助: 20 (小学生)
- 金币: 2460.3
- 红花: 5
- 帖子: 431
- 在线: 541.4小时
- 虫号: 488931
- 注册: 2007-12-31
- 专业: 宇宙学
8楼2010-09-26 12:30:27
solarman
新虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 1210.9
- 散金: 492
- 红花: 2
- 帖子: 458
- 在线: 210.5小时
- 虫号: 599917
- 注册: 2008-09-11
- 专业: 凝聚态物性 II ：电子结构
9楼2010-09-26 13:15:02
ghcacj
荣誉版主 (著名写手)
阿超
- QC强帖: 2
- 应助: 26 (小学生)
- 贵宾: 2.764
- 金币: 2543.5
- 散金: 8899
- 红花: 35
- 沙发: 3
- 帖子: 1376
- 在线: 1199小时
- 虫号: 513281
- 注册: 2008-02-27
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
10楼2010-09-29 15:34:23