Dear Prof. Hutter,<div><br /></div><div>Thank you very much for your reply and for the example.</div><div>Let us investigate it with my colleagues and see if we can resolve our problem.</div><div><br /></div><div>Best regards,</div><div>Aleksandros<br /><br /></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, March 10, 2023 at 3:18:28 PM UTC+1 Jürg Hutter wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hi
<br>
<br>are you assuming normalized or un-normalized Gaussians?
<br>The basis set input in CP2K uses (like all QC codes) normalized Gaussians.
<br>Internally, CP2K works with un-normalized Cartesian Gaussians, i.e. the
<br>coefficients are adapted at the beginning of the calculation.
<br>
<br>I have attached a simple example where you can play with the basis set
<br>in the input and the overlap matrix is printed.
<br>
<br>regards
<br>JH
<br>
<br>________________________________________
<br>From: <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of Aleksandros Sobczyk <<a href data-email-masked rel="nofollow">sobcz...@gmail.com</a>>
<br>Sent: Friday, March 10, 2023 2:26 PM
<br>To: cp2k
<br>Subject: [CP2K:18530] Re: Evaluation of basis functions over a grid and computation of overlap integrals
<br>
<br>Update: I have also calculated several overlap integrals analytically (for s-orbitals which are simpler), they still don't match the values of the S matrix.
<br>Any feedback would be greatly appreciated.
<br>
<br>On Wednesday, March 8, 2023 at 12:57:02 PM UTC+1 Aleksandros Sobczyk wrote:
<br>Hello,
<br>
<br>I have a set of atoms in real-space and the corresponding SZV basis sets.
<br>I want to evaluate each basis function over a grid of points in the cell.
<br>E.g., I have a grid of 3d points [r1, r2, ..., rk] and I want to evaluate each
<br>Φj(r1), Φj(r2), ... Φj(rk)
<br>As a test, I tried to numerically integrate Φj * conj(Φj) over the grid that it was evaluated, and compare the result with the corresponding entry S[j, j] of the overlap matrix that
<br>is returned by CP2K.
<br>Unfortunately my integral differs substantially from the element S[j, j], so I am doing something wrong.
<br>
<br>Can we find somewhere more detailed documentation on the precise mathematical formulation of the basis sets, and also on the specific algorithms that are used by CP2K to compute the overlap integrals?
<br>(So far I have followed as precisely as possible the following page: <a href="https://www.cp2k.org/basis_sets" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://www.cp2k.org/basis_sets&source=gmail&ust=1678546869252000&usg=AOvVaw0KT0eEBUCt1OaGLIcrQrZe">https://www.cp2k.org/basis_sets</a>
<br>but it is still missing information, e.g. are the coefficients normalized? do we assume that the spherical harmonics include the phase factor? etc.)
<br>
<br>Thanks a lot in advance!
<br>Aleksandros
<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 data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a><mailto:<a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>>.
<br>To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/fa61a2bc-2c9c-4771-8606-445d46af9148n%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/fa61a2bc-2c9c-4771-8606-445d46af9148n%2540googlegroups.com&source=gmail&ust=1678546869252000&usg=AOvVaw2eCuvUL5ARxk_h_eVFTtf8">https://groups.google.com/d/msgid/cp2k/fa61a2bc-2c9c-4771-8606-445d46af9148n%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/fa61a2bc-2c9c-4771-8606-445d46af9148n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/fa61a2bc-2c9c-4771-8606-445d46af9148n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1678546869252000&usg=AOvVaw2JoeH_FcLSy_Wc0IsSt9-F">https://groups.google.com/d/msgid/cp2k/fa61a2bc-2c9c-4771-8606-445d46af9148n%40googlegroups.com?utm_medium=email&utm_source=footer</a>>.
<br></blockquote></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:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/d971cfb2-5704-492b-9b08-aace0513065an%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/d971cfb2-5704-492b-9b08-aace0513065an%40googlegroups.com</a>.<br />