In this paper, we mainly employed the idea of the previous paper [36] to study the sharp uniform W 1,p estimates with 1 < p ≤ ∞ for more general elliptic systems with the Neumann boundary condition on a bounded C 1,η domain, arising in homogenization theory. Based on the skills developed by Z. Shen in [29] and by T. Suslina in [33,34], we also established the L 2 convergence rates on a bounded C 1,1 domain and a Lipschitz domain, respectively. Here we found a "rough" version of the first order correctors (see (1.12)), which can unify the proof in [29] and [34]. It allows us to skip the corresponding convergence results on R d that are the preconditions in [33,34]. Our results can be regarded as an extension of [24] developed by C. Kenig, F. Lin, Z. Shen, as well as of [34] investigated by respectively. Then the result of [24, Theorem 1.2] can be directly applied to (D 2 ). Due to the same reason as explained in [36], we need to derivewhich follows from the decay estimates of Neumann matrixes defined in [24, pp.916] (see Theorem 4.2), as well as ∇Ψ ε,0 C 0,σ 1 (Ω) = O(ε −σ 1 ) and ∇u ε C 0,σ 1 (Ω) = O(ε −σ 2 ) as ε → 0, which are the main conclusions of Corollary 4.6 and Lemma 4.8, where 0 < σ 1 < σ 2 < 1 are independent of ε. The above two estimates together with (1.6) guarantee that the right hand side of (D 2 ) can be uniformly bounded by the given data in Theorem 1.2. Note that the right hand side of (D 2 ) which involves div(f ) is, as a matter of fact, more general than that in [24, Theorem 1.2]. We find a simple way inspired by [31] to derive the uniform Lipschitz estimate for the weak solution u ε to L ε (u ε ) = div(f ) in Ω and ∂u ε /∂ν ε = −n · f on ∂Ω. The key ingredient is to construct the auxiliary functions {v ε,k } d k=1 , which satisfy L ε (v ε,k ) = 0 in Ω, ∂v ε,k /∂ν ε = −n k I on ∂Ω