计算化学公社

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

[GROMACS] 模拟IL和水这种二元混合物气液界面体系,为什么水分子不能跑到真空区,如何处理

[复制链接 Copy URL]

52

帖子

0

威望

789

eV
积分
841

Level 4 (黑子)

模拟了IL和水二元混合物体系,想研究相平衡和界面的性质。
通过预平衡后,为什么用gmx editconf -f npt1.gro -o nptslab.gro -box 5 5 25 -center 2.5 2.5 12.5这种拉盒子后,再373.15K下跑了15000ps的NVT后水分子仍然不会跑到气相区,怎么才能让水分子跑到气相区呢,不然统计z方向的密度,气相的密度不是一直都会是0吗?然后我想算逸度需要知道气相的体积。。。。。怎末获得呢?求各位大神帮助

280

帖子

0

威望

2294

eV
积分
2574

Level 5 (御坂)

20#
发表于 Post on 2022-5-12 20:46:19 | 只看该作者 Only view this author
qiyan 发表于 2022-5-2 09:36
感谢老师的回答,是的,在不产生汽化的情况下z轴密度分布在真空区是严格为0的。如果汽化就不为0了。老师 ...

我知道有用gibbs ensemble MC模拟蒸汽压的,没见过跑md模拟的。。。

52

帖子

0

威望

789

eV
积分
841

Level 4 (黑子)

19#
 楼主 Author| 发表于 Post on 2022-5-3 17:40:20 | 只看该作者 Only view this author
本帖最后由 qiyan 于 2022-5-3 17:42 编辑
sobereva 发表于 2022-4-29 08:44
你先搜搜文献看别人蒸气压怎么模拟的

十分不好意思再次打扰老师,老师我想获得像这种图片上温度的蒸汽压,但是通过模拟我发现slab后这个温度下我做的IL和水根本不会汽化,而实验是降压降到很低产生汽化得到蒸汽压,那我想问一下老师,我这种slab后nvt不能改变压力啊,他就不会汽化,那真空区的Pzz能代表气相压力吗,不能我该怎么处理呢,希望老师有时间回答一下,非常感谢老师

202205031742063203..png (126.26 KB, 下载次数 Times of downloads: 24)

202205031742063203..png

52

帖子

0

威望

789

eV
积分
841

Level 4 (黑子)

18#
 楼主 Author| 发表于 Post on 2022-5-2 09:36:26 | 只看该作者 Only view this author
本帖最后由 qiyan 于 2022-5-3 08:25 编辑
bobosiji 发表于 2022-5-1 20:51
"真空区没有水分子的话": 你算过密度分布,是严格为零么?还是随便看了几帧snapshot,没看到?

感谢老师的回答,是的,在不产生汽化的情况下z轴密度分布在真空区是严格为0的。如果汽化就不为0了。老师您能解答一下,像我做离子液体和水体系,正常400k大概0.8bar达到气液相平衡,那我做slab后nvt模拟400K,体系是不会汽化的,z轴的气相区密度为0,这时的Pzz可以代表蒸汽压嘛。或者我该怎么处理让他汽化呢(sob老师说提高温度才会汽化,但是我想要400K下的蒸汽压),希望老师可以解答一下

280

帖子

0

威望

2294

eV
积分
2574

Level 5 (御坂)

17#
发表于 Post on 2022-5-1 20:51:08 | 只看该作者 Only view this author
qiyan 发表于 2022-4-29 09:12
感谢老师回复,我看了一些文献,说计算获得体系的Pzz就是蒸汽压。我的疑问是,没有汽化的时候,真空区没 ...

"真空区没有水分子的话": 你算过密度分布,是严格为零么?还是随便看了几帧snapshot,没看到?

52

帖子

0

威望

789

eV
积分
841

Level 4 (黑子)

16#
 楼主 Author| 发表于 Post on 2022-4-29 09:12:24 | 只看该作者 Only view this author
sobereva 发表于 2022-4-29 08:44
你先搜搜文献看别人蒸气压怎么模拟的

感谢老师回复,我看了一些文献,说计算获得体系的Pzz就是蒸汽压。我的疑问是,没有汽化的时候,真空区没有水分子的话,Pzz也能代表蒸汽压嘛

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

15#
发表于 Post on 2022-4-29 08:44:17 | 只看该作者 Only view this author
qiyan 发表于 2022-4-28 16:04
感谢sob老师的回复.就是我想获得水各种温度下的汽化压力和NIST对比,来看一下我设置的模拟误差有多大。像 ...

你先搜搜文献看别人蒸气压怎么模拟的
北京科音自然科学研究中心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

52

帖子

0

威望

789

eV
积分
841

Level 4 (黑子)

14#
 楼主 Author| 发表于 Post on 2022-4-28 16:04:57 | 只看该作者 Only view this author
sobereva 发表于 2022-4-28 05:09
气体区域密度相对于液体极低。用极大的盒子模拟,真空区气体总分子数才显得比较明显

感谢sob老师的回复.就是我想获得水各种温度下的汽化压力和NIST对比,来看一下我设置的模拟误差有多大。像NIST里278.12K的蒸汽压为0.87kpa,那我做nvtslab这个温度下,体系也不会汽化啊,我获得的真空压力可以和NIST对比嘛

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

13#
发表于 Post on 2022-4-28 05:09:06 | 只看该作者 Only view this author
qiyan 发表于 2022-4-27 08:56
sob老师,请问一下,我如何能让体系在NVT下得到体系各种温度下汽化压力啊。例如纯水在308.18K下的蒸汽压 ...

气体区域密度相对于液体极低。用极大的盒子模拟,真空区气体总分子数才显得比较明显
北京科音自然科学研究中心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

52

帖子

0

威望

789

eV
积分
841

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 2022-4-27 08:56:41 | 只看该作者 Only view this author
sobereva 发表于 2022-4-17 03:13
继续提升温度再看
也可以尝试先模拟少量离子液体+大量水的情况
为了确认模拟设置得当,也最好先模拟一 ...

sob老师,请问一下,我如何能让体系在NVT下得到体系各种温度下汽化压力啊。例如纯水在308.18K下的蒸汽压为5.638kpa,那再MD中nvt下这个温度下水分子也不会汽化啊,我该怎么处理呢,希望老师有时间可以回答一下,感谢老师

52

帖子

0

威望

789

eV
积分
841

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2022-4-17 11:41:07 | 只看该作者 Only view this author
sobereva 发表于 2022-4-17 03:13
继续提升温度再看
也可以尝试先模拟少量离子液体+大量水的情况
为了确认模拟设置得当,也最好先模拟一 ...

非常老师的回复,我会去尝试一下的

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

10#
发表于 Post on 2022-4-17 03:13:15 | 只看该作者 Only view this author
qiyan 发表于 2022-4-16 12:32
老师看了您写的怎么判断体系达到平衡的帖子里水分子蒸发汽化的模拟,为什么气相区有那么多水分子啊,我做 ...

继续提升温度再看
也可以尝试先模拟少量离子液体+大量水的情况
为了确认模拟设置得当,也最好先模拟一个纯水沸腾过程
北京科音自然科学研究中心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

52

帖子

0

威望

789

eV
积分
841

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2022-4-16 12:32:57 | 只看该作者 Only view this author
sobereva 发表于 2022-4-13 17:42
温度不够高本来就没法显著气化,离子液体与水之间的相互作用可能导致水比纯水更难气化
当前情况气相体积 ...

老师看了您写的怎么判断体系达到平衡的帖子里水分子蒸发汽化的模拟,为什么气相区有那么多水分子啊,我做了375.15k-473.15k还都只有一个水分子在气相区,这是为什么呢(是5000ps时间太少了嘛),希望老师有时间可以指导一下学生,谢谢老师

5.png (110.11 KB, 下载次数 Times of downloads: 32)

5.png

52

帖子

0

威望

789

eV
积分
841

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2022-4-13 19:23:16 | 只看该作者 Only view this author
本帖最后由 qiyan 于 2022-4-16 11:32 编辑
sobereva 发表于 2022-4-13 17:42
温度不够高本来就没法显著气化,离子液体与水之间的相互作用可能导致水比纯水更难气化
当前情况气相体积 ...

非常感谢老师回答,老师我尝试了400.15k-473.15K的温度,为什么仍然只有一个水分子在气相区啊,这样的结果有意义嘛,是做错了嘛。

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

7#
发表于 Post on 2022-4-13 17:42:36 | 只看该作者 Only view this author
qiyan 发表于 2022-4-13 09:19
非常感谢sob老师的回复,昨天修改后并在400.15K跑完的体系如下,只有一个水分子跑到了真空区。我想问一下 ...

温度不够高本来就没法显著气化,离子液体与水之间的相互作用可能导致水比纯水更难气化
当前情况气相体积就是真空区体积
北京科音自然科学研究中心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

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

GMT+8, 2026-2-24 05:58 , Processed in 0.179909 second(s), 24 queries , Gzip On.

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