<div dir="ltr">Dear Professor Hutter,<div><br></div><div>Thank you very much for your advice! </div><div><br></div><div>I will read Carla Roeti's paper, and his book on periodic Hartree Fock: Hartree-Fock Ab Initio Treatment of Crystalline Systems, Cesare Pisani, Roberto Dovesi, Carla Roetti. </div><div><br></div><div>Best wishes,</div><div><br></div><div>Fangyong</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, May 25, 2020 at 7:35 AM <<a href="mailto:hut...@chem.uzh.ch">hut...@chem.uzh.ch</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi<br> <br> integrals are calculated in real space and then we use Fourier <br> transformation to get the matrix elements of the KS matrix in<br> reciprocal space. See the many papers on PBC implementations <br> from GAUSSIAN, Turbomole, FHI-AIMS, and of course CRYSTAL groups.<br> <br> For example<br> <br> Carla Roetti "The Crystal Code", in C.Pisani, Lecture Notes in Chemistry 67 (1996)<br> <br> regards<br> <br> Juerg Hutter<br> --------------------------------------------------------------<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 href="mailto:hut...@chem.uzh.ch" target="_blank">hut...@chem.uzh.ch</a><br> Winterthurerstrasse 190<br> CH-8057 Zürich, Switzerland<br> ---------------------------------------------------------------<br> <br> -----<a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a> wrote: -----<br> To: <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a><br> From: "Lucas Lodeiro" <br> Sent by: <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a><br> Date: 05/25/2020 10:57AM<br> Subject: Re: [CP2K:13366] a question on the Kohn-Sham orbital expressed as Gaussian basis function, is the Kohn-Sham orbital translational invariant (obey periodicity)?<br> <br> Hi Juerg,<br> <br> Reading your answer, I had a question. I understand the case at gamma point, the orbitals are traslational invariant as the local functions, due to the modulation function of Bloch function is equal to unity. But, what happened with non-gamma points? In this case the orbitals are Bloch functions, where the modulation function is a complex non-invariant one... Do you construct the invariant part of orbital with local functions and use a PW for modulation function?<br> <br> El lun., 25 may. 2020 a las 4:11, <<a href="mailto:hut...@chem.uzh.ch" target="_blank">hut...@chem.uzh.ch</a>> escribió:<br> Hi<br> <br> there is nothing special about local basis functions and PBC<br> in the GPW method. You can find the same reasonings in any <br> approach using local basis functions (Gaussians, STO, numerical).<br> <br> At the gamma point the basis functions can be considered infinite<br> sums of local functions replicated in all cells.<br> Energy is defined per unit cell and this means integrals are over<br> the unit cell only. As this is inconvenient, one tries to get them<br> converted into all space integrals over single functions only and then<br> summed over all combinations.<br> <br> Check the many papers for details.<br> <br> regards<br> <br> Juerg Hutter<br> --------------------------------------------------------------<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 href="mailto:hut...@chem.uzh.ch" target="_blank">hut...@chem.uzh.ch</a><br> Winterthurerstrasse 190<br> CH-8057 Zürich, Switzerland<br> ---------------------------------------------------------------<br> <br> -----<a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a> wrote: -----<br> To: "cp2k" <<a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>><br> From: "Fangyong Yan" <br> Sent by: <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a><br> Date: 05/23/2020 11:25AM<br> Subject: [CP2K:13360] a question on the Kohn-Sham orbital expressed as Gaussian basis function, is the Kohn-Sham orbital translational invariant (obey periodicity)?<br> <br> Dear CP2K developers,<br> <br> I have a question regarding the translational invariant of Kohn-Sham orbital in the Gaussian plane-wave method. <br> <br> The Kohn-Sham orbital should be translational invariant, and if we express the orbital using Gaussian basis function, these Gaussian basis function needs also be translational invariant, in the paper "A hybrid Gaussian and plane wave density functional scheme", by LIPPERT, HUTTER and PARRINELLO, MOLECULAR PHYSICS, 1997, VOL. 92, NO. 3, 477-487. The Kohn-Sham orbital has been expanded by Gaussian basis function which includes the periodicity. (Eq. 4 and 5 in the paper).<br> <br> However, in the recent paper of CP2K, "QUICKSTEP: Fast and accurate density functional calculations using a mixed Gaussian and plane waves approach", Joost VandeVondele, Matthias Krack , Fawzi Mohamed , Michele Parrinello , Thomas Chassaing , Jürg Hutter, Computer Physics Communications 167 (2005) 103–128, from Eq. 1, the definition of electron density, has been expanded based on Gaussian basis function, but these Gaussian basis function has not mentioned to take account of periodicity. But if the Gaussian basis function has not taken account the periodicity, the electron density is not translational invariant. <br> <br> Could you please help me with Eq. 1 regarding the periodicity of electron density in this recent paper of C2PK? <br> <br> Thank you very much!<br> <br> Fangyong<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 href="mailto:cp2k%2Bun...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/48567acf-57f7-4751-86f6-2b5b128796d6%40googlegroups.com" rel="noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/48567acf-57f7-4751-86f6-2b5b128796d6%40googlegroups.com</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 href="mailto:cp2k%2Bun...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/OF366936E3.28162C3A-ONC1258573.002CFDC1-C1258573.002CFDC2%40lotus.uzh.ch" rel="noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/OF366936E3.28162C3A-ONC1258573.002CFDC1-C1258573.002CFDC2%40lotus.uzh.ch</a>.<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 href="mailto:cp2k%2Bun...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAOFT4PLKMV8og%3D-fFP1ib9TW8Qi7MVzi4omOds90Xyo2oC5Evw%40mail.gmail.com" rel="noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/CAOFT4PLKMV8og%3D-fFP1ib9TW8Qi7MVzi4omOds90Xyo2oC5Evw%40mail.gmail.com</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 href="mailto:cp2k%2Bun...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/OFC668D3DA.627DC07D-ONC1258573.003FA34C-C1258573.003FA34D%40lotus.uzh.ch" rel="noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/OFC668D3DA.627DC07D-ONC1258573.003FA34C-C1258573.003FA34D%40lotus.uzh.ch</a>.<br> </blockquote></div>