<div dir="ltr"><div>Hi all,</div><div>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.</div><div>To obtain the charges I set:</div><div>                &PRINT<br>                        &HIRSHFELD<br>                                FILENAME<br>                                SELF_CONSISTENT .TRUE.<br>                                SHAPE_FUNCTION  DENSITY<br>                        &END HIRSHFELD<br>                &END PRINT</div><div>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. <br></div><div>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.</div><div><br></div><div>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.</div><div><br></div><div><span class="gmail-tlid-translation gmail-translation" lang="en"><span title="" class="gmail-">Is what I am getting correct?</span></span></div><div><span class="gmail-tlid-translation gmail-translation" lang="en"><span title="" class="gmail-"><br></span></span></div><div><span class="gmail-tlid-translation gmail-translation" lang="en"><span title="" class="gmail-">Regards - Lucas Lodeiro<br></span></span></div></div>