We consider harmonic Toeplitz operators
where
is the orthogonal projection onto
,
,
, is a bounded domain with boundary
, and
is an appropriate multiplier. First, we complement the known criteria which guarantee that
is in the pth Schatten–von Neumann class
, by simple sufficient conditions which imply
, the weak counterpart of
. Next, we consider symbols
which have a regular power-like decay of rate
at
, and we show that
is unitarily equivalent to a classical pseudo-differential operator of order
, self-adjoint in
. Utilizing this unitary equivalence, we obtain the main asymptotic term of the eigenvalue counting function for
, and establish a sharp remainder estimate. Further, we assume that Ω is the unit ball in
, and
is compactly supported in Ω, and investigate the eigenvalue asymptotics of the Toeplitz operator
. Finally, we introduce the Krein Laplacian K, self-adjoint in
, perturb it by a multiplier
, and show that
. Assuming that
and
, we study the asymptotic distribution of the discrete spectrum of
near the origin, and find that the effective Hamiltonian which governs this distribution is the Toeplitz operator
.