Bosonic squeezed-vacuum contractions ($⟨a_k^{†}{a}_q⟩={\sinh }^2{\theta }_k{\delta }_{kq},⟨a_k{a}_q⟩=-\sinh {\theta }_k\cosh {\theta }_k{\delta }_{k,-q}$)

Let $b_k=\cosh {\theta }_k{a}_k+\sinh {\theta }_k{a}_{-k}^{†}$ for real even ${\theta }_k$, and let every $b_k$ annihilate a normalized vacuum state. Inverting the Bogoliubov transformation gives $a_k=\cosh {\theta }_k{b}_k-\sinh {\theta }_k{b}_{-k}^{†}$. Only $⟨b_k{b}_q^{†}⟩={\delta }_{kq}$ survives, proving the displayed vacuum expectation values. Taking the Hermitian adjoint gives the same real negative pair amplitude for two creation operators. The normal occupation number is quadratic at small ${\theta }_k$, whereas the pair amplitude is linear. For one opposite-momentum pair the two-mode squeezed vacuum is $\mathrm{sech}\theta {\sum }_{n\ge 0}(-\tanh \theta {)}^n|n,n⟩$, fixing the sign convention directly.