Variational quantum algorithm for Gaussian discrete solitons and their boson sampling