The question below relates to the paper "Jensen Polynomials for the Riemann Zeta Function and Other Sequences" of Griffin, Ono, Rolen and Zagier. I'm asking it here because I am sure the answer is well-known to a specialist, and was unsuccessful getting an answer at MSE. All definitions included below can be found in the above paper.
Recall that George Pólya proved that the Riemann Hypothesis (RH) is equivalent to the hyperbolicity of the Jensen polynomials associated to the sequence of Taylor coefficients $\{\gamma(n)\}$ defined by \begin{align} (-1+4z^{2})\Lambda(\frac{1}{2}+z)=\sum_{n=0}^{\infty}\frac{\gamma(n)}{n!}z^{2n}, \nonumber \end{align} where $\Lambda(s):=\pi^{-s/2}\Gamma(\frac{s}{2})\zeta(s)$.
There is a table in the above paper containing some values of $\gamma(n)$ for large $n$, and all of these values are positive (albeit only slightly so).
Question: Is it known whether $\gamma(n)>0$ for all $n\in\mathbb{N}$?
In case it is helpful, there are formulas for these beasts, when $n$ is positive:
$$\displaystyle\gamma(n)=\frac{n!}{(2n)!2^{2n-1}}\left(32{2n \choose 2}F(2n-2)-F(2n)\right)$$
where $$\displaystyle F(n)=\int_{1}^{\infty}\ln(t)^{n}t^{-3/4}\theta_{0}(t)\mathrm{d}t$$
and the theta function $\displaystyle \theta_{0}(t)=\sum_{k=1}^{\infty}e^{-\pi k^{2}t}$.
(I expect the answer is something like "Of course they aren't all positive, by the result of the paper you link to...it's because they are modeled eventually by Hermite Polynomials which themselves don't have this property", but I can't see this clearly, so am hoping for expert clarification. Thanks, in advance!)