After solving and converging a fragment embedding calculation, the total energy can be obtained from the individual 1-and 2-PDMs of the fragments. There are two options for energy evaluation. Here, we use the cumulant PDMs, which are reported to give more accurate expectation values in energy. The cumulant form of the 2-PDM is related to the original 2-PDM via: