Computing an orthonormal basis of symmetric or antisymmetric hyperspherical harmonics
A numerical method to build an orthonormal basis of properly symmetrized hyperspherical harmonic functions is developed. As a part of it, refined algorithms for calculating the transformation coefficients between hyperspherical harmonics constructed from different sets of Jacobi vectors are derived and discussed. Moreover, an algorithm to directly determine the numbers of independent symmetric hyperspherical states (in case of bosonic systems) and antisymmetric hyperspherical-spinisospin states (in case of fermionic systems) entering the expansion of the A-body wave functions is presented. Numerical implementations for systems made with up to five bodies are reported.