Bottom-up quantification of energy flexibility in cluster of residential heat pump systems