计算化学公社

标题: 求助:使用semiisotropic对聚合物压缩,被冻结的水墙也被压缩后模拟崩溃 [打印本页]

作者
Author:
karmie    时间: 2023-10-24 10:37
标题: 求助:使用semiisotropic对聚合物压缩,被冻结的水墙也被压缩后模拟崩溃
力场使用的是AMBER14sb-OL15 + beta-glucose GLYCAM06h,gromacs版本为2019.6。
想通过上下两层结晶纤维素CC压缩中间层的非晶纤维素AC,获得CC/AC/CC的一个三层结构,试过不加水墙,无法获得想要的结构。
将水冻结后,在压缩时,水也被压缩,且AC层密度也远不能达到我想要的,再继续压缩就出现报错。
1.请问如何能做到只让上下两层结晶纤维素CC向下压缩,而冻结的水不动呢?以下为文献原文:
the forces applied on the atoms of water and the CC regions were set to 0, with atoms only given velocities in the CC layer.
2.施加在某一部分上的力设置为0,将电荷修改为0可以做到嘛?仅给一部分施加初速度是手动在.gro文件中添加速度信息嘛?
mdp文件:
pcoupl            = berendsen
pcoupltype        = semiisotropic
tau-p             = 1.0         
compressibility   = 1.0e-5 0   
ref-p             = 1.01325 1.01325   
deform            =  0 0 -0.08 0 0 0  
nstpcouple        = -1     
refcoord-scaling  =  com
文献结构压缩快照:
(, 下载次数 Times of downloads: 15)
我的结构压缩快照:
(, 下载次数 Times of downloads: 11)





欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3