Valence density cube file

I would like to produce a density cube file, but only with valence electrons. A bit like the cube file produced when ECP is used for core electrons.

How do I do this easiest? Can I remove core electrons from wfn Da and DB before making cubefile, or should I modify in (how)?

It would be good to comment on the following PR which may do what you need:

1 Like

And follow up question: How to pickle wfn Da and Db for later use?

That has been on a wishlist for a long time unfortunately: