A beam splitter transforms incoming mode operators $\hat{a}_{i}, \hat{b}_{i}$ to outgoing operators $\hat{a}_{o}, \hat{b}_{o}$ as $$ \begin{aligned} &\hat{a}_{o}=\sqrt{\eta} \hat{a}_{i}+i \sqrt{1-\eta} \hat{b}_{i} \\ &\hat{b}_{o}=i \sqrt{1-\eta} \hat{a}_{i}+\sqrt{\eta} \hat{b}_{i} \end{aligned} $$
$\text { where } \eta=\cos ^{2} \theta$ I want to find for the incoming state $\left|\alpha_{i}\right\rangle \otimes\left|\beta_{i}\right\rangle$, which is a product of two coherent states what will be the outgoing states?
First I have showed that: $ \hat{T}=\exp \left[-i \theta\left(\hat{a}^{\dagger} \hat{b}+\hat{a} \hat{b}^{\dagger}\right)\right]$
is equivalent to the transformation matrix of the beam splitter using the below expansion:
$$ \exp \left(\mathrm{i} \theta \sigma_{1}\right)=\cos \theta+\mathrm{i} \sigma_{1} \sin \theta $$
Then expanded the coherent states as follows:
$$ |\alpha\rangle=e^{-\frac{|\alpha|^{2}}{2}} \sum_{n=0}^{\infty} \frac{|\alpha|^{n}}{\sqrt{n !}}|n\rangle $$
and obtained: $\left(\hat{I} \cos \theta+i\left(\hat{a} \hat{b}^{+}+\hat{b} \hat{a}^{+}\right) \sin \theta\right)e^{-|\alpha|^{2} / 2} \cdot e^{-|\beta|^{2} / 2} \sum_{m,n} \frac{|\beta|^{m}|\alpha|^{n}}{\sqrt{ nm !}}(|m\rangle \otimes|n\rangle)$
and further obtained:
$e^{-|\alpha|^{2} / 2} \cdot e^{-|\beta|^{2} / 2} \sum_{m,n} \frac{|\beta|^{m}|\alpha|^{n}}{\sqrt{ nm !}}(|m, n\rangle \cos \theta+i(|m-1, n+1\rangle+|m+1, n-1\rangle) \sin \theta)$
But couldn't continue further at this point. Thanks in advance.