<div dir="ltr">Dear Mr. Bálint:<br>Thank you for your reply. In my opinion, I think there are two factors which might make the partial charges different for DFT and dftb+:<br><br><br>1. The population analysis method.<br><br>I have read some articles addressing that Mulliken charge can't capture the anisotropic feature of the charge. In some special systems, such as gold, it is manifested with 5d-6s hybridization and an anisotropic feature in charge may play an important rule in gold clusters, especially for predicting lowest energy structures. In order to capture this anisotropic feature, a different population method is necessary, such as Lowdin or Bader.<div><br></div><div>2. higher pole terms are neglected in charge fluctuation term.</div><div><br>    Currently, dipole, quadruple and higher order terms are not included and those direct and cross interaction terms are therefore not accounted in the charge fluctuation term, which might affect the total energy as well as the final partial charge.<br><br>It seems to me that since dftb+ always neglected some terms in DFT, the partial charge can't never be the same with DFT. As a matter of fact, some papers proposed to use a so called CM3 method to remedy such differences in partial charge predicted by DFT and dftb+. I would like to know that whether the two factors i listed above will affect the global structures predicted by dftb+ ? Thank you sinceriously.</div><div><br></div><div><br></div><div>                                                            With Best Regards</div><div>                                                             Peter               </div><div><br></div><div><br><br><br><div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">2016-10-18 15:13 GMT+08:00 Bálint Aradi <span dir="ltr"><<a href="mailto:aradi@uni-bremen.de" target="_blank">aradi@uni-bremen.de</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Peter,<br>
<span class=""><br>
> Dear DFTB+ users and developers Recently, we have encountered some<br>
> issues in calculating partial charge in small gold clusters. The<br>
> charge we obtained seem to deviate systematically with respect to<br>
> those calculated by DFT. I am wondering if a numerical scheme<br>
> proposed in this paper ( Phys. Status. solidi B, No.2 ,259-269<br>
> (2012)) is available, perhaps in the newest version of dftb+, or not<br>
> ? I would like to learn whether adding higher order terms in addition<br>
> to a mono-pole charge term can remedy the systematic deviance of<br>
> partial charge predicted by dftb+ and DFT. Any suggestions or advice<br>
> are much appreciated.Thanks.<br>
<br>
</span>Yes, we are working on the multipoles and hope to come up with a<br>
reliably working within a few months. (First results are promising, but<br>
a few peaces, like forces are still missing.)<br>
<br>
As for the partial charges, you should be careful on comparing DFT and<br>
DFTB charges. Depending on the charge partitioning scheme you use, you<br>
would have different values. In DFTB+ (and in DFTB in general) we use<br>
Mulliken-charges, which are, however, strongly basis dependent, so<br>
comparing to DFT results could be tricky. One could try to plot total<br>
charges on a grid and make Bader charge analysis for both DFT and DFTB,<br>
that would be probably more comparable, as at least you would use the<br>
same (basis-independent) charge partition scheme in both methods.<br>
<br>
  Best regards,<br>
<br>
  Bálint<br>
<span class="HOEnZb"><font color="#888888"><br>
--<br>
Dr. Bálint Aradi<br>
Bremen Center for Computational Materials Science, University of Bremen<br>
<a href="http://www.bccms.uni-bremen.de/cms/people/b_aradi/" rel="noreferrer" target="_blank">http://www.bccms.uni-bremen.<wbr>de/cms/people/b_aradi/</a><br>
<br>
<br>
</font></span><br>______________________________<wbr>_________________<br>
DFTB-Plus-User mailing list<br>
<a href="mailto:DFTB-Plus-User@mailman.zfn.uni-bremen.de">DFTB-Plus-User@mailman.zfn.<wbr>uni-bremen.de</a><br>
<a href="https://mailman.zfn.uni-bremen.de/cgi-bin/mailman/listinfo/dftb-plus-user" rel="noreferrer" target="_blank">https://mailman.zfn.uni-<wbr>bremen.de/cgi-bin/mailman/<wbr>listinfo/dftb-plus-user</a><br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">Research Assistant,Physics<br>Department,NCU,Taiwan</div>
</div>