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

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

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

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

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

似乎应该建一个有很多层的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.
这是因为节点数设置的不对吗？

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