code: def get_density(self):
          
          strc = self.get_output_structure()
          density = sum(strc.get_masses()) / strc.get_volume() * 1.660539040
          return Property(scalars=[Scalar(value=density)], units="g/(cm^3)")
docs: Compute the density from the output structure