[DFTB-Plus-User] question on Ewald summation in coulomb.F90

Wei Lai wei.weilai at gmail.com
Mon Jan 7 19:53:29 CET 2019


Dear Balint,

Thanks for your reply.  If the calculation is correct, I think those
factors of 2 are for the matrix form where there is double counting.  They
were removed later when the total coulmobic energy was calculated.

Thanks, Wei

On Mon, Jan 7, 2019 at 6:12 AM Bálint Aradi <aradi at uni-bremen.de> wrote:

> Dear Wei,
>
> the calculation of the Ewald is for sure correct, as you obtain the
> correct Madelung energy with it. (I just checked it for simple cubic
> lattice).
>
> You should take into account, that ususally we only sum over the half of
> the Brillouin-zone, which must be compensated by a factor of two in the
> reciprocal part. I am not sure, where the other factor 2's are comming
> from, it could be the energy units? (Hartree vs. Rydberg?)
>
>   Best regards,
>
>   Bálint
>
> --
> Dr. Bálint Aradi
> Bremen Center for Computational Materials Science, University of Bremen
> http://www.bccms.uni-bremen.de/cms/people/b-aradi/
>
>
> _______________________________________________
> DFTB-Plus-User mailing list
> DFTB-Plus-User at mailman.zfn.uni-bremen.de
> https://mailman.zfn.uni-bremen.de/cgi-bin/mailman/listinfo/dftb-plus-user
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.zfn.uni-bremen.de/pipermail/dftb-plus-user/attachments/20190107/703b1160/attachment.html>


More information about the DFTB-Plus-User mailing list