This paper is concerned with the critical quasilinear Schrödinger systems in R N : −∆w + (λa(x) + 1)w − (∆|w| 2)w = p p+q |w| p−2 w|z| q + α α+β |w| α−2 w|z| β −∆z + (λb(x) + 1)z − (∆|z| 2)z = q p+q |w| p |z| q−2 z + β α+β |w| α |z| β−2 z, where λ > 0 is a parameter, p > 2, q > 2, α > 2, β > 2, 2 • (2 * − 1) < p + q < 2 • 2 * and α + β = 2 • 2 *. By using variational method , we prove the existence of positive ground state solutions which localize near the set Ω = int a −1 (0) ∩ int b −1 (0) for λ large enough.