I compared results of HF calculation with that of another quantum chemistry cord. They yield different density matrix although their energy are same.This difference occurs to any basis set. I suspect that the definition of the atomic orbital(normalization factor etc.) could be different. Is there any description of the definition?

I got the density matrix represented with AO via this code.

wfn.Da_subset(“AO”)