Matrix phase-space representations in quantum optics

We introduce matrix quantum phase-space distributions. These extend the idea of a quantum phase-space representation via projections onto a density matrix of global symmetry variables. The method is applied to verification of low-loss Gaussian boson sampling (GBS) quantum computational advantage experiments with up to 10,000 modes, where classically generating photon-number counts is exponentially hard. We demonstrate improvements in sampling error by a factor of 1000 or more compared to unprojected methods, which are infeasible for such cases.