The eigenvalue problem for the classical elliptic operator, for example, the Laplace operator $-\Delta$, was well-studied. It is well-known that the eigenvalue of the Dirichlet Laplacian is monotone with respect to the domain. In this thesis, we will generalize the ideas to the spectral fractional elliptic operators. We consider the weighted eigenvalue problem involving the spectral fractional elliptic operator $L^{\gamma}$ for $\gamma \in (0,1)$ on a bounded Lipschitz domain $\Omega \subset \mathbb{R}^{n}$. Here we remark that the operator $L^{\gamma}$ is well-defined for all $\gamma \in \mathbb{R}\setminus\mathbb{N}$. Similar to the classical case, we want to show that the eigenvalue of $L^{\gamma}$ is monotone with respect to the domain. Moreover, the principal eigenvalue is continuous when it is viewed as a function of domains. For our purpose, we will establish a corresponding Rayleigh quotient. Last but not least, we remark that the domain monotonicity of the eigenvalues does not hold for Neumann Laplacian. We will elaborate a counter-example in [13] and [14].