24小时热门版块排行榜

返回列表

本帖产生 1 个 QC强帖，点击这里进行查看

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

ggdh

专家顾问 (正式写手)

专家经验: +69
QC强帖: 2
应助: 52 (初中生)
金币: 1747.4
散金: 10
红花: 18
帖子: 697
在线: 206.1小时
虫号: 96968
注册: 2005-11-08
性别: GG
专业: 化学生物学与生物有机化学
管辖: 量子化学

[交流] 【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

» 猜你喜欢

求助立方晶体In(OH)3的CIF卡片已经有2人回复
大湾区大学刘天辉课题组招收2026级博士生（与中山大学物理学院联合培养）已经有22人回复
物理化学论文润色/翻译怎么收费? 已经有276人回复
PbS量子点紫外吸收已经有0人回复
比利时鲁汶大学与国家留学基金委共同资助博士研究生CSC-KU Leuven PhD Scholarship 已经有0人回复
请问四氢呋喃溶解的聚合物用甲醇沉淀时，如何使沉淀过程加速？已经有2人回复
[电子教材]天津大学物理化学（第七版）上、下册已经有1人回复
阴离子交换膜电解二氧化碳还原已经有0人回复
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生已经有11人回复

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:

2012湿地农业可持续发展——洪山论坛——“湿地农业与涝渍灾害防治”国际研讨会已经有19人回复
J Opt. A:Pure Appl. Opt 已经停刊了吗？已经有13人回复
Gaussian计算得到的前线分子轨道图能说明什么？已经有7人回复
从刀笔小吏到第一重臣——李鸿章全传【转载】已经有215人回复
问个基础问题——树脂是热固性还是热塑性与是否结晶有关吗？已经有8人回复
Gaussian09中文用户参考手册（PDF，DOC）已经有724人回复
【求助】同一个结构，G09和G03 opt计算的能量差别大吗？已经有7人回复

1楼 2010-09-26 03:54:36

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

solarman

新虫 (正式写手)

应助: 5 (幼儿园)
金币: 1210.9
散金: 492
红花: 2
帖子: 458
在线: 210.5小时
虫号: 599917
注册: 2008-09-11
专业: 凝聚态物性 II ：电子结构

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

你不妨试一个大的体系，基函数超过2000的，这时RIJCOSX的优势会体现出来。
要确保基组，SCF和OPT的收敛精度，对称性均一致
此外，你在orca中可以用一下soscf。高斯虽没有rijcosx，但里面的densityfit也相当于用了RI
这种情况下，我还是相信orca会完胜

[ Last edited by solarman on 2010-9-26 at 11:51 ]