<div dir="ltr">Hello,<div><br></div><div>I'd like to decouple the the Coulomb interaction between the electron of one specific state, say HOMO, </div><div>and its periodic images, for a fully periodic DFT calculation. The interested charge density is localized.</div><div>I have tried to use different poisson solvers, say MT or WAVELET, to achieve my goal. So first I extracted </div><div>the the charge density from mo_coeff. Then called the poisson solver.</div><div><span style="color: rgb(36, 41, 46); font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, Courier, monospace; font-size: 12px; white-space: pre;"><br></span></div><div><span style="color: rgb(36, 41, 46); font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, Courier, monospace; font-size: 12px; white-space: pre;">pw_poisson_solve(poisson_env, orb_rho_g</span><span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73); font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, Courier, monospace; font-size: 12px; white-space: pre;">%</span><span style="color: rgb(36, 41, 46); font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, Courier, monospace; font-size: 12px; white-space: pre;">pw, ener1, v_gspace1</span><span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73); font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, Courier, monospace; font-size: 12px; white-space: pre;">%</span><span style="color: rgb(36, 41, 46); font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, Courier, monospace; font-size: 12px; white-space: pre;">pw) </span><br></div><div><br></div><div>with poisson environment PERIODIC3D. Next I changed the poisson_env to MT0D, then called  poisson</div><div>solver once more.</div><div><br></div><div><span style="color: rgb(36, 41, 46); font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, Courier, monospace; font-size: 12px; white-space: pre;">pw_poisson_solve(poisson_env, orb_rho_g</span><span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73); font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, Courier, monospace; font-size: 12px; white-space: pre;">%</span><span style="color: rgb(36, 41, 46); font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, Courier, monospace; font-size: 12px; white-space: pre;">pw, ener2, v_gspace2</span><span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73); font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, Courier, monospace; font-size: 12px; white-space: pre;">%</span><span style="color: rgb(36, 41, 46); font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, Courier, monospace; font-size: 12px; white-space: pre;">pw)</span><br></div><div><span style="color: rgb(36, 41, 46); font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, Courier, monospace; font-size: 12px; white-space: pre;"><br></span></div><div>Finally, the decoupling energy is deltaE = ener1 - ener2. I thought deltaE should be a very small</div><div>number, because the charge density of that state is quite localized and my unit cell is big enough for</div><div>the MT solver. However, I got a very large deltaE 0.05 Ha. Also the value is negative, which means the </div><div>Hartree energy is higher for the decoupled case. I cannot understand this, because I think the image </div><div>interaction would increase the energy. So can anyone give some advice?</div><div><br></div><div>Best,</div><div>Xiaoming Wang</div></div>