Hi,
I have used the scripts as showed below to take advantage of Reax force field to run molecular dynamic simulations, but the run end up with an error like this "ia(i1,2) = 21 reax_defs.h::MBONDDEF = 20 after 709517 of 1634650 pairs completed.STOP Too many bonds on atom. Increase MBONDDEF ". I have no idea about how this comes up and how to fix it. Can anybody give me some advice? Shoud it be neccessary for me to combine commands like delete_bonds or some other commands related to my script? Or could there be something else I should pay attention to?Any sugesstion will be appreciated! Thank you in advance!
The scipt I have used is as follows:
I'll let Aidan answer this. It's been discussed on the mail list
before. It's an internal ReaxFF library error that has to
do with the memory it allocates.
I'll let Aidan answer this. It's been discussed on the mail list
before. It's an internal ReaxFF library error that has to
do with the memory it allocates.
This error occurred because at that point in the simulation, one of the
atoms had 21 bonds, whereas memory was only allocated 20 bonds. You can
increase the memory allocation by increasing MBONDDEF (as the message says).
See lib/reax/README for more info.
Thanks, Steve and Adian. But how could this happen? Why didn’t any bonds on that specific atom break? Although this is a stage for the system to reach equilibrium, it can never be physics that an atom could be connected with such a number of bonds. What is the point of such an operation? Could you please give me some spark to get through this?
Thanks again.
Sometimes a large value of mbond is due to an unphysical geometry. Other
times, it is just a spontaneous occurrence of a very low probability
geometry that will quickly dissappear. Finally, it could be how ReaxFF
handles certain chemical environments, for example, some metals (where
second and third shell neighbors are treated as bonded). How do you tell the
difference?
"ia(i1,2) = 21 reax_defs.h::MBONDDEF = 20 after
709517
of 1634650 pairs completed.STOP Too many bonds on atom. Increase
MBONDDEF ".
The anser is hinted at by the fact that you processed more than half the
bonds before encountering a problem. This suggest it was just a very low
probability geometry. Increase MBONDDEF to 30 and you should be okay.