In Sean Carroll's spacetime and geometry appendix B he derives the action of the Lie derivative on 1-forms. He finds that $\mathcal{L}_X Y^\mu = [X, Y]^\mu$, which I believe is meant as $\mathcal{L}_X(Y)^\mu$ since other books and wikipedia quote $\mathcal{L}_X(Y) = [X, Y]$ for some vector fields $X$ and $Y$. He then acts on a contraction $Y^\mu \omega_\mu$. Since this is equivalent to simply applying the vector field to the scalar he finds
$$
\mathcal{L}_X(Y^\mu \omega_\mu) = X^\lambda (\partial_\lambda \omega_\mu) Y^\mu + X^\lambda \omega_\mu (\partial_\lambda Y^\mu)
$$
He then compares this to the expression one finds by instead applying the Leibniz rule. Here he then goes on to write
\begin{align*}
\mathcal{L}_X(Y^\mu \omega_\mu) = \mathcal{L}_X(\omega)_\mu Y^\mu + \omega_\mu \mathcal{L}_X(Y)^\mu
\end{align*}
But I don't understand why the indices appear outside all of the sudden. I am unable to justify this step, but without further knowledge I am unable to derive the Lie derivative of a 1-form.
Here's my attempt: starting only with linearity, Leibniz rule, aswell as $\mathcal{L}_X(f) = X(f)$ and $\mathcal{L}_X(Y) = [X, Y]$. By expanding the commutator defintion I find \begin{align*} \mathcal{L}_X(Y) &= \mathcal{L}_X(Y^\mu) \partial_\mu + Y^\mu \mathcal{L}_X(\partial_\mu) \\ \left(X^\lambda \partial_\lambda Y^\mu - Y^\lambda \partial_\lambda X^\mu \right)\partial_\mu &= \mathcal{L}_X(Y^\mu)\partial_\mu + Y^\mu \left(-(\partial_\mu x^\lambda) \partial_\lambda\right) \end{align*} Where on the right hand side I treat $\partial_\mu$ as a vector field and then compute the commutator according to the defintion of $\mathcal{L}_X(\partial_\mu)$ as is done here. Comparing both sides I then find $\mathcal{L}_X(Y^\mu) = X^\lambda\partial_\lambda Y^\mu$, which makes sense to me since the components of a vector are just functions again. But this is at odds with the derivation of Carroll since clearly $\mathcal{L}_X(Y^\mu) \neq \mathcal{L}_X(Y)^\mu$ in this case. I don't know how to continue then however, I can show again that $\mathcal{L}_X(\omega_\mu) = X^\lambda \partial_\lambda \omega_\mu$, but it seems to me that I need further information such as the action on a basis 1-form which is quoted on wikipedia as $\mathcal{L}_X(\mathrm{d}x^\mu) = (\partial_\lambda X^\mu)\mathrm{d}x^\lambda$. But I do not know enough about differential geometry to understand the derivation of that. This question originates from an old exam of mine so I believe that it should be possible without any further information.