Entanglement and Wigner function negativity of multimode non-Gaussian states

Non-Gaussian operations are essential to exploit the quantum advantages in optical continuous variable quantum information protocols. We focus on mode-selective photon addition and subtraction as experimentally promising processes to create multimode non-Gaussian states. Our approach is based on correlation functions, as is common in quantum statistical mechanics and condensed matter physics, mixed with quantum optics tools. We formulate an analytical expression of the Wigner function after subtraction or addition of a single photon, for arbitrarily many modes. It is used to demonstrate entanglement properties specific to non-Gaussian states, and also leads to a practical and elegant condition for Wigner function negativity. Finally, we analyse the potential of photon addition and subtraction for an experimentally generated multimode Gaussian state.