[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.htm>
More information about the DFTB-Plus-User
mailing list