A Gaussian moment method and its augmentation via LSTM recurrent neural networks for the statistics of cavitating bubble populations