In this note, we introduce a family of “power sum” kernels and the corresponding Gaussian processes on symmetric groups Sn. Such processes are bi-invariant: the action of Sn on itself from both sides does not change their finite-dimensional distributions. We show that the values of power sum kernels can be efficiently calculated, and we also propose a method enabling approximate sampling of the corresponding Gaussian processes with polynomial computational complexity. By doing this we provide the tools that are required to use the introduced family of kernels and the respective processes for statistical modeling and machine learning. Bibliography: 20 titles. © The Author(s), under exclusive licence to Springer Nature Switzerland AG 2025.