The finite field method using ab initio hybrid functional approach was applied to estimate nonlinear properties of the cluster based on paratellurite crystal structure. The good agreement was found with previously reported data.The method demonstrate high efficiency in prediction of nonlinear properties for a wide ensemble of sufficiently large polyatomic clusters.