  We investigate the rational approximation of fractional powers of unbounded
positive operators attainable with a specific integral representation of the
operator function. We provide accurate error bounds by exploiting classical
results in approximation theory involving Pad\'{e} approximants. The analysis
improves some existing results and the numerical experiments proves its
accuracy.
