In this paper we study the existence of bound state solutions for stationary Schrodinger systems of the form {-Delta u + V (x) u = K(x) F-u(u, v) in R-N, -Delta u + V (x) v = K(x) F-v(u, v) in R-N, where N >= 3, V and K are bounded continuous nonnegative functions, and F(u, v) is a C-1 and p-homogeneous function with 2 < p < 2N/(N - 2). We give a special attention to the case when V may eventually vanishes. Our arguments are based on penalization techniques, variational methods and Moser iteration scheme.