[CP2K-user] [CP2K:21886] Full charge density cube is zero everywhere

[Jürg Hutter]

Hi

I see you are doing a GPW calculation. The DENSITY_INCLUDE keyword only applies to GAPW calculations.
see
DENSITY_INCLUDE: enum= TOTAL_HARD_APPROX

Which parts of the density to include. In GAPW the electronic density is divided into a hard and a soft component, and the default (TOTAL_HARD_APPROX) is to approximate the hard density as a spherical gaussian and to print the smooth density accurately. This avoids potential artefacts originating from the hard density. If the TOTAL_DENSITY keyword is used the hard density will be computed more accurately but may introduce non-physical features. The SOFT_DENSITY keyword will lead to only the soft density being printed. In GPW these options have no effect and the cube file will only contain the valence electron density.

If you want to have all electron densities, you need to do GAPW with POTENTIAL ALL.

regards
JH

________________________________________
From: cp2k at googlegroups.com <cp2k at googlegroups.com> on behalf of hfruchtl <hfruchtl at gmail.com>
Sent: Monday, October 6, 2025 7:06 PM
To: cp2k
Subject: [CP2K:21881] Full charge density cube is zero everywhere

Folks,
I have a 2-dimensionally periodic system and want to create a cube file for Bader analysis (with the Henkelman program). I am using the TOTAL_HARD_APPROX keyword to take care of the core charges. But the cube file I get contains no charges. The first lines are:

-Quickstep-
ELECTRON DENSITY
  217    0.000000    0.000000    0.000000
  122    0.234963    0.000000    0.000000
  108    0.000000   -0.228920    0.000000
  200    0.000000    0.000000    0.236216
    8    0.000000   21.546059  -14.495211   21.850861
    1    0.000000   22.840666  -13.745938   23.284968
    6    0.000000    2.379002   -6.501607   19.060353
...

and when it comes to the values on the grid after the atom information, all vales are zero. Am I doing something wrong, or does this not work for periodic systems? I attach the input file for a smallish model system. The version number is CP2K version 2025.1.

Thanks in advance,

  Herbert

--
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com<mailto:cp2k+unsubscribe at googlegroups.com>.
To view this discussion visit https://groups.google.com/d/msgid/cp2k/dfaf463e-5f7a-46d5-ba50-5e3cc22b424fn%40googlegroups.com<https://groups.google.com/d/msgid/cp2k/dfaf463e-5f7a-46d5-ba50-5e3cc22b424fn%40googlegroups.com?utm_medium=email&utm_source=footer>.

-- 
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/cp2k/ZR0P278MB0759E64CFF69040B98F980E49FEFA%40ZR0P278MB0759.CHEP278.PROD.OUTLOOK.COM.