Let $M$ be a closed manifold, equipped with a volume form $\omega$. I understand that the vector space
$$ \{ X \in \Gamma(TM) \, | \, L_X\omega=0 \}$$
is always infinite-dimensional.
Is there an elementary argument showing this?
We can assume of course $\omega$ comes from a Riemannian metric $g$, and then $L_X\omega=0$ if and only if $\text{trace}(\nabla X)=0$.
Let us show first that the space of compactly supported, divergence-free vector fields on $\mathbb{R}^n$ is infinite dimensional. Let $\rho:\mathbb{R}^+\to\mathbb{R}^+$ be a smooth, compactly supported function, such that $\rho\equiv1$ on a neighborhood of $0$. Let $A=(a_i^j)\in M_{n\times n}(\mathbb{R})$ be anti-symmetric. In particular, all the diagonal entries of $A$ vanish. Define a vector field $X_{\rho,A}$ on $\mathbb{R}^n$ by $$X_{\rho,A}(x)=\rho\left(|x|\right)A(x).$$(Here we think of a vector field as a function $\mathbb{R}^n\to\mathbb{R}^n$). This is compactly supported, as $\rho$ is. We compute the divergence: $$\begin{align}\mathrm{div}X_{\rho,A}&=\frac{\partial}{\partial x^1}\left(\rho(|x|)a_1^ix_i\right)+\ldots+\frac{\partial}{\partial x^n}\left(\rho(|x|)a_n^ix_i\right)\\&=\dot{\rho}(|x|)\frac{x_1}{|x|}a_1^ix_i+\ldots+\dot{\rho}(|x|)\frac{x_n}{|x|}a_n^ix_i\\&=0.\end{align}$$ This clearly provides an infinite-dimensional family of vector fields, due to the freedom in the choice of $\rho$.
We return to the closed manifold $M$. Recall that a volume form is always integrable. This means that there is a coordinate neighborhood $U$, in which we have $\omega=dx^1\wedge\ldots\wedge dx^n.$ Now, we may choose some $X_{\rho,A}$ in $U$ and extend it by $0$ to a global vector field.