I am currently interested in using pySCF to generate density matrices (first order single particle density matrices) and then using them in Molpro. In principle, this should be possible if the matrix coefficients refer to the same basis set. However, the matrop program of Molpro has some conventi...