¿$ f $ también es una autofunción de $ L_- $?
En general, no. Introducimos la notación $$ L^2 |\ell;m\rangle=\ell(\ell+1)|\ell;m\rangle \tag{1} $$ y $$ L_z|\ell;m\rangle=m |\ell;m\rangle \tag{2} $$
En general$^1$, $$ L_\pm|\ell;m\rangle=\sqrt{\ell(\ell+1)-m(m\pm1)}\;|\ell;m\pm 1\rangle \tag{3} $$ lo que significa que $|\ell;m\rangle$ no es un autovector de $L_\pm$, porque al actuar con $L_\pm$ en $|\ell;m\rangle$, no se obtiene un múltiplo escalar de $|\ell;m\rangle$, sino un vector diferente.
En el caso específico $|\ell;\ell\rangle$ (es decir, cuando $m=\ell$), tenemos $$ L_+|\ell;\ell\rangle=0 \tag{4} $$ lo que significa que $|\ell;\ell\rangle$ es un autovector de $L_+$, con valor propio 0. Podemos decir lo mismo acerca de $L_-$: $|\ell;-\ell\rangle$ es un autovector de $L_-$, con valor propio $0$.
En general, ¿cuáles son las autofunciones de $L_\pm$?
Para empezar, $L_\pm$ no son hermitianos, por lo que no hay garantía de que sean diagonalizables, y si lo son, los autovalores no serán reales (es decir, $L_\pm$ no es un observable). En el párrafo anterior, argumentamos que $|\ell;\ell\rangle$ es un autovector de $L_+$, por lo que al menos existe un autovector. Ahora vamos a demostrar que $|\ell;\ell\rangle$ es el único autovector de $L_+$.
Supongamos que existen un conjunto de vectores $|\alpha;\ell\rangle$ tales que $$ L_+|\alpha;\ell\rangle=\alpha|\alpha;\ell\rangle \tag{5} $$ y $$ L^2|\alpha;\ell\rangle=\ell(\ell+1)|\alpha;\ell\rangle \tag{6} $$ (nótese que podemos hacer esto porque $[L_+,L^2]=0$).
Como el conjunto $\{|\ell;m\rangle\}$ es una base, podemos escribir $|\alpha;\ell\rangle$ como una combinación lineal de estos vectores: $$ |\alpha;\ell\rangle=\sum_{m=-\ell}^{+\ell} c^m_\ell|\ell;m\rangle \tag{7} $$
Es bastante fácil verificar que $(5)$ solo se puede cumplir si todos los coeficientes $c^m_\ell=0$ excepto por $c^\ell_\ell$, de forma que $$ |\alpha;\ell\rangle=|\ell;\ell\rangle \tag{8} $$ se deduce fácilmente.
En el caso del oscilador armónico simple, donde el álgebra es similar, la situación es diferente: allí, tenemos una expresión similar a $(7)$, pero donde la suma es sobre $n=0,1,\cdots,\infty$. En este caso, la conclusión es diferente, porque hay un número infinito de términos en la suma. Ahora, es fácil demostrar que existen un conjunto de coeficientes no nulos $c_n$, lo que significa que existen autovectores no nulos de los operadores de creación y destrucción. Estos se llaman estados coherentes, y son bastante interesantes.
$^1$ La demostración de esta expresión es bastante estándar y se puede encontrar en línea y en cualquier libro de Mecánica Cuántica. Voy a reproducer la demostración aquí para hacer la publicación más autocontenido.
El álgebra de los operadores del momento angular es $$ [L_x,L_y]=iL_z\qquad [L_y,L_z]=iL_x \qquad [L_z,L_x]=iL_y \tag{9} $$
Si definimos $L_\pm=L_x\pm i L_y$, entonces es fácil ver que $(9)$ es equivalente a $$ [L_z,L_\pm]=\pm L_\pm \tag{10} $$
Por ejemplo, $[L_z,L_+]=[L_z,L_x+iL_y]=[L_z,L_x]+i[L_z,Ly]$, que en virtud de $(9)$ es igual a $=iL_y+L_x=L_+$.
Con esto, podemos demostrar $(3)$. Sea $$ |\varphi\rangle\equiv L_+|\ell;m\rangle \tag{11} $$ por definición. Si actuamos a la izquierda con $L_z$, obtenemos $$ L_z|\varphi\rangle=L_z L_+|\ell;m\rangle \tag{12} $$
Luego, escribimos $L_zL_+=L_+L_z+[L_z,L_+]$ (esto debería ser obviamente cierto: simplemente expande el conmutador y verifica que funcione). Dado que sabemos que $[L_z,L_+]=L_+$, obtenemos $$ (12)=(L_+ L_z+L_+)|\ell;m\rangle \tag{13} $$ que, usando $L_z|\ell;m\rangle=m|\ell;m\rangle$, es igual a $$ (12)=(1+m)L_+|\ell;m\rangle \tag{14} $$
Finalmente, notemos que $L_+|\ell;m\rangle$ es, por definición, $|\varphi\rangle$, lo que significa que $$ L_z|\varphi\rangle=(m+1)|\varphi\rangle \tag{15} $$
¡Esta relación es muy importante! Intenta pensar en ella por un minuto. Mírala detenidamente. Esta relación significa que $|\varphi\rangle$ es un autovector de $L_z$, y su valor propio es $m+1$. Por lo tanto, debemos tener $|\varphi\rangle\propto |\ell;m+1\rangle$, ya que $|\ell;m+1\rangle$ está definido como el autovector de $L_z$ con valor propio $m+1$.
Por lo tanto, podemos escribir $$ L_+|\ell;m\rangle=c|\ell;m+1\rangle \tag{16} $$ donde $c$ es una constante de normalización, que es fácil de encontrar, porque sabemos que $L_- L_+=L^2-L_z^2-L_z$: $$ |c|^2=\langle\ell;m|L_-L_+|\ell;m\rangle=\langle\ell;m|L^2-L_z^2-L_z|\ell;m\rangle=\ell(\ell+1)-m^2-m \tag{17} $$ donde utilicé $L^2|\ell;m\rangle=\ell(\ell+1)|\ell;m\rangle$ y $L_z|\ell;m\rangle=m|\ell;m\rangle. Esto completa la prueba de $(3)$.