请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 890|回复: 2

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

[复制链接]

1万

帖子

25

威望

1万

eV
积分
32416

管理员

公社社长

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

文/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, 下载次数: 5)

评分

参与人数 4eV +17 收起 理由
Jack + 5 精品内容
三石草祭 + 5 谢谢
函数与激情 + 2 谢谢分享
vigaryang + 5 精品内容

查看全部评分

北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

116

帖子

0

威望

514

eV
积分
630

Level 4 (黑子)

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

1万

帖子

25

威望

1万

eV
积分
32416

管理员

公社社长

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


是。

得到本文的最终包裹结构后,应当根据要用的计算级别(DFT还是半经验)、计算程序(诸如同等计算量ORCA能算得动的体系远比Gaussian大得多)、具体情况(那些水可能比较重要、有多少水是必须保留的),来手动修改,去掉没太大意义的水分子,进一步减小体系来降低后续计算耗时
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2018-7-22 16:17 , Processed in 0.206026 second(s), 28 queries .

快速回复 返回顶部 返回列表