模型 这是我用moltemplate做出来的一个聚合度为50的PE链,共有72条,盒子体积64x64x350,ppf边界。初始结构虽然经过最小化,但是由于间距很小,在lammps中的弛豫过程需要很长的时间甚至会出错。 解决方案 采用soft势 先用soft势替换原先设定的lj势,达到推开原子,避免重叠的情况,同时也可以做出一个较为良好的初始结构。 pair_style soft 2.0 pair_coeff * * 50.0 推开后的效果图: 2. 定义随时间变化的wall variable z equal 311-0.001*elapsed fix 1 all nvt t