We prove compactness of the embeddings in Sobolev spaces for fractional super and sub harmonic functions with radial symmetry. The main tool is a pointwise decay for radially symmetric functions belonging to a function space defined by finite homogeneous Sobolev norm together with finite L2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$L^2$$\\end{document} norm of the Riesz potentials. As a byproduct we prove also existence of maximizers for the interpolation inequalities in Sobolev spaces for radially symmetric fractional super and sub harmonic functions.