 Unit-speed cubic time-derivative bispectrum

For $H_{\mathrm{int}}=-\int d^{3}x{a}^3ϵ\mathcal{F}{\dot{\zeta }}^3/H$, $a=-1/(H\tau )$ and stipulated modes $u_k=H(1+ik\tau )e^{-ik\tau }/\sqrt{4ϵk^3}$, six connected Wick contractions and ${\int }_{-\infty (1-i0)}^0{\tau }^2{e}^{iK\tau }d\tau =2i/K^3$ yield $B=3H^4\mathcal{F}/(8{ϵ}^2{k}_1{k}_2{k}_3{K}^3)$. Here $⟨{\zeta }^3{⟩}_c=(2\pi {)}^3\delta (\sum {\mathbf{k}}_i)B$. Using unstarred vertex modes with the external-before-vertex prefactor incorrectly reverses the sign. Non-unit sound speed also changes the quadratic modes, so this result is conditional on the specified unit-speed modes.