Dear all,

Is there a function which returns one of the properties accessible through `cubeprop()`

(e.g. electron density), but at any arbitrary 3D coordinate, given a wavefunction? I assume the information itself could be easily computed, given that `cubeprop()`

can return this property if you choose the grid correctly, but this of course comes along with the entire cubefile. I’m thinking of something like this:

```
```E, wfn = psi4.energy(level_of_theory, return_wfn=True)
result = pointprop(wfn, [x, y, z], tasks="density")```
```

`oeprop()`

seems to provide some of this functionality, but not for instance electron density.

Thanks and best regards,

Clemens