<div dir="ltr">Dear Mr. Balint<br>Thank you for your kind reply. You have mention that the repulsive energy (RP) is fitted to bulk system of gold and since crystal gold will not have charge transfer or nearly zero charge transfer, the RP fitted will be similar to the non-scc case. However, in the case of gold clusters, there are in fact non zero charge transfer among atoms. And i believe the Rp in this case will be quite different to the bulk case since  charge transfer now occurs in the former. In these two extreme cases ( bulk gold and gold clusters Au_n, where n is less than 10 ), do you think two separate RPs is necessary ? Because the Rp in Au-Au.skf is actually fitted to system with nearly zero charge transfer, yet in small clusters case, the charge transfer can't be neglected. I was speculating maybe  the Rps will be quite different in these two cases. Would you provide any comments or advice ? Many thanks.<br><br><br>                                                      With Best Regards<div>                                                      Peter<br><br><br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2017-03-24 10:39 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>
> I would like to carry out a non-scc dftb calculation. But i have one<br>
> part not clear and need your advice. Say for example, the gold case.<br>
> The SK parameter for gold in auorg folder is Au-Au.skf. Do i just<br>
> simply set the SCC=no option  in the  script for my non-scc<br>
> calculation? One part to pay attention to is that the repulsive<br>
> energy in Au-Au.skf is fitted by assuming "band structure" and<br>
> "charge fluctuation energy" existed and both terms are subtracted<br>
> with respect to DFT energy during the fitting process of the<br>
> repulsive energy. If now i just turn off " charge fluctuation energy<br>
> " by setting "SCC=no", will the repulsive energy i use in Au-Au.skf<br>
> not compatible to my non-scc calculation ? Shouldn't it be more<br>
> reasonable to have a separate repulsive energy profile for the<br>
> non-scc case ? Does the admin have this non-scc repulsive energy<br>
> profile that i can request from ?  Thank for your suggestions!<br>
<br>
</span>You are correct, if you set SCC=no, the repulsives in the Auorg set<br>
(which where made with SCC = Yes) will not be correct. But depending on<br>
the amount of charge transfer between the atoms in your system, you may<br>
still get qualitatively correct forces this way, so you could prerelax<br>
structures in non-scc mode and refine them with scc at the end. (This is<br>
mostly what non-scc is used for.)<br>
<br>
We usually do not bother to create a non-scc repulsive, as it is not<br>
really useful. Either your system has not much charge transfer (e.g. a<br>
slightly distorted Au supercell), then scc and non-scc repulsives would<br>
be very similar. Or your system has a substantial charge transfer (e.g.<br>
some defect in gold, or some molecule containing one Au atom), but then<br>
making a non-scc calculation would give a wrong Hamiltonian (one not<br>
considering this charge transfer), so fitting a repulsive with that<br>
unphysical non-scc Hamiltonain would yield a non-physical and therefore<br>
non-transferable repulsive.<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"><div dir="ltr"><div>Postdoctoral fellow, Physics department,<br>National Central University, Taiwan(R.O.C)</div></div></div>
</div>