Photon number resolving detectors are the ultimate measurement of quantum optics, which is the reason why developing the technology is getting significant attention in recent years. With this arises the question of how to evaluate the performance of the detectors. We suggest that performance of a photon number detector can be evaluated by comparing it to a multiplex of on-off detectors in a practical scenario: conditional preparation of a photon number state. Here, both the quality of the prepared state and the probability of the preparation are limited by the number of on-off detectors in the multiplex, which allows us to set benchmarks that can be achieved or surpassed by the photon number resolving detectors.

We prove experimentally the predicted existence of a three-qubit quantum state with genuine multipartite entanglement which can be certified solely from its separable two-qubit reduced density matrices. The qubits are encoded into different degrees of freedom of a pair of correlated photons and the state is prepared by letting the photons to propagate through a linear optical circuit. The presence of genuine multipartite entanglement is verified by finding numerically a fully decomposable entanglement witness acting nontrivially only on the reductions of the global state. Our result confirms viability of detection of emerging global properties of composite quantum systems from their parts which lack the properties.

Manipulating light by adding and subtracting individual photons is a powerful approach with a principal drawback: the operations are fundamentally probabilistic and the probability is often small. This limits not only the fundamental scalability but also the number of operations that can be applied in realistic experimental settings. We propose and analyze an adaptive technique which can significantly increase the probability of success while preserving the quality of the photon subtraction. We show the improvement both in single mode preparation and manipulation of non-Gaussian states with negative Wigner functions and in two-mode entanglement distillation protocol with Gaussian states of light.

Genuine multipartite entanglement underlies correlation experiments corroborating quantum mechanics and it is an expedient empowering many quantum technologies. One of many counterintuitive facets of genuine multipartite entanglement is its ability to exhibit an emergent character. That is, one can infer its presence in some multipartite states merely from a set of its separable marginals. Here we show that the effect can also be found in the context of Gaussian states of bosonic systems. Specifically, we construct examples of multimode Gaussian states carrying genuine multipartite entanglement which can be verified solely from separable nearestneighbor two-mode marginals. The key tool of our construction is an entanglement witness acting only on some two-mode reductions of the global covariance matrix, which we find by a numerical solution of a semidefinite program. We also propose an experimental scheme for preparation of the simplest three-mode state, which requires interference of three correlatively displaced squeezed beams on two beam splitters. Besides revealing the concept of emergent genuine multipartite entanglement in the Gaussian scenario and bringing it closer to experimentally testable form, our results pave the way to effective diagnostics methods of global properties of multipartite states without complete tomography.

Numerical simulation of continuous variable quantum state preparation is a necessary tool for optimization of existing quantum information processing protocols. A powerful instrument for such simulation is the numerical computation in the Fock state representation. It unavoidably uses an approximation of the infinite-dimensional Fock space by finite complex vector spaces implementable with classical digital computers. In this approximation we analyze the accuracy of several currently available methods for computation of the truncated coherent displacement operator. To overcome their limitations we propose an alternative with improved accuracy based on the standard matrix exponential. We then employ the method in analysis of non-Gaussian state preparation scheme based on coherent displacement of a two mode squeezed vacuum with subsequent photon counting measurement. We compare different detection mechanisms, including avalanche photodiodes, their cascades, and photon number resolving detectors in the context of engineering non-linearly squeezed cubic states and construction of qubit-like superpositions between vacuum and single photon states.