加上隐式溶剂，并给泛函加上D3校正，如果计算出来还是这样的，说明当前结构下，质子在亚胺上才是合理的质子化配置。用GFN2-xTB级别跑了一下，确实是在亚胺上的。

截取优化后的部分片段，酚O与锌离子结合后本来结合质子能力就会降低，质子在亚胺上能与酚O和肟的N形成更稳定的双重氢键，所以这个结论是没啥问题的。

我在想，如果你对锌原子用了个诸如cc-pVQZ-pp之类的赝势基组，但是你还是把genecp写出gen，能不能正常结束任务并且保持住基本正常的结构呢？然后你就发现不了这个错误了。

gen改成genecp
6-31g(d)后头别写个莫名其妙的0

检查Zn是不是只指定了赝势基组没指定赝势。