Green Kubo calculation of Ploymer

Dear Lammps users,

I used the below command to calculate the local stress of polymer.

fix SS all ave/correlate $s $p $d &

v_pxy v_pyz v_pxz type auto file profile.gk ave running

variable scale equal 1.0/t/{kb}vol$s*dt

variable vxy equal trap(f_SS[3])*${scale}

variable vyz equal trap(f_SS[4])*${scale}

variable vxz equal trap(f_SS[5])*${scale}

thermo_style custom step temp press pxy pyz pxz v_vxy v_vyz v_vxz

run 10000000

variable v equal (v_vxy+v_vyz+v_vxz)/3.0

print “running average viscosity: $v”

How can I calculate the normal stress. Would you give me some suggestion? Thank you for your help.

Best Regards,

Gerry Wong

发送自 Windows 10 版邮件应用