Final answer:
The reasoning that the canonical expression for pressure should be evaluated at the mean number of particles <N> is correct, since this mean is representative in the Grand Canonical Ensemble. Density fluctuations are derived by differentiating the expression of pressure with respect to the chemical potential. Textbooks sometimes assume the reader will interpret N as its mean value without explicitly stating so.
Step-by-step explanation:
Your teacher's reasoning appears to be correct; when dealing with density fluctuations in the context of the Grand Canonical Ensemble, it is standard to evaluate expressions such as pressure at the mean particle number, <N>, especially in the thermodynamic limit where this mean value becomes representative of the system.
In statistical mechanics, particularly for a system in equilibrium with a reservoir that allows for exchange of heat and particles, the Grand Canonical Ensemble is the appropriate framework. Here, fluctuations about this mean number are central to the physics of the system.
A canonical expression for pressure is given by P = -(∂A/∂V)_N,T, where A is the Helmholtz free energy, and this expression should be differentiated twice with respect to the chemical potential μ to derive the density fluctuations. However, throughout these derivations, the number of particles N appears variably and must eventually be replaced by its average <N> to correctly represent the equilibrium properties of the ensemble.
As for literature like Huang's, it is not uncommon for textbooks in statistical mechanics to gloss over some of the subtleties in the derivations, assuming that the reader will interpret N as the mean number of particles <N> in the context of ensemble averages. A rigorous approach would indeed highlight when approximations are made or when variables are to be taken as mean values and not as fluctuating quantities.