I has a system in which I am optimizing a molecule with methanol and water molecules at B2PLYPD/6-311++(2p,d) level of theory with maxstep =10, tight convergence criteria and ultra-fine integration grid. In this system I am getting a flat potential energy surface as shown in the attached figure. I have tried scf=qc but result is same. How to rectify the problem?
I am using Gaussian 09 revision D.01