lammps结构优化中遇到结构变形严重问题如何解决

lzl913821827

请求各位老师指导，我想通过lammps计算体系的热性能，首先我构建了只有两层氮化硼(BN)片的模型，设置了对应的LJ势，能够顺利的跑出结果并正常。
当我构建复合体系时，即在两层BN中间加入了聚3-己基噻吩链（五个重复单元），根据文献加入对应LJ势。发现当进行结构优化时，在低精度时（1.0e-4）结构仍然和建模时保持一样；在高精度时（1.0e-6），结构变形严重：有部分BN原子进行堆叠，并且有的BN像断开一样。根据一样（与上文提到的只有两层BN的in文件）的步骤同样可以运行完成，但是输出结果不理想。










这种情况，请问各位老师这种情况是因为哪些因素导致的，我应该从哪些参数中进行修改呢，我目前已经试过改动参数是：S元素的LJ势参数、控温系数、延长NVT时间，最终输出结果均不理想。