<div dir="ltr">Dear Thomas,<div><br></div><div>Thanks for your explanation! Yes, periodicity has to be kept for condensed phase calculations.</div><div><br></div><div>Best,</div><div><br></div><div>Fangyong</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Dec 14, 2020 at 5:12 AM Thomas Kühne <<a href="mailto:tku...@gmail.com">tku...@gmail.com</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"><div style="overflow-wrap: break-word;">Dear Fangyong, <div><br></div><div>within the GPW scheme, the orbitals are described by contracted Gaussian functions, </div><div>whereas electron density is expanded in plane-waves. You are right that the extent of </div><div>the Gaussian functions is finite, which is controlled by eps values that are eventually </div><div>related to EPS_DEFAULT. You are also right wrt to the plane-waves, whose truncation </div><div>is controlled by CUTOFF. However, contrary to the cutoff of conventional PW codes, </div><div>the latter differs by a factor of 4, which due to the fact that the density is the sum of the </div><div>orbitals squared combined with the application of the Shannon-Nyquist sampling theorem. </div><div>Lastly, within the GAPW approach, the Gaussian functions are indeed also periodic </div><div>Gaussians. </div><div><br></div><div>Best, </div><div>Thomas</div><div><div><br><blockquote type="cite"><div>Am 08.12.2020 um 06:22 schrieb Fangyong Yan <<a href="mailto:fyya...@gmail.com" target="_blank">fyya...@gmail.com</a>>:</div><br><div><div dir="ltr" style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none">Hi, Professor Hutter,<div><br></div><div>I read your paper, and I think when you calculate the electron density based on periodic Gaussian, you need to truncate the periodic Gaussian because it is an infinite series, but it decays fast since the periodic Gaussian is a localized basis set. Am I understanding correctly? </div><div><br></div><div>For electron density represented by planewaves, it also has the same problem, it is an infinite series, and has to be truncated, which is the planewave cutoff. </div><div><br></div><div>Another question, in the Gaussian augmented planewave method, the atomic hard and soft electron density, they are atom-based electron density, and are represented by Gaussian functions, they are also periodic Gaussian, in order to fulfill the periodicity, am I correct? Thanks!</div><div><br></div><div>Fangyong</div><div><br></div><div><br></div><div><br></div></div><br style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><div class="gmail_quote" style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><div dir="ltr" class="gmail_attr">On Fri, Dec 4, 2020 at 9:49 AM Fangyong Yan <<a href="mailto:fyya...@gmail.com" target="_blank">fyya...@gmail.com</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"><div dir="ltr"><div>Dear Professor Hutter, </div><div><br></div>I will read the paper. Thanks!<div><br></div><div>Fangyong</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Dec 4, 2020 at 4:37 AM <<a href="mailto:hut...@chem.uzh.ch" target="_blank">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>yes, it is periodic Gaussians. For details of the calculation<br>read the publications:<br><br> Lippert, G; Hutter, J; Parrinello, M.<br> MOLECULAR PHYSICS, 92 (3), 477-487 (1997).<br> A hybrid Gaussian and plane wave density functional scheme.<br> <a href="https://doi.org/10.1080/002689797170220" rel="noreferrer" target="_blank">https://doi.org/10.1080/002689797170220</a><br><br> VandeVondele, J; Krack, M; Mohamed, F; Parrinello, M; Chassaing, T;<br> Hutter, J. COMPUTER PHYSICS COMMUNICATIONS, 167 (2), 103-128 (2005).<br> QUICKSTEP: Fast and accurate density functional calculations using a<br> mixed Gaussian and plane waves approach.<br> <a href="https://doi.org/10.1016/j.cpc.2004.12.014" rel="noreferrer" target="_blank">https://doi.org/10.1016/j.cpc.2004.12.014</a><br><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:<span> </span><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><span> </span>wrote: -----<br>To: "cp2k" <<a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>><br>From: "<a href="mailto:fy...@gmail.com" target="_blank">fy...@gmail.com</a>"<span> </span><br>Sent by:<span> </span><a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a><br>Date: 12/04/2020 08:42AM<br>Subject: [CP2K:14345] how does CP2K construct electron density in real-space FFT mesh point? Is the Gaussian product a periodic Gaussian product?<br><br>Hi,<span> </span><br><br>The electron density in real-space FFT mesh point is defined as:<br><br>n(Ri) = summation of m, n {P_mn * Gaussian_mn(Ri),<span> </span><br><br>Ri is the grid point, m and n are the index for the Gaussian basis set, P_mn is the density matrix, Gaussian_mn(Ri) is the product of Gaussian basis set m and Gaussian basis set n,<br><br>my question is:<br>1)<span> </span><br>is Gaussian_mn(Ri) a periodic Gaussian? (Product Gaussian is also a Gaussian).<span> </span><br>That is,<br><br>Gaussian_mn(Ri) = Gaussian_mn(Ri+L_unitcell), where L_unitcell is the length for the unit cell in either 1,2 or 3 dimension.<span> </span><br>2) if it is the periodic Gaussian, and it should be in order to fulfill the periodicity of electron density, how do you do the calculation?<span> </span><br>Because the periodic Gaussian is actually an infinite series, that is,<br>Gaussian(Ri) = summation of L (Gaussian(Ri+L), where L = n*L_unitcell, and n is an integer and span from minus infinite to positive infinite.<br><br>Thanks!<br><br>Fangyong<br><br><br><br><br><br><br> <span> </span>--<span> </span><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<span> </span><a href="mailto:cp2k%2Bun...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br> To view this discussion on the web visit<span> </span><a href="https://groups.google.com/d/msgid/cp2k/e023de64-1040-49bd-a228-2aa4cf0732a5n%40googlegroups.com" rel="noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/e023de64-1040-49bd-a228-2aa4cf0732a5n%40googlegroups.com</a>.<br><br><br>--<span> </span><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<span> </span><a href="mailto:cp2k%2Bun...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>To view this discussion on the web visit<span> </span><a href="https://groups.google.com/d/msgid/cp2k/OFE157B452.F6F180E1-ONC1258634.0034CF74-C1258634.0034DAA2%40lotus.uzh.ch" rel="noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/OFE157B452.F6F180E1-ONC1258634.0034CF74-C1258634.0034DAA2%40lotus.uzh.ch</a>.<br></blockquote></div></blockquote></div><div style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><br></div><span style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;float:none;display:inline">--<span> </span></span><br style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><span style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;float:none;display:inline">You received this message because you are subscribed to the Google Groups "cp2k" group.</span><br style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><span style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;float:none;display:inline">To unsubscribe from this group and stop receiving emails from it, send an email to<span> </span></span><a href="mailto:cp...@googlegroups.com" style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px" target="_blank">cp...@googlegroups.com</a><span style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;float:none;display:inline">.</span><br style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><span style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;float:none;display:inline">To view this discussion on the web visit<span> </span></span><a href="https://groups.google.com/d/msgid/cp2k/CAJjmc%2BJyp%3DM1Wmu6TResTf%3Dw3oNkcZB%2B%2BPL10Z6aURJ_x%2BoKfg%40mail.gmail.com?utm_medium=email&utm_source=footer" style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px" target="_blank">https://groups.google.com/d/msgid/cp2k/CAJjmc%2BJyp%3DM1Wmu6TResTf%3Dw3oNkcZB%2B%2BPL10Z6aURJ_x%2BoKfg%40mail.gmail.com</a><span style="font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;float:none;display:inline">.</span></div></blockquote></div><br></div><br><br><div>
<div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">==============================</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Thomas D. Kühne</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Dynamics of Condensed Matter</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Chair of Theoretical Chemistry</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">University of Paderborn</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Warburger Str. 100</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">D-33098 Paderborn</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Germany</div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><a href="mailto:tdku...@mail.upb.de" target="_blank">tdku...@mail.upb.de</a></div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-east-asian:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">+49/(0)5251/60-5726</div></div></div></div>
</div>
<br></div>
<p></p>
-- <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:cp...@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/83C07016-CE0A-4E58-9D3B-DDD18511A546%40gmail.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/83C07016-CE0A-4E58-9D3B-DDD18511A546%40gmail.com</a>.<br>
</blockquote></div>