def METHOD_NAME(self) -> FloatQuantity:
    nonbonded_energy = 0.0 * kj_mol
    for key in ("Nonbonded", "vdW", "Electrostatics"):
        if key in self.energies is not None:
            nonbonded_energy += self.energies[key]
    return nonbonded_energy