版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

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

» 猜你喜欢

CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优已经有6人回复
德国Karlsruhe Institute of Technology招收电化学储能及联合培养CSC博士等已经有20人回复
物理化学论文润色/翻译怎么收费? 已经有285人回复
德国Karlsruhe Institute of Technology招收电化学储能及联合培养CSC博士已经有23人回复
德国Karlsruhe Institute of Technology招收电化学储能及联合培养CSC博士已经有1人回复
固态锂金属电池迁移数计算求助已经有1人回复
荷兰奈梅亨大学赵文博与Galimberti课题组招收理论计算CSC博士 2000 欧元/月+房补已经有0人回复
26储能博士申请自荐已经有22人回复
cif求助下载已经有1人回复
高效微纳米储能材料与产业化推进实验室苏州招收硕士研究生（含调剂）已经有2人回复
上海交通大学--宁波东方理工大学电池方向博士招生已经有10人回复

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

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

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

ikea1984

金虫 (小有名气)

应助: 1 (幼儿园)
金币: 1074.8
散金: 132
红花: 2
帖子: 130
在线: 152.6小时
虫号: 1083533
注册: 2010-08-27
专业: 凝聚态物性 II ：电子结构

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
小红豆(金币+2):感谢交流 2010-09-26 08:45:45

基组可能太小。一般来说RI技术在原基组较大时占优势，例如使用TZ或者QZ基组。

赞一下(2人)

回复此楼

2楼2010-09-26 06:29:17

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

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上

QC强帖: 2
应助: 132 (高中生)
金币: 7495.4
散金: 9378
红花: 74
帖子: 3107
在线: 1177.3小时
虫号: 448403
注册: 2007-11-01
性别: GG
专业: 理论和计算化学

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
小红豆(金币+2):辛苦了 2010-09-26 08:54:39

有两点一问，一时ORCA加了这个TightSCF ，不知道自洽场收敛是否一致？
二，基组不一致，不知道楼主是否比较了基函数的数目？
其实楼主可以去基组库找到def2-SVP之后，在高斯中也采用这个基组。
还有，二者的拟合基也不一样的。

赞一下(2人)

回复此楼

前途光明，出路难觅！

3楼2010-09-26 08:52:05

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

sculhf

荣誉版主 (知名作家)

断翼天使

专家经验: +1
应助: 243 (大学生)
贵宾: 4.018
金币: 43334.1
散金: 1574
红花: 64
沙发: 16
帖子: 7659
在线: 2006.6小时
虫号: 497659
注册: 2008-02-02
性别: GG
专业: 凝聚态物性 II ：电子结构
管辖: 基金申请

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
小红豆(金币+2):辛苦了 2010-09-26 09:32:05

选择的基组大小不同
计算结果的可比性有待提高

赞一下(2人)

回复此楼

谋事在人成事在天

4楼2010-09-26 09:23:17

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

shalene

至尊木虫 (著名写手)

应助: 4 (幼儿园)
金币: 10199.1
散金: 2
帖子: 2953
在线: 459.9小时
虫号: 3657
注册: 2003-03-18
性别: GG
专业: 化学动力学

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

楼主，把你的初始结构发出来，我也来算算。

赞一下(1人)

回复此楼

5楼2010-09-26 09:32:29

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

beefly

专家顾问 (职业作家)

地沟油冶炼专家

专家经验: +458
QC强帖: 12
应助: 88 (初中生)
金币: 12931.9
散金: 27307
红花: 300
帖子: 4670
在线: 926.5小时
虫号: 408372
注册: 2007-06-21
性别: GG
专业: 理论和计算化学
管辖: 计算模拟

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

lz忽视了一个至关重要的因素：对称性。

orca目前并未充分利用对称性，充其量也就是分析轨道、指定轨道的时候方便一点，计算过程并不节省时间。如果和Gaussian做比较，应该在Gaussian输入文件里指定Nosymm。

赞一下(1人)

回复此楼

beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地

6楼2010-09-26 11:18:11

已阅回复此楼关注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 ]

赞一下(1人)

回复此楼

7楼2010-09-26 11:50:35

已阅回复此楼关注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 ]

赞一下(1人)

回复此楼

未知必须借由更深的未知去获悉，隐晦必须借由更深的隐晦来明了

8楼2010-09-26 12:30:27

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