<div dir="ltr"><div>Hello,</div><div><br></div>In the lines 613 and 614 of the subroutine qs_ks_utils.F, key operations of the SIC correction are shown:<div><table class="highlight tab-size js-file-line-container" data-tab-size="8" style="border-collapse: collapse; box-sizing: border-box; tab-size: 8; color: rgb(36, 41, 46); font-family: -apple-system, BlinkMacSystemFont, 'Segoe UI', Helvetica, Arial, sans-serif, 'Apple Color Emoji', 'Segoe UI Emoji', 'Segoe UI Symbol'; font-size: 14px;"><tbody style="box-sizing: border-box;"><tr style="box-sizing: border-box;"><td id="LC613" class="blob-code blob-code-inner js-file-line" style="padding-right: 10px; padding-left: 10px; box-sizing: border-box; position: relative; line-height: 20px; vertical-align: top; overflow: visible; word-wrap: normal;"><div><span class="pl-k" style="font-family: SFMono-Regular, Consolas, 'Liberation Mono', Menlo, Courier, monospace; font-size: 12px; white-space: pre; box-sizing: border-box; color: rgb(215, 58, 73);">CALL</span><span style="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="font-family: SFMono-Regular, Consolas, 'Liberation Mono', Menlo, Courier, monospace; font-size: 12px; white-space: pre; box-sizing: border-box; color: rgb(215, 58, 73);">%</span><span style="font-family: SFMono-Regular, Consolas, 'Liberation Mono', Menlo, Courier, monospace; font-size: 12px; white-space: pre;">pw, ener, work_v_gspace</span><span class="pl-k" style="font-family: SFMono-Regular, Consolas, 'Liberation Mono', Menlo, Courier, monospace; font-size: 12px; white-space: pre; box-sizing: border-box; color: rgb(215, 58, 73);">%</span><span style="font-family: SFMono-Regular, Consolas, 'Liberation Mono', Menlo, Courier, monospace; font-size: 12px; white-space: pre;">pw) </span> </div><div><span style="font-family: SFMono-Regular, Consolas, 'Liberation Mono', Menlo, Courier, monospace; font-size: 12px; white-space: pre;">energy</span><span class="pl-k" style="font-family: SFMono-Regular, Consolas, 'Liberation Mono', Menlo, Courier, monospace; font-size: 12px; white-space: pre; box-sizing: border-box; color: rgb(215, 58, 73);">%</span><span style="font-family: SFMono-Regular, Consolas, 'Liberation Mono', Menlo, Courier, monospace; font-size: 12px; white-space: pre;">hartree </span><span class="pl-k" style="font-family: SFMono-Regular, Consolas, 'Liberation Mono', Menlo, Courier, monospace; font-size: 12px; white-space: pre; box-sizing: border-box; color: rgb(215, 58, 73);">=</span><span style="font-family: SFMono-Regular, Consolas, 'Liberation Mono', Menlo, Courier, monospace; font-size: 12px; white-space: pre;"> energy</span><span class="pl-k" style="font-family: SFMono-Regular, Consolas, 'Liberation Mono', Menlo, Courier, monospace; font-size: 12px; white-space: pre; box-sizing: border-box; color: rgb(215, 58, 73);">%</span><span style="font-family: SFMono-Regular, Consolas, 'Liberation Mono', Menlo, Courier, monospace; font-size: 12px; white-space: pre;">hartree</span><span class="pl-k" style="font-family: SFMono-Regular, Consolas, 'Liberation Mono', Menlo, Courier, monospace; font-size: 12px; white-space: pre; box-sizing: border-box; color: rgb(215, 58, 73);">-</span><span style="font-family: SFMono-Regular, Consolas, 'Liberation Mono', Menlo, Courier, monospace; font-size: 12px; white-space: pre;">dft_control</span><span class="pl-k" style="font-family: SFMono-Regular, Consolas, 'Liberation Mono', Menlo, Courier, monospace; font-size: 12px; white-space: pre; box-sizing: border-box; color: rgb(215, 58, 73);">%</span><span style="font-family: SFMono-Regular, Consolas, 'Liberation Mono', Menlo, Courier, monospace; font-size: 12px; white-space: pre;">sic_scaling_a</span><span class="pl-k" style="font-family: SFMono-Regular, Consolas, 'Liberation Mono', Menlo, Courier, monospace; font-size: 12px; white-space: pre; box-sizing: border-box; color: rgb(215, 58, 73);">*</span><span style="font-family: SFMono-Regular, Consolas, 'Liberation Mono', Menlo, Courier, monospace; font-size: 12px; white-space: pre;">ener</span> </div></td></tr></tbody></table></div><div>The first line calculates the electrostatic energy of the unpaired electron (self-interaction) and the second line subtracts this energy from the total Hartree energy. For open boundary calculations, this is correct, since the electrostatic energy calculated this way is the self-interaction energy. However, for periodic boundary conditions, the electrostatic energy calculated by the first line is not the self-interaction energy. Instead, the interaction energy of the unpaired electron and its periodic images is also included!! (Please correct me if I am wrong.)</div><div> So it seems that to do a correct SIC calculation for PBC systems, one needs to either add back the interaction energy of the unpaired electron and its periodic images using perhaps Ewald summation method (not sure) or correct the self-interaction term by changing a poisson solver. Is it possible to calculate the Hartree energy using non-periodic poisson solver, say multipole, for a periodic system?</div><div><br></div><div>Best,</div><div>Xiaoming<br><div><br><div><br><table class="highlight tab-size js-file-line-container" data-tab-size="8" style="box-sizing: border-box; border-collapse: collapse; tab-size: 8; color: rgb(36, 41, 46); font-family: -apple-system, BlinkMacSystemFont, 'Segoe UI', Helvetica, Arial, sans-serif, 'Apple Color Emoji', 'Segoe UI Emoji', 'Segoe UI Symbol'; font-size: 14px;"><tbody style="box-sizing: border-box;"></tbody></table></div></div></div></div>