I am trying to simulate polymers between two surfaces. For a test case, I simulated a few polymer particles between two reflecting surfaces (fix wall/lj93). Polymers are homopolymers with LJ-FENE interactions.
Q: Even when I start with zero momentum initial condition, when I check the result after some timesteps (ran using Nose Hoover thermostat - which is supposed to have the momentum conserving property) the momentum in the direction perpendicular to the surface is not zero. I could guess that the momentum in the direction perpendicular is not zero since there is an unbalanced force from the surface. Question is whether fix momentum will be a good fix (or it will just cure the symptoms and not the problem?) for this or there is something else that needs to be taken care of LAMMPS which is equivalent to the force holding the surface in place?
PS: If required I can send a small datafile and an input file to generate the scenario. The datafile I am having is big to load here. Sorry about that
Vaidyanathan M S
University of Texas at Austin