NOTE:
I've seen this line of reasoning and those calculations somewhere, but I don't remember where. Here is my reconstruction of them.
Type I:
Let's have in mind the sketch of fusion type I from the original paper, Fig. 2a). Let us call the left port A, the down port B (these are the input ports), the right port C and the upper port D (these are the port where the detectors are located). In the middle there is a PBS, and before port D there is a waveplate which rotate the incoming signal by polarization.
In type I, you only detect the incoming photons from port D, but you can distinguish the polarization, that it to say you can either get no photons at all, a single H or V photons, or 2H or 2V photons.
Now, how does the input change after the transformation of the PBS + rotation? Let us assume that you have a general incoming wave, which is already connected to other qubits (for example, you have two large cluster states and you want to connect them by fusing together one photon from each). So, we have from port A: $f_{1}a_H^{\dagger} + f_{2}a_V^{\dagger}$, where $f_{1}$ and $f_{2}$ are the other modes, $a$ stands for the port and $H$ or $V$ is the polarization. From port B we have the same: $f_{3}b_H^{\dagger} + f_{4}b_V^{\dagger}$.
The H mode from port A doesn't feel the PBS because the photon comes in parallel to it, so it continues unharmed to a horizontal mode in port C. On the contrary, the vertical mode will turn to port D and then will be rotated. The opposite happens for the photons from port B. Finally, what will happen is:
- $a_{H}^{\dagger} \to c_{H}^{\dagger}$
- $a_{V}^{\dagger} \to \dfrac{d_{H}^{\dagger} - d_{V}^{\dagger}}{\sqrt{2}}$
- $b_{V}^{\dagger} \to c_{V}^{\dagger}$
- $b_{H}^{\dagger} \to \dfrac{d_{H}^{\dagger} + d_{V}^{\dagger}}{\sqrt{2}}$
(The reason why port B doesn't get a minus sign, is because it comes with H polarization, so it go through unharmed, and than get a rotation).
Finally, one can write:
$(f_{1}a_{H}^{\dagger} + f_{2}a_{V}^{\dagger})(f_{3}b_{H}^{\dagger}f_{4}b_{V}^{\dagger}) \to ...$, replace $a$ and $b$ with $c$ and $d$ using relations 1-4, do some algebra, and you should get:
$\to \dfrac{1}{\sqrt{2}}f_{1}f_{3}c_{H}^{\dagger}(d_{H}^{\dagger} + d_{V}^{\dagger}) + f_{1}f_{4}c_{H}^{\dagger}c_{V}^{\dagger} + \dfrac{1}{2}f_{2}f_{3}((d_{H}^{\dagger})^2 - (d_{V}^{\dagger})^2) + \dfrac{1}{\sqrt{2}}f_{2}f_{4}c_{V}^{\dagger}(d_{H}^{\dagger} - d_{V}^{\dagger})$.
For simplicity, and to be aligned with the paper, let us assume that you started with two bell states, such that $f_{1}a_{H}^{\dagger} = |H, H\rangle$, $f_{2}a_{V}^{\dagger} = |V, V\rangle$, and the same for port B ($f_{3}$ is H and $f_{4}$ is V). Also, the final 3-qubit state (after post-selecting the output from the detector) is written in the order $|A, B, C\rangle$.
Now, when do you get a GHZ? If you obtained a single photon in detector D, no matter if it is in the H or V polarization. Because, in that case, you only left with the first and last terms, which are $\dfrac{1}{\sqrt{2}}(|H, H, H\rangle \pm |V, V, V\rangle)$, depending if you detected (H, 0) or (0, V), respectively.
You can also check what are the results from the other options, and see that this is not the case.
Type II:
If you understand this, type II is very similar, in the sense that you just need to write the general input signal (which is exactly the same) and do the transformation for the output ports. In that case, the sketch is somewhat different (see Fig. 2b) in the paper), because each port first has a rotation. So first of all, the waveplate transforms the input as follows:
- $a_{H}^{\dagger} \to \dfrac{a_{H}^{\dagger} + a_{V}^{\dagger}}{\sqrt{2}}$
- $a_{V}^{\dagger} \to \dfrac{a_{H}^{\dagger} - a_{V}^{\dagger}}{\sqrt{2}}$
- $b_{H}^{\dagger} \to \dfrac{b_{H}^{\dagger} + b_{V}^{\dagger}}{\sqrt{2}}$
- $b_{V}^{\dagger} \to \dfrac{b_{H}^{\dagger} - b_{V}^{\dagger}}{\sqrt{2}}$
Then there is a PBS and another rotation, so each individual polarization goes to:
- $a_{H}^{\dagger} \to \dfrac{c_{H}^{\dagger} + c_{V}^{\dagger}}{\sqrt{2}}$
- $a_{V}^{\dagger} \to \dfrac{d_{H}^{\dagger} - d_{V}^{\dagger}}{\sqrt{2}}$
- $b_{H}^{\dagger} \to \dfrac{d_{H}^{\dagger} + d_{V}^{\dagger}}{\sqrt{2}}$
- $b_{V}^{\dagger} \to \dfrac{c_{H}^{\dagger} - c_{V}^{\dagger}}{\sqrt{2}}$
Now what is left is to do the same algebra, to obtain that when one gets a single detection in each detector it acts like a parity measurement that projects the "f"'s into the even/odd parity of the input, depends if one gets the same or different polarization.
Notice that, in fact, in type II fusion the two qubits are gone. It seems wired, but because (when you have a "successful" fusion, in the aforementioned sense) you know exactly to which subspace you projected the state, you can do a computation, similar to a single-qubit measurement in measurement-based quantum computation.