24小时热门版块排行榜    

查看: 2254  |  回复: 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 的主题更新
信息提示
请填处理意见