  We study the behavior of eigenvalues for magnetic Aharonov-Bohm operators
with half-integer circulation and Dirichlet boundary conditions in a planar
domain. We analyse the leading term in the Taylor expansion of the eigenvalue
function as the pole moves in the interior of the domain, proving that it is a
harmonic homogeneous polynomial and detecting its exact coefficients.
