Inconsistency in molecular dipole calculation

Daniele Passerone dpass... at gmail.com
Fri Jun 5 14:10:31 UTC 2009


Just for the sake of completeness:

exactly the same input file as before, but with the following cell:

    A  24.22 13.43   00
    B  0.0  25.80  0.0
    C  0.0  00.00  30.0


gives rise to the following DIPOLE

MOLECULAR=TOTAL

 17.948087   -2.330780    0.923091


So we have a problem with non-orthorombic cells


On 5 Giu, 15:58, Daniele Passerone <dpass... at gmail.com> wrote:
> On 5 Giu, 14:08, Teodoro Laino <teodor... at gmail.com> wrote:
>
> > The way it is computed is similar to TOTAL_DIPOLE. This means that you
> > should get big problems as well
> > with the TOTAL_DIPOLE..
>
> > There could be a bug somewhere both in TOTAL_DIPOLE and
> > MOLECULAR_DIPOLES for non-orthorombic cells..
> > Is this the case?
>
> Exactly.


More information about the CP2K-user mailing list