ESP charge visualization

7 views
Skip to first unread message

Andreas Döll

unread,
Jan 19, 2026, 9:42:12 AM (2 days ago) Jan 19
to cp2k
Dear cp2k community,

I want to visualize my RESP charges of my quasi 2d periodic system using the electrostatic potential. For that I the following print section

 &E_DENSITY_CUBE
                        STRIDE 1 1 1
                        FILENAME ./eden
                &END E_DENSITY_CUBE
                &V_HARTREE_CUBE
                        STRIDE 1 1 1
                        FILENAME ./hart
                &END V_HARTREE_CUBE

For visualization I use VESTA. When inspecting the resulting surface I see that potential does not look good. I want to change two things: 

- First I want to change how much the electrostatic potential around a single atom is displayed, so the range around a singular atom. I think I somehow "cut off" my considered ESP in cp2k, but I am unsure where. How can I control that?

- The potential looks as if certain atoms would belong to a positive RESP charge, but in the results, they are negative. I use the following RESP section with additional constraints for physical reasons. Especially the Carbon molecules which are bonded to the triazine molecule have a negative RESP charge, but you would not guess that if you just look at the potential. Im unsure if I fail at the visualization or if I expect too much, especially since the surface colouring goes from -9 to =0.07 which I do not understand...

  &RESP
     USE_REPEAT_METHOD
        &SPHERE_SAMPLING
     AUTO_VDW_RADII_TABLE UFF
   &END

I attach both my resp.input and output below. I cant attach the cube or eden file due to the size. If you would recommend another programm for easier visualization (or a totally different approach) please go ahead.

Sincerely, 
Andreas

epot.png
resp.inp
output.dat
Reply all
Reply to author
Forward
0 new messages