计算化学公社

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

[GROMACS] 真空条件下模拟分子发生聚集是否正常?

[复制链接 Copy URL]

44

帖子

0

威望

87

eV
积分
131

Level 2 能力者

charmm36力场,在上一次模拟当中我选取了一帧构相,加入了配体小分子,然后尝试在真空体系下,先进行nvt平衡后进行模拟,但是结果是所有物质都聚集到了一块去,请问这是正常现象吗,还是有什么办法可以解决?谢谢

1687753187645.png (61.18 KB, 下载次数 Times of downloads: 22)

初始构象

初始构象

1687753131504.png (55.3 KB, 下载次数 Times of downloads: 23)

模拟过程中的构象

模拟过程中的构象

323

帖子

0

威望

1980

eV
积分
2303

Level 5 (御坂)

2#
发表于 Post on 2023-6-26 14:10:47 | 只看该作者 Only view this author
真空里只受到吸引力的话不聚集还能干啥

44

帖子

0

威望

87

eV
积分
131

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2023-6-26 14:15:13 | 只看该作者 Only view this author
pal 发表于 2023-6-26 14:10
真空里只受到吸引力的话不聚集还能干啥

那请问如果我希望模拟溶剂很少的一个过程,有什么办法吗,我这是一个包埋体系,我希望模拟它在水含量很低情况下的跟小分子配体结合的过程

1万

帖子

0

威望

9888

eV
积分
22133

Level 6 (一方通行)

4#
发表于 Post on 2023-6-26 14:33:41 | 只看该作者 Only view this author
sd6512066 发表于 2023-6-26 07:15
那请问如果我希望模拟溶剂很少的一个过程,有什么办法吗,我这是一个包埋体系,我希望模拟它在水含量很低 ...

溶剂很少有两个可能的意思:(1)这个体系是气相里的一个几乎蒸干了的微液滴;(2)这个体系是几乎完全干燥的固相。先确定你是哪个意思
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

44

帖子

0

威望

87

eV
积分
131

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2023-6-27 09:58:50 | 只看该作者 Only view this author
wzkchem5 发表于 2023-6-26 14:33
溶剂很少有两个可能的意思:(1)这个体系是气相里的一个几乎蒸干了的微液滴;(2)这个体系是几乎完全干 ...

谢谢您的回复,我的主体是固态的,配体小分子是液态的,反应过程我把两者放在了高压反应釜当中进行反应,并没有加入其它溶剂,我是希望能够模拟这个过程来探究包合方式的

1万

帖子

0

威望

9888

eV
积分
22133

Level 6 (一方通行)

6#
发表于 Post on 2023-6-27 14:09:26 | 只看该作者 Only view this author
sd6512066 发表于 2023-6-27 02:58
谢谢您的回复,我的主体是固态的,配体小分子是液态的,反应过程我把两者放在了高压反应釜当中进行反应, ...

那显然不能在真空里模拟啊,得建立一个周期性体系做NPT模拟,使得体系里有无数多个主体和配体,把空间都填满
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

44

帖子

0

威望

87

eV
积分
131

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2023-6-27 20:49:50 | 只看该作者 Only view this author
wzkchem5 发表于 2023-6-27 14:09
那显然不能在真空里模拟啊,得建立一个周期性体系做NPT模拟,使得体系里有无数多个主体和配体,把空间都 ...

不好意思我不太了解具体要怎么操作,是指构建一个含有无数多个主体配体的盒子吗?我接触MD没多久,只熟悉教程上的操作,谢谢您

44

帖子

0

威望

87

eV
积分
131

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2023-6-27 20:52:30 | 只看该作者 Only view this author
wzkchem5 发表于 2023-6-27 14:09
那显然不能在真空里模拟啊,得建立一个周期性体系做NPT模拟,使得体系里有无数多个主体和配体,把空间都 ...

如果方便的话,不知道能不能私下跟老师您交流

1万

帖子

0

威望

9888

eV
积分
22133

Level 6 (一方通行)

9#
发表于 Post on 2023-6-27 21:29:27 | 只看该作者 Only view this author
sd6512066 发表于 2023-6-27 13:49
不好意思我不太了解具体要怎么操作,是指构建一个含有无数多个主体配体的盒子吗?我接触MD没多久,只熟悉 ...

随便一个MD的教程都会教怎么建周期性盒子(或者叫做考虑周期性边界条件)
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

44

帖子

0

威望

87

eV
积分
131

Level 2 能力者

10#
 楼主 Author| 发表于 Post on 2023-6-27 21:45:35 | 只看该作者 Only view this author
wzkchem5 发表于 2023-6-27 21:29
随便一个MD的教程都会教怎么建周期性盒子(或者叫做考虑周期性边界条件)

周期性的话,我所使用的NPT.mdp里面是有pbc的,但我还是不太理解如果不能在真空里模拟的话,我应该怎么构建我这个体系呢?

1万

帖子

0

威望

9888

eV
积分
22133

Level 6 (一方通行)

11#
发表于 Post on 2023-6-27 22:17:22 | 只看该作者 Only view this author
sd6512066 发表于 2023-6-27 14:45
周期性的话,我所使用的NPT.mdp里面是有pbc的,但我还是不太理解如果不能在真空里模拟的话,我应该怎么构 ...

你已经用了PBC的话,就在周期性盒子里放足够多的主体和客体(应该可以用packmol等工具),使得把盒子差不多填满,然后跑一个短的NVT预平衡(可能可以省略),再跑一个NPT到完全平衡就行了
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

44

帖子

0

威望

87

eV
积分
131

Level 2 能力者

12#
 楼主 Author| 发表于 Post on 2023-6-27 22:23:06 | 只看该作者 Only view this author
wzkchem5 发表于 2023-6-27 22:17
你已经用了PBC的话,就在周期性盒子里放足够多的主体和客体(应该可以用packmol等工具),使得把盒子差不 ...

①这个时候体系里面不加溶剂的话,不算是真空体系吗?②我接触的例子少,不太清楚这样做能够得到些怎样的结果,也就是这样的目的是什么。可能有些问题比较浅薄但我确实不太理解,非常感谢老师您的耐心回复!

1万

帖子

0

威望

9888

eV
积分
22133

Level 6 (一方通行)

13#
发表于 Post on 2023-6-27 23:02:29 | 只看该作者 Only view this author
sd6512066 发表于 2023-6-27 15:23
①这个时候体系里面不加溶剂的话,不算是真空体系吗?②我接触的例子少,不太清楚这样做能够得到些怎样的 ...

(1)不算,因为跑NPT的时候真空部分会自动塌掉,变成凝聚相体系
(2)你的实际体系是无穷多个主体和客体组成的凝聚相体系,这个理解吧?但是无穷多个分子一般是没法算的,计算量是无穷大,所以必须要做个近似,就是近似认为体系具有周期性,这样就可以简化计算,让计算量变成有限大。你需要做的就是在近似认为体系有周期性的情况下,模拟这个凝聚相体系
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

44

帖子

0

威望

87

eV
积分
131

Level 2 能力者

14#
 楼主 Author| 发表于 Post on 2023-6-27 23:26:56 | 只看该作者 Only view this author
wzkchem5 发表于 2023-6-27 23:02
(1)不算,因为跑NPT的时候真空部分会自动塌掉,变成凝聚相体系
(2)你的实际体系是无穷多个主体和客 ...

好的,非常感谢老师您,但因为我的主体它具有不同的空腔,也就是配体分子能够进入的位置不止一个,所以我希望能够通过模拟来探究配体分子的一个倾向,进入哪个空腔的可能性更大、以何种方式进入,这个时候我看文献中一般都是用一个主体分子及多个配体分子来探究的,请问老师我应该怎么模拟这样的情况呢?

1万

帖子

0

威望

9888

eV
积分
22133

Level 6 (一方通行)

15#
发表于 Post on 2023-6-28 00:07:14 | 只看该作者 Only view this author
sd6512066 发表于 2023-6-27 16:26
好的,非常感谢老师您,但因为我的主体它具有不同的空腔,也就是配体分子能够进入的位置不止一个,所以我 ...

一样的,你在每个周期性盒子里放至少一个主体分子以及足够多的配体分子就行了
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

本版积分规则 Credits rule

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

GMT+8, 2026-2-24 09:52 , Processed in 0.235666 second(s), 23 queries , Gzip On.

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