Given a star-shaped bounded Lipschitz domain Ω ⊂ R d , we consider the Schrödinger operator L G = −∆ + V on Ω and its restrictions L Ωt G on the subdomains Ωt, t ∈ [0, 1], obtained by shrinking Ω towards its center. We impose either the Dirichlet or quite general Robin-type boundary conditions determined by a subspace G of the boundary space H 1/2 (∂Ω) × H −1/2 (∂Ω), and assume that the potential is smooth and takes values in the set of symmetric (N × N ) matrices. Two main results are proved: First, for any t 0 ∈ (0, 1] we give an asymptotic formula for the eigenvalues λ(t) of the operator L Ωt G as t → t 0 up to quadratic terms, that is, we explicitly compute the first and second t-derivatives of the eigenvalues. This includes the case of the eigenvalues with arbitrary multiplicities. Second, we compute the first derivative of the eigenvalues via the (Maslov) crossing form utilized in symplectic topology to define the Arnold-Maslov-Keller index of a path in the set of Lagrangian subspaces of the boundary space. The path is obtained by taking the Dirichlet and Neumann traces of the weak solutions of the eigenvalue problems for L Ωt G .