https://github.com/lanl/qmd-progress/blob/917f1242e9c20a4cdc61b3354b0ca4389a603271/src/prg_densitymatrix_mod.F90#L631C57-L631C65