Sub-Bath Cooling in Bosonic Systems: Gaussian Constraints and Non-Gaussian Enhancements

Cooling quantum systems with finite resources is a central task in quantum technologies and has been extensively explored in discrete-variable settings. As continuous-variable (CV) platforms play an increasingly important role in quantum information processing, it becomes crucial to understand the fundamental limitations of cooling bosonic systems. In this work, we develop a general framework for cooling CV systems, identifying both the constraints imposed by Gaussianity and the advantages enabled by non-Gaussian interactions. We derive a reachable bound on the cooling performance of Gaussian operations that applies to arbitrary cooling architectures. By optimizing over all protocols saturating this bound, we further identify the most efficient scheme, which minimizes dissipated energy for a given number of ancilla modes. Beyond Gaussian operations, we show that $p$-excitation exchange exploits non-Gaussian resources to achieve a $p$-fold enhancement of the cooling limit. Our results establish the fundamental limits of CV heat-bath algorithmic cooling and reveal the crucial role of non-Gaussianity in surpassing Gaussian cooling barriers.