Hi Amiri,
You can solve it in two ways.
- In the ‘ffield.reax’ force field file, reduce off-diagonal terms deleting the unnecessary ones.
15 ! Nr of off-diagonal terms; Ediss;Ro;gamma;rsigma;rpi;rpi2
1 2 0.1239 1.4004 9.8467 1.1210 -1.0000 -1.0000
2 3 0.0283 1.2885 10.9190 0.9215 -1.0000 -1.0000
- Increase the number of off-diagonal terms in the Reax source code and then compile it.
Hope this will work.
Best,
Sanjib