The Bose-Einstein distribution n_B(E) = 1/(exp((E-μ)/kT) - 1) allows unlimited occupancy of a single-particle state. Unlike fermions, μ must remain less than the ground-state energy, creating a maximum particle density at fixed T. When particle density exceeds this limit, the chemical potential hits zero and a finite fraction of particles condenses into the ground state.
From quantum statistics, you know that identical particles come in two types: fermions (half-integer spin) obeying the Pauli exclusion principle, and bosons (integer spin) that can occupy the same state without restriction. The Bose-Einstein distribution n_B(E) = 1/(exp((E−μ)/kT) − 1) gives the average number of bosons occupying a single-particle state of energy E. The minus sign in the denominator — compared to the +1 for fermions — is what makes all the difference: it means the occupation number can be arbitrarily large when E is close to μ.
The chemical potential μ plays a controlling role. For the distribution to be positive at all energies, the denominator must be positive, which requires E − μ > 0 for all states. If the ground state has energy E₀, then μ must satisfy μ < E₀ at all times. As you add more particles to a fixed-volume system at fixed temperature, μ must increase to accommodate them — but it is bounded above by E₀. At high temperatures, particles spread across many excited states and the constraint is easily satisfied. As temperature drops (or density increases), μ approaches E₀ from below.
The critical point is when μ reaches E₀ exactly: n_B(E₀) diverges. Physically, this signals Bose-Einstein condensation. The thermal occupation of excited states has a maximum value — there is a maximum number of particles that can "fit" into excited states at a given temperature. Any particles above this limit have nowhere to go except the ground state, which they flood with macroscopic occupation. The condensate fraction — the fraction of all particles sitting in the ground state — grows as temperature falls below the critical temperature T_c. Above T_c, no macroscopic occupation exists; below it, a finite fraction occupies a single quantum state.
This condensation is a purely quantum statistical effect with no classical analogue. It does not require interactions — an ideal gas of bosons condenses purely because of quantum indistinguishability and the structure of the Bose-Einstein distribution. The grand partition function, which you used to derive n_B in the first place, captures this transition through the behavior of the fugacity z = exp(μ/kT): z approaches 1 at the condensation point, and the sum over excited states saturates, leaving the ground state to absorb the overflow.