Hi LAMMPS developers,

There may be an error in pair_gromacs. I advise you check it. Because I push a rigid to a bilayer and use compute_group/group command to calculate the interactions of rigid and bilayer. There exists a strange resultes of energy and force.

- I do not know the c_fff[0] is right, because manual points that " If
*c_ID[N]*is used, then N must be an index from 1-M where M is the length of the vector calculated by the compute." , But if we use c_fff[1]…C_fff[4], It give an error " the c_ID is not that long",.

2.The energy is decreased but all the forces of x/y/z is positive when pushing two group closely along z-axis. I mean if the forces is positive, the energy will increase.

- Especially, The result of our simple simulation using LAMMPS is counterintuitive.

compute fff bilayer group/group rigid

thermo_style custom step temp etotal press f_rigd[3] c_fff[0] c_fff[1] c_fff[2] c_fff[3]