 Solution

For a multiplication operator $A_j(x)$, $[p_i,A_j]=-i\hslash {\partial }_i{A}_j$. Therefore the mechanical momentum satisfiesThe displayed ladder operators implicitly require $qB>0$, so that their denominator is real and the second is the adjoint of the first. Then $[a,a^{†}]=1$ andThe free longitudinal motion gives a continuous spectrum in $k_z$, and the transverse Landau levels have guiding-centre degeneracy of a Landau level. For arbitrary sign of $qB$, use $|qB|$ in the energy and $a=({\pi }_x+i\mathrm{sgn}(qB){\pi }_y)/\sqrt{2\hslash |qB|}$. The question's normalizable wavefunctions likewise assume $qB>0$.

For the stated symmetric-gauge wavefunction, $|\psi {|}^2\propto r^{2M}{e}^{-qB{r}^2/(2\hslash )}$, whose peak for $M>0$ is at $r_M^2=2\hslash M/(qB)$. Thus the largest useful angular index is about $qB{R}^2/(2\hslash )$, givingThe $M=0$ state should also be included in the usual lowest Landau level; its omission by “positive integer” changes only an order-one term. Including the radial measure shifts the most probable radius to $r^2=(2M+1)\hslash /(qB)$, again only an order-one counting correction. Boundary conditions and the tails of the wavefunctions make this a large-flux estimate, not an exact finite-disc count.