I want to simulate a polymer in shear flow. There may be three examples to create a shear flow in LAMMPS, they are in nemd,shear, flow sub-file respectively. But I don’t know which method is suitabe for my case.
What is the difference of nemd and flow except simulation method?
In the folder of \Lammps-22May08\examples, there are three folders named nemd, shear and flow. That’s my meaning. Which method is suitable for me to simulate a polyelectrolyte in shear flow, which should compute the electrostatic interaction in simulation system.
That will give you a better idea about what the simulations do and which one best matches what you’re looking for.
None of these simple examples include electrostatic interactions. You’ll have to add that. I’d recommend that you take a look at the kspace_modify command with the slab option: http://lammps.sandia.gov/doc/kspace_modify.html
One additional note. If you want long-range electrostatics, then
you can't do that in a NEMD (tilted) box. So you should
do that with something like fix viscosity which is the Muller-Plathe
method.
See the howto section of the LAMMPS manual for more discussion
on these methods.