A. Neumaier said:
But the Boltzmann equation is valid only for very weak interactions.
... where the approach to equilibrium is governed by a different extremal principle (the grand potential is minimized), not by max entropy. Only the final equilibrium result is the same.
This shows that the max entropy principle cannot be regarded as being fundamental.
I don't understand what you mean. The grand canonical ensemble is given by the stat. Op.
$$\hat{\rho}=\frac{1}{Z} \exp(-\beta \hat{H}+\mu \beta \hat{Q}),$$
for ##\hat{H}## the Hamiltonian and ##\hat{Q}## some conserved charge (you can also have more conserved charges and more chemical potentials, but that's of course not much different). This stat. Op. follows from the maximum entropy principle with ##U## and ##Q## as the constraints. The independent variables are ##T=1/\beta## and ##\mu## (and ##V## which you can introduce as a "quantization volume" with periodic boundary conditions).
The associated potential usually used is the grand potential
$$\Omega(V,T,\mu)=-T \ln Z,$$
fulfilling
$$\mathrm{d} \Omega = -p \mathrm{d} V + S \mathrm{d} T -N \mathrm{d} \mu.$$
Equilibrium is characterized by a minimum of ##\Omega## (corresponding to fixed ##V##, ##T##, and ##\mu##).
You can Legendre transform to other thermodynamic potentials with other "natural" independent thermodynamic quantities, which characterize equilibrium as minima (or maxima) keeping these quantities fixed.
Nevertheless the grand-canonical stat. op. is determined by the maximum-entropy principle under the given constraints, ##\langle H \rangle=U## and ##\langle Q \rangle=N## fixed.