<p>Hi all</p>

<p>I’m trying to construct two diabatic
states: NO/Au(111) state and NO-/Au(111)+ state. Firstly,
I used Becke constraint, but this method can’t get convergence during SCF
calculations. I know the Becke method has poor scaling with system size. Can I
use Becke constraint to metal surface/molecule syatem?</p>

<p>After that, I tried Hirshfeld constraint.
This method can get convergence, but the results go against the physical intuition,
especially for NO-/Au(111)+ state. Is Hirshfeld method
suitable for metal surface/molecule syatem?</p><p>I have attached the relevant files.<br></p>

<p>To compare these two methods, I calculated
a simple system: potential curve for NO-. fig1 is the CCSD(T) result, fig2 is the Becke result, fig3 is the Hirshfeld result. Becke method qualitatively reproduce the result of CCSD(T), and Hirshfeld method get the wrong result.</p><p><br></p><p>Any help in understanding this would be greatly appreciated!<br></p><div>Regards,</div><div> Gang Meng</div>