Accelerating analysis of Boltzmann equations using Gaussian mixture models: Application to quantum Bose-Fermi mixtures

The Boltzmann equation is a powerful theoretical tool for modeling the collective dynamics of quantum many-body systems subject to external perturbations. Analysis of such a model gives access to linear response properties such as collective modes and transport coefficients that characterize a system, but often proves intractable due to computational costs associated with multidimensional integrals describing collision processes. Here, we present a method to resolve this bottleneck, enabling efficient study of the linear response of a broad class of quantum many-body systems whose behavior can be described using a Boltzmann equation. Specifically, we demonstrate that a Gaussian mixture model can accurately represent equilibrium distribution functions and, when combined with the variational method of moments framework, allows efficient computation of collision integrals. We apply this method to investigate the collective behavior of a quantum Bose-Fermi mixture of cold atoms in a cigar-shaped trap, focusing on monopole and quadrupole collective modes above the Bose-Einstein transition temperature. We find a rich phenomenology that spans interference effects between bosonic and fermionic collective modes, dampening of these modes, and the emergence of hydrodynamics in various parameter regimes. Typical spectral functions exhibit Fano interference profiles, systems with dilute fermions manifest behavior reminiscent of Bose polarons, and systems with comparable bosonic and fermionic densities display hallmarks of hydrodynamics such as mode-locking. These effects are readily verifiable with modern cold-atom experiments, and the method developed here opens the door to understanding the collective behavior of many fundamental and technologically-relevant systems.