Mixing time of a kinetically constrained spin model on trees: power law scaling at criticality