计算化学公社

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

[Gaussian/gview] 金球表面含有多苯环的π-π作用和氢键的体系的建模问题

[复制链接 Copy URL]

27

帖子

0

威望

675

eV
积分
702

Level 4 (黑子)

sob老师好,这个一个分子印迹体系,想计算一下模板分子和基质是以何种形式相作用的,比如π-π堆积,氢键等。
因为之前只做过有机体系的,前两天上的基础班也没怎么讲过渡金属和周期性的方面,所以问题有很多。
金球表征过,是3.5-5纳米。金球表面,以金硫键连接着4-氨基苯硫酚,文献里说金硫键是这样连接的。

用高斯pm6d3预优化过,当然也没收敛。

模板分子是个多巴胺。
因为之前只做过有机体系的,前两天上的基础班也没怎么讲过渡金属和周期性的方面,所以问题有很多。
1,多巴胺分子四周大概都有4-氨基苯硫酚与之连接,所以该放多少4-氨基苯硫酚呢。
2,其实包括金硫键在内的金球都不是要描述的重点,但是另一方面,金球的存在使4-氨基苯硫酚在多巴胺周围有个密度的问题(如1),即多巴胺周围是否有可能接触到那么多的4-氨基苯硫酚。那么是否需要画出金球呢,还是忽略掉。
3,如果金球忽略掉,之前用b3lyp-D3/6-31+G**和lanl2dz基组算过,24核超级慢,也不收敛。之后考虑用ORCA算是否可行,老师说开了RI超级快的?
4,所以最重要的还是如何合理构建这个体系,以及用啥程序啥方法基组的问题啊摔。
5,啊咧,老师哪个培训班会好好讲讲ORCA,好感兴趣的说。
谢谢sob老师




6万

帖子

99

威望

6万

eV
积分
125152

管理员

公社社长

2#
发表于 Post on 2018-3-29 16:37:44 | 只看该作者 Only view this author
过渡金属体系在培训中已经涉及不少了,很多实例和练习都是有机金属化合物。在没来得及讲但是在讲义上详细写了的周期性计算、赝势与相对论部分都有大量涉及。一定要注意仔细阅读讲义上没来及讲的部分并仔细复习讲授过的内容。计算当前体系所需知识绝对没有超越量化基础班的内容范畴。

由于对你的体系不了解,前两个问题我不清楚

3 可以用ORCA通过BP86-D3-gCP/def2-SVP计算。或者用B97-3c,性价比更高。
只要计算资源不是很拮据,用高斯算也完全可以。

4 什么量化程序都能算。适合的泛函,就是按照课上讲DFT时候对泛函的推荐来选择。基组也没有特殊性

5 ORCA是高级班的范畴,基础知识通过阅读本版ORCA分类的帖子以及手册很容易就能学会。
北京科音自然科学研究中心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

909

帖子

37

威望

5527

eV
积分
7176

Level 6 (一方通行)

3#
发表于 Post on 2018-3-30 03:49:44 | 只看该作者 Only view this author
本帖最后由 ggdh 于 2018-3-30 03:50 编辑

水溶液中多巴胺哪个胺基质子化,然后和苯胺上的氨基形成H键作用苯环之间基本没有stacking, 有点C-H-π作用

27

帖子

0

威望

675

eV
积分
702

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2018-4-2 14:41:05 | 只看该作者 Only view this author
sobereva 发表于 2018-3-29 16:37
过渡金属体系在培训中已经涉及不少了,很多实例和练习都是有机金属化合物。在没来得及讲但是在讲义上详细写 ...

谢谢sob老师,我用orca或者pm6d3跑之后,都会使不同分子的s原子和au原子成键,苯环也堆积了起来,然后scf就不收敛了。请问要如何解决呢, 限制性优化吗?我觉得还是这个体系构建的有问题,但不知如何解决



6万

帖子

99

威望

6万

eV
积分
125152

管理员

公社社长

5#
发表于 Post on 2018-4-2 15:40:27 | 只看该作者 Only view this author
cottondog 发表于 2018-4-2 14:41
谢谢sob老师,我用orca或者pm6d3跑之后,都会使不同分子的s原子和au原子成键,苯环也堆积了起来,然后scf ...

成键是正常的,硫很亲金。必要时候该冻结就做冻结
北京科音自然科学研究中心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

27

帖子

0

威望

675

eV
积分
702

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2018-4-4 12:05:30 | 只看该作者 Only view this author
sobereva 发表于 2018-4-2 15:40
成键是正常的,硫很亲金。必要时候该冻结就做冻结

请教老师,用IGM做分子间弱相互作用的时候,载入IGM_inter.vmd,VMD中金和硫没有显示键,请问怎么手动连接两个原子成键呢,或者如何载入包含成键信息的文件?

27

帖子

0

威望

675

eV
积分
702

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2018-4-4 12:23:51 | 只看该作者 Only view this author
sobereva 发表于 2018-4-2 15:40
成键是正常的,硫很亲金。必要时候该冻结就做冻结

我先去读下    使VMD根据pdb文件中的CONECT字段设定原子连接关系 这篇博文

6万

帖子

99

威望

6万

eV
积分
125152

管理员

公社社长

8#
发表于 Post on 2018-4-4 13:31:18 | 只看该作者 Only view this author
cottondog 发表于 2018-4-4 12:23
我先去读下    使VMD根据pdb文件中的CONECT字段设定原子连接关系 这篇博文

用mouse-add/remove bonds,手动点击两个原子也可以直接给连上
北京科音自然科学研究中心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

27

帖子

0

威望

675

eV
积分
702

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2018-4-4 14:04:56 | 只看该作者 Only view this author
sobereva 发表于 2018-4-4 13:31
用mouse-add/remove bonds,手动点击两个原子也可以直接给连上

sob老师,是不是因为金硫键是手动连接的,所以设置所有fragment后,金原子还是无法被设置,是否可以直接设置某个原子如atom 6?

6万

帖子

99

威望

6万

eV
积分
125152

管理员

公社社长

10#
发表于 Post on 2018-4-4 14:09:28 | 只看该作者 Only view this author
cottondog 发表于 2018-4-4 14:04
sob老师,是不是因为金硫键是手动连接的,所以设置所有fragment后,金原子还是无法被设置,是否可以直接 ...

不是已经连上了?被设置指什么?
北京科音自然科学研究中心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

27

帖子

0

威望

675

eV
积分
702

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2018-4-4 14:13:06 | 只看该作者 Only view this author
sobereva 发表于 2018-4-4 14:09
不是已经连上了?被设置指什么?

抱歉,我没说清楚,我是在模仿老师这个图但是金原子不能使用fragment设置成主体的那个显示模式,其他片段都可以,唯独剩下了金。

6万

帖子

99

威望

6万

eV
积分
125152

管理员

公社社长

12#
发表于 Post on 2018-4-4 17:48:46 | 只看该作者 Only view this author
cottondog 发表于 2018-4-4 14:13
抱歉,我没说清楚,我是在模仿老师这个图但是金原子不能使用fragment设置成主体的那个显示模式,其他片段 ...

还是不清楚你的情况...
再鼓捣鼓捣吧
北京科音自然科学研究中心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

27

帖子

0

威望

675

eV
积分
702

Level 4 (黑子)

13#
 楼主 Author| 发表于 Post on 2018-4-4 19:15:11 | 只看该作者 Only view this author
sobereva 发表于 2018-4-4 17:48
还是不清楚你的情况...
再鼓捣鼓捣吧

我是做IGM嘛,载入vmd时候,金和硫是不连着的,手动给连上了,但是vmd还是不能判断是同一个分子片段。 后面做主客体分别设置显示模式(Representation)时候,用fragment N也设置不上,大概因为不是片段吧,所以问问老师,单个原子怎么设置,或者其他方法给连上?

6万

帖子

99

威望

6万

eV
积分
125152

管理员

公社社长

14#
发表于 Post on 2018-4-4 23:23:54 | 只看该作者 Only view this author
cottondog 发表于 2018-4-4 19:15
我是做IGM嘛,载入vmd时候,金和硫是不连着的,手动给连上了,但是vmd还是不能判断是同一个分子片段。 后 ...


VMD是在刚载入结构的时候根据判断的成键关系确定fragment是怎么划分的
如果是想选定某些原子区域,而这些区域又不正好对应某个自动确定的fragment,就用诸如serial 3 6 10 to 20这样的语句按原子序号选择就可以了,不用利用fragment
北京科音自然科学研究中心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

27

帖子

0

威望

675

eV
积分
702

Level 4 (黑子)

15#
 楼主 Author| 发表于 Post on 2018-4-11 09:52:47 | 只看该作者 Only view this author
sobereva 发表于 2018-4-4 23:23
VMD是在刚载入结构的时候根据判断的成键关系确定fragment是怎么划分的
如果是想选定某些原子区域,而 ...

老师,算相互作用能的时候还需要单独优化片段的结构吗,以及ORCA有像高斯Counterpoise的方法吗?gCP可以做BSSE矫正,但是输出相互作用能吗

本版积分规则 Credits rule

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

GMT+8, 2026-2-23 18:00 , Processed in 0.268198 second(s), 23 queries , Gzip On.

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