Global geometrical optics method for vector-valued Schrödinger problems