Let the input probe state be the pure state $\rho_s = |0\rangle\langle0|$ (unitary covariance of the depolarizing channel $\mathcal{E}_f(\rho_s) = (1-f)\rho_s + \frac{f}{2}I$ allows us this choice without loss of generality, see reference [15] in the paper you linked for a similar argument). The output state of the depolarizing channel is
$\rho_f = \mathcal{E}_f(|0\rangle\langle0|) = (1-f)|0\rangle\langle0| + \frac{f}{2}I.$
The derivative with respect to $f$ is
$\partial_f \rho_f = \begin{pmatrix} -1/2 & 0 \\ 0 & 1/2 \end{pmatrix}$.
Solving for the elements of $L_f$ from $2\partial_f \rho_f=L_f\rho_f+\rho_f L _f$ gives $(L_f)_{00} = \frac{-1}{2-f}$ and $(L_f)_{11} = \frac{1}{f}$.
Thus $$L_f = \begin{pmatrix} -1/(2-f) & 0 \\ 0 & 1/f \end{pmatrix}$$
and the QFI is
\begin{align*} J(f) = \mathrm{Tr}(\rho_f L_f^2) &= \mathrm{Tr}\left(\begin{pmatrix} 1-f/2 & 0 \\ 0 & f/2 \end{pmatrix} \begin{pmatrix} 1/(2-f)^2 & 0 \\ 0 & 1/f^2 \end{pmatrix}\right) \\[1em] &= \frac{1}{f(2-f)}. \end{align*}
Let the input probe be a maximally entangled state $\rho_{AS} = |\Phi^+\rangle\langle\Phi^+|$, where $|\Phi^+\rangle = \frac{1}{\sqrt{2}}(|00\rangle + |11\rangle)$. The output state $\rho_f = (I_A \otimes \mathcal{E}_f)(\rho_{AS})$ is diagonal in the Bell basis $\{|\Phi^+\rangle, |\Phi^-\rangle, |\Psi^+\rangle, |\Psi^-\rangle\}$ with eigenvalues:
\begin{align*} \lambda_{\Phi^+}(f) = 1-3f/4, \quad \lambda_{\Phi^-}(f) = f/4, \quad \lambda_{\Psi^+}(f) = f/4, \quad \lambda_{\Psi^-}(f) = f/4. \end{align*}
The operator $2\partial_f \rho_f$ has diagonal elements $
-3/2, 1/2,1/2, 1/2.$
Simplifying $2\partial_f \rho_f=L_f\rho_f+\rho_f L _f$ we obtain diagonal elements of $L_f$:
\begin{align*} (L_f)_{\Phi^+\Phi^+} &= \frac{-3}{4-3f}, \\ (L_f)_{\Phi^-\Phi^-} &= (L_f)_{\Psi^+\Psi^+} = (L_f)_{\Psi^-\Psi^-} = \frac{1}{f}. \end{align*}
$L_f^2$ has diagonal entries $\frac{9}{(4-3f)^2}$ and $\frac{1}{f^2}$ (with multiplicity 3).
The QFI is
\begin{align*} J(f) &= \mathrm{Tr}(\rho_f L_f^2) = (1-3f/4) \frac{9}{(4-3f)^2} + 3 \cdot (f/4) \frac{1}{f^2} = \frac{1}{f(4/3-f)}. \end{align*}