In this paper where the authors derive the formula for coincidence probability in a Hong-Ou Mandel (HOM) interference effect as a function of time delay $\tau$, they arrive at an equation (15) with the following matrix element:
$$\langle\text{vac}|~\hat{c}(\omega)\hat{d}(\omega')~[\hat{c}^{\dagger}(\omega_1)\hat{d}^{\dagger}(\omega_2) - \hat{c}^{\dagger}(\omega_1)\hat{d}^{\dagger}(\omega_2)+\hat{d}^{\dagger}(\omega_1)\hat{c}^{\dagger}(\omega_2)-\hat{d}^{\dagger}(\omega_1)\hat{d}^{\dagger}(\omega_2)]~|\text{vac}\rangle\tag{1}$$
They then use a commutation relation (equation 2) to arrive at a set of delta functions. The commutation relation they exploit is:
$$[\hat{a}(\omega),\hat{a}^{\dagger}(\omega')]=\delta(\omega-\omega')\tag{2}$$
However, if you expand out equation (1), it is clear that all $\hat{c}$'s appear to the left of all $\hat{c}^{\dagger}$'s, so it is completely unclear to me where the commutation relation is coming in!
Can anyone help?