计算化学公社

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

[Material Studio] 求助两个小分子在溶液中的自组装过程研究

[复制链接 Copy URL]

97

帖子

0

威望

317

eV
积分
414

Level 3 能力者

各位老师好,我是做纳米聚合物仿生亲和材料这一块的,其中实验最重要的一部分是在氯仿中让两个小分子结合到一起,小分子的结构式我在下图给出来了,所以我想通过分子动力学的手段分析一下两个小分子是否能够结合在一起,同时添加比例为多少的情况下,这两个分子形成复合物的比例高,我用的是material studio的forcite模块进行的动力学分析,通过优化、退火、动力学运行了2ns之后,结果两个小分子根本没有相互靠近的趋势,所以想问一下大家分子动力学是不是不能模拟这种分子之间弱相互作用比较弱的自组装过程,或者说MS 不可以的话有其他软件可以做相应的研究吗。



94-75-7.png (9.34 KB, 下载次数 Times of downloads: 44)

94-75-7.png

100-43-6.png (7.98 KB, 下载次数 Times of downloads: 53)

100-43-6.png

326

帖子

0

威望

5312

eV
积分
5638

Level 6 (一方通行)

2#
发表于 Post on 2020-6-26 10:03:37 | 只看该作者 Only view this author
2ns时间太短了,感觉至少得跑几十ns才能说事,lz可以多放点分子,这样碰到一起的几率大点

97

帖子

0

威望

317

eV
积分
414

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2020-6-26 11:46:57 | 只看该作者 Only view this author
mol 发表于 2020-6-26 10:03
2ns时间太短了,感觉至少得跑几十ns才能说事,lz可以多放点分子,这样碰到一起的几率大点

我是按照实验过程中这两个分子和溶剂分子的摩尔比加的,如果多放这两个分子又害怕和实验不相符呢。
请问常规来说分子动力学可以做两个结合力特别弱分子的自组装过程吗,谢谢

326

帖子

0

威望

5312

eV
积分
5638

Level 6 (一方通行)

4#
发表于 Post on 2020-6-26 16:28:00 | 只看该作者 Only view this author
华理小妮 发表于 2020-6-26 11:46
我是按照实验过程中这两个分子和溶剂分子的摩尔比加的,如果多放这两个分子又害怕和实验不相符呢。
请问 ...

我认为没关系,完全可以多加,没必要严格和实验对应,实际上到微观水平,分子也不是均匀分散的了,完全可以认为我的模拟格子里瞬间有很多分子。lz可以试试把两个小分子摆到一块,慢慢升温,看看他俩是否稳定。

评分 Rate

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

查看全部评分 View all ratings

97

帖子

0

威望

317

eV
积分
414

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2020-6-26 16:41:26 | 只看该作者 Only view this author
mol 发表于 2020-6-26 16:28
我认为没关系,完全可以多加,没必要严格和实验对应,实际上到微观水平,分子也不是均匀分散的了,完全可 ...

这个我做过了,我是先在gaussian里面做的它们两个形成复合物的结构优化,然后用gromacs测试它们两个形成复合物的稳定性,结果发现在60摄氏度的时候它们复合物就非常不稳定了。

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

6#
发表于 Post on 2020-6-27 00:34:07 | 只看该作者 Only view this author
本身这俩分子间的相互作用就不强,实际情况中,又有热运动又有溶剂的碰撞,在不是特别低的温度下能稳定结合才显得奇怪。
也应当说清楚用的是什么溶剂。

算这种问题gromacs比M$强得多得多,速度快一个数量级以上,用起来还灵活,轨迹分析还方便,还不要钱。
北京科音自然科学研究中心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

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

7#
发表于 Post on 2020-6-27 01:40:43 | 只看该作者 Only view this author
如果是在氯仿中,作为介电常数不大的溶剂,分离的电荷并不会被溶剂稳定,基本来说不会是盐桥的作用,应该都是分子形态。
如果考虑氢键,唯一的氢键供体就是羧酸的OH,而吡啶分子中唯一的氢键受体即为吡啶N,不过这种优化结构的时候很容易会变成酸碱反应的过渡态,这个你貌似已经在水里优化出来过了。
如果考虑色散或者说pi-pi作用,这里两个芳环间的pi-pi作用应该是很弱的,尤其在室温,溶剂下,应该几乎体现不出来,即使是教科书里典型的比较强的pi-pi作用,苯与六氟苯的堆积也只有在晶体中才能体现。
所以得不到结果很正常

97

帖子

0

威望

317

eV
积分
414

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2020-6-27 14:00:40 | 只看该作者 Only view this author
sobereva 发表于 2020-6-27 00:34
本身这俩分子间的相互作用就不强,实际情况中,又有热运动又有溶剂的碰撞,在不是特别低的温度下能稳定结合 ...

sober老师我用的溶剂是氯仿,如果将温度设定在4摄氏度下,并且保证溶剂分子占的比重小一点的话,增加这两个分子的量能形成复合物吗。
我现在正在学着用garomacs,因为MS 的forcite模块用起来比较简单,我想先用这个简单操作的看一下具不具备做的可能性,然后在用forcite.

97

帖子

0

威望

317

eV
积分
414

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2020-6-27 14:07:26 | 只看该作者 Only view this author
喵星大佬 发表于 2020-6-27 01:40
如果是在氯仿中,作为介电常数不大的溶剂,分离的电荷并不会被溶剂稳定,基本来说不会是盐桥的作用,应该都 ...

谢谢您的回复,我的溶剂是氯仿,您说的很有道理,因为我用Gaussian算过这两个分子结合的结合自由能ΔG,发现这个ΔG虽然是负值但是负值很大为-0.12849Kcal.mol, 如果将Gaussian的结果和分子动力学的结果相结合可不可以说明实际情况下这两个分子结合不成复合物。
我之所这么关注这两个分子是否形成复合物,是因为我做的非共价分子印迹聚合物在聚合的时候有一个条件必须要得到就是用到的模板分子和功能单体一定要先形成复合物,之前我们没有考虑用分子模拟的手段去分析这个复合物形成的合理性,现在一分析不要紧,发现形成复合物的几率很小,这样真的可以推翻很多这样的文章

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

10#
发表于 Post on 2020-6-27 18:04:39 | 只看该作者 Only view this author
华理小妮 发表于 2020-6-27 14:00
sober老师我用的溶剂是氯仿,如果将温度设定在4摄氏度下,并且保证溶剂分子占的比重小一点的话,增加这两 ...

不是增加分子的能量,而是低温下热运动、溶剂碰撞的影响会相对较弱,相对较强的分子间作用对二聚体的稳定化效果更容易被展现出来。

在我来看M$那种东西算这个很不适合,除了又贵又慢又不灵活以外,由于过程太傻瓜化,很多对结果有关键性影响的东西(比如原子电荷的合理性等)容易被用户忽视,导致模拟结果其实不合理但用户可能都意识不到。
北京科音自然科学研究中心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

97

帖子

0

威望

317

eV
积分
414

Level 3 能力者

11#
 楼主 Author| 发表于 Post on 2020-6-27 18:39:25 | 只看该作者 Only view this author
sobereva 发表于 2020-6-27 18:04
不是增加分子的能量,而是低温下热运动、溶剂碰撞的影响会相对较弱,相对较强的分子间作用对二聚体的稳定 ...

好的,谢谢sober老师,您说的很有道理。
另外有个问题想听听您的建议,我做的是非共价分子印迹聚合物方面的工作,在聚合的时候有一个条件必须要得到就是用到的模板分子和功能单体一定要先通过非共价键的方式(主要是氢键)结合形成复合物,之前我们没有考虑用分子模拟的手段去分析这个复合物形成的合理性,现在一分析不要紧,发现通过非共价键形成复合物的几率很小,并且就算形成复合物复合物的稳定性也不好,如果这个结果准确的话,会推翻非常多的现在已经发表过的非共价分子印迹方面的文章呢。

97

帖子

0

威望

317

eV
积分
414

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2020-6-28 12:12:27 | 只看该作者 Only view this author
sobereva 发表于 2020-6-27 18:04
不是增加分子的能量,而是低温下热运动、溶剂碰撞的影响会相对较弱,相对较强的分子间作用对二聚体的稳定 ...

sober老师,如果我想知道A B两个分子在氯仿溶液中能否自发的形成复合物的话,是不是可以分两步走,从热力学角度考虑,就是算一下结合能,看结合能是否是负值,如果是负值则说明能够自发的形成复合物,如果不是负值则没办法自发的形成复合物
从动力学角度,就是跑一下动力学,看看实际情况下有溶剂参与的情况下能否自发的形成复合物

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

13#
发表于 Post on 2020-6-29 01:11:11 | 只看该作者 Only view this author
华理小妮 发表于 2020-6-28 12:12
sober老师,如果我想知道A B两个分子在氯仿溶液中能否自发的形成复合物的话,是不是可以分两步走,从热力 ...

北京科音自然科学研究中心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

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

14#
发表于 Post on 2020-6-29 01:12:58 | 只看该作者 Only view this author
华理小妮 发表于 2020-6-27 18:39
好的,谢谢sober老师,您说的很有道理。
另外有个问题想听听您的建议,我做的是非共价分子印迹聚合物方 ...

只要发生双分子化学反应(除非是无势垒的),肯定先要形成弱相互作用复合物,这是众所周知的事情。如果连形成复合物的几率都非常小,反应的速率就必然更是非常小了。
但也应当注意微观和宏观的差异性,模拟微观情况时可能很难出现的过程,在宏观尺度上由于体系尺寸、时间尺度都大得多得多,并非一定宏观上难以发生。比如冰结冰就是。

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
辣椒丛中水果糖 + 2 我很赞同

查看全部评分 View all ratings

北京科音自然科学研究中心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

144

帖子

0

威望

2820

eV
积分
2964

Level 5 (御坂)

15#
发表于 Post on 2020-7-16 15:27:04 | 只看该作者 Only view this author
试试Gromacs + Plumed用Well-tempered Metadynamics跑。这个模拟的CV很好定义,一个分子间质心距离,一个是分子间的取向

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 04:00 , Processed in 0.184151 second(s), 24 queries , Gzip On.

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