Derivative expansion for computing critical exponents of O(N) symmetric models at NNLO

We apply the derivative expansion of the effective action in the exact renormalization group equation up to fourth order to the $Z_2$ and $O(N)$ symmetric scalar models in $d=3$ Euclidean dimensions. We compute the critical exponents $\nu$, $\eta$ and $\omega$ using polynomial expansion in the field. We obtain our predictions for the exponents employing two regulators widely used in ERG computations. We apply Wynn's epsilon algorithm to improve the predictions for the critical exponents, extrapolating beyond the next-to-next-to-leading order prediction of the derivative expansion.