<div dir="ltr">Thanks Matt.<div><br></div><div>Best,</div><div>Xiaoming Wang<br><br>On Friday, August 17, 2018 at 1:42:00 PM UTC-4, Matt W wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><Shameless advertisement><div><br></div><div><a href="https://aip.scitation.org/doi/abs/10.1063/1.5029818" target="_blank" rel="nofollow" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Faip.scitation.org%2Fdoi%2Fabs%2F10.1063%2F1.5029818\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNGDQ3wVCVH96IPXd1tnOXA8mNaB5w';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Faip.scitation.org%2Fdoi%2Fabs%2F10.1063%2F1.5029818\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNGDQ3wVCVH96IPXd1tnOXA8mNaB5w';return true;">https://aip.scitation.org/doi/<wbr>abs/10.1063/1.5029818</a></div><div><br></div><div>Certainly contains a fairly good literature on the subject.</div><div><br></div><div><div><\Shameless advertisement></div></div><div><br></div><div>Matt</div><div><br><br>On Friday, August 17, 2018 at 3:53:56 PM UTC+1, Xiaoming Wang wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi,<div><br></div><div>Thanks again.</div><div><br></div><div>Best,</div><div>Xiaoming Wang<br><br>On Friday, August 17, 2018 at 3:34:52 AM UTC-4, jgh wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex">Hi <br> <br>I am not sure, but think DDAPC suffers from the same problem. <br> <br>regards <br> <br>Juerg <br>------------------------------<wbr>------------------------------<wbr>-- <br>Juerg Hutter Phone : ++41 44 635 4491 <br>Institut für Chemie C FAX : ++41 44 635 6838 <br>Universität Zürich E-mail: <a rel="nofollow">hut...@chem.uzh.ch</a> <br>Winterthurerstrasse 190 <br>CH-8057 Zürich, Switzerland <br>------------------------------<wbr>------------------------------<wbr>--- <br> <br>-----<a rel="nofollow">cp...@googlegroups.com</a> wrote: ----- <br>To: "cp2k" <<a rel="nofollow">cp...@googlegroups.com</a>> <br>From: "Xiaoming Wang" <br>Sent by: <a rel="nofollow">cp...@googlegroups.com</a> <br>Date: 08/16/2018 04:48PM <br>Subject: Re: [CP2K:10640] electrostatic decoupling <br> <br>Dear Prof. Hutter, <br> <br>Thanks for your explanations and suggestion. I will read the <br>corresponding papers. Btw, does the ddapc method solve <br>this problem? There is a 'e_decpl' decoupling energy in the <br>cp_ddapc_apply_CD. Or it still suffers the same <br>problem? <br> <br>Best, <br>Xiaoming <br> <br> <br>On Thursday, August 16, 2018 at 8:12:32 AM UTC-4, jgh wrote:Hi <br> <br>your problem is related to the ill-defined energy of a charge <br>in a periodic system. <br>Your energy 1 is calculated with a background charge to nuetralize <br>the charge of your orbital. In energy 2, the isolated system, <br>no such background charge is needed. <br> <br>If you want to get some idea how to attack this problem, I would <br>suggest to read the vast literature on the calculation of <br>charged defects in solids. Start with the recent work of Pasquarello. <br> <br>regards <br> <br>Juerg <br>------------------------------<wbr>------------------------------<wbr>-- <br>Juerg Hutter Phone : ++41 44 635 4491 <br>Institut für Chemie C FAX : ++41 44 635 6838 <br>Universität Zürich E-mail: <a>hut...@chem.uzh.ch</a> <br>Winterthurerstrasse 190 <br>CH-8057 Zürich, Switzerland <br>------------------------------<wbr>------------------------------<wbr>--- <br> <br>-----<a>cp...@googlegroups.com</a> wrote: ----- <br>To: "cp2k" <<a>cp...@googlegroups.com</a>> <br>From: "Xiaoming Wang" <br>Sent by: <a>cp...@googlegroups.com</a> <br>Date: 08/15/2018 10:05PM <br>Subject: [CP2K:10631] electrostatic decoupling <br> <br>Hello, <br> <br>I'd like to decouple the the Coulomb interaction between the electron of one specific state, say HOMO, <br>and its periodic images, for a fully periodic DFT calculation. The interested charge density is localized. <br>I have tried to use different poisson solvers, say MT or WAVELET, to achieve my goal. So first I extracted <br>the the charge density from mo_coeff. Then called the poisson solver. <br> <br>pw_poisson_solve(poisson_env, orb_rho_g%pw, ener1, v_gspace1%pw) <br> <br>with poisson environment PERIODIC3D. Next I changed the poisson_env to MT0D, then called poisson <br>solver once more. <br> <br>pw_poisson_solve(poisson_env, orb_rho_g%pw, ener2, v_gspace2%pw) <br> <br>Finally, the decoupling energy is deltaE = ener1 - ener2. I thought deltaE should be a very small <br>number, because the charge density of that state is quite localized and my unit cell is big enough for <br>the MT solver. However, I got a very large deltaE 0.05 Ha. Also the value is negative, which means the <br>Hartree energy is higher for the decoupled case. I cannot understand this, because I think the image <br>interaction would increase the energy. So can anyone give some advice? <br> <br>Best, <br>Xiaoming Wang <br> -- <br> You received this message because you are subscribed to the Google Groups "cp2k" group. <br> To unsubscribe from this group and stop receiving emails from it, send an email to <a>cp2k+...@googlegroups.com</a>. <br> To post to this group, send email to <a>cp...@googlegroups.com</a>. <br> Visit this group at <a href="https://groups.google.com/group/cp2k" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/group/cp2k';return true;" onclick="this.href='https://groups.google.com/group/cp2k';return true;">https://groups.google.com/<wbr>group/cp2k</a>. <br> For more options, visit <a href="https://groups.google.com/d/optout" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>. <br> <br> <br> -- <br> You received this message because you are subscribed to the Google Groups "cp2k" group. <br> To unsubscribe from this group and stop receiving emails from it, send an email to <a rel="nofollow">cp2k+...@googlegroups.com</a>. <br> To post to this group, send email to <a rel="nofollow">cp...@googlegroups.com</a>. <br> Visit this group at <a href="https://groups.google.com/group/cp2k" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/group/cp2k';return true;" onclick="this.href='https://groups.google.com/group/cp2k';return true;">https://groups.google.com/<wbr>group/cp2k</a>. <br> For more options, visit <a href="https://groups.google.com/d/optout" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>. <br> <br></blockquote></div></div></blockquote></div></div></blockquote></div></div>