Eigenvalue counting function for Robin Laplacians on conical domains

We study the discrete spectrum of the Robin Laplacian $Q^{Omega}_alpha$ in $L^2(Omega)$, [ umapsto -Delta u, quad dfrac{partial u}{partial n}=alpha u text{ on }partialOmega, ] where $Omegasubset mathbb{R}^{3}$ is a conical domain with a regular cross-section $Thetasubset mathbb{S}^2$, $n$ is the outer unit normal, and $alpha>0$ is a fixed constant. It is known from previous papers that the bottom of the essential spectrum of $Q^{Omega}_alpha$ is $-alpha^2$ and that the finiteness of the discrete spectrum depends on the geometry of the cross-section. We show that the accumulation of the discrete spectrum of $Q^Omega_alpha$ is determined by the discrete spectrum of an effective Hamiltonian defined on the boundary and far from the origin. By studying this model operator, we prove that the number of eigenvalues of $Q^{Omega}_alpha$ in $(-infty,-alpha^2-lambda)$, with $lambda>0$, behaves for $lambdato0$ as [ dfrac{alpha^2}{8pi lambda} int_{partialTheta} kappa_+(s)^2d s +oleft(frac{1}{lambda}right), ] where $kappa_+$ is the positive part of the geodesic curvature of the cross-section boundary.