High energy asymptotics for eigenvalues of the Schrödinger operator with a matrix potential

We consider a Schrodinger operator with a matrix potential defined in L-2(m) (Q) by the differential expression Lu = -Delta u + Vu and the Neumann boundary condition, where Q is a d-dimensional parallelepiped and V a matrix potential, d >= 2, m >= 2. We obtain the high energy asymptotics of arbitrary order for a rich set of eigenvalues.