24小时热门版块排行榜    

查看: 2857  |  回复: 16

ghcacj

荣誉版主 (著名写手)

阿超

优秀版主

[求助] 为何用同样的方法,软件之间差异这么大?

用高斯和Orca ,均采用DFT B3LYP 6-31G* + Lanl2dz计算一个片段的chelpg电荷,为何结果差别这么大?尤其是Zn的电荷,我后来用gamess也算过,和高斯也不一样,主要是体现在Zn上,如果说是ecp的差别,那么不采用ecp,还是这样,这是为何?请指教,附上输入输出文件。
=======================
前两个附件分别是Gaussian的输入输出文件,后两个是Orca的。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

coolrainbow

木虫 (著名写手)

未来国家冻凉

【答案】应助回帖


ghcacj(金币+4): 谢谢 2011-05-18 19:02:15
gmy1990(金币+1): 非常感谢你对论坛的贡献,奖励一下! 2011-05-18 21:37:42
ghcacj(金币+5): 谢谢 2011-05-20 10:13:36
所谓的B3LYP,在不同软件中的定义是不一样的

要复现结果,最好用些单纯的泛函,或者手动设置泛函,这样结果才有可比性
技术博客:http://hi.baidu.com/coolrainbow/blog
2楼2011-05-18 18:41:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chrinide

木虫 (正式写手)

【答案】应助回帖


ghcacj(金币+5): 谢谢 2011-05-18 19:02:29
gmy1990(金币+1): 非常感谢你对论坛的贡献,奖励一下! 2011-05-18 21:37:48
ghcacj(金币+6): 谢谢 2011-05-20 10:12:00
首先,ORCA的6-31G*如果直接定义的话,与Gaussian有很大的不同。
二,ORCA只使用5D 7F 球谐基组, Gaussian 在默认情况下使用笛卡尔基。
三,在ORCA中正确的使用6-31G*要用 %basis 6_31G_d end。
四,ORCA的COSMO与Gaussian 的执行不大相同。

个人浅见,仅供参考
未知必须借由更深的未知去获悉,隐晦必须借由更深的隐晦来明了
3楼2011-05-18 19:00:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghcacj

荣誉版主 (著名写手)

阿超

优秀版主

引用回帖:
Originally posted by chrinide at 2011-05-18 19:00:08:
首先,ORCA的6-31G*如果直接定义的话,与Gaussian有很大的不同。
二,ORCA只使用5D 7F 球谐基组, Gaussian 在默认情况下使用笛卡尔基。
三,在ORCA中正确的使用6-31G*要用 %basis 6_31G_d end。
四,ORCA的CO ...

那么我应该如果使得两者尽量一致呢?
差别别那么大就好。
4楼2011-05-18 19:03:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chrinide

木虫 (正式写手)

【答案】应助回帖

ghcacj(金币+5): 谢谢 2011-05-18 19:17:22
ghcacj(金币+5): 谢谢 2011-05-20 10:13:58
引用回帖:
Originally posted by ghcacj at 2011-05-18 19:03:15:
那么我应该如果使得两者尽量一致呢?
差别别那么大就好。

基组 按照%basis 6_31G_d end 设置结果应该会有改进,但是应该不会太大。

至于COSMO的执行,ORCA将在下一个版本进行升级,到时候应该和Gaussian相同或者更好...

当然,结果到底咋样,要算了后才知道的,哈哈
未知必须借由更深的未知去获悉,隐晦必须借由更深的隐晦来明了
5楼2011-05-18 19:08:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghcacj

荣誉版主 (著名写手)

阿超

优秀版主

引用回帖:
Originally posted by chrinide at 2011-05-18 19:08:52:
基组 按照%basis 6_31G_d end 设置结果应该会有改进,但是应该不会太大。

至于COSMO的执行,ORCA将在下一个版本进行升级,到时候应该和Gaussian相同或者更好...

当然,结果到底咋样,要算了后才知道的, ...

非常感谢,我不是量化专业选手,没想到不同的软件之间水这么深。
那么我还想问的是,如果结果和高斯不一致,能说我的是错的吗?
6楼2011-05-18 19:13:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

coolrainbow

木虫 (著名写手)

未来国家冻凉

引用回帖:
Originally posted by ghcacj at 2011-05-18 19:13:07:
非常感谢,我不是量化专业选手,没想到不同的软件之间水这么深。
那么我还想问的是,如果结果和高斯不一致,能说我的是错的吗?

凭什么非得Gaussian才是对的?
技术博客:http://hi.baidu.com/coolrainbow/blog
7楼2011-05-18 19:15:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihb734

铁杆木虫 (职业作家)

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


gmy1990(金币+1): 非常感谢你对论坛的贡献,奖励一下! 2011-05-18 21:37:59
引用回帖:
Originally posted by chrinide at 2011-05-18 19:00:08:
首先,ORCA的6-31G*如果直接定义的话,与Gaussian有很大的不同。
二,ORCA只使用5D 7F 球谐基组, Gaussian 在默认情况下使用笛卡尔基。
三,在ORCA中正确的使用6-31G*要用 %basis 6_31G_d end。
四,ORCA的CO ...

您很厉害的,
COSMO在高斯中没有执行,只能说与CPCM比较类似。
Neese是反对Gaussian的,应该不会执行意大利人的CPCM。只会去优化执行Klaut的COSMO
前途光明,出路难觅!
8楼2011-05-18 19:18:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghcacj

荣誉版主 (著名写手)

阿超

优秀版主

引用回帖:
Originally posted by lihb734 at 2011-05-18 19:18:59:
您很厉害的,
COSMO在高斯中没有执行,只能说与CPCM比较类似。
Neese是反对Gaussian的,应该不会执行意大利人的CPCM。只会去优化执行Klaut的COSMO

Neese反对Gaussian我能看出来,因为早期Orca是可以通过外部接口调用Gaussian中的函数,后来Gaussian公司把这个给掐了,不反对就怪了,手册中原话在这:
“Earlier versions of ORCA have supported the Gau_External option of the Gaussian optimizer. Since the ORCA team got banned by Gaussian in January 2007 we can no longer support this option. It has not been removed from the code and may or may not work. If there is trouble with it we can–unfortunately–not offer any help since we do not have access to the Gaussian code any longer. If you feel uncomfortable with this situation you have to complain to Gaussian Inc.!”
9楼2011-05-18 20:32:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihb734

铁杆木虫 (职业作家)

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

引用回帖:
Originally posted by ghcacj at 2011-05-18 20:32:53:
Neese反对Gaussian我能看出来,因为早期Orca是可以通过外部接口调用Gaussian中的函数,后来Gaussian公司把这个给掐了,不反对就怪了,手册中原话在这:
“Earlier versions of ORCA have supported the Gau_E ...

所以为了比较的目的,最好都在气相中进行,排除溶剂化的干扰,再比较一下。
前途光明,出路难觅!
10楼2011-05-19 07:47:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ghcacj 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +10 yufeiwaner 2026-08-09 10/500 2026-08-10 13:44 by cdutwml
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +10 archvillain 2026-08-04 12/600 2026-08-10 12:45 by yudaoqian88
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +5 家与远方 2026-08-10 8/400 2026-08-10 12:45 by 家与远方
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 关于代码变化问题,想知道的进来 +16 且听虎啸 2026-08-07 23/1150 2026-08-10 12:35 by majunge000
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 5/250 2026-08-08 21:54 by 13586093586
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见