On the direct diagonalization method for a few particles trapped in harmonic potentials

We describe a procedure to systematically improve direct diagonalization results for few-particle systems trapped in one-dimensional harmonic potentials interacting by contact interactions. We start from the two-body problem to define a renormalization method for the interparticle interactions. The procedure is benchmarked with state-of-the-art numerical results for three and four symmetric fermions.