Hello

I’m pretty new in QM methods and I’m trying to perform all my calculation using psi4 from python. Becouse it’s my intention to integrate the procedures in our group molecular discovery pipeline.

I would like to generate spin density maps. It seems pretty straight forward as it’s said in the documentation. However I don’t know how to calculate ρα(r)−ρβ(r).

```
DENSITY
This task can be used to obtain the alpha and beta electron densities, ρα(r) and ρβ(r), together with the total density ρ(r)=ρα(r)+ρβ(r), and the spin density ρ(r)=ρα(r)−ρβ(r).
```

Can anyone help me with this?

Best,

E~