请问这个问题最后解决了吗？我也是开始结冰速度比较快（10ns内），但到后面剩下一半的水一直都无法结冰，从10ns-100ns一直没有明显变化。

博主这个纯冰水结冰模拟文件能分享一下嘛，万分感谢！

sobereva 发表于 2023-10-28 03:52
1 semiisotropic不是给膜体系唯一专用的

2、3 半各项同性控压时你只需设两个值，分别负责XY和Z（你的 ...

好的，谢谢老师

牧生 发表于 2023-10-27 19:46
社长这一页ppt好像是这次培训的新增内容。有个疑问，为什么一定要“应当用半各项同性控压，在垂直于和平行 ...

原理上不是最好的
在完全结冰之前，体系是固-液界面体系，平行于界面方向和垂直于界面方向的特征很不同，不应各方向约束盒子尺寸同步等比例变化

非好汉 发表于 2023-10-27 18:49
好的，老师。还有几个疑问，麻烦老师看一下。
1.体系中没有位置限制也没有冻结，也不是膜，可以直接用se ...

1 semiisotropic不是给膜体系唯一专用的

2、3 半各项同性控压时你只需设两个值，分别负责XY和Z（你的模型得让冰种表面平行于XY）。可压缩系数都设为水的就行，1 atm


Pcoupl     = Berendsen
pcoupltype = semiisotropic
tau_p = 0.5
ref_p = 1.01325 1.01325
compressibility = 4.5e-5 4.5e-5

非好汉 发表于 2023-10-27 18:49
好的，老师。还有几个疑问，麻烦老师看一下。
1.体系中没有位置限制也没有冻结，也不是膜，可以直接用se ...

社长这一页ppt好像是这次培训的新增内容。有个疑问，为什么一定要“应当用半各项同性控压，在垂直于和平行于冰种界面的方向要分别控压”
我做结冰，控压直接用的这个，感觉也没觉得有什么问题啊？

Pcoupl     = parrinello-rahman
pcoupltype = isotropic
tau_p = 2.0
ref_p = 20.0
compressibility = 4.5e-5

sobereva 发表于 2023-10-27 08:17
当前不正常。参考北京科音分子动力学与GROMACS培训班（http://www.keinsci.com/workshop/KGMX_content.html ...

好的，老师。还有几个疑问，麻烦老师看一下。
1.体系中没有位置限制也没有冻结，也不是膜，可以直接用semiisotropic吗？
2.然后三个方向上的ref_p都设为1bar，平行冰种方向yz方向上的压缩系数设为4.5e-5，垂直冰种方向x方向设为0吗？
3.以下写法可以吗？
compressibility = 0 4e-5 4e-5
ref_p = 1 1 1

当前不正常。参考北京科音分子动力学与GROMACS培训班（http://www.keinsci.com/workshop/KGMX_content.html）里我在260K下用TIP4P/Ice做结冰的模拟：



注意控压方式，应当用半各项同性控压，在垂直于和平行于冰种界面的方向要分别控压