In the limit ℏ→0ℏ→0, we analyze a class of Schrödinger operators Hℏ=ℏ2L+ℏW+V⋅idℰ acting on sections of a vector bundle ℰ over a Riemannian manifold M where L is a Laplace type operator, W is an endomorphism field and the potential energy V has a non-degenerate minimum at some point p∈M. We construct quasimodes of WKB-type near p for eigenfunctions associated with the low-lying eigenvalues of Hℏ. These are obtained from eigenfunctions of the associated harmonic oscillator Hp,ℏ at p, acting on smooth functions on the tangent space.