Is it possible to apply a constraint in the charge equilibration process in ReaxFF simulation so that the charge on a particular particle type stays constant? This is very easily done in the serial ReaxFF code (using the MOLCHARGE keyword) but I don’t think it is available as part of pair_style reaxff command, or is there a way around that?
A similar questions was asked 11 years ago and the answer at that time was no.
The functionality of fix qeq/reaxff honoring the group name of the fix command and applying the charge equilibration only to atoms within the fix group has since been implemented and thus you can have charges unchanged by defining a group that doe not include those atom types and then use that group instead of “all” when defining you qeq/reaxff fix.
this does work but I noticed that when a type is excluded from qeq/reaxff it is not integrated in time at all. The atoms stay where they are.
I attached my input files to this message. I ran it once with qeq/reaxff operating on all atoms and once where it was only acting on molecular species. The atoms in the second case did not move at all. I can’t think of a physical reason for that.