GROMACS no longer supports cutoffs set to zero, as it required the group nonbonded scheme. The only way to do these calculations in current versions is to set a very large box, use very large cutoffs (sufficient to encompass all nonbonded pairs) and apply an NVT ensemble.