Entanglement in N-harmonium: bosons and fermions

The entanglement of the N-harmonium system (i.e., a completely-integrable model of N particles where both the confinement and the two-particle interaction are har-monic) is shown to be analytically determined in terms of N and the relative interaction strength. First, the entanglement of the N-boson system via the von Neumann entropy is computed by use of the corresponding natural occupation numbers. Below a critical value Nc of N it is found that, for positive values of the coupling constant the entan-glement grows when the number of particles is increasing; opposite situation occurs for N> Nc. Then the entanglement of the N-fermion system is investigated by means of the linear entropy, which is a linearized form of the von Neumann entropy of the one-body reduced density matrix. We have shown that the entanglement grows when the number of particles is increasing in the region of negative coupling. For positively increasing values of the coupling the entanglement as a function of N varies according to various regimes. Moreover, in the strong coupling regime the entropy of the system decreases when N is growing. For fixed N, the entanglement is found (a) to decrease (increase) for negatively (positively) increasing values of the coupling constant, and (b) to grow when the energy is increasing. Finally, let us point out that the spatial and spin contributions to the total entanglement are of comparable size