Hi, dear LAMMPS colleges:
I wish to use LAMMPS to simulate a single Carbon or Silicon system.
My plan is to make the velocity of center of mass constant and in -z direction. Periodically adding perfect crystal atoms(with some initial speed and initial temperature) from the right side and if atoms reach the left side of box, they will be deleted.
However, there’re some difficulties.
I understand I can use create_atoms with the help of region command to add atoms. But I need to do this periodically, say, every 100 timesteps, I need to add 8 C-atoms. Is there a way of doing this?
To delete atoms, is there a way to delete atoms whenever they’re dropping out of the left side of the box?
The newly added atoms must undergo a langevin thermostat for a fixed amount of time before they eventually join the rest. Is this possible to define such behavior?
Thanks in advance for your help.