求助：模拟电极电解液固液界面遇到的问题

不明嚼栗子

老师们好，我想要模拟Li金属和含Li盐的电解液之间的界面。建的模型如下图，其中最长轴为z轴，Li金属层之间的空间比平衡状态的电解液体积更大，因为建模的时候和体积相同大小的空间反而塞不下电解液。把xy压缩系数设成0。一开始没有固定Li金属，电解液的体积能缩回到平衡状态附近，但是Li金属层有些变形。第二次我把Li金属层冻住再跑，结果一直崩溃。第三次我想先平衡了电解液相再在两边加上Li金属板，结果发现周期性边界上有断裂的长链的溶剂分子，没办法在两边加上金属板。现在我想请教一下：(1)第二次一直崩溃是因为冻住金属板和压缩系数设成0冲突导致的吗？我该如何在固定金属板的同时能让体积过大的电解液缩回平衡状态？(2)对于第三次，我该如何得到一个分子都完整的较平整的表面，方便在两侧加金属板，还是说这是行不通的呢？万分感谢_(:з」∠)_

xjw

您好，我最近也想做锂金属和电解液的界面模拟，但是优化的时候一直报错，请问您可以把您成功的文件发我学习下嘛

sobereva

不明嚼栗子 发表于 2021-1-28 23:48
谢谢社长！
我按第一种方法再给锂加上xy方向的限制势，但是出现报错
Fatal error:

先用-ntmpi 1不用域分解跑跑，跑平衡再说

不明嚼栗子

flwboka 发表于 2021-1-29 16:10
似乎应该建一个有很多层的Li平板模型，然后把平板的中间的一些层固定，让表面的几层可以弛豫，在平板间添加 ...

感谢回复！因为我现在想要z方向控压，而金属层横跨了z边界，想请问您，冻住金属层会不会导致z方向无法控压啊？主要是冻住和控压一起用会崩溃，但是我既不想让金属板扭曲变形也不想固定住溶液的体积，而是希望金属板在固定住形状的同时溶液能缩到平衡状态的体积

flwboka

似乎应该建一个有很多层的Li平板模型，然后把平板的中间的一些层固定，让表面的几层可以弛豫，在平板间添加电解液，然后做模拟

不明嚼栗子

sobereva 发表于 2021-1-28 17:32
还是按第一种跑，给锂都加上xy方向的限制势。还有奇怪的现象的话z两侧都加上真空区再试，锂层应该再厚一些 ...

谢谢社长！
我按第一种方法再给锂加上xy方向的限制势，但是出现报错
Fatal error:
33 particles communicated to PME rank 2 are more than 2/3 times the cut-off
out of the domain decomposition cell of their charge group in dimension x.
This usually means that your system is not well equilibrated.
这是因为节点数设置的不对吗？

sobereva

还是按第一种跑，给锂都加上xy方向的限制势。还有奇怪的现象的话z两侧都加上真空区再试，锂层应该再厚一些。