Let
- $\Omega\subseteq\mathbb{R}^n$ be a bounded domain
- $\lambda_1$ be the first weak eigenvalue of $-\Delta$ in $\Omega$
- $\varphi_1$ be the weak eigenfunction associated with $\lambda_1$
- $H:=W_0^{1,2}(\Omega)$ be the Sobolev space
- $|\;\cdot\;|_p$ be the seminorm $$|u|_p^p:=\int_\Omega|\nabla u|^p\;d\lambda^n\;\;\;\text{for }u:\Omega\to\mathbb{R}\;\text{weakly differentiable}$$
- $\left\|\;\cdot\;\right\|_p$ be the $L^p(\Omega)$-norm and $\langle\;\cdot\;,\;\cdot\;\rangle_2$ be the scalar product in $L^2(\Omega)$
From the basic theory of the eigenvalue problem of the Laplacian, one knows that $$\lambda_1=\inf_{H\setminus\left\{0\right\}}R\;\;\;\;\;\text{with }R(u):=\frac{|u|_2^2}{\left\|u\right\|_2^2}\;\;\;\text{for }u\in H\setminus\left\{0\right\}$$ and $\varphi_1$ can be constructed such that
- $\left\|\varphi_1\right\|_2=1$
- $R(\varphi_1)=\lambda_1$
Since different eigenvalues correspond to $L^2$-orthogonal eigenfunctions (with $L^2$-orthogonal gradients), it's plausible to seek for the next eigenfunction $\varphi_2$ in $$U_1:=\left\{u\in H:\langle u,\varphi_1\rangle_2=0 \right\}\;.$$ We may observe, that $U_1$ is the null space of the continuous operator $\left.\langle\;\cdot\;,\varphi\rangle_2\right|_H$. Thus, $H$ is a closed subspace of $H$ and thereby a Hilbert space itself.
So, we can apply the same process used to find $(\lambda_1,\varphi_2)$ in order to find $$\lambda_2:=\inf_{U_1\setminus\left\{0\right\}}R$$ and $\varphi_2$.
Now, using the same arguments as discussed in a similar question, it's easy to see, that $\lambda_2\ge \lambda_1$ (Please don't get confused with the naming of the $\lambda_i$ in the other question).
However, is it possible that we have $\lambda_1=\lambda_2$?