I'm trying to understand the algebraic method of formulating the Landau level problem better. I'm referring to David Tong's notes on the Quantum Hall effect for this (but not exactly following his notation). The relevant section to this discussion is Chapter 1, pages 14-17 and 22-24.
Working in natural units (I set the magnetic length $l_B=1$ and the cyclotron frequency $w_c=2$) the Landau level Hamiltonian $ \hat{H} = (\mathbf{\hat{p}}-\mathbf{\hat{A}})^2$ is rewritten as a quantum harmonic oscillator (QHO) Hamiltonian: $ \hat{H} = 2\left(\hat{b}^\dagger \hat{b}+\frac{1}{2}\right) $ where the bosonic operators are given by $\hat{b}=\frac{1}{\sqrt{2}}(\hat{\pi}_x-i\hat{\pi}_y)$ and $\hat{b}^\dagger=\frac{1}{\sqrt{2}}(\hat{\pi}_x+i\hat{\pi}_y)$ and $\mathbf{\hat{\pi}}=\mathbf{\hat{p}}-\mathbf{\hat{A}}$ is the canonical momentum.
The usual Fock space ladder of states and energies of the QHO however does not account for the degeneracy of the Landau levels. To account for this, Tong constructs another set of non gauge-invariant momenta $\mathbf{\hat{\Pi}}=\mathbf{\hat{p}}+\mathbf{\hat{A}}$ (see Eq. 1.28 in the notes) and shows that another set of bosonic operators (let me call these $\hat{d}, \hat{d}^\dagger$) constructed using these momenta commutes with the Hamiltonian (see Section 1.4.3). Specifically, the Landau level wavefunctions are then given by: $$ |n,m> = \frac{(\hat{b}^\dagger)^n(\hat{d}^\dagger)^m}{\sqrt{n!m!}}|0,0> $$ where $n,m$ are the eigenvalues of the number operators $\hat{b}^\dagger\hat{b}$ and $\hat{d}^\dagger\hat{d}$ respectively. The levels are degenerate because the energies of the Landau levels are still labeled only by $n$ (see Eq. 1.16 in the notes): $$ E_{n,m} = 2\left(n+\frac{1}{2}\right) $$
Crucially, however, (and as he remarks) this analysis only holds for the symmetric gauge. The $\hat{\Pi}$ operator only commutes with the Hamiltonian in this gauge. This is somewhat troubling to me.
How do we show the degeneracy of the Landau levels in other gauges---let's say the Landau gauge within the algebraic approach? What is the operator there that commutes with the bosonic operators of the QHO Hamiltonian?