This paper concerns the validity of bootstrap in high dimensional principal component analysis. We show that the bootstrap estimator for the distribution of the largest eigenvalue converges to Tracy Widom law. Numerical simulations reveal a substiantial improvement in approximation. To assess the accuracy of the approximation, we investigate the high order refinements in the asymptotic expansion.