[lammps-users] deposition simulation

Greetings. First, I'd like to say I'm new to LAMMPS but that I've already found it immensely useful. I'm a bit stuck now, though.

I'm trying to model a deposition process using the fix deposit command. However, the atoms created by the deposit command don't move. This is true regardless of what I specify for vx, vy, and vz, as well as any external forces I try to apply with the fix addforce command. I see this problem when I'm looking directly at my dumpfile (.atom) so I believe the problem is with my input script and not with my visualization procedures (lmp2cfg + atomeye + Mencoder).

I'm pasting a sample input script to the end of this email. I've played with different values for timestep, boundary, vx, vy, vz, global, local, etc, all with the same result. Is there something I'm missing here? Any advice would be greatly appreciated.

Daniel Bowring

Attached is an in.deposit test script that I know works.


in.deposit (847 Bytes)