 Solution

Use SI units and normalize the condensate wavefunction so $|\psi {|}^2=n$. With mechanical momentum operator $-i\hslash \nabla -qA$, the charge current density isFor uniform $n$ and $\psi =\sqrt{n}{e}^{i\vartheta }$, this becomes $j=(nq/m)(\hslash \nabla \vartheta -qA)$.

Under the gauge transformation $A\mapsto A+\nabla \chi$, the wavefunction transforms as $\psi \mapsto e^{iq\chi /\hslash }\psi$, so $\vartheta \mapsto \vartheta +q\chi /\hslash$. The two added gradients cancel in $\hslash \nabla \vartheta -qA$, while $|\psi {|}^2$ stays fixed. Thus the superconducting current is gauge invariant. The same cancellation in the full displayed current proves it without requiring a constant amplitude.

For uniform carrier density in a region with smooth single-valued phase, taking a curl gives the London equation $\nabla \times j=-n{q}^{2}B/m$. In the time-independent setting Ampère's circuital law gives $\nabla \times B={\mu }_{0}j$, and Gauss's law for magnetism gives $\nabla \cdot B=0$. Thereforeso the Helmholtz equation and London penetration depth areThe equation has the modified, exponentially decaying Helmholtz sign. For a planar surface at $x=0$ with material in $x>0$, the bounded solution is $B(x)=B(0)e^{-x/\ell }$; the exponentially growing solution is excluded in the bulk. Hence a macroscopic region far from the surface has negligible field: magnetic flux is expelled from the bulk, the Meissner effect. The literal assertion of zero flux everywhere is the bulk limit: finite penetration remains within distance $\ell$ of a surface. Vortices or flux trapped through a hole require different phase or topology assumptions and are not ruled out by this smooth-phase derivation.