NEB calculation don't run with ReaxFF force field

Thanks for the nice small example. This is not a ReaxFF problem. Rather, you
used NEB without an atom map. The documentation tells you not to do this.
You can get the same behavior in the NEB examples by the following edit to

#atom_modify map array sort 0 0.0
atom_modify sort 0 0.0


I'm confused. NEB has error checks for not using an atom map
or having atom sorting enabled. How was the user able to run
if those settings were not correct?


The error occurs before the error check.

Just posted a patch that should flag this error instead of crashing.