I am going to do md simulation of 10 receptor-ligand complex. Is the gmx_energy function is enough to decide the binding strength or I have to do gmx_mmpbsa to find free binding energy ? Is interaction energy or binding free energy are different and what they convey, is they are positively correlated?