24小时热门版块排行榜

>论坛更新日志 (2920)
>虫友互识 (231)
>导师招生 (172)
>文献求助 (171)
>论文投稿 (68)
>休闲灌水 (52)
>硕博家园 (50)
>公派出国 (48)
>考博 (40)
>招聘信息布告栏 (32)
>论文道贺祈福 (27)
>催化 (25)
>教师之家 (24)
>找工作 (24)
>考研 (22)
>博后之家 (21)

返回列表

本帖产生 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

» 猜你喜欢

香港理工大学计算材料与物理课题组招聘博后已经有15人回复
西班牙Institut Català de Nanociència i Nanotecnologia, CSIC-BIST招收公派博士生已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有152人回复
英国利物浦大学屠昕和雍雪课题组招收CSC博士生已经有27人回复
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生（还有1个名额）已经有0人回复
华北电力大学（北京）（第一性原理计算）博士招生——学博，专博各1人已经有45人回复
英国贝尔法斯特女王大学招收 2026 年 CSC/BBSRC phd 机器学习催化剂设计方向已经有0人回复
法国里昂大学招生CSC 资助博士生已经有0人回复
法国里昂大学接收CSC 资助访问学者已经有0人回复
北京工业大学材料科学与工程学院《院士团队》招收工程专业博士研究生（计算优先）已经有1人回复

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

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

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的回帖

chrinide

木虫 (正式写手)

QC强帖: 1
应助: 20 (小学生)
金币: 2460.3
红花: 5
帖子: 431
在线: 541.4小时
虫号: 488931
注册: 2007-12-31
专业: 宇宙学

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

引用回帖:

Originally posted by solarman at 2010-09-26 11:50:35:
你不妨试一个大的体系，基函数超过2000的，这时RIJCOSX的优势会体现出来。
要确保基组，SCF和OPT的收敛精度，对称性均一致
此外，你在orca中可以用一下soscf。高斯虽没有rijcosx，但里面的densityfit也相当于用 ...

这位是明白人啊 O(∩_∩)O哈哈哈~

[ Last edited by chrinide on 2010-9-26 at 12:42 ]