My apologies for the long delay. (Also, please disregard the other email I sent you today privately. This email is more complete.) Here are some suggestions:
- It sounds like you want to move one plate at a constant velocity in one direction, while keeping the other plate stationary in order to exert a shear force on the liquid. There are several other ways to do this. My favorite is “fix move”:
several other fixes come to mind, although I think “fix move” is probably the best choice (because you can leave the force-field between carbon atoms unspecified):
Comment: I don’t think you should use “addforce” or “momentum” in this case, unless you have
defined a suitable force-field to keep the carbon atoms in place. (Moltemplate’s graphene example does not bother to define a force-field between the carbon atoms, but it does not allow them to move, so it does not matter in that case. I’m not sure if this is helpful, but the AIREBO force-field has been used in LAMMPS to simulate bond-breakage in carbon nanotubes, although it might require using a smaller timestep.)
Additional comment: If you decide to use “fix move”, then you should get rid of the
line which applies “fix nvt” to these boundary atoms:
fix 1 boundary nvt temp 300.0 300.0 100.0
(From the documentation: “The atoms affected by this fix should not normally be time integrated by other fixes (e.g. fix nve, fix nvt), since that will change their positions and velocities twice.”)
Also keep in mind that you should probably pick the spacing between the carbon layers, and the number of water molecules between them carefully so that the density of water is realistic at this pressure. (Probably the best way to do this is by using an NPT simulation, although we warned that this can tricky to do correctly in LAMMPS.)
You need to think more carefully about the units you are using and the timescales of your system. (Axel has a point.)
Why aren’t the water molecules moving? How long will it take for the water molecules to reach uniform shear-flow? It definitely won’t happen immediately. You only ran your simulation for 9000 time steps. (On top of that you were using a timestep of 0.01fs (which is 200 times smaller than the recommended timestep). I would not expect anything to happen on a 90fs timescale. If I remember correctly, aqueous hydrogen bonds last longer than that. Perhaps you can use a more reasonable timestep if you slow down the movement of the walls. 10 Angstroms per femto-second is ridiculously fast, compared to thermal movement at 300K.
Use “fix move” instead of “velocity”, use a more sensible velocity (I don’t know, how about 0.002 Angstroms per fs), a larger timestep (1fs or 2fs), and then run a long simulation (say 200000 timesteps) and visualize the movement of the atoms using OVITO or VMD.
I hope that helps a little bit.
It’s probably a bad idea to post a trajectory/dump file to the LAMMPS mailing list. Thankfully yours was small.