计算化学公社

 找回密码 Forget password
 注册 Register
Views: 22987|回复 Reply: 17
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 用ONIOM进行QMMM计算时对QM区添加dispersion出错

[复制链接 Copy URL]

3

帖子

0

威望

41

eV
积分
44

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
我试图用高斯的oniom优化一个金属配合物和DNA的相互作用。对于QM区域,因为有弱相互作用,所以希望加上dispersion关键字,但是这样一来就会报l301的错。
输入信息如下:
%chk=test2.chk
%nproc=12
%mem=12GB
#p oniom(PBE1PBE/genecp empiricaldispersion=gd3:amber=Hardfirst) geom=connectivity opt SCRF=(solvent=water)

Title Card Required

-4 1 2 1 2 1
  C-CT--0.006900(PDBName=C5S,ResName=DG,ResNum=4)                  0    8.86045015    4.82512912    7.03733107 L
  O-OH--0.631800(PDBName=O5S,ResName=DG,ResNum=4)                  0    7.77345015    5.60612912    6.55133107 L
  C-CT-0.162900(PDBName=C4S,ResName=DG,ResNum=4)                   0    8.91445015    3.49512912    6.31933107 L
  O-OS--0.369100(PDBName=O4S,ResName=DG,ResNum=4)                  0    7.77045015    2.66312912    6.66633107 L
  C-CT-0.071300(PDBName=C3S,ResName=DG,ResNum=4)                   0    8.91945015    3.57812912    4.80033107 L
  O-OS--0.523200(PDBName=O3S,ResName=DG,ResNum=4)                  0    9.86345015    2.61412912    4.36033107 L
  C-CT--0.085400(PDBName=C2S,ResName=DG,ResNum=4)                  0    7.49245015    3.21312912    4.43933107 L
  C-CT-0.035800(PDBName=C1S,ResName=DG,ResNum=4)                   0    7.12645015    2.20912912    5.51233107 L
  N-NA--0.505300(PDBName=N1,ResName=DG,ResNum=4)                   0    3.19645015   -0.72987088    6.79333107 L
  C-CA-0.743200(PDBName=C2,ResName=DG,ResNum=4)                    0    4.47945015   -1.08287088    6.48733107 L
  N-N2--0.923000(PDBName=N2,ResName=DG,ResNum=4)                   0    4.73945015   -2.39187088    6.54333107 L
  N-NC--0.663600(PDBName=N3,ResName=DG,ResNum=4)                   0    5.42945015   -0.23887088    6.15933107 L
  C-CB-0.181400(PDBName=C4,ResName=DG,ResNum=4)                    0    4.98745015    1.03712912    6.15933107 L

  ...
  ...
  ...
  H-HC(PDBName=H21B,ResName=UNK,ResNum=1)                  0    7.56500000   -4.28800000    0.41700000 H
  H-HC(PDBName=H21C,ResName=UNK,ResNum=1)                  0    2.82600000   -7.46300000    1.43200000 H
  H-HC(PDBName=H21C,ResName=UNK,ResNum=1)                  0    6.08800000   -4.36600000   -0.10500000 H

1 3 1.0 2 1.0 29 1.0 30 1.0
2 31 1.0
3 4 1.0 5 1.0 28 1.0
  ...
  ...
  ...

340
341
342


HrmStr1 Pt Pt  100.0 1.00
HrmStr1 Pt CB  100.0 1.00
HrmStr1 Pt NC  100.0 1.00
HrmStr1 Pt CY  100.0 1.00
HrmStr1 CB CM  100.0 1.00
HrmStr1 NC CY  100.0 1.00
HrmStr1 NC CT  100.0 1.00
HrmBnd1 Pt Pt CB  100.0 180.00
HrmBnd1 Pt Pt NC  100.0 180.00
HrmBnd1 Pt Pt CY  100.0 180.00
HrmBnd1 Pt CB CM  100.0 180.00
HrmBnd1 Pt CB CB  100.0 180.00
HrmBnd1 Pt NC CB  100.0 180.00
HrmBnd1 Pt NC CA  100.0 180.00
HrmBnd1 Pt CY NC  100.0 180.00
HrmBnd1 CB Pt NC  100.0 180.00
HrmBnd1 CB Pt CY  100.0 180.00
HrmBnd1 CB CM CM  100.0 180.00
HrmBnd1 CB CM HA  100.0 180.00
HrmBnd1 CB CB CM  100.0 180.00
HrmBnd1 CB CB CB  100.0 180.00
HrmBnd1 NC Pt NC  100.0 180.00
HrmBnd1 NC Pt CY  100.0 180.00
HrmBnd1 NC Pt CM  100.0 180.00
HrmBnd1 NC CB CM  100.0 180.00
HrmBnd1 NC CA CA  100.0 180.00
HrmBnd1 NC CA H4  100.0 180.00
HrmBnd1 CM CM CM  100.0 180.00
HrmBnd1 NC CT CT  100.0 180.00
HrmBnd1 CB NC CB  100.0 180.00
HrmBnd1 CM CA HA  100.0 180.00
HrmBnd1 CM CA CA  100.0 180.00
HrmBnd1 CY NC CT  100.0 180.00

C H N -O 0
6-31G*
****
Pt 0
lanL2DZ
****

Pt 0
lanL2DZ

由于体系太大,输入文件中省略了关于坐标和成键的一部分信息。

报错信息如下:
Integral buffers will be    131072 words long.
Raffenetti 1 integral format.
Two-electron integral symmetry is turned on.
R6DS8: Unable to choose the S8 parameter, IExCor=    0 IXCFnc=  0 ScaHFX=  1.000000 IDFTD=3
Error termination via Lnk1e in /share1/g09/g09d01/l301.exe at Fri Nov  6 14:24:08 2015.
Job cpu time:       0 days  0 hours  0 minutes  1.3 seconds.
File lengths (MBytes):  RWF=     96 Int=      0 D2E=      0 Chk=      1 Scr=      1


我很困惑到底应该如何在oniom中加入dispersion关键字?多谢指教。



6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2015-11-6 15:29:02 | 只看该作者 Only view this author
你就直接用本身自带色散校正的泛函,诸如wB97XD、B97D,或者用M062X这样本身描述色散作用好的泛函,这样绝对没问题。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

40

帖子

0

威望

97

eV
积分
137

Level 2 能力者

3#
发表于 Post on 2015-11-6 18:42:36 | 只看该作者 Only view this author
sobereva 发表于 2015-11-6 15:29
你就直接用本身自带色散校正的泛函,诸如wB97XD、B97D,或者用M062X这样本身描述色散作用好的泛函,这样绝 ...

请假sob老师,我用oniom计算分层时层与层之间没有链接原子,输入文件是oniom=embedcharge,我计算能量的应该怎么算

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

4#
发表于 Post on 2015-11-6 19:34:21 | 只看该作者 Only view this author
那年·时光飞 发表于 2015-11-6 18:42
请假sob老师,我用oniom计算分层时层与层之间没有链接原子,输入文件是oniom=embedcharge,我计算能量的 ...

如果你在gv里,原子之间都有键连着,直接保存出来的输入文件里那些成键原子间应该有链接信息。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

40

帖子

0

威望

97

eV
积分
137

Level 2 能力者

5#
发表于 Post on 2015-11-6 20:36:55 | 只看该作者 Only view this author
sobereva 发表于 2015-11-6 19:34
如果你在gv里,原子之间都有键连着,直接保存出来的输入文件里那些成键原子间应该有链接信息。

我分层是分三层,比如:内层是两个分子,中层是两个分子,外层是一个分子,这样分子之间是不存在链接原子的,我不知道在输入文件中写oniom=embedcharge 对不对

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

6#
发表于 Post on 2015-11-6 21:34:07 | 只看该作者 Only view this author
这你就不用管链接原子了。这根写不写embedcharge完全没关系。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

40

帖子

0

威望

97

eV
积分
137

Level 2 能力者

7#
发表于 Post on 2015-11-7 15:06:19 | 只看该作者 Only view this author
sobereva 发表于 2015-11-6 21:34
这你就不用管链接原子了。这根写不写embedcharge完全没关系。

sob老师,我是自己摸索的oniom,有很多不懂的地方,在写输入文件的时候,不写oniom=ebedcharge,直接是101错误,加上的话还能算出来,但是我不知道这个embedcharge加上有什么影响,会不会影响到计算的结果

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

8#
发表于 Post on 2015-11-7 15:30:29 | 只看该作者 Only view this author
那年·时光飞 发表于 2015-11-7 15:06
sob老师,我是自己摸索的oniom,有很多不懂的地方,在写输入文件的时候,不写oniom=ebedcharge,直接是10 ...


应该看具体错误提示判断。

embedcharge这只不过是让程序使用电子嵌入方法,也就是允许MM区域的原子电荷极化QM部分,使描述更合理。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

54

帖子

0

威望

591

eV
积分
645

Level 4 (黑子)

9#
发表于 Post on 2015-11-7 20:26:03 | 只看该作者 Only view this author
貌似你后面的配体没有电荷,建议你先用antechamber处理小分子配体,gaussian里面的amber力场也比较旧,自己添加吧

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
sobereva + 2

查看全部评分 View all ratings

3

帖子

0

威望

41

eV
积分
44

Level 2 能力者

10#
 楼主 Author| 发表于 Post on 2015-11-8 14:13:49 | 只看该作者 Only view this author
sobereva 发表于 2015-11-6 15:29
你就直接用本身自带色散校正的泛函,诸如wB97XD、B97D,或者用M062X这样本身描述色散作用好的泛函,这样绝 ...

多谢回复。我测试了一下B97D和M062X,感觉和PBE0+G3比起来,高估了弱相互作用,和实验结果相比,还是偏差有些大,不知道还有没有可以取代PBE0+G3的泛函。

3

帖子

0

威望

41

eV
积分
44

Level 2 能力者

11#
 楼主 Author| 发表于 Post on 2015-11-8 14:16:54 | 只看该作者 Only view this author
planet5460 发表于 2015-11-7 20:26
貌似你后面的配体没有电荷,建议你先用antechamber处理小分子配体,gaussian里面的amber力场也比较旧,自己 ...

配体是用QM处理的,没有加电荷应该没有关系吧。根据我的其它作业的经验,不对QM区域加电荷应该不会影响结果的。我的MM区域是DNA,amber力场应该是现成的,不需要另行添加。

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

12#
发表于 Post on 2015-11-8 15:45:50 | 只看该作者 Only view this author
sherry3116 发表于 2015-11-8 14:13
多谢回复。我测试了一下B97D和M062X,感觉和PBE0+G3比起来,高估了弱相互作用,和实验结果相比,还是偏差 ...

PBE0本身没法描述色散作用。
G3不知道你具体指什么,热力学组合方法原理上是没法混合进ONIOM来用的。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

40

帖子

0

威望

97

eV
积分
137

Level 2 能力者

13#
发表于 Post on 2015-11-8 16:19:00 | 只看该作者 Only view this author
sobereva 发表于 2015-11-7 15:30
应该看具体错误提示判断。

embedcharge这只不过是让程序使用电子嵌入方法,也就是允许MM区域的原子 ...

如果不加embedcharge的话,直接是101错误,输入文件的问题,我觉得是我分层没有链接原子,但是加入embedcharge就能继续往下算,有时候构型好的话还能normal,愈加不理解为什么会出现这种问题了,麻烦sob老师了,非常感谢

40

帖子

0

威望

97

eV
积分
137

Level 2 能力者

14#
发表于 Post on 2015-11-8 16:24:31 | 只看该作者 Only view this author
sobereva 发表于 2015-11-7 15:30
应该看具体错误提示判断。

embedcharge这只不过是让程序使用电子嵌入方法,也就是允许MM区域的原子 ...

sob老师  您那有关于oniom分层,没有链接原子的输入文件吗?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

15#
发表于 Post on 2015-11-8 19:46:33 | 只看该作者 Only view this author
那年·时光飞 发表于 2015-11-8 16:24
sob老师  您那有关于oniom分层,没有链接原子的输入文件吗?


自己弄一个就行了。没有连接原子的话两层之间肯定得是非共价相互作用,比如你算个分子包在几个水分子中央,水用低层,溶质分子用高层。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-22 04:06 , Processed in 0.189441 second(s), 27 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list