It seems there is mismatch between pqr and *param.dat file. Do you really need WCA in the calculation. According to the g_mmpbsa publication, it does not improves the correlation and significantly increases mean absolute error. I have seen this error for some system in forum, but I do not have test cases to reproduce them. Furthermore, we authors are not currently in a situation to test and modify the code due to the time constraint.