This tutorial introduces the physics of many-boson and many-fermion interference required for the description of current experiments and for the understanding of novel approaches to quantum computing.
This tutorial presents an illustrative introduction to the physics of many-boson and many-fermion interference, which supplies the reader the necessary tools to describe photonic entanglement manipulation and novel approaches to optical quantum computing.
Any incoming particle is reflected with probability R and transmitted with probability . An event with and particles in the first and second output mode, respectively, occurs with probability . In particular, for a balanced setup, , the distribution of distinguishable particles in the output modes follows a binomial, and
The creation operator of a particle in input mode is denoted by , such that the initial state reads
where denotes the vacuum state. Throughout this tutorial,the letter denotes Fock-states, i.e. is interpreted as mode occupation , etc. The letters and are reserved for qudits controlled by observing parties; then denotes the state in which is observed by the first party, by the second, and so on.
The time-evolution of particle creation operators in the Heisenberg picture reads accordingly the single-particle wavefunction acquires a phase-shift of upon reflection. The single-particle time-evolution can be inserted into the initial state, and, for a balanced beam splitter , the final state becomes Fermions for fermionic creation operators. The final state becomes from which we immediately read off
The Pauli principle prohibits the double population of any output mode, which is compatible with the constructive collective interference that enhances the probability to find the particles in different modes.
Bosons Bosonic bunching and the Pauli principle are not the consequence of coherently superimposed many-particle paths, but they are rather kinematic boundary conditions on the state-space of identical particles.
Collective many-particle interference can only affect many-particle observables
As a realistic example for such transition, we consider photons and choose the time of arrival at the beam splitter as the degree of freedom through which the indistinguishability of the light quanta may be affected. In quantum optics experiments , the relative delay in the arrival time is manipulated by the spatial displacement of the wave-packets to the beam splitter.
The resulting event probability for partially distinguish- able particles becomes the average of the probabilities assigned to bosons and to distinguishable particles, weighted by and , respectively: To model typical experiments with photons , a Gaussian frequency distribution around a central frequency with width is assumed, The temporal overlap of two photons then becomes
The visibility of the HOM-dip constitutes today the defacto standard for ensuring the indistinguishability of two bosonic particles in the experiment.
By post-selecting events with a well-defined total number of photons N, i.e. by neglecting all final events in which the total number of detected photons does not match the desired N, one can effectively project with high fidelity onto the component with a total of N light quanta in the initial state.
the outcomes of the measurements are determined by the distribution of the objects among the observers, i.e. before the detection actually takes place.
By post-selecting those events with exactly one particle per spatial output mode, i.e. coincident events, we can effectively project the state onto where for bosons and takes into account the anti-commutation relation of fermionic creation operators.
If the particles prepared in the input modes do not differ only in the degree of freedom in which entanglement is to be created, but also in another physical property, the correlations measured at the output modes will eventually degrade to purely classical ones. In practice, this deteriorating degree of freedom is often the time of arrival.
The indistinguishability directly quantifies the many-body coherences of the emerging state and, thus, the quantum nature of the induced correlations. Indeed, a quantitative indicator for entanglement, the concurrence directly reflects the indistinguishability
Relying on propagation and detection for entanglement manipulation also leads to natural constraints on operations. For example, a measurement that fully discriminates all four Bell states is impossible, which inhibits deterministic teleportation, for which nonlinear interactions are necessary.
For indistinguishable particles, the argument can be repeated. Now, however, the exchange of two particles in the output modes does not change the quantum state (up to a global sign change for fermions), such that all possibilities to distribute the particles among the output modes contribute coherently to the final state. As a consequence, many-body amplitudes of the form need to be summed
where takes into account the phase acquired upon exchange of two fermions in the output modes.
For Pauli arrangements, the difference between the average probability for bosons and fermions dies out when the density of particles N/n is decreased: in the dilute limit N/n → 0, bunched events with several particles per mode are a priori very unlikely, and their influence on the statistics vanishes
For fully bunched events for which one mode receives all particles, , we can formulate an exact relation between the probabilities for distinguishable particles and bosons On the other hand, if all bosons are prepared in the same mode, the probabilities for bosons and distinguishable particles do not differ at all.
In order to find a more compact expression for the probabilities, a characteristic -matrix can be defined as i.e. M contains those rows and columns of U that correspond to (possibly multiply) occupied in- and output modes, respectively. With M defined as such, we can write
where the function is the permanent of the matrix , for which the absolute-square is taken component-wise. For bosons, we find i.e. the permanent of the complex matrix needs to be computed before taking the absolute-square to obtain the probability. For fermions, the sign-change allows us to identify the total amplitude as a determinant, The evaluation of the determinant can be performed in time polynomial in the matrix size .
For distinguishable particles, by sampling over many randomly chosen permutations , the permanent can be estimated.
It is tempting to state that quantum mechanics makes it possible to 'compute the permanent of a complex matrix with exponential speedup with respect to classical computers' using bosons that propagate through a multimode setup. However, even though the amplitudes of final states characterized by correspond to the permanent of a certain matrix , these amplitudes are very small, which makes them difficult to extract. Although it is possible to define an observable whose expectation value coincides with the permanent, its large variance requires an exponential number of measurements to extract the permanent experimentally . This caveat jeopardizes any strategy to directly compute the permanent via the scattering of bosons, although there are current attempts in this direction.
An alternative approach to certifying boson-sampling devices is to leave the space of computationally complex problems: instead of choosing the scattering matrix randomly according to the unitary Haar-measure (which is a crucial assumption for the very hardness of the problem ), one can artificially design test-cases for which pertinent properties of the resulting probability distribution are obtained efficiently on a classical computer. In other words, the boson-sampling device is assessed by a task whose exact or approximate solution is known beforehand.
One approach is based on processes with well-understood average behaviour, such as many-particle quantum walks.In other words, although bosonic statistical effects witness features that go beyond the uncorrelated behaviour of distinguishable particles, the observation of these effects is nevertheless insufficient to strictly guarantee a behaviour in accordance to equation (PB).
The superior computational complexity of many-boson scattering in comparison to the analogous problem with fermions is rooted in the absence of symmetries of the permanent compared to the more benevolent determinant.
Our strategy to devise a solvable instance of boson-sampling therefore relies on imposing artificial symmetries on the scattering matrix , which then allows us to exactly compute the permanent of the pertinent submatrices in a large number of cases.