LAUSR

Advanced quantum technologies require scalable and controllable quantum resources 1 , 2 . Gaussian states of multimode light, such as squeezed states and cluster states, are scalable quantum systems 3 – 5 , which can be generated on demand. However, non-Gaussian features are indispensable in many quantum protocols, especially to reach a quantum computational advantage 6 . Embodying non-Gaussianity in a multimode quantum state remains a challenge as non-Gaussian operations generally cannot maintain coherence among multiple modes. Here, we generate non-Gaussian quantum states of a multimode light field by removing a single photon in a mode-selective manner from a Gaussian state 7 . To highlight the potential for continuous-variable quantum technologies, we first demonstrated the capability to generate negativity of the Wigner function in a controlled mode. Subsequently, we explored the interplay between non-Gaussianity and quantum entanglement and verify a theoretical prediction 8 about the propagation of non-Gaussianity along the nodes of photon-subtracted cluster states. Our results demonstrate large-scale non-Gaussianity with great flexibility along with an ensured compatibility with quantum information protocols. This range of features makes our approach ideal to explore the physics of non-Gaussian entanglement 9 , 10 and to develop quantum protocols, which range across quantum computing 11 , 12 , entanglement distillation 13 and quantum simulations 14 . Continuous-variables quantum information processing requires non-Gaussian states and operations. The generation of non-Gaussian quantum states of a multimode field is now reported through a mode-selective photon-subtraction scheme