Hi Don,
the BCS state of superconductivity is not an eigenstate of the
particle
number operator. Thus there is a uncertainty DeltaN in particle
number.
In the thermodynamic limit (BCS theory for infinite number of
electrons)
DeltaN -> \infty but DeltaN/N -> 0. As a consequence of number-phase
uncertainty relation the phase \phi becomes sharply defined in the
N-> \infty limit (This leads to a macrocopic quantum state).
It is probably true that a different statistical framework (other
than
grand-canonical in standard BCS theory) will - in the thermodynamical
limit - lead to identical results. I have never seen a concrete work
using
this alternative formulation (but it might exist).
The symmetry that is broken (typically) in the superconducting state
is
gauge symmetry, as you wrote. It is not that intuitively accessible
as
broken rotational symmetry in the ferromagnet but it is analogous.
The simplest formulation to see this is Ginzburg-Landau theory. Above
T_c the order parameter vanishes, the phase is undefined and
consequently one has perfect gauge symmetry. Below Tc a finite
order parameter exists and consequently the phase must take a
definite
value (even if it is arbitrary). Gauge symmetry is broken.
Best
Ulf