如何正确进行聚合物的应力松弛模拟？

weixi

1.体系：粗粒化粒子的聚合物交联网络，约7个粒子一个交联点，常用的Fene成键势(rmax=1.5,k=30)和LJ（1,1,1)；

2.软件：pygamd
2.目的：想计算松弛时间，看文献和教程常使用应力松弛，拉伸（约5%）后固定，拟合应力与时间的公式 δ=δ0exp（-t/τ)来计算松弛时间。
3.问题：搜索看到常用拉伸方式有
        a,改变盒子尺寸实现单轴拉伸（如lammps的deform）
        b,给盒子某一侧粒子施加力然后拉伸
        c,用NPT系综但压力向异性，（如Px=-1，Py=Pz=0）模拟x方向进行拉伸
       
Q1：请问对于我这个聚合物网络体系，应该用哪一种更适合？
Q2：我用(a,c)两种方法拉伸结束后收集应力信息时，在NVT系综下，得到的Press-x起步下降的非常快（如图），无法用 δ/δ0 = exp（-t/τ)来拟合，
       求教我的模拟和分析可能哪一步出了问题？


（图横坐标是 time，纵坐标是应力比值 δ/δ0，黑点是模拟收集的Press-x数据，红线是拟合的指数曲线 ）

admin

Q1应该都可以abc都可以，你只要跑出你需要的结果就行。但是个人更建议用a和c，b应该需要一些特定的体系下才能做，因为就像你发力在一个点上进行测试，我也不太熟，个人理解
Q2感觉差不多，看那个符合你的实际效果吧