Preserving the Hermiticity of the one-body density matrix for a non-interacting Fermi gas

The one-body density matrix (ODM) for a zero temperature non-interacting Fermi gas can be approximately obtained in the semiclassical regime through different h-expansion techniques. One would expect that each method of approximating the ODM should yield equivalent density matrices which are both Hermitian and idempotent to any order in h. However, the Kirzhnits and Wigner-Kirkwood methods do not yield these properties, while the Grammaticos-Voros method does. Here we show explicitly, for arbitrary d-dimensions through an appropriate change into symmetric coordinates, that each method is indeed identical, Hermitian, and idempotent. This change of variables resolves the inconsistencies between the various methods, showing that the non-Hermitian and non-idempotent behavior of the Kirzhnits and Wigner-Kirkwood methods is an artifact of performing a non-symmetric truncation to the semiclassical h-expansions. Our work also provides the first explicit derivation of the d-dimensional Grammaticos-Voros ODM, originally proposed by Redjati et al (2019 J. Phys. Chem. Solids 134 313-8) based on their d=1,2,3,4 expressions.