24小时热门版块排行榜    

查看: 825  |  回复: 0

bluejh

金虫 (小有名气)

[求助] 如何在液/固混合相的npt系综中固定固相原子?

最近需要模拟一个液体/固体混合相体系,其中固体材料要保持不动,试了几个方法都不行,急求请高手指点?

首先做一个nvt系综收缩原子,模拟1ns左右后发现,边界上出现真空层,于是考虑加一个npt系综使模拟盒子收缩,但是实际情况时体系要么膨胀要么崩溃。

1. fix            fnpt1 solv npt temp 300 300 100.0 z 1.0 1.0 1000 couple none dilate solv drag 1.0

   使用上述命令,体系直接蹦了,solv是液相group

2. fix            fnpt1 solv npt temp 300 300 100.0 z 1.0 1.0 1000 couple none  drag 1.0
   去掉dilate solv,体系倒是能算,但是体积翻倍了,同时固相材料也跟着膨胀了,不知道什么原因。

这里因为是二维固体材料,所以couple none,希望固定固相材料和模拟盒子x与y两个维度,只让盒子在z方向上伸缩。请问哪里不对呢?

辛苦高手们指点一二,谢谢啦
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bluejh 的主题更新
信息提示
请填处理意见