It is often said in the physics papers or textbooks that in equilibrium if there is no external potential then number density is uniform $n = N/V$, where $N$ is number of particles, $V$ is volume. I would like to see if this follows from statistical mechanics (or maybe contradicts statistical mechanics).
Let us assume that we are given $N$ particles under Hamiltonian that has only pairwise interaction.
$$ H = \sum \limits_{i=1}^{N} \frac{\vec{p}_i^2}{2m} + \sum \limits_{i=1}^{N} \sum \limits_{j=i+1}^N \Phi(|\vec{x}_i - \vec{x}_j|)$$
Assume that distribution in the phase space is given by Boltzmann distribution.
$$ f(\vec{x}_1, ..., \vec{x}_N, \vec{p}_1, ..., \vec{p}_N) = \begin{cases} \frac{1}{Z} e^{-\beta H(\vec{x}_1, ..., \vec{x}_N, \vec{p}_1, ..., \vec{p}_N)} \quad \textrm{ if all } \vec{x}_i \in V \\ 0 \textrm{,} \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \textrm{else.} \end{cases} $$
Define number density as the following quantity.
$$ n(\vec{x}) = N \int \delta(\vec{x}-\vec{x}_1) f(\vec{x}_1, ..., \vec{x}_N, \vec{p}_1, ..., \vec{p}_N) \prod \limits_{i=1}^N \mathrm{d}\vec{x}_i \mathrm{d}\vec{p}_i$$
Because of properties of Boltzmann distribution, it is possible to integrate out all momentum coordinates. In that case, we have the following simplification, where $Q_N$ is defined in the following equation.
$$ n(\vec{x}) = \frac{N}{Q_N} \int \limits_{V^N} \delta(\vec{x} - \vec{x}_1) e^{-\beta \sum \limits_{i=1}^{N} \sum \limits_{j=i+1}^N \Phi(|\vec{x}_i - \vec{x}_j|)} \prod \limits_{i=1}^N \mathrm{d}\vec{x}_i$$
$$ Q_N = \int \limits_{V^N} e^{-\beta \sum \limits_{i=1}^{N} \sum \limits_{j=i+1}^N \Phi(|\vec{x}_i - \vec{x}_j|)} \prod \limits_{i=1}^N \mathrm{d}\vec{x}_i$$
$\textbf{Question}$: Is it possible using these definitions and assumptions to show that in this Boltzmann equilibrium model number density is uniform (or maybe that it is false that it is uniform)?