计算化学公社

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

[GROMACS] 如何在平行板狭缝间施加压差驱动

[复制链接 Copy URL]

43

帖子

0

威望

368

eV
积分
411

Level 3 能力者

各位老师好,我想要研究在不同压差下狭缝间流体流动特性,有一些小问题想请问一下:
为了尽可能接近实验中通过压力驱动液体在狭缝间流动特性,是采用构建实体板+pull(constant force)的方式,还是采用accelerate呢?

202507051701566679..png (771.99 KB, 下载次数 Times of downloads: 106)

202507051701566679..png

6万

帖子

99

威望

5万

eV
积分
120081

管理员

公社社长

2#
发表于 Post on 2025-7-6 08:22:19 | 只看该作者 Only view this author
常外力pull和设置加速度组可以实现相同的效果,只不过目前最新版gromacs已不再支持加速度组了,只能用pull实现
北京科音自然科学研究中心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

43

帖子

0

威望

368

eV
积分
411

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2025-7-6 13:48:29 | 只看该作者 Only view this author
sobereva 发表于 2025-7-6 08:22
常外力pull和设置加速度组可以实现相同的效果,只不过目前最新版gromacs已不再支持加速度组了,只能用pull ...

sob老师,我gmx是2025.2的版本的,我在手册里面还能看到这条指令,不知道还能不能用呢

43

帖子

0

威望

368

eV
积分
411

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2025-7-6 23:54:56 | 只看该作者 Only view this author
sobereva 发表于 2025-7-6 08:22
常外力pull和设置加速度组可以实现相同的效果,只不过目前最新版gromacs已不再支持加速度组了,只能用pull ...

sob老师,我在进行了accelerate的模拟后,出现了一个不能接受的问题,我明明把氧化铝骨架部分冻结了,采用的是NVT系综(nose-hoover),我把速度调大一些,结果冻结的氧化铝骨架随离子液体流体的流动,也发生了平移【之前NVT和NPT的平衡过程中都很正常】。而后我在想是不是控温的问题,于是把氧化铝和流体部分分别控温,但是依然出现了这个问题
我的MDP中涉及到冻结和accelerate的代码如下
  1. ;Freeze
  2. freezegrps         = bone
  3. freezedim         = Y Y Y
  4. ; Acceleration
  5. acc_n_grps     = 1          ;
  6. acc_grps       = flow        ;
  7. accelerate     = 0 0.01 0    ;
复制代码

【上下两层是氧化铝】

加压.mp4

1.71 MB, 下载次数 Times of downloads:

6万

帖子

99

威望

5万

eV
积分
120081

管理员

公社社长

5#
发表于 Post on 2025-7-9 02:45:54 | 只看该作者 Only view this author
方方方 发表于 2025-7-6 23:54
sob老师,我在进行了accelerate的模拟后,出现了一个不能接受的问题,我明明把氧化铝骨架部分冻结了,采 ...

检查冻结组的原子序号确保定义正确
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2025-8-13 05:26 , Processed in 0.156016 second(s), 23 queries , Gzip On.

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