Consider the optimization problem pmin,Q:=minx∈Qp(x), where p is a degree m multivariate polynomial and Q:=[0,1]n is the hypercube. We provide explicit degree and error bounds for the sums of squares approximations of pmin,Q corresponding to the Positivstellensatz of Putinar. Our approach uses Bernstein multivariate approximation of polynomials, following the methodology of De Klerk and Laurent to provide error bounds for Schmüdgen type positivity certificates over the hypercube. We give new bounds for Putinar type representations by relating the quadratic module and the preordering associated with the polynomials gi:=xi(1−xi),i=1,…,n, describing the hypercube Q.