<div dir="ltr"><div>Thanks!</div><div>Changing to GPW the reference charge and population is the same, but the problem on the total charge persists... Reducing the EPS_DEFAULT value resolves the problem.</div><div> </div><div>Regards! </div></div> <div class="gmail_quote"><div dir="ltr" class="gmail_attr">El mié., 1 jul. 2020 a las 5:16, <<a href="mailto:hut...@chem.uzh.ch">hut...@chem.uzh.ch</a>> escribió: </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi the problem is in the use of the GAPW method. The Hirshfeld code only works with the (soft) density on the realspace grid. There should be a warning or a stop, my mistake. Fixing this is probably not too hard. If I find some time I will do it. No promises. regards Juerg Hutter -------------------------------------------------------------- Juerg Hutter Phone : ++41 44 635 4491 Institut für Chemie C FAX : ++41 44 635 6838 Universität Zürich E-mail: <a href="mailto:hut...@chem.uzh.ch" target="_blank">hut...@chem.uzh.ch</a> Winterthurerstrasse 190 CH-8057 Zürich, Switzerland --------------------------------------------------------------- -----<a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a> wrote: ----- To: <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a> From: "Lucas Lodeiro" Sent by: <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a> Date: 07/01/2020 03:03AM Subject: Re: [CP2K:13577] Hirshfeld charges, neutral system with no zero total Hirshfeld charge Sorry Now I attach the files (input+outputs) for isolated molecule of Urea, and the hirshfeld file for the bulk liquid system. For Bulk system calculation I use the same input, just activating the periodicity and cell size. El lun., 29 jun. 2020 a las 5:52, <<a href="mailto:hut...@chem.uzh.ch" target="_blank">hut...@chem.uzh.ch</a>> escribió: Hi there is not enough information to investigate your problem and to suggest a solution. In all my tests I get combined Hirshfeld charges < 0.0005. regards Juerg Hutter -------------------------------------------------------------- Juerg Hutter Phone : ++41 44 635 4491 Institut für Chemie C FAX : ++41 44 635 6838 Universität Zürich E-mail: <a href="mailto:hut...@chem.uzh.ch" target="_blank">hut...@chem.uzh.ch</a> Winterthurerstrasse 190 CH-8057 Zürich, Switzerland --------------------------------------------------------------- -----<a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a> wrote: ----- To: <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a> From: "Lucas Lodeiro" Sent by: <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a> Date: 06/28/2020 08:08PM Subject: [CP2K:13555] Hirshfeld charges, neutral system with no zero total Hirshfeld charge Hi all, I want to obtain Hirshfeld charges for a neutral bulk liquid system (372 atoms), but when I compute the Hirshfeld charges, the total system charge is different to zero. To obtain the charges I set: &PRINT &HIRSHFELD FILENAME SELF_CONSISTENT .TRUE. SHAPE_FUNCTION DENSITY &END HIRSHFELD &END PRINT I explore the four combinations between nonSCf / SCF and density / gaussian shape functions. The SCF does not change the behavior of total charges (just correct the atomic charges as the reference paper explains). The density shape function corrects the total charge from -1.1 to -0.167, but it still remains different from zero. I compute the total charge of an isolated molecule (the charge is different of zero too, but little) and multiply this value for the number of molecules in the bulk liquid system, and I get -0.084, this is a half the bulk total charge. I understand there is numerical "noise" and approximations (just 3 decimals), but I am not pretty sure if this is "normal" or I am missing something. I think if I use more accurate SCF convergence for Hirshfeld charges it could be corrected, but I do not find any parameter that chose it. Is what I am getting correct? Regards - Lucas Lodeiro -- You received this message because you are subscribed to the Google Groups "cp2k" group. 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>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAOFT4P%2BsvSemLPx_gL%3DL70TxoCQsS6npVdnbeqwpmHS8xkhzag%40mail.gmail.com" rel="noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/CAOFT4P%2BsvSemLPx_gL%3DL70TxoCQsS6npVdnbeqwpmHS8xkhzag%40mail.gmail.com</a>. -- You received this message because you are subscribed to the Google Groups "cp2k" group. 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>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/OFA385F851.8772786F-ONC1258596.00363672-C1258596.00363673%40lotus.uzh.ch" rel="noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/OFA385F851.8772786F-ONC1258596.00363672-C1258596.00363673%40lotus.uzh.ch</a>. -- You received this message because you are subscribed to the Google Groups "cp2k" group. 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>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAOFT4PLzXt7ByCLauu-Ztk-O7okyiifPLhSORUSnFqvOK7E2_Q%40mail.gmail.com" rel="noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/CAOFT4PLzXt7ByCLauu-Ztk-O7okyiifPLhSORUSnFqvOK7E2_Q%40mail.gmail.com</a>. [attachment "Bulk-1.hirshfeld" removed by Jürg Hutter/at/UZH] [attachment "Urea.xyz" removed by Jürg Hutter/at/UZH] [attachment "Urea-1.hirshfeld" removed by Jürg Hutter/at/UZH] [attachment "system.out" removed by Jürg Hutter/at/UZH] [attachment "input.inp" removed by Jürg Hutter/at/UZH] -- You received this message because you are subscribed to the Google Groups "cp2k" group. 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>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/OFD76EC35C.E126A33D-ONC1258598.0032F5E2-C1258598.0032F5E4%40lotus.uzh.ch" rel="noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/OFD76EC35C.E126A33D-ONC1258598.0032F5E2-C1258598.0032F5E4%40lotus.uzh.ch</a>. </blockquote></div>