Momentum distribution for a one-dimensional trapped gas of hard-core bosons

Using the exact $N$-particle ground state wave function for a one-dimensional gas of hard-core bosons in a harmonic trap we develop an algorithm to compute the reduced single-particle density matrix and corresponding momentum distribution. Accurate numerical results are presented for up to N=8 particles, and the momentum distributions are compared to a recent analytic approximation.