Proving Darboux's theorem, Hofer-Zehnder try to find, given $\omega$ a closed nondegenerate 2-form and $\omega_0$ the canonical symplectic form, a family of diffeomorphisms $\phi^t$ such that for all $t$, if $\omega_t=\omega_0+t(\omega-\omega_0)$, then $(\phi^t)^\ast\omega_t=\omega_0$, with $0\leq t\leq1$. To find them, the book tries to construct a time-dependent field $X_t$ such that $\phi^t$ is the flow of $X_t$. Differentiating the relation of $\phi^t$ and $\omega_t$ with $\omega_0$ given above, we get, according to the book:
$$0=\frac{d}{dt}(\varphi^t)^* \omega_t=(\varphi^t)^*\left\{L_{X_t}\omega_t+\frac{d}{dt}\omega_t\right\}.$$
The first equals sign is just the differentiation. It is the second equals that is obscure to me. $L_{X_t}\omega_t$ should be just the LHS, AFAIK. But then I think, $\omega_t$ also depends on $t$, so there must be a second term. I can guess the second term is $(\phi^t)^\ast\frac{d}{dt}\omega_t$, i.e. the derivative passed on the right of the pushforwards, whereas the first term should be differentiating "w.r.t. the flow". BUt still I can't see how I get that brace. How do I prove this formula really holds?
When the vector field and the $2$-form are time independent I hope it is clear that:
$$(\phi^t)^*\circ\mathcal{L}_X(\omega)=\frac{\mathrm{d}}{\mathrm{d}s}(\phi^s)^*\omega\big|_{s=t}.$$
A simple change of variables $t-s$ proves this.
Now, it is important to understand what it means to take the derivative of something which is not a function, as $(\phi^s)^*\omega$. The symbol $\frac{\mathrm{d}}{\mathrm{d}t}(\phi^t)^*\omega$ stands for the derivative of the following function: for each point $p\in M$ and pair of vector fields $Y,Z\in\mathfrak{X}(M;\mathbb{R})$ one has a smooth function:
$$\mathbb{R}\ni t\mapsto (\phi^t)^*\omega(Y,Z)(p)\in\mathbb{R},$$
and $\frac{\mathrm{d}}{\mathrm{d}t}(\phi^t)^*\omega\big|_{t=t_0}$ is the $2$-form defined by:
$$\frac{\mathrm{d}}{\mathrm{d}t}(\phi^t)^*\omega\big|_{t=t_0}(Y,Z)(p):=\frac{\mathrm{d}}{\mathrm{d}t}(\phi^t)^*\omega(Y,Z)(p)\big|_{t=t_0}.$$
Once this is clarified, the time dependent case follows.
Update:
Let me show here how the computation works for time dependent functions. If $f_t\in C^\infty(M;\mathbb{R})$ is a time dependent function (i.e. $f(t,p):=f_t(p)$ defines a smooth function on $\mathbb{R}\times M$), then $((\phi^t)^*f_t)(p)=f(t,\phi^t(p))$ and $\frac{\mathrm{d}}{\mathrm{d}t}(\phi^t)^*f_t$ is obtained by computing:
\begin{align*} \frac{\mathrm{d}}{\mathrm{d}t}((\phi^t)^*f_t)(p)\big|_{t=\tau}={}&\frac{\partial}{\partial x}f(x,y)\big|_{(x,y)=(\tau,\phi^\tau(p))}+\frac{\partial}{\partial y}f(x,y)\big|_{(x,y)=(\tau,\phi^\tau(p))}\frac{\mathrm{d}}{\mathrm{d}t}\phi^t(p)\big|_{t=\tau}={} \\ {}={}&\left(\frac{\mathrm{d}}{\mathrm{d}t}(f_t)\big|_{t=\tau}\right)\circ\phi^\tau(p)+X_\tau(f_\tau)\circ\phi^\tau(p). \end{align*}
Therefore:
$$\frac{\mathrm{d}}{\mathrm{d}t}(\phi^t)^*f_t=(\phi^t)^*\left\{\mathcal{L}_{X_t}(f_t)+\frac{\mathrm{d}}{\mathrm{d}t}(f_t)\right\}.$$