I want to find the root space decomposition of the symplectic lie algebra $\mathfrak{sp}(2n,F)=C_n$.
I use the notation from Humphreys. The root space decomposition of a semisimple lie algebra $L$ is $L=H\oplus \bigoplus_{\alpha \in \Phi} L_\alpha$. Where $H$ is a maximal total subalgebra (this is more typically called cartan subalgebra). The $L_\alpha$ are the root spaces, and $\Phi$ is the root system.
First we must determine a suitable $H$. For this it seems we can pick the diagonal matrices in $C_n$. I think that this is toral since all its elements are diagonal and hence semisimple? To see that it is maximal, suppose not i.e. $H\subset H'$ where $H'$ is the maximal toral subalgebra. There there must be some $a\in H'$ that commutes with every $ha=ah$ for every $h\in H$. But I think by picking some $h$'s cleverly this implies that $a$ must also be diagonal.
My main confusion is about trying to find the roots, and then the root spaces. The roots are the $\alpha$ such that $L_\alpha$ is non zero. How are we supposed to find which $L_\alpha$ are non zero before finding the $\alpha$'s?
If we try to work directly from the definition we have $L_\alpha=\{x\in L \,|\, [h,x]=\alpha(h)x \quad \forall h \in H \}$, we are left with quite a complicated eigen value equation to solve. I think if we had an intuition for for the spaces should look like, we could use the fact that root spaces are one dimensional.
I have done this calculation for $\mathfrak{sl}(n,F)$ but that feel too prototypical to help get a feel for doing these. I would like to complete this calculation for $\mathfrak{sp}(2n,F)$ and then try again myself to do the other classical lie algebras.
You're right, $H$, the subalgebra consisting of the diagonal matrices in ${sp}(2n, \mathbb C)$ do form a Cartan subalgebra. Clearly, $H$ is an abelian subalgebra consisting of diagonalisable elements.
But how can we see there is no bigger abelian subalgebra consisting of diagonalisable elements? I suggest you don't worry about showing that this straight away. Instead, I suggest you go ahead and show that the $L = sp(2n, \mathbb C)$ has a decomposition $$ L = H \oplus \bigoplus_{\alpha} L_\alpha,$$ where each $\alpha \in H^\star$ is non-zero, and where $[h , x] = \alpha(h) x$ for $h \in H $ and $x \in L_\alpha$. Once you know that the $\alpha$'s are non-zero, then you know that no element outside of $H$ commutes with all the elements in $H$, which means that you can't make $H$ bigger and have it remain abelian.
So how do we find generators for these $L_\alpha$'s? I agree that it's a hard eigenvalue equation to solve. But actually, it's not that hard to guess a set of generators for the $L_\alpha$'s. In fact, the most natural vector-space basis for $L$ that you can think of gives you a set of generators for these $L_\alpha$'s!
I'll do the case of $sp(4, \mathbb C)$. A convenient basis for $H$ is $$ H_1 := \begin{bmatrix} 1 & 0 & 0 &0 \\ 0 & 0 & 0 & 0 \\ 0 &0 &-1 & 0 \\ 0 &0 &0 &0\end{bmatrix}, \ \ H_2 := \begin{bmatrix} 0 & 0 & 0 &0 \\ 0 & 1 & 0 & 0 \\ 0 &0 & 0 & 0 \\ 0 &0 &0 &-1\end{bmatrix}$$
And then, the $L_\alpha$'s are generated by $$ X_{1, 2} = \begin{bmatrix} 0 & 1 & 0 &0 \\ 0 & 0 & 0 & 0 \\ 0 &0 &0 & 0 \\ 0 &0 &-1 &0\end{bmatrix} \ \ \ \ ({\rm with \ } \alpha(H_1) = 1, \ \alpha(H_2) = -1)$$ $$ X_{2, 1} = \begin{bmatrix} 0 & 0 & 0 &0 \\ 1 & 0 & 0 & 0 \\ 0 &0 &0 & -1 \\ 0 &0 &0 &0\end{bmatrix} \ \ \ \ ({\rm with \ } \alpha(H_1) = -1, \ \alpha(H_2) = 1)$$ $$ Y_{1, 2} = \begin{bmatrix} 0 & 0 & 0 &1 \\ 0 & 0 & 1 & 0 \\ 0 &0 &0 & 0 \\ 0 &0 &0 &0\end{bmatrix} \ \ \ \ ({\rm with \ } \alpha(H_1) = 1, \ \alpha(H_2) = 1)$$ $$ Z_{1, 2} = \begin{bmatrix} 0 & 0 & 0 &0 \\ 0 & 0 & 0 & 0 \\ 0 &1 &0 & 0 \\ 1 &0 &0 &0\end{bmatrix} \ \ \ \ ({\rm with \ } \alpha(H_1) = -1, \ \alpha(H_2) = -1)$$ $$ U_{1} = \begin{bmatrix} 0 & 0 &1 & 0\\ 0 & 0 & 0 & 0 \\ 0 &0 &0 & 0 \\ 0 &0 &0 &0\end{bmatrix} \ \ \ \ ({\rm with \ } \alpha(H_1) = 2, \ \alpha(H_2) = 0)$$ $$ U_{2} = \begin{bmatrix} 0 & 0 & 0 &0 \\ 0 & 0 & 0 & 1 \\ 0 &0 &0 & 0 \\ 0 &0 &0 &0\end{bmatrix} \ \ \ \ ({\rm with \ } \alpha(H_1) = 0, \ \alpha(H_2) = 2)$$ $$ V_{1} = \begin{bmatrix} 0 & 0 & 0 &0 \\ 0 & 0 & 0 & 0 \\ 1 & 0 & 0 & 0 \\ 0 &0 &0 &0\end{bmatrix} \ \ \ \ ({\rm with \ } \alpha(H_1) = -2, \ \alpha(H_2) = 0)$$ $$ V_{2} = \begin{bmatrix} 0 & 0 & 0 &0 \\ 0 & 0 & 0 & 0 \\ 0 &0 &0 & 0 \\ 0 &1 &0 & 0\end{bmatrix} \ \ \ \ ({\rm with \ } \alpha(H_1) = 0, \ \alpha(H_2) = -2)$$
It shouldn't be too hard to generalise this to higher $n$!