One of the ways to determine the energy is the family of ab initio methods: the approximate solution of the Schrödinger equation for a system of atomic nuclei and electrons. One approach to such a solution is to approximate the wave function with the Slater determinant 1 and calculate the self-consistent field for this approximate equation (Hartree-Fock method). The computation involves two-electron integrals corresponding to the electron-electron interaction and has a high computational complexity. Due to that, and also a large number of local minima, global optimisation of such a function is not practical.