In this paper we present the asymptotic analysis of the realised quadratic variation for multivariate symmetric $β$-stable Lévy processes, $β\in (0,2)$, and certain pure jump semimartingales. The main focus is on derivation of functional limit theorems for the realised quadratic variation and its spectrum. We will show that the limiting process is a matrix-valued $β$-stable Lévy process when the original process is symmetric $β$-stable, while the limit is conditionally $β$-stable in case of integrals with respect to symmetric $β$-stable motions. These asymptotic results are mostly related to the work [5], which investigates the univariate version of the problem. Furthermore, we will show the implications for estimation of eigenvalues and eigenvectors of the quadratic variation matrix, which is a useful result for the principle component analysis. Finally, we propose a consistent subsampling procedure in the Lévy setting to obtain confidence regions.