I want to have a simple shear flow in my simulation box which means zero velocity in the middle of the simulation box.
I am using :
fix 1 all nvt/sllod temp 1.35 1.35 0.5 tchain 1
fix 2 all deform 1 xz erate 0.03 remap v
but this is producing a linear flow with the zero and the maximum velocity in the bottom and top of the simulation box respectively.
Is there any idea?