  What if Big Bang was hot from its very inception? This is possible in a
bimetric theory where the source of fluctuations is thermal, requiring the
model to live on a critical boundary in the space of parameters and can be
realized when an anti-DBI brane moves within an $EAdS_2 \times E_3$ geometry.
This setup renders the model unique, with sharp predictions for the scalar
spectral index and its running. We investigate the non-Gaussian signatures of
this thermal bimetric model, or "bi-thermal" for short. We adapt the standard
calculation of non-Gaussianities for $P(X,\phi)$ models to the thermal nature
of the model, emphasising how the bi-thermal peculiarities affect the
calculation and alter results. This leads to precise predictions for the shape
and amplitude of the three-point function of the bi-thermal model (at
tree-level): $f^{\rm local} _{\rm NL} = -3/2$ and $f^{\rm equil} _{\rm NL} = -2
+ 4 \sqrt{3}\pi/9 \simeq 0.4$. We also discover a new shape of flattened
non-gaussianity $\propto (k_1+k_2-k_3)^{-3/2} +$ permutations, which is
expected due to the excited thermal initial conditions. These results, along
with our earlier predictions for the scalar power spectrum, provide sharp
targets for the future generation of cosmological surveys.
