[CP2K:3673] All electron calculations

[Peter Mamonov]

Hi there!

1. Try finer grid (larger cutoff). Also try adding "commensurate"
keyword into "mgrid" input section.

3. Have you multiplied the sum of electronic density over grid points
by the product of three grid step sizes? Also take in account, that
inner electrons density is pretty sharp, and you need fairly small
grid step size to get some reasonable numbers.

Regards,
Peter

On Tue, Dec 13, 2011 at 7:25 PM, Eva Pluharova <pluha... at gmail.com> wrote:
> Hi,
>
> I am trying to perform all electron calculation of water molecule with
> print out of all electron density using CP2K version 2.2.366. Input
> files see in attachment.
>
> I obtained following results that I wouldn't expect:
>
> 1. Total electronic charge (G-space):      -10.0349568912
>     Total electronic charge (R-space):      -10.0349568912
>
> Why it isn't exactly -10.000000? In other systems that I have tried
> the discrepancy is significantly higher, for example CH3CONHCH3
> -40.0000 vs. -41.758.
>
> 2. Mulliken population analysis gives the following results:
> # Atom  Element  Kind  Atomic population                Net charge
>       1     O        1           7.943086                  0.056914
>       2     H        2           1.028457                 -0.028457
>       3     H        2           1.028457                 -0.028457
>  # Total charge                 10.000000                  0.000000
>
> Which aren't in agreement with BLYP/6-31g* (identical basis set)
> calculation in Gaussian:
>      1  O   -0.742711
>      2  H    0.371355
>      3  H    0.371355
> that make more sense.
>
> 3. Total electronic density printed to cube file doesn't integrate to
> 10, but something orders of magnitude larger.
>
> Any suggestions? Thanks.
>
> Best regards,
> Eva
>
> --
> You received this message because you are subscribed to the Google Groups "cp2k" group.
> To post to this group, send email to cp... at googlegroups.com.
> To unsubscribe from this group, send email to cp2k+uns... at googlegroups.com.
> For more options, visit this group at http://groups.google.com/group/cp2k?hl=en.
>