计算化学公社

 找回密码 Forget password
 注册 Register

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

查看数: 19713 | 评论数: 20 | 收藏 Add to favorites 3
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2020-6-25 23:18

正文摘要:

各位老师好,我是做纳米聚合物仿生亲和材料这一块的,其中实验最重要的一部分是在氯仿中让两个小分子结合到一起,小分子的结构式我在下图给出来了,所以我想通过分子动力学的手段分析一下两个小分子是否能够结合在一 ...

回复 Reply

华理小妮 发表于 Post on 2020-7-25 09:10:30
sobereva 发表于 2020-7-25 00:04
不代表ORCA就一定适合大体系
ORCA的振动分析不仅慢,而且容易出虚频,而且溶剂模型做得没有Gaussian成熟 ...

好的sober老师
sobereva 发表于 Post on 2020-7-25 00:04:07
华理小妮 发表于 2020-7-24 22:26
sober老师,小分子的话可以用Gaussian来计算结合自由能,但是体系大的话是不是需要用ORCA这种软件来计算 ...

不代表ORCA就一定适合大体系
ORCA的振动分析不仅慢,而且容易出虚频,而且溶剂模型做得没有Gaussian成熟稳健

计算步骤和使用Gaussian没有什么本质差异
华理小妮 发表于 Post on 2020-7-24 22:26:23
sobereva 发表于 2020-7-24 07:49
从精度上说,用量子化学程序(不仅限于Gaussian)来做更理想,也更容易。
二聚体的最稳定构型可以通过本 ...

sober老师,小分子的话可以用Gaussian来计算结合自由能,但是体系大的话是不是需要用ORCA这种软件来计算,老师这边有用orca软件计算结合自由能的帖子吗,想学一下。
华理小妮 发表于 Post on 2020-7-24 10:14:03
sobereva 发表于 2020-7-24 07:49
从精度上说,用量子化学程序(不仅限于Gaussian)来做更理想,也更容易。
二聚体的最稳定构型可以通过本 ...

好的,谢谢sober老师。
sobereva 发表于 Post on 2020-7-24 07:49:38
华理小妮 发表于 2020-7-23 11:01
sober老师,这个热力学角度算结合能的时候,为了更加贴近真实情况你建议用Gaussian来做结合能呢,还是gro ...

从精度上说,用量子化学程序(不仅限于Gaussian)来做更理想,也更容易。
二聚体的最稳定构型可以通过本文方式很容易地搜索:
genmer:生成团簇初始构型结合molclus做团簇结构搜索的超便捷工具
http://bbs.keinsci.com/thread-2369-1-1.html

使用Molclus结合xtb做的动力学模拟对瑞德西韦(Remdesivir)做构象搜索
http://bbs.keinsci.com/thread-16255-1-1.htmlhttp://sobereva.com/532

之后按照下文计算二聚体和单体自由能,求差即可
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327http://bbs.keinsci.com/thread-3345-1-1.html
华理小妮 发表于 Post on 2020-7-23 11:01:08

sober老师,这个热力学角度算结合能的时候,为了更加贴近真实情况你建议用Gaussian来做结合能呢,还是gromacs来做,有必要这两种都做一下嘛,另外从分子动力学角度的话用的gromacs, GAFF力场RESP2电荷。想问问您的建议。谢谢
yjmaxpayne 发表于 Post on 2020-7-16 15:27:04
试试Gromacs + Plumed用Well-tempered Metadynamics跑。这个模拟的CV很好定义,一个分子间质心距离,一个是分子间的取向
sobereva 发表于 Post on 2020-6-29 01:12:58
华理小妮 发表于 2020-6-27 18:39
好的,谢谢sober老师,您说的很有道理。
另外有个问题想听听您的建议,我做的是非共价分子印迹聚合物方 ...

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

评分 Rate

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

查看全部评分 View all ratings

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

华理小妮 发表于 Post on 2020-6-28 12:12:27
sobereva 发表于 2020-6-27 18:04
不是增加分子的能量,而是低温下热运动、溶剂碰撞的影响会相对较弱,相对较强的分子间作用对二聚体的稳定 ...

sober老师,如果我想知道A B两个分子在氯仿溶液中能否自发的形成复合物的话,是不是可以分两步走,从热力学角度考虑,就是算一下结合能,看结合能是否是负值,如果是负值则说明能够自发的形成复合物,如果不是负值则没办法自发的形成复合物
从动力学角度,就是跑一下动力学,看看实际情况下有溶剂参与的情况下能否自发的形成复合物
华理小妮 发表于 Post on 2020-6-27 18:39:25
sobereva 发表于 2020-6-27 18:04
不是增加分子的能量,而是低温下热运动、溶剂碰撞的影响会相对较弱,相对较强的分子间作用对二聚体的稳定 ...

好的,谢谢sober老师,您说的很有道理。
另外有个问题想听听您的建议,我做的是非共价分子印迹聚合物方面的工作,在聚合的时候有一个条件必须要得到就是用到的模板分子和功能单体一定要先通过非共价键的方式(主要是氢键)结合形成复合物,之前我们没有考虑用分子模拟的手段去分析这个复合物形成的合理性,现在一分析不要紧,发现通过非共价键形成复合物的几率很小,并且就算形成复合物复合物的稳定性也不好,如果这个结果准确的话,会推翻非常多的现在已经发表过的非共价分子印迹方面的文章呢。
sobereva 发表于 Post on 2020-6-27 18:04:39
华理小妮 发表于 2020-6-27 14:00
sober老师我用的溶剂是氯仿,如果将温度设定在4摄氏度下,并且保证溶剂分子占的比重小一点的话,增加这两 ...

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

在我来看M$那种东西算这个很不适合,除了又贵又慢又不灵活以外,由于过程太傻瓜化,很多对结果有关键性影响的东西(比如原子电荷的合理性等)容易被用户忽视,导致模拟结果其实不合理但用户可能都意识不到。
华理小妮 发表于 Post on 2020-6-27 14:07:26
喵星大佬 发表于 2020-6-27 01:40
如果是在氯仿中,作为介电常数不大的溶剂,分离的电荷并不会被溶剂稳定,基本来说不会是盐桥的作用,应该都 ...

谢谢您的回复,我的溶剂是氯仿,您说的很有道理,因为我用Gaussian算过这两个分子结合的结合自由能ΔG,发现这个ΔG虽然是负值但是负值很大为-0.12849Kcal.mol, 如果将Gaussian的结果和分子动力学的结果相结合可不可以说明实际情况下这两个分子结合不成复合物。
我之所这么关注这两个分子是否形成复合物,是因为我做的非共价分子印迹聚合物在聚合的时候有一个条件必须要得到就是用到的模板分子和功能单体一定要先形成复合物,之前我们没有考虑用分子模拟的手段去分析这个复合物形成的合理性,现在一分析不要紧,发现形成复合物的几率很小,这样真的可以推翻很多这样的文章
华理小妮 发表于 Post on 2020-6-27 14:00:40
sobereva 发表于 2020-6-27 00:34
本身这俩分子间的相互作用就不强,实际情况中,又有热运动又有溶剂的碰撞,在不是特别低的温度下能稳定结合 ...

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

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

GMT+8, 2026-2-21 08:25 , Processed in 0.183981 second(s), 28 queries , Gzip On.

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