I am attempting to perform a QM/MM calculation using the interface between Gromacs and Gaussian. I was only able to get a hold of the g16 binary. It has been stated that if you use the g16 binary in conjunction with Gromacs you must use the "gau" script provided at "http://wwwuser.gwdg.de/~ggroenh/qmmm.html". How do you use this properly? Has anyone successfully used it?