Concentration of eigenfunctions on singular Riemannian manifolds

We consider a compact Riemannian manifold with boundary and a metric that is singular at the boundary. The associated Laplace-Beltrami operator is of the form of a Grushin operator plus a singular potential. In a supercritical parameter regime, we identify the rate of concentration and profile of the high-frequency eigenfunctions that accumulate at the boundary. We give an application to acoustic modes on gas planets.