<div dir="ltr"><pre style="color:rgb(0,0,0)">Dear pwscf users,</pre><pre style=""><font color="#000000">I want to calculate charged system binding energy.</font></pre><pre style=""><font color="#000000">I have been trying to calculate slab + molecule binding energy. For non vdw(van der waals) and vdw interaction</font><span style="font-family:arial,sans-serif"><font color="#000000">(input_dft='vdw-df2')</font></span><span style="font-family:arial,sans-serif"><font color="#000000">, I got good results. The binding energy with vdw is about 0.62 ev . But </font></span><span style="font-family:arial,sans-serif"><font color="#000000">When I gave </font></span><font color="#000000" style="font-family:arial,sans-serif">one additional electron</font><span style="font-family:arial,sans-serif"><font color="#000000"> charge to the system with tot_charge=-1 , the binding energy with vdw in charged system became about 8 eV.  ( By the way, I applied </font></span><span style="font-family:arial,sans-serif"><font color="#000000">tot_charge=-1 code to slab and molecule </font></span><font color="#000000" face="arial, sans-serif">individually)</font></pre><pre style=""><font color="#000000" face="arial, sans-serif"><br></font></pre><pre style=""><font color="#000000" face="arial, sans-serif">Is this energy wrong? Because, also in one electron missing system, the result is the same. </font><span style="color:rgb(0,0,0);font-family:arial,sans-serif"> </span></pre><pre style=""><span style="color:rgb(0,0,0);font-family:arial,sans-serif"><br></span></pre><pre style=""><span style="color:rgb(0,0,0);font-family:arial,sans-serif">Is there any mistake which I made? or What can I do when I want to calculate charged system binding energy with vdw? </span></pre><pre style="color:rgb(0,0,0)"><br></pre><pre style="color:rgb(0,0,0)"><pre>Best regards,
Bahadir</pre></pre></div>