Currently we don't overwrite __rtruediv__ for FloatWithUnit, so when do float / FloatWithUnit the unit would be lost:

from pymatgen.core.units import FloatWithUnit


x = FloatWithUnit(5, "MPa")
y = 5 / x

print(y)
print(type(y))

1.0
<class 'float'>