I was idly thinking about why one might naïvely expect a discrete spectrum of eigenvalues for a linear operator $L$ when I dreamt up the following argument (which I expect isn't new instead - references?). Can it be made rigorous, with any necessary added restrictions, or does it at least offer some genuine informal insight? (Or did I just mess something up?)
Suppose $L$ is a linear, self-adjoint operator acting on (sufficiently) continuously differentiable functions $\phi:[0,1]\to\mathbb{R}$ with boundary conditions $\phi(0)=\phi(1)=0$. An eigenfunction $\phi_\lambda(x)$ is a non-zero such function such that
$$\left[L \phi_\lambda\right](x) = \lambda \phi_\lambda(x)$$
We expect that such functions occur only for 'discrete' $\lambda$, in particular not for an interval $(a,b)$ of eigenvalues. Clarification: I'm thinking of $L$ being a 'nice' integral/differential operator in my head. What I'm interested in is what further restrictions are necessary to make the spectrum simple.
Here's a reason to suspect this is the case: suppose $\phi\equiv \phi(x;\lambda)$ is a (sufficiently) continuously differentiable function $[0,1]\times(a,b)\to\mathbb{R}$ such that $L \phi = \lambda\phi$. Suppose further that $L,\phi$ is sufficiently nicely behaved that $L (\partial_\lambda \phi) \equiv \partial_\lambda(L\phi)$. Note that we expect $\partial_\lambda \phi|_{x=0,1}=0$ since all $\phi$ obey the boundary conditions, so $L$ should be self-adjoint on this function too.
Then, with a simple $L^2$ inner product $\left<\cdot,\cdot\right>$ say,
$$\boxed{\lambda \left<\phi,\partial_\lambda \phi\right> = \left<L\phi,\partial_\lambda \phi\right> = \left<\phi,L\partial_\lambda \phi\right> = \left<\phi,\partial_\lambda L\phi\right> = \left<\phi,\partial_\lambda (\lambda\phi)\right> = \left<\phi,\phi\right> + \lambda \left<\phi,\partial_\lambda \phi\right>}$$
and hence $$\left<\phi,\phi\right> = 0$$
Hence there is no continuously differentiable family of solutions which vary the eigenvalue throughout some interval.
Example By contrast, if we dropped our boundary conditions and worked on the fully infinite domain $x\in(-\infty,\infty)$ the operator $Lf=f''(x)$ acquires all real eigenvalues on the smooth families $\cos kx,\exp kx$. But for s finite interval with homogeneous BCs, we only get e.g. $0,\pi,2\pi,\ldots$
Any interesting thoughts are welcome! I haven't thought about this in great detail, so it may be that this just leads to a standard spectral theory argument when one tries to make it rigorous, but I'm curious as to whether it has an interesting interpretation, and whether you think it would be a good, informal justification to teach someone first coming to Sturm-Liouville theory.
Mainly, I am interested in whether this can be tweaked to a proof without losing all the simple spirit of the boxed equation above.