I would like to use "fix deposit" with ReaxFF to add atoms to a simulation in which I deposit molecules on a surface.

I have built working scripts for my system with other force fields. However, there is an issue with the ReaxFF implementation in LAMMPS where depositing new molecules into a system will cause an error (too many ghost atoms). As per the link, the error most likely arises because the USER-REAXC package pre-allocates memory for atoms based on the initial number of atoms in the simulation box.

I do know the maximum number of atoms needed before the simulation starts. I am looking for a workaround for this issue, which would allow me to pre-allocate memory for this number of atoms, which I can later add to the simulation box.

  • 1
    $\begingroup$ +1. Welcome to our site and thank you for asking your question here! We hope to see much more of you !!! $\endgroup$ – Nike Dattani Sep 4 '20 at 18:21

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.