The following error was obtained for the system of protein with two ligands when running the minimization step of the QM/MM.
***** LEVEL -1 WARNING FROM ***** ***** shake already active and may affect fixed atoms ****************************************** BOMLEV ( -1) IS REACHED - TERMINATING. WRNLEV IS 5
The respective .inp and .out files are attached here.Appreciate if you can give a help on this.
One possibility is that the PSF was saved with CONS FIX active for some atoms. If so, adding the following line after reading the .psf file should cure the problem--