Thank you for your kind reply. I was wondering if we can have a big enough simulation cell and perform Bader charges for some atoms in the middle, would that work?
It may work, you can have a try.
]]>PS: Bader charge is not a good choice, it is not only very time-consuming but also poor. You may consider calculating 1.2*CM5 charge using Multiwfn based on .molden file of CP2K.
]]>