We consider an elliptic system of equations on the torus − L 2 , L 2 d with random coefficients A, that are assumed to be coercive and stationary. Using two different approaches we obtain moment bounds on the gradient of the corrector, independent of the domain size L. In the first approach we use Green function representation. For that we require A to be locally Hölder continuous and distribution of A to satisfy Logarithmic Sobolev inequality. The second method works for non-smooth (possibly discontinuous) coefficients, and it requires that statistics of A satisfies Spectral Gap estimate. 1 Both of these articles state their results only for equations, but their methods extends also to the systems case.2 Here we assumed that A is elliptic in some sense (see discussion on different notions of ellipticity below). If this was not the case, to obtain the formula for A hom one would need to consider the equation for the corrector over multiple of cells (see work of Müller [24] for a similar result in a more general setting).