[CP2K-user] CDFT applied to Au(111)/NO system

[menggang meng]

Hi all

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?

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?

I have attached the relevant files.

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.


Any help in understanding this would be greatly appreciated!
Regards,
 Gang Meng
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20210526/506fdc32/attachment.htm>