计算化学公社

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

[GROMACS] 求助:如何在盐溶液浓度有要求的情况下加离子使体系呈电中性

[复制链接 Copy URL]

88

帖子

0

威望

1057

eV
积分
1145

Level 4 (黑子)

本帖最后由 confidence 于 2022-12-23 17:47 编辑

各位老师好,我想研究凝胶在高盐环境下对离子的吸附情况, 设定盐浓度为5 mol/L,目前尝试后建模如图 发现两个问题,想请问老师们该怎么解决呢?

1、理论上盐浓度为5 mol/L,即100个水分子和9个Na/9个Cl所适配。但是所构建的凝胶分子带负电,那么在跑md的过程中会提示电荷不为0,若我加Na离子中和体系,则会导致体系中NaCl浓度与预期设定的不符(Na会多于Cl)。请问我能通过加其他阳离子使体系呈电中性吗?或者是否可以忽略这个note?还是有其他办法呢

2、是否应该在盐溶液上端加真空层来防止周期性镜像作用使一部分盐被吸引到上端呢?(还是说通过加墙的方式)
3、对于此模拟应该也是用npt就行吧?那跑npt的时间有要求么,我用如图模型跑了1ns的npt后,再跑10ns的产生相模拟感觉两个结果就基本一样



6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2022-12-23 21:26:58 | 只看该作者 Only view this author
1 先搞清楚实际凝胶中是什么物质中和了凝胶分子的负电荷。模拟时该加什么加什么,和实际情况一致

2 如果只希望有一个界面表现吸附情况,可以加真空层

3 足够达到平衡就够了,怎么判断看
谈谈怎么判断分子动力学模拟是否达到了平衡
http://sobereva.com/627http://bbs.keinsci.com/thread-27122-1-1.html

如果Z方向加了真空层,至少Z方向不能控压,或者Z方向可压缩系数设成0
北京科音自然科学研究中心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

88

帖子

0

威望

1057

eV
积分
1145

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2022-12-25 09:53:14 | 只看该作者 Only view this author
sobereva 发表于 2022-12-23 21:26
1 先搞清楚实际凝胶中是什么物质中和了凝胶分子的负电荷。模拟时该加什么加什么,和实际情况一致

2 如果 ...

谢谢sob老师的回答,我认为凝胶分子带负电是因为带有磺酸,最后会变成SO3-,该凝胶也只与盐溶液接触,所以应该是na离子中和了该凝胶所带的负电荷吧。如果是na的话,那么像我上述所说的用na中和导致na多于cl就是正常的吗?(即此情况下我能否说盐溶液浓度是符合我初始设定值的)

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

4#
发表于 Post on 2022-12-26 00:07:19 | 只看该作者 Only view this author
confidence 发表于 2022-12-25 09:53
谢谢sob老师的回答,我认为凝胶分子带负电是因为带有磺酸,最后会变成SO3-,该凝胶也只与盐溶液接触,所 ...

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

88

帖子

0

威望

1057

eV
积分
1145

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2022-12-31 17:20:24 | 只看该作者 Only view this author

谢谢sob老师的回复,经过尝试后有了一些新的疑惑。

1、本意上是想建立一层凝胶层和盐溶液层,但是凝胶层由于凝胶分子不是很规整,无法把盒子指定部分完全填充,即我所构建的凝胶层会有一部分真空区,这种情况下我能否还是把这一层当作凝胶层呢?

2、凝胶如果不加固定的话会发生移动,我现在是直接将其冻结,这种做法应该合理吧?

3、预期结果是盐溶液中的NA离子往凝胶层方向移动而CL离子留在原地,但是现在模拟结果为NA\CL\水全都往凝胶聚集,这种情况我应该如何分析呢?(取更短时间内的变化来说明其移动趋势?)

请各位老师解答一下呢,谢谢老师们

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

6#
发表于 Post on 2022-12-31 17:31:52 | 只看该作者 Only view this author
confidence 发表于 2022-12-31 17:20
谢谢sob老师的回复,经过尝试后有了一些新的疑惑。

1、本意上是想建立一层凝胶层和盐溶液层,但是凝胶 ...

1 你先对凝胶体系跑NPT的动力学,得到纯凝胶的gro文件,再用gmx solvate对盐溶液体系扩出来的真空层部分填充凝胶,这样填得能比较满

2 没有固定的必要,该怎样就怎么样。冻结明显不合理

3 先把模型弄得当了再说别的。几个凝胶分子根本无法描述凝胶层,否则相当于把界面体系的问题弄成了溶质-溶剂体系
北京科音自然科学研究中心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

88

帖子

0

威望

1057

eV
积分
1145

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2023-1-5 19:21:58 | 只看该作者 Only view this author
sobereva 发表于 2022-12-31 17:31
1 你先对凝胶体系跑NPT的动力学,得到纯凝胶的gro文件,再用gmx solvate对盐溶液体系扩出来的真空层部分 ...

谢谢sob老师的回答,我按照您说的先对纯凝胶跑了5 ns的动力学后取其gro作为初始结构,在8*8*4 nm3的盒子中用packmol填充了26个这样的凝胶分子(再多packmol就会无法顺利完成),后构建了如图所示的模型 ,总共包含了4888个水分子,26个凝胶分子,784个NA离子以及442个CL离子。
后经过em\npt预平衡和2ns的产生相模拟,图为em后的结果 ,而2ns的产生相模拟结果与此图也几乎一致。

但我本意是想模拟出NA离子会往凝胶层移动而CL离子不动,请问老师我这里是还有什么地方改进吗?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

8#
发表于 Post on 2023-1-7 02:18:48 | 只看该作者 Only view this author
confidence 发表于 2023-1-5 19:21
谢谢sob老师的回答,我按照您说的先对纯凝胶跑了5 ns的动力学后取其gro作为初始结构,在8*8*4 nm3的盒子 ...

有经过NPT的纯物质的gro文件后,应当用gmx solvate来填空真空区而不是用packmol,前者填充得会明显更为致密

氯离子跑下去并没什么意外的,本来Na+下去了,靠静电作用,也必然会把Cl-拉下去一部分,而且Cl-自身也有扩散运动。不可能Cl-一点都不动。关注两层区域里Cl-相对浓度,如果浓度的差异显著到能说明你想说明的问题就够了。
北京科音自然科学研究中心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

88

帖子

0

威望

1057

eV
积分
1145

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2023-1-7 20:23:30 | 只看该作者 Only view this author
本帖最后由 confidence 于 2023-1-9 17:47 编辑
sobereva 发表于 2023-1-7 02:18
有经过NPT的纯物质的gro文件后,应当用gmx solvate来填空真空区而不是用packmol,前者填充得会明显更为致 ...

谢谢Sob老师的回答,1、我原本也是打算用solvate命令的,但是我在盐溶液扩充出的真空区中填充跑过NPT的纯物质gro文件时只能填充1个(命令为solvate -cp yan.gro -cs snj.gro ,其中snj.gro为纯物质的gro,yan.gro中有8*8*4nm3的真空区),所以我才改用了packmol。不知道是我哪里出错了呢?

2、如上图,我的凝胶层在跑完后会偏上移,那我还能说这是我原来设定那个凝胶层么?(一开始凝胶层设定在Z<4nm,盐在4<Z<8nm处)

3、如老师您说通过两端的离子浓度来对比,那我在计算凝胶层钠离子数时是否需要减去初始建模时为了电中性加入的na离子数呢?(按照设定浓度盐溶液有442个na和442个cl,保持电中性又多加了312个na)

4、老师我最终在凝胶下层加了低盐(此时因为浓度太低,低盐层就没加NA/CL离子),目的是想让凝胶一端处于高盐、一端处于低盐,跑了em后如图 ,发现低盐层也有部分NA离子(但理论上此时不应该有),我该怎么才能不让低盐层出现离子呢?

88

帖子

0

威望

1057

eV
积分
1145

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2023-1-10 14:48:57 | 只看该作者 Only view this author
再想请教一下老师,我跑了20ns的动力学后,提取轨迹发现,模拟过程中NA/CL离子在凝胶层/高盐环境/低盐环境中的个数基本没有发生改变,与跑了npt后的情况几乎一致。那我都看不出来这个过程中的离子迁移情况,这是哪里出了问题吗?

1万

帖子

0

威望

9867

eV
积分
22107

Level 6 (一方通行)

11#
发表于 Post on 2023-1-10 20:45:20 | 只看该作者 Only view this author
confidence 发表于 2023-1-10 07:48
再想请教一下老师,我跑了20ns的动力学后,提取轨迹发现,模拟过程中NA/CL离子在凝胶层/高盐环境/低盐环境 ...

实验上已知这个时间尺度下会发生明显的迁移吗?
可能实验只证明了在毫秒时间尺度下有迁移,你计算只能跑纳秒级的动力学,那看不出来迁移就是正常的
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
本团队长期招收研究生,有意者可私信联系

46

帖子

0

威望

197

eV
积分
243

Level 3 能力者

12#
发表于 Post on 2023-1-31 09:05:02 | 只看该作者 Only view this author
试问一下楼主,考虑过用过AIMD跑类似这样的结构嘛?

88

帖子

0

威望

1057

eV
积分
1145

Level 4 (黑子)

13#
 楼主 Author| 发表于 Post on 2023-2-21 19:29:16 | 只看该作者 Only view this author
ShuaiCui 发表于 2023-1-31 09:05
试问一下楼主,考虑过用过AIMD跑类似这样的结构嘛?

没有考虑哦

88

帖子

0

威望

1057

eV
积分
1145

Level 4 (黑子)

14#
 楼主 Author| 发表于 Post on 2023-2-21 19:37:39 | 只看该作者 Only view this author
wzkchem5 发表于 2023-1-10 20:45
实验上已知这个时间尺度下会发生明显的迁移吗?
可能实验只证明了在毫秒时间尺度下有迁移,你计算只能跑 ...

是的老师,之后我又延长了模拟时间,就观察到少部分NA/CL离子会迁移到低盐部分。在此基础上我又创建了带电荷量较少以及不带电的凝胶层。发现了一些问题:

1、跑md过程中盒子会向xy方向收缩,使得最后整个盒子都充满了水,因此无法分清NA/CL是通过凝胶层渗透到低盐部分,还是在盒子顶端跑过去的(此时z方向可压缩系数为0,知道会往xy方向收缩,但是这也缩小太多了),请问老师有什么办法解决么?
2、对于不带电的凝胶层来说,跑完发现几乎没有离子处在凝胶层中,也没有离子到达低盐部分。如果考虑渗透压的作用,那么即使凝胶层不带电,也是会有离子透过现象的出现的。所以想请教一下各位老师,gmx中是存在渗透压的么?跑出这种结果是否合理呢?

1万

帖子

0

威望

9867

eV
积分
22107

Level 6 (一方通行)

15#
发表于 Post on 2023-2-21 20:05:34 | 只看该作者 Only view this author
confidence 发表于 2023-2-21 12:37
是的老师,之后我又延长了模拟时间,就观察到少部分NA/CL离子会迁移到低盐部分。在此基础上我又创建了带 ...

1. 一般解决方法是在上层水和下层水的周期镜像之间加一层不可穿透的膜(如石墨烯),或者冻结一层水
2. 不带电的凝胶层可能碰巧透过率低,导致模拟的时间内不足以观察到离子透过的现象。渗透压纯粹是统计学的现象,因此任何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
本团队长期招收研究生,有意者可私信联系

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 12:27 , Processed in 0.278771 second(s), 23 queries , Gzip On.

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