In theoretical studies of Bose-Einstein condensates, it's common to look at the one-body density matrix: $$ n_{ab}=\langle\hat{c}_a^\dagger\hat{c}_b\rangle $$ where the $\hat{c}_j$ are annihilation operators for some single-particle modes. (Often, it's written in the position basis as $n(\mathbf{x},\mathbf{y})=\langle\hat{\psi}^\dagger(\mathbf{x})\hat{\psi}(\mathbf{y})\rangle$, but I'm working in a discrete basis.)
The von Neumann entropy of a system is defined as: $$S=-\text{Tr}\left\{\hat{\rho}\log\left(\hat{\rho}\right)\right\}$$ where $\hat{\rho}$ is the full many-body density matrix. I'm working with a method in which I can't tractably obtain the full many-body density matrix, but frequently use the one-body density matrix. Consider the quantity: $$S^{(n)}=-\text{Tr}\left\{\mathbf{n}\log\left(\mathbf{n}\right)\right\}$$ where $\mathbf{n}$ is the one-body density matrix, with matrix elements $n_{ab}$ as defined in my first equation. In the case of a single particle, the one-body density matrix is the full density matrix, and this would give us the von Neumann entropy. What about for a many-body system? Is it meaningful in any way then? Intuition makes me suspect that it's something like a lowest-order approximation to the full von Neumann entropy that discards many-body effects (and my needs make me hope so!), but I'm not sure how I'd go about formally showing this, or deducing its meaning.