Use policyThe full-text may be used and/or reproduced, and given to third parties in any format or medium, without prior permission or charge, for personal research or study, educational, or not-for-prot purposes provided that:• a full bibliographic reference is made to the original source • a link is made to the metadata record in DRO • the full-text is not changed in any way The full-text must not be sold in any format or medium without the formal permission of the copyright holders.Please consult the full DRO policy for further details. We investigate finite-number effects in collisions between two states of an initially well-known number of identical bosons with contact interactions, oscillating in the presence of harmonic confinement in one dimension. We investigate two N/2 (interacting) ground states, which are initially displaced from the trap center, and the effects of varying interaction strength. The numerics focus on the simplest case of N = 4. In the noninteracting case, such a system would display periodic oscillation with a half harmonic oscillator period (due to the left-right symmetry). With the addition of contact interactions between the bosons, collisions generate entanglement between each of the states and distribute energy into other modes of the oscillator. We study the system numerically via an exact diagonalization of the Hamiltonian with a finite basis, investigating left-right number uncertainty as our primary measure of entanglement. Additionally, we study the time evolution and equilibration of the single-body von Neumann entropy for both the attractive and repulsive cases. We identify parameter regimes for which attractive interactions create behavior qualitatively different from that of repulsive interactions, due to the presence of bound states (quantum solitons), and explain the processes behind this.