@felix00 calling properties(‘ccsd’, properties=[‘dipole’], return_wfn=True) will indeed
update the DM. This way you don’t need to calculate gradients. However, if you require
the relaxed 1PDM, then you do need to use the gradients approach instead. Also, you can’t
use frozen core for calculating relaxed densities as it has not been implemented in PSI4 .