We study equidistribution properties of translations on nilmanifolds along functions of polynomial growth from a Hardy field. More precisely, if is a nilmanifold, are commuting nilrotations, and are functions of polynomial growth from a Hardy field then we show that the distribution of the sequence is governed by its projection onto the maximal factor torus, which extends Leibman's Equidistribution Criterion form polynomials to a much wider range of functions; and the orbit closure of is always a finite union of sub-nilmanifolds, which extends some of the previous work of Leibman and Frantzikinakis on this topic.