It is well known that if $M, \Omega$ is a symplectic manifold then the Poisson bracket gives $C^\infty(M)$ the structure of a Lie algebra. The only way I have seen this proven is via a calculation in canonical coordinates, which I found rather unsatisfying. So I decided to try to prove it just by playing around with differential forms. I got quite far, but something isn't working out and I am hoping someone can help. Forgive me in advance for all the symbols.

Here is the setup. Given $f \in C^\infty(M)$, let $X_f$ denote the unique vector field which satisfies $\Omega(X_f, Y) = df(Y) = Y(f)$ for every vector field $Y$. We define the Poisson bracket of two functions $f$ and $g$ to be the smooth function $\{f, g \} = \Omega(X_f, X_g)$. I can show that the Poisson bracket is alternating and bilinear, but the Jacobi identity is giving me trouble. Here is what I have.

To start, let's try to get a handle on $\{ \{f, g \}, h\}$. Applying the definition, this is given by $d(\Omega(X_f, X_g))X_h$. So let's try to find an expression for $d(\Omega(X,Y))Z$ for arbitrary vector fields $X, Y, Z$.

Write $\Omega(X,Y) = i(Y)i(X)\Omega$ where $i(V)$ is the interior product by the vector field $V$. Applying Cartan's formula twice and using the fact that $\Omega$ is closed, we obtain the formula

$$d(\Omega(X,Y)) = (L_Y i(X) - i(Y) L_X) \Omega$$

where $L_V$ is the Lie derivative with respect to the vector field $V$. Using the identity $L_V i(W) - i(W) L_V = i([V,W])$, we get:

$$(L_Y i(X) - i(Y) L_X) = L_Y i(X) - L_X i(Y) + i([X,Y])$$

Now we plug in the vector field $Z$. We get $(L_Y i(X) \Omega)(Z) = Y(\Omega(X,Z)) - \Omega(X,[Y,Z])$ by the definition of the Lie derivative, and clearly $(i([X,Y])\Omega)(Z) = \Omega([X,Y],Z)$. Putting it all together:

$$d(\Omega(X,Y))Z = Y(\Omega(X,Z)) - X(\Omega(Y,Z)) + \Omega(Y, [X,Z]) - \Omega(X, [Y,Z]) + \Omega([X,Y], Z)$$

This simplifies dramatically in the case $X = X_f, Y = X_g, Z = X_h$. The difference of the first two terms simplifies to $[X_f, X_g](h)$, and we get:

$$ \begin{align} \{\{f, g\}, h\} &= [ X_f, X_g ](h) + [ X_f, X_h ](g) - [ X_g, X_h ](f) - [ X_f, X_g ](h)\\ &= [ X_f, X_h ](g) - [ X_g, X_h ](f) \end{align}$$

However, this final expression does not satisfy the Jacobi identity. It looks at first glance as though I just made a sign error somewhere; if the minus sign in the last expression were a plus sign, then the Jacobi identity would follow immediately. I have checked all of my signs as thoroughly as I can, and additionally I included all of my steps to demonstrate that if a different sign is inserted at any point in the argument then one obtains an equation in which the left hand side is alternating in two of its variables but the right hand side is not. Can anybody help?