Dear Prof. Hutter,<div><br /></div><div>I was wondering if you always need to do an all-electron calculation to reproduce the all-electron density or if you could reconstruct the all-electron density from a pseudopotential calculation using GAPW? I ask this because I used to use Quantum ESPRESSO and there you could use PAW pseudopotentials during your calculation and reconstruct the all-electron density through the projector functions, so I am wondering if this is also possible in CP2K when using the GAPW method (which is based on the PAW method) with pseudopotentials?</div><div><br /></div><div>Thanks in advance,</div><div>Léon<br /><br /></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, 10 October 2025 at 11:01:21 UTC+2 Jürg Hutter wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hi
<br>
<br>I see you are doing a GPW calculation. The DENSITY_INCLUDE keyword only applies to GAPW calculations.
<br>see
<br>DENSITY_INCLUDE: enum= TOTAL_HARD_APPROX
<br>
<br>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.
<br>
<br>If you want to have all electron densities, you need to do GAPW with POTENTIAL ALL.
<br>
<br>regards
<br>JH
<br>
<br>________________________________________
<br>From: <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of hfruchtl <<a href data-email-masked rel="nofollow">hfru...@gmail.com</a>>
<br>Sent: Monday, October 6, 2025 7:06 PM
<br>To: cp2k
<br>Subject: [CP2K:21881] Full charge density cube is zero everywhere
<br>
<br>Folks,
<br>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:
<br>
<br>-Quickstep-
<br>ELECTRON DENSITY
<br> 217 0.000000 0.000000 0.000000
<br> 122 0.234963 0.000000 0.000000
<br> 108 0.000000 -0.228920 0.000000
<br> 200 0.000000 0.000000 0.236216
<br> 8 0.000000 21.546059 -14.495211 21.850861
<br> 1 0.000000 22.840666 -13.745938 23.284968
<br> 6 0.000000 2.379002 -6.501607 19.060353
<br>...
<br>
<br>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.
<br>
<br>Thanks in advance,
<br>
<br> Herbert
<br>
<br>--
<br>You received this message because you are subscribed to the Google Groups "cp2k" group.
<br>To unsubscribe from this group and stop receiving emails from it, send an email to <a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a><mailto:<a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>>.
<br>To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/dfaf463e-5f7a-46d5-ba50-5e3cc22b424fn%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/dfaf463e-5f7a-46d5-ba50-5e3cc22b424fn%2540googlegroups.com&source=gmail&ust=1760175476379000&usg=AOvVaw2MbSVGldBRol0Utd0Ero7W">https://groups.google.com/d/msgid/cp2k/dfaf463e-5f7a-46d5-ba50-5e3cc22b424fn%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/dfaf463e-5f7a-46d5-ba50-5e3cc22b424fn%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/dfaf463e-5f7a-46d5-ba50-5e3cc22b424fn%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1760175476379000&usg=AOvVaw0disEqkd4dM-fy2f38_MpN">https://groups.google.com/d/msgid/cp2k/dfaf463e-5f7a-46d5-ba50-5e3cc22b424fn%40googlegroups.com?utm_medium=email&utm_source=footer</a>>.
<br></blockquote></div>
<p></p>
-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/98b9b314-13e2-4c62-bb4e-5a96aa5d9963n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/98b9b314-13e2-4c62-bb4e-5a96aa5d9963n%40googlegroups.com</a>.<br />