计算化学公社

 找回密码 Forget password
 注册 Register
Views: 14381|回复 Reply: 8

[综合交流] 构建被显式溶剂层包裹的分子的简单方法

[复制链接 Copy URL]

4万

帖子

99

威望

4万

eV
积分
89863

管理员

公社社长+计算化学玩家

发表于 Post on 2018-1-21 22:58:43 | 显示全部楼层 Show all |阅读模式 Reading model
构建被显式溶剂层包裹的分子的简单方法

文/Sobereva@北京科音  2018-Jan-20

0 前言

溶剂模型分两大类,显式溶剂模型、隐式溶剂模型。在量子化学研究中,一般为了省事、方便以及为了考虑溶剂平均效应而用隐式溶剂模型,但是有时候溶剂效应没法靠连续介质的考虑方式描述,非得用显式溶剂模型不可。有些情况只需要很少的显式溶剂即可,比如水催化的反应,只需要把参与反应的水分子纳入体系即可。但研究许多问题的时候,显式溶剂放在哪、怎么放,哪些位置的显式溶剂会起到关键作用并不明确(比如考虑显式溶剂对UV-Vis光谱影响的问题),那么为了保险就得给整个分子加上一个溶剂层。通常显式溶剂和隐式溶剂是一起使用的,后者描述显式溶剂层外部的体相溶剂效应,这称为杂化溶剂模型。

加溶剂层比较严谨、可靠的做法是对体系跑分子动力学,然后从中提取一些帧(可以均匀间隔采样,也可以做周期性退火,取对应0K的帧),然后把靠近溶质的溶剂分子连同溶质提取出来,再用量化方法优化。但这样做对于不懂MD的人门槛略高,有时候处理的体系还不好生成拓扑文件、弄到适合的力场参数。

本文目的是介绍比较简单的给分子加上溶剂层的做法,不需要跑动力学。虽然没有MD严格,但对于搞量化的人来说多数情况也够用了,本文的方法也比较普适。本文使用的程序为Gaussian16 A.03、Packmol 16.344、VMD 1.9.3、MOPAC2016。本文以构建靛青(indigo)分子被显式水包围的体系作为例子。

本文涉及的所有文件可在此处下载: 406.7z (2.87 MB, 下载次数 Times of downloads: 189)

评分 Rate

参与人数
Participants 7
eV +28 收起 理由
Reason
yangyangyang + 3 谢谢
cyx98 + 4 赞!
keiskei + 4 好物!
Jack + 5 精品内容
三石草祭 + 5 谢谢
函数与激情 + 2 谢谢分享
vigaryang + 5 精品内容

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班,内容介绍以及往届资料购买请点击链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的最佳途径。培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人,讨论范畴相同
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

292

帖子

0

威望

1813

eV
积分
2105

Level 5 (御坂)

发表于 Post on 2018-1-22 00:56:39 来自手机 | 显示全部楼层 Show all
好 :鼓掌:鼓掌  只是不知道加上显式水分子之后,计算量如何? 比如50个原子的体系,加水之后共有300个原子;那这个体系的计算量是不是跟平常算300原子体系的计算量一样大?  还是说加上去的显式溶剂有简化的计算方法?

4万

帖子

99

威望

4万

eV
积分
89863

管理员

公社社长+计算化学玩家

 楼主 Author| 发表于 Post on 2018-1-22 10:55:28 | 显示全部楼层 Show all
霜晨月 发表于 2018-1-22 00:56
好 :鼓掌:鼓掌  只是不知道加上显式水分子之后,计算量如何? 比如50个原子的体系,加水之后共有300个原子 ...


是。

得到本文的最终包裹结构后,应当根据要用的计算级别(DFT还是半经验)、计算程序(诸如同等计算量ORCA能算得动的体系远比Gaussian大得多)、具体情况(那些水可能比较重要、有多少水是必须保留的),来手动修改,去掉没太大意义的水分子,进一步减小体系来降低后续计算耗时
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班,内容介绍以及往届资料购买请点击链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的最佳途径。培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人,讨论范畴相同
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

85

帖子

0

威望

1577

eV
积分
1662

Level 5 (御坂)

发表于 Post on 2020-1-13 12:40:46 | 显示全部楼层 Show all
本帖最后由 wanlichuan 于 2020-1-13 12:41 编辑

我做了一个练习。到了“4 用MOPAC2016优化溶质+溶剂层体系”这一步发现一个问题,不固定溶质分子的坐标,可以正常优化,但是一旦给溶质分子的坐标后面都加上0,就不优化了,报错信息如下:
           Faulty atom:  307
            Faulty line: " C        -0.948000        2.360000        0.116000 0"
          Unless MINI is used, optimization flags must be 1, 0, or -1
后来把需要优化溶剂分子的坐标后面加上1,不需要优化的溶质分子的坐标后面加上0,还是不能优化,报错信息如下:
          Faulty atom:    5
          Faulty line: " H        -6.112000        1.627000        4.565000 1"
          Unless MINI is used, optimization flags must be 1, 0, or -1
如果坐标后面什么都不加,所有分子都优化,则一切正常。
这是为什么呢?
我把第二种情况的输入输出文件上传了(整个体系带一个正电荷,所以加上了“CHARGE=1”),请卢老师帮忙看看。
谢谢。

B-UFF-select.out

885 Bytes, 下载次数 Times of downloads: 0

B-UFF-select.mop

16.89 KB, 下载次数 Times of downloads: 0

4万

帖子

99

威望

4万

eV
积分
89863

管理员

公社社长+计算化学玩家

 楼主 Author| 发表于 Post on 2020-1-13 20:54:24 | 显示全部楼层 Show all
wanlichuan 发表于 2020-1-13 12:40
我做了一个练习。到了“4 用MOPAC2016优化溶质+溶剂层体系”这一步发现一个问题,不固定溶质分子的坐标,可 ...

把含有结构信息的文件(如xyz、pdb等)载入Multiwfn,进入主功能100的子功能2,选择产生MOPAC输入文件,在界面里提供了选项来设定对哪些原子进行冻结,设好之后选择计算级别,就可以得到你需要的mop文件。这样产生mop文件又方便又不会弄错。
(必须是目前官网上最新版本Multiwfn才有此功能)
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班,内容介绍以及往届资料购买请点击链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的最佳途径。培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人,讨论范畴相同
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

85

帖子

0

威望

1577

eV
积分
1662

Level 5 (御坂)

发表于 Post on 2020-1-14 08:49:36 | 显示全部楼层 Show all
好的,谢谢卢老师。

595

帖子

11

威望

9375

eV
积分
10190

Level 6 (一方通行)

发表于 Post on 2020-2-22 21:56:49 | 显示全部楼层 Show all
本帖最后由 captain 于 2020-2-22 21:59 编辑

请问大神,关于这句话,“然后把靠近溶质的溶剂分子连同溶质提取出来,再用量化方法优化”,
用量化方法优化时,应该最好带着隐式溶剂模型优化吧?(即采用的是杂化溶剂模型)
计算化学与分子模拟

4万

帖子

99

威望

4万

eV
积分
89863

管理员

公社社长+计算化学玩家

 楼主 Author| 发表于 Post on 2020-2-22 22:36:14 | 显示全部楼层 Show all
captain 发表于 2020-2-22 21:56
请问大神,关于这句话,“然后把靠近溶质的溶剂分子连同溶质提取出来,再用量化方法优化”,
用量化方法优 ...


通常没有任何理由只用显式而不用杂化溶剂模型
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班,内容介绍以及往届资料购买请点击链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的最佳途径。培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人,讨论范畴相同
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

595

帖子

11

威望

9375

eV
积分
10190

Level 6 (一方通行)

发表于 Post on 2020-2-22 23:13:45 | 显示全部楼层 Show all
sobereva 发表于 2020-2-22 22:36

通常没有任何理由只用显式而不用杂化溶剂模型

明白了,谢大神指点!
计算化学与分子模拟

本版积分规则 Credits rule

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

GMT+8, 2023-2-2 00:22 , Processed in 0.211558 second(s), 25 queries .

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