如下经验上的建议：
1、ALGO换成N。
2、LDIPOL=T和IDIPOL=3。
这两个会导致难收敛，额外适当调整AMIX和BMIX（通常是降低减小震荡）。如果不行真空下面先做个波函数然后读取进去。额外可能需要考虑PREC=A或适当提升ENAUG，PB方程对这个会有点敏感。

可以给出完整的OSZICAR吗，C3N4和负载了单原子的都要

因为看起来还是有收敛的趋势的，所以可以通过调整NELM参数来多跑几步看看能否收敛。
或者算一个不加溶剂化的单点写一个WAVECAR，然后加溶剂化计算的时候读取这个WAVECAR，一般也会加快SCF收敛。