计算化学公社

 找回密码 Forget password
 注册 Register

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

查看数: 5680 | 评论数: 19 | 收藏 Add to favorites 1
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2022-4-9 13:37

正文摘要:

模拟了IL和水二元混合物体系,想研究相平衡和界面的性质。 通过预平衡后,为什么用gmx editconf -f npt1.gro -o nptslab.gro -box 5 5 25 -center 2.5 2.5 12.5这种拉盒子后,再373.15K下跑了15000ps的NVT后水分子 ...

回复 Reply

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

我知道有用gibbs ensemble MC模拟蒸汽压的,没见过跑md模拟的。。。
qiyan 发表于 Post on 2022-5-2 09:36:26
本帖最后由 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下的蒸汽压),希望老师可以解答一下
bobosiji 发表于 Post on 2022-5-1 20:51:08
qiyan 发表于 2022-4-29 09:12
感谢老师回复,我看了一些文献,说计算获得体系的Pzz就是蒸汽压。我的疑问是,没有汽化的时候,真空区没 ...

"真空区没有水分子的话": 你算过密度分布,是严格为零么?还是随便看了几帧snapshot,没看到?
qiyan 发表于 Post on 2022-4-29 09:12:24
sobereva 发表于 2022-4-29 08:44
你先搜搜文献看别人蒸气压怎么模拟的

感谢老师回复,我看了一些文献,说计算获得体系的Pzz就是蒸汽压。我的疑问是,没有汽化的时候,真空区没有水分子的话,Pzz也能代表蒸汽压嘛
sobereva 发表于 Post on 2022-4-29 08:44:17
qiyan 发表于 2022-4-28 16:04
感谢sob老师的回复.就是我想获得水各种温度下的汽化压力和NIST对比,来看一下我设置的模拟误差有多大。像 ...

你先搜搜文献看别人蒸气压怎么模拟的
qiyan 发表于 Post on 2022-4-28 16:04:57
sobereva 发表于 2022-4-28 05:09
气体区域密度相对于液体极低。用极大的盒子模拟,真空区气体总分子数才显得比较明显

感谢sob老师的回复.就是我想获得水各种温度下的汽化压力和NIST对比,来看一下我设置的模拟误差有多大。像NIST里278.12K的蒸汽压为0.87kpa,那我做nvtslab这个温度下,体系也不会汽化啊,我获得的真空压力可以和NIST对比嘛
sobereva 发表于 Post on 2022-4-28 05:09:06
qiyan 发表于 2022-4-27 08:56
sob老师,请问一下,我如何能让体系在NVT下得到体系各种温度下汽化压力啊。例如纯水在308.18K下的蒸汽压 ...

气体区域密度相对于液体极低。用极大的盒子模拟,真空区气体总分子数才显得比较明显
qiyan 发表于 Post on 2022-4-27 08:56:41
sobereva 发表于 2022-4-17 03:13
继续提升温度再看
也可以尝试先模拟少量离子液体+大量水的情况
为了确认模拟设置得当,也最好先模拟一 ...

sob老师,请问一下,我如何能让体系在NVT下得到体系各种温度下汽化压力啊。例如纯水在308.18K下的蒸汽压为5.638kpa,那再MD中nvt下这个温度下水分子也不会汽化啊,我该怎么处理呢,希望老师有时间可以回答一下,感谢老师
qiyan 发表于 Post on 2022-4-17 11:41:07
sobereva 发表于 2022-4-17 03:13
继续提升温度再看
也可以尝试先模拟少量离子液体+大量水的情况
为了确认模拟设置得当,也最好先模拟一 ...

非常老师的回复,我会去尝试一下的
sobereva 发表于 Post on 2022-4-17 03:13:15
qiyan 发表于 2022-4-16 12:32
老师看了您写的怎么判断体系达到平衡的帖子里水分子蒸发汽化的模拟,为什么气相区有那么多水分子啊,我做 ...

继续提升温度再看
也可以尝试先模拟少量离子液体+大量水的情况
为了确认模拟设置得当,也最好先模拟一个纯水沸腾过程
qiyan 发表于 Post on 2022-4-13 19:23:16
本帖最后由 qiyan 于 2022-4-16 11:32 编辑
sobereva 发表于 2022-4-13 17:42
温度不够高本来就没法显著气化,离子液体与水之间的相互作用可能导致水比纯水更难气化
当前情况气相体积 ...

非常感谢老师回答,老师我尝试了400.15k-473.15K的温度,为什么仍然只有一个水分子在气相区啊,这样的结果有意义嘛,是做错了嘛。
sobereva 发表于 Post on 2022-4-13 17:42:36
qiyan 发表于 2022-4-13 09:19
非常感谢sob老师的回复,昨天修改后并在400.15K跑完的体系如下,只有一个水分子跑到了真空区。我想问一下 ...

温度不够高本来就没法显著气化,离子液体与水之间的相互作用可能导致水比纯水更难气化
当前情况气相体积就是真空区体积

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

GMT+8, 2026-2-24 04:36 , Processed in 0.176711 second(s), 25 queries , Gzip On.

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