Boundary behaviour of eigenfunctions and superharmonic functions on harmonic manifolds of purely exponential volume growth

On $\mathbb{X}$, a non-positively curved harmonic manifold of purely exponential volume growth, of dimension $n \ge 3$, we study certain quantitative aspects of the boundary behaviour of eigenfunctions and superharmonic functions. We first focus on complex-valued eigenfunctions lying outside the $L^2$-spectrum of $\Delta$ and obtain the almost everywhere existence of weighted non-tangential limits, sharp Hausdorff dimension and Hausdorff measure estimates of the boundary exceptional sets for radial limits. Then in the second part, we shift our attention to non-tangential and tangential boundary behaviour of positive superharmonic functions. Most of our results are new even for the homogeneous setting of rank one Riemannian symmetric spaces of non-compact type and Damek-Ricci spaces. Our arguments are based on potential theory adapted to the intrinsic Gromov hyperbolic geometry of $\mathbb{X}$.