计算化学公社

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

[GROMACS] 求助如何分析界面扩散速率

[复制链接 Copy URL]

512

帖子

0

威望

4062

eV
积分
4574

Level 6 (一方通行)

truffle

跳转到指定楼层 Go to specific reply
楼主
各位好,我想请教一下如何用GROMACS反应不同双向体系萃取过程的快慢,具体就是想得到溶质从一相扩散到另一相的速率(界面是垂直Z轴的,不知道算组成分布随时间变化能不能反映问题)。请问有什么好的方法实现吗,谢谢各位啦~
No problem is insoluble in all conceivable circumstances.

6万

帖子

99

威望

5万

eV
积分
120205

管理员

公社社长

2#
发表于 Post on 2019-6-21 10:59:33 | 只看该作者 Only view this author
自己写个脚本,比如统计从一相里有一半分子转移到另一相所需的时间来衡量
北京科音自然科学研究中心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

553

帖子

0

威望

3324

eV
积分
3877

Level 5 (御坂)

3#
发表于 Post on 2019-6-21 11:38:04 | 只看该作者 Only view this author
算Z方向的扩散系数/算溶质从一相移动到另一相的PMF

512

帖子

0

威望

4062

eV
积分
4574

Level 6 (一方通行)

truffle

4#
 楼主 Author| 发表于 Post on 2019-6-21 13:46:19 | 只看该作者 Only view this author
sobereva 发表于 2019-6-21 10:59
自己写个脚本,比如统计从一相里有一半分子转移到另一相所需的时间来衡量

谢谢Sob老师回复,请问有没有类似的脚本供参考学习一下呢
No problem is insoluble in all conceivable circumstances.

512

帖子

0

威望

4062

eV
积分
4574

Level 6 (一方通行)

truffle

5#
 楼主 Author| 发表于 Post on 2019-6-21 13:50:33 | 只看该作者 Only view this author
k64_cc 发表于 2019-6-21 11:38
算Z方向的扩散系数/算溶质从一相移动到另一相的PMF

谢谢k64兄回复,我想请教下您,算Z方向的扩散系数是要在萃取结束达到平衡的状态下吗,得到的扩散系数对于萃取过程来说有怎样的物理意义呢?谢谢~
No problem is insoluble in all conceivable circumstances.

553

帖子

0

威望

3324

eV
积分
3877

Level 5 (御坂)

6#
发表于 Post on 2019-6-21 16:38:19 | 只看该作者 Only view this author
本帖最后由 k64_cc 于 2019-6-21 16:43 编辑
naoki 发表于 2019-6-21 13:50
谢谢k64兄回复,我想请教下您,算Z方向的扩散系数是要在萃取结束达到平衡的状态下吗,得到的扩散系数对于 ...

萃取过程就是你的solute从一相移动到另一相,然后你想研究的是这个过程的动力学和平衡态性质,那么对单个solute来说,它对两相的倾向性用PMF可以很好地表示,而它在Z轴每个位置的动力学性质可以用XY面上的扩散系数反映(是的我之前打错了……)。以上性质都是平衡态的。

https://pubs.rsc.org/en/Content/ ... 01542B#!divAbstract
我们组之前的文章,仅供参考

对于多个solute我暂时还没想出什么好方法来,毕竟体系平衡的时候萃取已经结束了,对于一堆solute分子你也不能想出个好的CV来……

512

帖子

0

威望

4062

eV
积分
4574

Level 6 (一方通行)

truffle

7#
 楼主 Author| 发表于 Post on 2019-6-21 19:47:59 | 只看该作者 Only view this author
k64_cc 发表于 2019-6-21 16:38
萃取过程就是你的solute从一相移动到另一相,然后你想研究的是这个过程的动力学和平衡态性质,那么对单个 ...

收到,非常感谢~我去研究研究PMF,有不懂的来骚扰老哥啊
No problem is insoluble in all conceivable circumstances.

6万

帖子

99

威望

5万

eV
积分
120205

管理员

公社社长

8#
发表于 Post on 2019-6-21 22:50:28 | 只看该作者 Only view this author
naoki 发表于 2019-6-21 13:46
谢谢Sob老师回复,请问有没有类似的脚本供参考学习一下呢

没有现成的
自己用VMD的tcl语言写一个脚本撑死了超不过20行

还可以通过写脚本导出“在另一相的百分比” vs “模拟时间”的数据然后作图,放到文章里更为直观
北京科音自然科学研究中心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

512

帖子

0

威望

4062

eV
积分
4574

Level 6 (一方通行)

truffle

9#
 楼主 Author| 发表于 Post on 2019-6-22 01:12:37 | 只看该作者 Only view this author
sobereva 发表于 2019-6-21 22:50
没有现成的
自己用VMD的tcl语言写一个脚本撑死了超不过20行

谢谢Sob老师的指点~感觉跑完动力学,各种信息都在里面了,怎么提取更是个技术活,还得不断学习
No problem is insoluble in all conceivable circumstances.

本版积分规则 Credits rule

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

GMT+8, 2025-8-18 06:58 , Processed in 0.367973 second(s), 20 queries , Gzip On.

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