Phase space sampling with Markov Chain Monte Carlo methods

The efficient exploration of the high-dimensional and multi-modal phase space of scattering events at high-energy particle colliders presents a severe computational challenge. We here discuss the application of Markov Chain Monte Carlo (MCMC) techniques, an approach that so far received rather little attention for this purpose. To this end, we combine the BAT.jl package that provides implementations of a variety of MCMC algorithms with the SHERPA event generator framework. We discuss technical aspects of the implementation and the resulting algorithm and present first results for the process of Z + 3 jets production at the LHC.